拟一维喷管流动的数值解法(MATLAB)代码(共3页).docx

拟一维喷管流动的数值解法（MATLAB）代码数值计算代码%拟一维喷管流动的数值解%亚声速-超声速，非守恒形式function main()clear;clc;r=1.4; %绝热指数N=1001; %时间步长i=31; %网格数目L=3; %喷管长度C=0.5; %柯朗数dx=L/(i-1); %空间步长dt(N)=0; %时间步长x=linspace(0,L,i); %网格点横坐标A=1+2.2*(x-1.5).2; %喷管面积%赋值M(N,i)=0;T(N,i)=0;V(N,i)=0;%初始条件M(1,:)=1-0.3146*x;T(1,:)=1-0.2314*x;V(1,:)=(0.1+1.09*x).*(1-0.2314*x).0.5;%按时间步长推进for k=1:N-1