曲线代码如下
syms wAngle dAngle M 
dAngle=0:5:35
k=0;
for k=1:8
dAngle(k)=dAngle(k)*(pi/180);
f(k)=2*cot(wAngle*(pi/180))*((M.^2*sin(wAngle*(pi/180)).^2-1)/((M.^2*(1.4+cos(2*wAngle*(pi/180)))+2)))-tan(dAngle(k));
ezplot(f(k),[1 4 0 90]),grid;
hold on;
end
hold off