偏微分方程解的几道算例（差分、有限元）-含matlab程序.pdfVIP

《偏微分方程数值解》 上机报告 实验内容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.