matlab 六点差分,偏微分方程解的几道算例(差分、有限元)+含matlab程序

《偏微分方程数值解》

上机报告

实验内容1：

分别用向前差分格式、向后差分格式及六点对称格式，求解下列问题：

u 2u=+2, 00, t x2

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时观察数值解与精确解u=e πsin(πx)+x(1 x)的误差.2

(一)算法描述：

(二)实验结果：

1.误差的数值解结果数值对比

(A)“向前差分格式”程序：

>>forward(0.1,0.01,0.25)

Currentplotheld

ans=

0.00000.00270.0051

0.00820.00700.0051

(B)“向后差分格式”程序：

>>back(0.1,0.01,0.25)

Currentplotheld

ans=

0.0000-0.0037-0.0071

-0.0114-0.0097-0.0071

(C)“六点差分格式”程序：

>>six(0.1,0.01,0.25)

Currentplotheld

ans=

0.0000-0.0005-0.0009

-0.0015-0.0013-0.00090.00700.00270.00820.00000.0087-0.0097-0.0037-0.01140.0000-0.0120-0.0013-0.0005-0.00150.0000-0.0016