【2019年整理】一维抛物线偏微分方程数值解法(4)(附图及matlab程序).docxVIP

一维抛物线偏微分方程数值解法(4) 上一篇参看 一维抛物线偏微分方程数值解法( 3)(附图及matlab程序) 解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法) Ut-Uxx=0, 0x1,0t=1 (Ut-aUxx=f(x,t),a0) U(x,0)=eAx, 0=x=1, U(0,t)=eAt,U(1,t)=eA(1+t), 0t=1 精确解为：U(x,t)=eA(x+t); 用紧差分格式： 此种方法精度为o(h1A2+h2A4), 无条件差分稳定； 一：用追赶法解线性方程组(还可以用迭代法解) Matlab程序为： function [u p e x t]=JCHGS(h1,h2,m,n) %紧差分格式解一维抛物线型偏微分方程 %t匕程序用的是追赶法解线性方程组 %h1为空间步长，h2为时间步长 %m,n分别为空间，时间网格数 %p为精确解，u为数值解,e为误差 x=(0:m)*h1+0; x0=(0:m)*h1; %定义x0,t0 是为了 f(x,t)~=0 的情况 % t=(0:n)*h2+0; t0=(0:n)*h2+1/2*h2; syms f; for (i=1:n+1) for (j=1:m+1) f(i,j)=0; %f(i,j)=f(x0(j),t0(i))==0% end end for (i=1:n+1) u(i,1)=exp(t(i));