8~数值分析.pptVIP

第8章 抛物型方程的差分格式 模型 长度为 L 的绝缘杆上的一维热流 热传导方程在时间 t 和位置 x 处的温度 u (x , t) 表示为 初始温度分布为 杆端点的边界值为 其中 k 是导热率系数，σ是热量，ρ是杆的密度。 8.1 差分格式建立的基础 一、网格 将求解区域Ω分割成 M?N个小矩形，长宽分 别为 ⊿x=h ⊿t=k 二、差商 一阶偏导数 的向前、向后、中心差商为 二阶偏导数 的中心差商为 用(xm , tn)处的一阶向前差商近似代替微商ux， 即 考虑 u(x , t) 的 Taylor 展式 截断误差 用(xm , tn)处的一阶中心差商近似代替微商ux，即 截断误差为 三、算子 为 x 方向偏导数算子 Tx为 x 方向位移算子 μx 为 x 方向平均算子 x 方向差分算子 前差算子 后差算子 中心差算子 （I为恒等算子） 8.2 显式差分格式 差分方法：用差商代替微商，在网格节点上 求出微分方程解的近似值的一种方法。 一、一维常系数热传导方程的古典显式格式 考虑一维热传导方程 (2.26) 在网格节点 的差分近似 取 r = k / h2 为步长比，得显式向前差分方程 (2.29) (2.29)也称为解热传导方程(2.26)的古典显式格 式。截断误差为 例1 用古典显式格式求解抛物型方程 初始条件为 边界条件为 取步长⊿x = h = 0.2 , ⊿t = k = 0.02 。 解 r = k / h2 = 0.02 / 0.22 = 0.5, 古典显式格式为 x=0.0 x=0.2 x=0.4 x=0.6 x=0.8 x=1.0 U = t=0.00 0 0.6400 0.9600 0.9600 0.6400 0 t=0.02 0 0.4800 0.8000 0.8000 0.4800 0 t=0.04 0 0.4000 0.6400 0.6400 0.4000 0 t=0.06 0 0.3200 0.5200 0.5200 0.3200 0 t=0.08 0 0.2600 0.4200 0.4200 0.2600 0 t=0.10 0 0.2100 0.3400 0.3400 0.2100 0 t=0.12 0 0.1700 0.2750 0.2750 0.1700 0 t=0.14 0 0.1375 0.2225 0.2225 0.1375 0 t=0.16 0 0.1113 0.1800 0.1800 0.1113 0 t=0.18 0 0.0900 0.1456 0.1456 0.0900 0 t=0.20 0 0.0728 0.1178 0.1178 0.0728 0 function u=gu_dian(f,a,b,c1,c2,m,n) %输入初值和U h=a/(m-1); k=b/(n-1); r=k/h^2; U=zeros(n,m); %赋边界条件 U(2:n,1)=c1; U(2:n,m)=c2; %赋初始条件 U(1,1:m)=fg(0:h:h*(m-1)); %计算内点