有如图所示方程组,以下是代码形式:
syms ma1 ma2 b
r = 1.4;
d = 0:5:35;
eq1 = tand(d) - (ma1^2*sind(b)^2-1)/((1+ma1^2*((r+1)/2-sind(b)^2))*tand(b));
eq2 = ma2^2 - (1+(r-1)/2*ma1^2*sind(b)^2)/((r*ma1^2*sind(b)^2-(r-1)/2)*sind(b-d)^2);
其中r为常数1.4,d为常数(d为从0开始以5为步长到35的数),两个方程消去b后绘制Ma2-Ma1的图像,结果图如图。
已经尝试过单独解出b再画但是不成功,有老师建议用等高线的方法,但是我不太清楚怎么确定那个矩阵Z
请各位老师指教!谢谢!
|