高斯求积代码matlab 计算物理 该存储库包含计算物理课程中讲授的概念，并已翻译为MATLAB代码。 话题 非线性方程组的牛顿拉夫森法 高斯消除 多项式插值，线性/二次/三次样条插值，Thomas算法 多项式曲线拟合（在某些文件中可以使用正交多项式，例如Legendre，Hermite，Chebyshev等） 数值积分：矩形，梯形，辛普森1/3，高斯求积 微分方程组的数值解： 欧拉法 龙格·库塔（Runge Kutta）二阶 龙格·库塔（Runge Kutta）4级 带边界条件的微分方程的射击方法 偏微分方程： 有限差分法：前进，后退，中央方案 Dirichlet边界条件，Newmann边界条件，重影边界点 蒙特卡洛技术 生成满足概率分布的数字 随机行走问题，布朗运动等 经典顺磁性的蒙特卡罗模拟 n维空间中的蒙特卡洛积分 伊辛模型

努梅罗夫 python脚本，用于解一维时间独立的Schrodinger方程的束缚态。 该脚本使用Numerov方法来求解微分方程，并显示所需的能级和带有这些能级中每一个的近似波动函数的图形。 跑步 要运行此代码，只需克隆此存储库并使用python运行Numerov.py脚本（需要numpy和matplotlib模块）： $ git clone https://github.com/FelixDesrochers/Numerov/ $ cd Numerov $ python Numerov.py 然后程序将要求您输入要显示的能级数和所需的电势（确保电势大约位于x = 0的中心）： $ >> Which first energy levels do you want (enter an integer) : 4 $ >> Potential (as a fonction of x):

双摆运动动画 使用Julia和Python的双摆的混沌运动动画。 julia版本是在Jupyter笔记本环境中实现的，而python版本是常规脚本。 可以使用此代码模拟任意长度，质量，初始位置和初始速度的双摆。 两种版本的代码均允许以mp4或gif格式保存动画。 Julia版本 该代码可以将动画另存为gif文件。 它还可以在Jupyter笔记本环境中播放视频。 以相等的时间步长录制动画，以保留视频中摆的相对速度。 Python版本 使用matplotlib的FuncAnimation（）完​​成动画。 该代码可以将动画另存为mp4文件。 python脚本已在spyder环境中进行了测试。 应该启用在外部窗口中进行绘制以查看动画。 Julia版本的包装要求 情节 微分方程 迪瑞克 Python版本的包装要求 麻木 scipy.integrate（odeint） pylab，matplotlib（动画） imagemagick（将动画另存为gif文件） 利用 摆的质量和长度，重力加速度常数，初始条件和模拟持续时间在代码的开头定义。 通过更改这些参数可以模拟不同的情况。 样品输出