第三章常微分方程的边值和本征值问题 .pptxVIP

第三章常微分方程的边值和本征值问题;本章内容;3、0 边值问题与本征值问题;本征值问题;物理学中许多重要的微分方程具有如下形式;例如泊松方程;作变换;量子力学中,中心势场 V(r) 中运动的粒子波函数满足定态薛定谔方程;关于该方程,我们感兴趣的是对哪些能量本征值E, 能够导出满足适当边界条件的物理上能够接受的非零解, 这个问题就是一个本征值问题;3、1 Numerov 算法;对 yn+1 或者 yn-1 解这个线性方程,就提供了一个对 x 向前或者向后积分的递推关系,其局部误差为 O(h6)、;注意这个算法比四阶Runge-Kutta算法高一个精度,而且Numerov算法更有效率,因为每一步只需要在一个格点上计算 k2 和 S。;例子——利用Numerov算法解初值问题;而 y1 是未知的,需要一个一步迭代格式来产生 y1 ,例如可选择 Euler 方法或 Taylor 级数展开,并利用初始条件;注意 Numerov 算法与前面所讲算法的区别;3、2 边值问题的直截了当积分;其中;为了求得 φ1,直截了当对方程积分;计算结果发现,当 r 增大时, φ的误差变大;这个渐进方程有两个线性独立的解 ;当 r 特别大时,位势 Φ→ r -1,因而 ? ~ 1 。;已知 Φ→ 1+ br, b 为未知常数,需要从解中减去 br 。;线性修正后得到的解;关于本例,我们也能够采纳向后积分的迭代格式来实施直截了当积分, 即从 r 特别大处(例如 r=20 )出发,取 ?n+1=?n=1 ,然后向后积分。;3、3 打靶法求边值问题;如此关于给定的参数 δ ,我们就能够通过积分这个初始问题得到 yδ(b) 、;一般来说,由于可调参数 δ 的随意选择, yδ(b) 和 yb 特别难相等。;yδ(b)= yb ;例子;关于边值条件 y’(0)=a, y(1)=b ,我们能够选择 y(0) 的值为可调参数 δ ,即 y(0) = δ ,如此就构成了一个含参数的初始问题,然后通过使用一个搜索算法去调整参数 δ ,使数值解在误差范围内等于 y(1) 。;3、4 打靶法求解本征值问题;策略:我们先推测一个试验本征值 k,同时任取一个非零数δ, 把微分方程变化为一个初始值问题;注意:试验本征值 k 是一个可调参数,而参数 δ 只是一个任意选定的辅助参数,它的任意性是由于解 的不唯一性引起的,并不影响本征值的求解,一般来说它能够由本征函数的归一化来确定。;3、5 一维薛定谔方程的定态解;定解问题;算法的思想为:由于我们要求的是一个束缚 态解,因此取一个负的试验本征值,从 xmin 出发向前直截了当积分,能够产生一个数值解 ψ。;因此比 较明智的做法是,在每一个试验本征值上,由 xmax 出发向后直截了当积分产生另一个数值解 ψ。;这个地方的ψ 是一个方便的标尺;先将Schro?dinger 方程做无量纲处理,能够写成;感谢您的聆听！