偏微分方程数值解教案.pdf

《偏微分方程数值解》 上机报告 实验内容1： 分别用向前差分格式、向后差分格式及六点对称格式，求解下列问题： ⎧ 2 ∂u ∂ u ⎪ = 2 + 2, 0 x1,t 0, ⎪∂t ∂x ⎨ ⎪u(0,t) = u(1,t) = 0, t 1, ⎪u(x, 0) = sin(πx) +x(1−x). ⎩ 方向 ， 方向 .在 时观察数值解与精确解 x h= 0.1 t τ = 0.01 t= 0.25 2 −π 的误差. u= e sin(πx) + x(1− x) （一）算法描述： （二）实验结果： 1.误差的数值解结果数值对比 (A)“向前差分格式”程序： forward(0.1, 0.01, 0.25) Current plot held ans = 0.0000 0.0027 0.0051 0.0070 0.0082 0.0087 0.0082 0.0070 0.0051 0.0027 0.0000 (B)“向后差分格式”程序： back(0.1, 0.01, 0.25) Current plot held ans = 0.0000 -0.0037 -0.0071 -0.0097 -0.0114 -0.0120 - 0.0114 -0.0097 -0.0071 -0.0037 0.0000 (C)“六点差分格式”程序： six(0.1, 0.01, 0.25) Current plot held ans = 0.0000 -0.0005 -0.0009 -0.0013 -0.0015 -0.0016 -0.0015 - 0.0013 -0.0009 -0.0005 0.0000 注：这里的误差=精确解-数值解. 2.精确解与数值解结果图像对比 “向前差分格式”： 注：曲线表示精确解，o表示数值解（t=0.25时）. “向后差分格式”： 注：曲线表示精确解，o表示数值解（t=0.25时）. “六点差分格式”： 注：曲线表示精确解，O表示数值解（t=0.25时）. (三)结果分析 通过（一），（二），我们检验了三种方法都能很好的求解此一维热传导方程，其 中明显能发现“六点对称格式”的误差更小。 （四）程序（附最后） 实验内容2： 用差分法求解如下自由振动问题的周期解: 2 2 ⎧∂ u ∂ u ⎪ 2 − 2 = 0, −∞ x ∞, t 0, ⎪∂t ∂x ⎨ ∂u u = 0, |t=0 = sin x, ⎪ t=0 ∂t ⎪ ⎩u(0,t) = u(2π,t). （一）算法描述： 1．网格剖分 取t∈[0, 2π], x∈[0, 2π] t t =t +ih, h = , i=0,1,..., n i 0 t t t n t x−x x = x + jh, h = o , j= 0,1,..., n j 0 x x x n x 2．差分格式 i 2 i−1 2 i−1 2 i−1 i− 2 ht u = r u + 2(1− r )u + r u −u ，r= ; j j+1 j j−1