南京大学计算物理课程课件5.pptVIP

建立差分格式： 流程图 开始 输入N,M,V,H 初始边界 条件赋值 K=1,2,…,M I=1,2,…,N 重新赋边值 输出 计算结果： 求解本征值问题 物理问题 1 ：两端固定的均匀弦自由振动 设 代入上述波动方程和边界条件得 x=0 u x=l x 用 遍除各项即得 上式成立要求： 即有： 和 本征值问题 特点：1 只对一些特定的k值才能找到满足 边界条件的解 2 方程是线性的和齐次的 数值求解：打靶法 先猜一个试验本征值 对微分方程作为初值问题求解 检验所得解是否满足边界条件 若满足，则该试验本征值为真实本征值，对应的解为本征函数，否则重复1,2,3步 对弦振动问题 1. 选取k值 2. 从x=0开始求微分方程，初始条件为 3. 求出x=l时的 值，并判断是否为0 4. 重新调整k，并再度求解微分方程，直到x=l 时 的值为0，这时就找到了本征值以及对 应的本征函数 IMPLICIT NONE REAL*8 L,H,K,TOLK,DK,PHIOLD,PHI INTEGER N COMMON /PHI_PSI/PHI,H,N,K PRINT *, input string length L READ *, L