倒向微分方程.docVIP

微分方程数值解 ———倒向微分方程 计算终边值问题： 利用显格式计算： 原式可化为： 为了简化计算，取： 由： 且由显格式稳定的条件： 化简得： 时该方法是稳定的 将带入计算，求得解如下: 编程代码： %赋初值 e=0.0001; h=0.02; t=0.001:e:0.1;%k s=0.2:h:2;%j tl=length(t);sl=length(s); %初始矩阵 F(1,1:tl)=4;%边值 F(sl,1:tl)=0;%边值 F(1:sl,tl)=max(4-s,0);%初值 %构造C R=1/2*(s).^2*0.25; S=(s)*0.025; a=R(2:sl-2)+0.5*S(2:sl-2); b=1-2*R(2:sl-1)-0.0025; c=R(2:sl-2)-0.5*S(2:sl-2); C=diag(a,1)+diag(b)+diag(c,-1); %构造H q=zeros(sl-4,tl); H=[(R(1)-0.5*S(1))*F(1,:);q;(R(sl)+0.5*S(sl))*F(sl,:)]; %算法 ： for n=tl:-1:2 F(2:sl-1,n-1)=C*F(2:sl-1,n)+H(:,n); end %画图 [T,S]=meshgrid(t,s(1:sl)); mesh(T,S,F) xlabel(T)