计算物理实验报告常微分方程.docxVIP

常微分方程的边值问题和本征值问题 一、问题描述 j2 -~~= ~k 卩(JT = 0) = 0* (p[X — 1) = 0 dy 利用搜索法和弦割法，得到该常微分方程的本征值，再利用打靶法计算多个本征值。 二、解决方法 (一)搜索法 1?先随便猜测k的一个试验值，程序中令k=1 2.由Numerov算法 =—(S^+i + 105^ + Sn—1 ) + O(h^) 根据本题的条件，kn+1=kn=kn-仁k,s=O,得到yn+2,yn+1,yn间的迭代公式 令 con=(k*h)A2/12 yn+2=2*(1 -5*con)*yn+1/(1+con) -yn 3自己给定$的初始条件，然后利用公式得到边界值$ (1) 然后以小的步长 dk增加k值，这里令dk=1,每当$ (1)改变符号时，就将步长减半后倒退回 来重复 当步长小于所要求的容许误差时终止程序，此时的 k值即为所求。 (二)弦割法 随便猜测两个 k值，这里令k0=1,k仁2 自己给定$的初始条件，对两个 k值分别利用上述公式进行迭代，得到边界值 y1(1)和 y2(1)。 比较y1(1)和y2(1)的绝对值大小。若绝对值大，说明对应的 k值距离本征值距离较远。 将(k0+k1) /2赋给k2,边界值绝对值小的对应的 k值保持不变，边界值绝对值大的对应 k值重新定位k2的值。 重复进行实验，当 y1(1)和y(2)的差