分类： 其他代码

一、定义 正格子基矢：。 倒格子基矢：。 需要满足关系：。 1. 三维的情况 倒格子基矢可以写为以下形式： 2… Continue reading 倒格子基矢的计算（附数值计算、符号计算Python代码）

施密特正交化（Schmidt orthogonalization）过程[1]： Python代码： 运行结果：… Continue reading 施密特正交化（附Python代码）

矩阵的特征向量是有朝向的，通常写为列向量。这里给出特性向量在矩阵中朝向的验证。 先给出结论： Python中n… Continue reading 验证特征向量在矩阵中的朝向（附Python、Matlab、Mathematica、Fortran代码）

这里给出结论[1-4]： 矩阵乘积：时间复杂度为 矩阵求逆：时间复杂度为 矩阵本征值：时间复杂度为 Pytho… Continue reading 矩阵一些运算的时间复杂度

Ising （伊辛）模型为： 这里要用到Metropolis采样，可看这篇文章：Metropolis采样 （附… Continue reading 蒙特卡洛模拟Ising模型（附Python代码）

采样的一个简单应用见这篇博文：使用蒙特卡洛计算定积分（附Python代码）。 一般的简单分布采样有：均匀分布、… Continue reading Metropolis采样 （附Python和Matlab代码）

网上有很多介绍蒙特卡洛计算定积分。因为这是蒙特卡洛方法的经典案例，所以我这里也摘抄整理下。 一般来说，数值计算… Continue reading 使用蒙特卡洛计算定积分（附Python代码）

三体一般来说没有解析解，只有几个特殊的初始条件才有解析解。数值解相对来说就比较简单了，只要套用万有引力公式即可… Continue reading 在二维平面模拟三体运动（附Python代码）