luyi07的笔记

摘要： numpy.loadtxt 可以载入 txt 数据。 matplotlib.pyplot.scatter 可以画散点图 matplotlib.pyplot.quiver 可以画箭头 matplotlib.pyplot.xlabel, matplot.pyplot.ylabel 可以增加 x, y 轴 阅读全文

摘要： numpy.loadtxt 可以从 txt 中读取数据。 numpy.mean 可以计算所有数据的平均数。 matplotlib.pyplot.hist 可以画柱状图，并且返还各个区间的统计数。 示例代码和效果图如下 import numpy as np import matplotlib.pypl 阅读全文

摘要： 1. patsubst $(patsubst PATTERN,REPLACEMENT,TEXT) 在 TEXT 中以空格分开的单词中，寻找 PATTERN，然后替换成 REPLACEMENT，然后返回替换后的单词。 2. 自动化变量 $@ ：表示规则中的目标文件集。 \(% ：当规则的目标文件是一个 阅读全文

摘要： 这是以前写的笔记，散在不同的地方，现在收集在这个随笔里，供以后自己和合作者翻阅。写了这个笔记以后，我还会在相应代码里加上注释，因为没有人能记住大串公式的每一个细节，所以只能把笔记、博客用作外部的记忆存储器。 1. 壳模型相互作用: pn格式 在 pn 格式下，单体 + 2 体相互作用写作 \[ \b 阅读全文

摘要： 我们写的 PVPC 的代码在 gitee 上可以下载：https://gitee.com/luyi07/pvpc.git 这里摘录一个脚本，用来自动做变分+投影，在最新的版本里会放在 example/run.sh 中，运行方法： bash run.sh #!/bin/bash nucleus=Cr4 阅读全文

摘要： BigStick 公开版本可在 https://github.com/cwjsdsu/BigstickPublick.git 下载。其中 docs/ 与 util/docs/ 下都有详细的文档说明，这里只是记下简单运用的步骤。 下面的脚本是根据 util/example/ 里面的脚本修改的，可以自动 阅读全文

摘要： 若有一个下三角阵，不包括对角元，其矩阵元为 \[ A_{ij}, ~~~~ i = 1, \cdots, n-1, ~~~~ j=0,\cdots, i-1. \] 那么，可以用一维数组储存这些矩阵元： \[ a[k] = A_{ij}, ~~~ k = i(i-1)/2 + j. \] 这样可以将 阅读全文

摘要： BCS 方程和 Bogoliubov 是处理对关联的经典内容。一半是学习，一半是为了教学，这里我写下了较为完整的推导过程。我欣赏 Bogoliubov 变换，它把 BCS 波函数变为准粒子真空，在此基础上，构造的准粒子激发可以看做是准粒子的自由气体，这很漂亮，虽然是在粒子数不守恒的 BCS 框架下， 阅读全文

摘要： 鼓面振动就是一种圆膜振动。知了的发声也是靠腹部的薄膜振动。 1. 圆膜振动的解析结果 圆膜振动的数学物理方程、及其求解如下。 2. 圆膜振动的可视化 已知上面的解，可以用计算机编程，把不同模式的圆膜振动画成动图，甚至还可以带上声音。 对于给定的$n,i$，已经确定一个振动模式， \[ u_{ni} 阅读全文

摘要： 这都是教科书上都有的内容，我只是整理整理，把一些约定也统一一下，方便以后写代码的时候参照。 1. 约定 1.1 约化矩阵约定 \[ \langle J'M' | s \sigma | J M \rangle = \langle J' || s || J \rangle (JM,s \sigma| J 阅读全文