时间步长与空间步长对稳定性影响.docVIP

对下面的一阶对流扩散方程 ?u+ a —u= 0 x €(- x,旳,t0 ?t ?t 八 (1.1)u(x,O)=U(x) x € (- x, x (1.1) 对求解区间[°，1] 10x+1 - 0.1 x0 U (x)= - 10x + 1 0x0.1 0 其余 Matlab程序 fun cti on u = yi ngfen g(a,dt, n, mi nx,maxx,M) %方程中的常数：a %时间步长：dt %空间节点个数：n %求解区间的左端：minx %求解区间的右端：naxx %时间步的个数：M %求解区间上的数值解：u format lo ng; h = (maxx-mi nx)/( n-1); if a0 for j=1:( n+M) u0(j) = In iU(mi nx+(j-M-1)*h); end else for j=1:( n+M) u0(j) = In iU(mi nx+(j-1)*h); end end u1 = u0; for k=1:M if a0 for i=(k+1): n+M u1(i) = -dt*a*(u0(i)-u0(i-1))/h+u0(i); end else for i=1: n+M-k u1(i) = -dt*a*(u0(i+1)-u0(i))/h+u0(i); end end u0 = u1; end if a0 u