模拟物理-05第三章.pptVIP

第三章 边值问题和本征值问题 线性二阶方程 物理学中许多重要的微分方程可以写成线性二阶方程的形式 Poisson方程 在球对称情况下为 做一个变换 边值问题 看起来，可以变换成两个一次方程，然后用第二章的解法来做积分。 但是，边值问题不同于初值问题。 初值问题是在初始时刻给定初条件，由它出发开始迭代。 边值问题是给定不同自变量上面对于因变量的一些约束。不是在一个自变量上的约束。 内容 推导一个特别适合于此线性微分方程的积分算法。 讨论边值问题。 3.1 Numerov算法 本节推导一个可以用来求解线性二阶微分方程的简单并且高效的方法。 这个方法一般被叫做Numerov方法或者Cowling方法。 由二阶导数的差分开始 由上述泰勒展开式形式进行组合得到 不忽略四次导数项，利用方程和差分得 把y”和y””它们带入微分方程。 递推公式 整理得到 其局部精度是O(h6)。比四阶Runge-Kutta方法要高一阶。 此方法的效率也更高，每步只需要计算k2和S。 例子3.1： 1. 对比标准方程： k2=4π2, S=0 2. 分析迭代需要的条件：y(0)=1 由于y’(0)=0，y(h)=y(0)=1 例子3.1： 开始 输入h 赋初值x=0, y0=1, y1=1 迭代 y2=… 更新格点：x=x+h y0=y1, y1=y2 X==1 输出y No 结束 Yes 流程图 注1：h必须保证[0.0, 1.0]被均分。 例子3.1： 精确解 作业： 比较Numerov方法与四阶Runge-Kutta方法的效率与精度 采用上一章的方法，必须引入一个新变量，把二阶方程转化为一阶方程组。 比较效率的一个方法是记录计算所需要的时间。 参见例子 记录程序运行时间 3.2 边值问题的直接积分 考虑电荷分布为 Poisson方程为 精确解 对方程的两边进行积分得 边条件为 r=0时， 时， 常数 精确解为 迭代关系 方程 其中 由 迭代的条件 在坐标原点 在r=h 此处有精确解，我们用精确解给出 否则，我们去计算 的值。利用Taylor展开得到 的近似值。 流程图 为了效率： 记住S在两步上的值。 把系数放在循环外。 开始 设定h 设定初值，由精确解给出phi(0), phi(1) 计算s0,s1 迭代phi(i)= 程序 误差 误差随着r增大 一般情况下的初值 如果没有精确解，需要Taylor展开。 由 得 通过数值积分计算 ， 从而可以求第一步的近似值。 叠加原理 球对称性 初值引入的误差 数值积分总包含误差。 设误差为5%。考虑它的影响。 修改程序，在phi（1）精确值基础上改变。 修正方法 理解误差起因。 齐次方程 的解可以加到方程的特解上得出另一个解。 齐次方程的两个线性独立的解为 第二个解是有物理意义的。 第二个解的误差不随着r增大。但是第一个解的误差随着r线性增大。 修正方法 进行修正的方法就是，从数值结果中减去齐次方程的坏的非物理的解的一个倍数，以保证渐进区域内的物理行为。 针对误差随着r线性增大的规律，把数值解的最后10个点拟合于直线 然后从数值解中减去mr。 向内积分 由于 当r很大时 从外向内积分也是可行的 避免了不得不进行一次求积以决定phi(1)的近似启动值。 作业 习题3.2 对于本节的电荷分布，采用Numerov方法从大r值向内积分。 有些同学迷信基础，其实这是不对的。与其说是基础重要不如说是方法重要。基础好的同学更多的是方法好。很有可能遇到基础知识不足的时候，他会快速地进行回顾。迷信vs坚韧 * * 有些同学迷信基础，其实这是不对的。与其说是基础重要不如说是方法重要。基础好的同学更多的是方法好。很有可能遇到基础知识不足的时候，他会快速地进行回顾。迷信vs坚韧 * *