《计算物理学》课件 第二章常微分方程.pptx

第二章常微分方程的数值求解方法;1、引言

2、数值求解的基本思想

3、欧拉方法

4、R-K方法

5、线性多步法

6、预估校正方法

7、步长选择

8、常微分方程组和高阶微分方程

9、边值问题

10、有限差分方法;1、引言;解析求解?;解析求解无能为力;根据常微分方程的定解条件，常微分方程可分为初值问题和边值问题。

如果定解条件是描述函数在初始点的状态，则称为初值问题。一个典型的常微分方程初值问题的形式如下：;常微分方程初值问题的数值求解;3、欧拉（Euler）方法;?;三种简单的数值解法;?;数值积分法;?;Talyor级数展开;欧拉公式及其几何意义;欧拉方法的稳定性分析;?;总体误差和局部误差的关系;?;?;收敛性和稳定性是两个不同的概念。

收敛性是指数值方法的截断误差对计算结果的影响

稳定性指的是某一步的计算误差对计算结果的影响

此外，稳定性与步长密切相关。;;;?;RUNGE-KUTTA方法;?;?;?;?;?;?;?;?;?;?;一步法在计算时只用到了前面一步的近似值，但要提高精度，需要增加中间函数值的计算，这就加大了计算量。;线性多步法;?;?;?;?;?;?;;;步长的自动选择;小步长高精度;;;;?;二分法选取步长;常微分方程组和高阶微分方程;边值问题;打靶法（shootingmethod）;?;?;?;?;?;?;;有限差分方法;在第一步中，通过所谓的网络分割法，将函数定义域分成大量相邻而不重合的子区域，通常采用的是规则的分割方式。

网络线划分的交点称为节点。

;差分形式;一、二阶微商;?;小结;;;;;;;;;;;;;;;;;例：BicycleRacing;;Considertheairresistance;Intheabovecalculations,theinitialvelocityisnotzero.

Whathappensifv0=0?;例：Thetrajectoryofacannonshell;Initialspeed=700m/s;Withairresistance;Equationofmotion:;300;ConsiderAirdensityvs.height;;2.量子力学中描述一维无限深势阱的微分方程为：