微分方程数值解法编程作业一(共1页).docx

精选优质文档-倾情为你奉上微分方程数值解法编程作业一考虑对流方程初值问题 其中，取，书中p26页显示在n=9个步长时的计算结果。不同的取值，结果不同。上述算例表明，差分格式的稳定性不仅与差分格式本身有关，而且还与网格比的大小有关。差分格式为：以上格式为二步显式格式，编程时将和分别存储。用Matlab编制程序duiliu(N,A)，输入时间步数N =9和网格比A=0.9,1.0,1.1后运行程序，画出图形。并分析下结果，谈谈自己的体会。function duiliu(N,A) %T是计算时间的步数，A是lamdax=-2:0.01:2;M=length(x);u=zeros(1,M);for i=1:M if x(i)(0.5+N*A*0.1) u(i)=1; %u是精确解 endenduh=YingFeng(N,A);plot(x,u,black,x,uh,r:,LineWidth,2);axis(0,1.6,