微分方程数值解法编程作业一

微分方程数值解法编程作业一考虑对流方程初值问题00,,0,,,0,uuxRtTtxuxux其中，11,2()10,2xgxx取0.1h，书中p26页显示0.9,1.0,1.1h在n=9个步长时的计算结果。不同的取值，结果不同。上述算例表明，差分格式的稳定性不仅与差分格式本身有关，而且还与网格比的大小有关。差分格式为：111nnnjjjuuu以上格式为二步显式格式，编程时将nu和1nu分别存储。用Matlab编制程序duiliu(N,A)，输入时间步数N=9和网格比A=0.9,1.0,1.1后运行程序，画出图形。并分析下结果，谈谈自己的体会。functionduiliu(N,A)%T是计算时间的步数，A是lamdax=-2:0.01:2;M=length(x);u=zeros(1,M);fori=1:Mifx(i)(0.5+N*A*0.1)u(i)=1;%u是精确解endend[uh]=YingFeng(N,A);plot(x,u,'black',x,uh,'r:','LineWidth',2);axis([0,1.6,0,3]);xlabel('x');ylabel('u');legend('精确解','数值解');function[uh]=YingFeng(N,A)%迎风格式……………….