抛物型方程的有限差分方法.docVIP

抛物型方程的有限差分方法 一，求解问题 考虑一维非齐次热传导方程的定解问题 其中为正长数，,,,为已知函数，, 式（2）为初值条件，（3）为边值条件。 二，网格剖分 取空间步长和时间步长，其中M、N都是整数。用两族平行直线和将矩形域分割成矩形网格，网格结点为。以表示网格内点集合，即位于开矩形G的网点集合；表示所有位于闭矩形的网点集合；是网格界点集合。 其次，用表示定义在网点的函数，。用适当的差商代替方程（1）中相应的偏微商。 三， 差分格式 向前差分 向前差分格式 以为网比。将上式改写为便于计算的形式，则得以下向量形式 上式表示第k层的值显示表示出来。已知第k层的值，则可以直接得到第k+1的值。将上式表示成矩阵形式： 向前Euler格式的截断误差 记 然后在进行泰勒展开从而得到 即 显然得到截断误差 在矩形中某点的值。 向前Euler格式的稳定性估计 向前Euler格式，其中 ，C的特征值。 为使或必须且只须 即。故向前差分格式当时稳定，当时不稳定。 算例：应用向前Euler格式计算定解问题 上述定解问题的精确解为。 解：通过编程可求得不同步长比的数值解。表1给出了取步长h=1/10和（步长比r=1/2）时计算得到的部分数值结果，数值解很好的逼近精确解。表2给出了取步长h=1/10和（步长比r=1）时计算得到的部分数值解。随着计算层数的增加，误差越来越大，数值结果无实用价值。从图1_1取时间步长和空间步长都为10可以看到后面数值解和精确解的误差非常大，前向差分格式不稳定。出表3给出了r=1/2时，取不同步长，数值解的最大误差 从表3可以看出当空间步长缩小到原来1/2,时间步长缩小到原来的1/4时，最大误差约缩小到原来的1/4. 表1给出了取步长h=1/10和时计算得到的部分数值结果 k (x,t) 数值解 精确解 |精确解-数值解| 20 (0.5,0.1) 1.822 1.8221 0.000166 40 (0.5,0.2) 2.0134 2.0138 0.000303 60 (0.5,0.3) 2.2251 2.2255 0.00041 80 (0.5,0.4) 2.4591 2.4596 0.000484 100 (0.5,0.5) 2.7178 2.7183 0.000521 120 (0.5,0.6) 3.0036 3.0042 0.000519 140 (0.5,0.7) 3.3196 3.3201 0.000472 160 (0.5,0.8) 3.6689 3.6693 0.000375 180 (0.5,0.9) 4.055 4.0552 0.00022 200 (0.5,1.0) 1.822 4.4817 0.000166 表2给出了取步长h=1/10和时计算得到的部分数值解 k (x,t) 数值解 精确解 |精确解-数值解| 1 (0.5,0.01) 1.6652 1.6653 6.90E-05 2 (0.5,0.02) 1.6819 1.6820 0.000139 3 (0.5,0.03) 1.6987 1.6989 0.000211 4 (0.5,0.04) 1.7157 1.7160 0.000284 5 (0.5,0.05) 1.7329 1.7333 0.000359 6 (0.5,0.06) 1.7504 1.7507 0.000279 7 (0.5,0.07) 1.7673 1.7683 0.000976 8 (0.5,0.08) 1.7875 1.7860 0.00142 9 (0.5,0.09) 1.7970 1.8040 0.007015 10 (0.5,0.10) 1.8423 1.8221