偏微分方程数值及matlab实验报告.docx

偏微分方程数值实验报告八 实验题目：利用有限差分法求解 u ( x) u(x) f (x), u( 1) 0, u(1) 0. 真解为 u( x) e x2 (1 x2 ) 实现算法：对于两点边值问题 d2u f , x l , dx 2 (1) u(a),u(b) , 其中 l ( a, b) (a b), f 为 l [ a,b] 上的连续函数， ， 为给定常数 . 其相应的有限差分法的算法如下： 1.对求解区域做网格剖分，得到计算网格 .在这里我们对区间 l 均匀剖分 n 段，每个剖分单元 b a 的剖分步长记为 h . n 2.对微分方程中的各阶导数进行差分离散，得到差分方程 .运用的离散方法有： 方法一 :用待定系数和泰勒展开进行离散 d 2u( xi ) i 1u( xi 1) i u( xi ) i 1u( xi 1) d( xi )2 方法二：利用差商逼近导数 d2u( xi ) u( xi 1 ) 2u( xi ) u( xi 1 ) d( xi )2 h2 将(2) 带入 (1)可以得到 u(xi 1) 2u(xi ) u(xi 1 ) ) Ri (u) , h2 f ( xi 其中 Ri (u) 为无穷小量，这时我们丢弃 Ri (u) ，则有在 xi 处满足的计算公式： u(xi 1) 2u( xi ) u( x