一维对流方程编程求解_一维对流方程编程求解_源码

function z=OneD_convection_Center(a,b,r,h,T) %空间中心差分 %参数设置 %区间左端点a=0 区间右端点b=1.2 时间步长r=0.01 空间步长h=0.01 时间迭代步数T=50 c=1; Alpha=r/h; N=round((b-a)/h); x=linspace(a,b,N+1); %t=0时，初始化 u(1,:)=Initialization(a,h,N); plot(x,u(1,:),'linewidth',1); set(0,'defaultfigurecolor','w'); pause(0.1); %开始迭代 for k=1:(T-1) u(k+1,1)=(1+c*Alpha)*u(k,1)-c*Alpha*u(k,2);%边界点处理 for n=2:N u(k+1,n)=(1+c*Alpha)*u(k,n)-c*Alpha*u(k,n+1); end u(k+1,N+1)=(1-c*Alpha)*u(k,N+1)+c*Alpha*u(k,N);%边界点处理 plot(x,u(k+1,:)); pause(0.1); end