如何用Matlab求题主给出问题的数值解?
第一步,将已知值分别赋值给变量,即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步,使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步,使用for双循环语句,得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步,显示A/Acr、 k、 Ma值,即
disp(A/Acr)
disp(k)
disp(Ma)