绘制简单的激波的图
%声速
c=340;
%物体速度
v=500;
%画的圆个数
n=3;
Ma=v/c;
miu=asin(1/Ma);
figure;
hold on
theta=0:0.01:2*pi;
for i=1:n
r=i*c;
x=v*i+r*cos(theta);
y=r*sin(theta);
plot(x,y,'b-')
hold on
end
x=0:(n+1)*v;
y1=tan(miu)*x;
y2=-tan(miu)*x;
plot(x,[y1;y2])
axis equal效果如图
