一维激波管数值模拟－修改

p1(1,i) + temp2(1,i));Р Р % 密度、速度乘积计算Р temp1(1,i) = (f2_Pos(1,i) - f2_Pos(1,i-1)) / delta_x;Р temp2(1,i) = (f2_Neg(1,i+1) - f2_Neg(1,i)) /delta_x;Р Den_u(1,i) = Den_u(1,i) - delta_t*(temp2(1,i) + temp1(1,i));Р Р % 速度计算Р u(1,i) = Den_u(1,i) / Den(1,i);Р Р % 能量计算Р temp1(1,i) = (f3_Pos(1,i) - f3_Pos(1,i-1)) /delta_x;Р temp2(1,i) = (f3_Neg(1,i+1) - f3_Neg(1,i)) /delta_x;Р E(1,i) = E(1,i) - delta_t*(temp2(1,i) + temp1(1,i));Р % 压强计算Р P(1,i) = (gama - 1)*(E(1,i) - 0.5*Den(1,i)*u(1,i)^2);Р endРendР%结果显示Рx=0:0.001:1;axis([0 1 0 1]);Рplot(x,u,'LineWidth',2);xlabel('\fontsize{14}x');ylabel('\fontsize{14}速度V');Рfigure(2);plot(x,P,'LineWidth',2);xlabel('\fontsize{14}x');ylabel('\fontsize{14}压力P')Рfigure(3);plot(x,Den,'LineWidth',2);xlabel('\fontsize{14}x');ylabel('\fontsize{14}密度ρ')РCalculate_time=toc