热传导方程的差分格式.docxVIP

一维抛物方程的初边值问题 分别用向前差分格式、向后差分格式、六点对称格式，求解下列问题: du~did du ~di d2u 0xl, u(x, 0) = sin 7Tx. 0 x 1 w(0,/) = w(l,/) = 0, / 0 在t = 0.05,0.1和0.2时刻的数值解，并与解析解u(x, t)=护sin(^)进行比较。 1差分格式形式 设空I可步长h = ( N,时I可步长r0, T = Mt,网比r = T/h2. (1)向前差分格式 该问题是第二类初边值问题(混合问题)，我们要求出所需次数的偏微商的函数 如,满足方碍M器, 0x 1,和初始条件u(x,0) = sinxr, 0x 1, 及边值条件以0,/) = 心)= 0, / 0 o 已知sin/rx在相应区域光滑，并且在x = 0J与边值相容，使问题有唯一充分光滑的 解。 取空间步长h = \IN ,和时间步长t = T/M ,其中N,M都是正整数。用两族平行 直 线 x = Xj = jh(j = 0,\,???，N) 和 t = tk=lcr(k 」,???，M) 将 矩形域 G = {0xl,z0}分割成矩形网络，网络格节点为W 以表示网格内点集合, 即位于矩形G的网点集合；Gh表示闭矩形石的网格集合；rh =Gh-Gh是网格界点的集 合。 向前差分格式，即 fi = /U,)，u(j = (pj =(p(Xj), Wq =U； = 0其中,) = 1,2,???，N_1,E = 1,2,???，M_1?以尸=M/胪表示网比。(1)式可改写成如下:R+1 fi = /U,)， u(j = (pj =(p(Xj), Wq =U； = 0 其中,) = 1,2,???，N_1,E = 1,2,???，M_1?以尸=M/胪表示网比。(1) 式可改写成如下: R+1 Wy+I = rukj+i + (1 - 2厂)iij + rukH + Tfj 此格式为显格式。 其矩阵表达式如下： l-2r r l-2r \ (j、 NJ ? Ui ■ ■ ? ■ W,V-I ? ■ “N-l 7 出丿 [) r l-2r r r l-2r （2）向后差分格式 向后差分格式，即 uk 〃k+l 。知1 I -/+1 UJ ~UJ _/屮一弘 + 勺T , r -^― 一 +力, h2 Wj = (pj =(P(Xj),Uq =Wn = °， 其中丿?=1,2,…，N —= …,M—1. （2）式可改写成 一皿;:;+ (1 + 2厂)Uj+] - = Uj + Tfj 此种差分格式被称为隐格式。 其矩阵表达式如下： l + 2r -r 、 缶田、 (j、 -r 1+2厂 ? ? ^2+， ■ ? ? -r 1 + 2厂 -r ■ ■ ■ ■ 必一】 一厂1 + 2乙 k Wiv （3）六点对称格式 六点差分格式: %— T ~ 2 喑-2咁+屹;+厂2与+屹打j h2 h2 心=Pj =(p(x),u；=说=0. 将（3）式改写成 _ 彳 W./+l +（1 + r）W./+1 - f W；-l =彳 0+1 +（1 _ 厂）UJ + 彳 0_1 + 戏j 其矩阵表达式如下： l + r -r/2 \ l-r r/2 \ (j y -r/2 l + r ? ■ ■ r/2 1-r ? ? ■ ■ -r/2 l + r -r/2 ■ ■ ? r/2 1-r r/2 ■ ■ UN-l -r 14- 2广丿 〔才丿 k r/2 1一2乙 k UN ) 2利用MATLAB求解问题的过程 对每种差分格式依次取N = 40?，el/1600, r=l/3200, r=l/6400,用 MATLAB 求解并图形比较数值解与精确解，用表格列岀不同剖分时的厶彳误差。 向前差分格式： t = 0.05: r=l/1600： x 10X 10 x 10 31 i I i I I I I i I ?数值解 精确解 2 ? ? ? - ? ? -■ ■ ? ■ ? ? ? ? ? ? ? ? 0 1 ? ? ? ■ ? ? -1 - . ? - ? ? -2 - - ? ? ? ? ? ?二； I I I I I I I I I 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 r = 1/3200： 0.7 0.6 0.5 0.4 0.3 0.2 ?数值解 精确解 0 彳 I 1 [ I II ] I ] [ 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 r = 1/6400： 0.1: r = 1/3600： I I? 数值解 I I ? 数值解 精确解 0.5 0.6 0.7 0.8 0.9 1 x 10