Matlab实现斜激波受偏向角、冲击角和马赫数的影响关系

1 简介

Matlab实现斜激波受偏向角、冲击角和马赫数的影响关系

2 完整代码

登录后复制

          
          

           
           close all; clear; clc;
           
           
g=1.4;
           
           
beta=0:0.01:(pi/2);
           
           
m=0;
           
           
for M1=[1:0.1:1.4,1.6,1.8,2,2.2,2.5,3,4]
           
           
    m=m+1;
           
           
    theta=atan(2*cot(beta).*(((M1^2)*((sin(beta)).^2))-1) ...
           
           
        ./(((g+(cos(2*beta)))*M1^2)+2));
           
           
    a(m)=max(theta);
           
           
    b(m)=beta(theta==a(m));
           
           
    plot(theta,beta)
           
           
    hold on
           
           
end
           
           
plot(a,b,'-r')
           
           
xlabel('Deflection angle, $\theta$ (rad)','interpreter','latex')
           
           
ylabel("Shock wave angle, $\beta$ (rad)",'interpreter','latex')
           
           
axis([0 0.7 0 pi/2])
           
           
legend('$M_1=1$','$M_1=1.1$','$M_1=1.2$','$M_1=1.3$','$M_1=1.4$',...
           
           
    '$M_1=1.6$','$M_1=1.8$','$M_1=2$','$M_1=2.2$','$M_1=2.5$','$M_1=3$',...
           
           
    '$M_1=4$','interpreter','latex','location','southeast')
           
           
text(0.3,0.6,'weak wave','BackgroundColor','white','interpreter','latex')
           
           
text(0.4,1.45,'strong wave','BackgroundColor','white','interpreter','latex')
           
           
text(0.55,1.21,'$\theta=\theta_{max}$','color','red','BackgroundColor',...
           
           
    'white','interpreter','latex')
           
           
text(0.35,1,'$M_2>1$','BackgroundColor','white','interpreter','latex')
           
           
text(0.35,1.25,'$M_2<1$','BackgroundColor','white','interpreter','latex')
           
           
% print('-dpdf', 'obliqueSW.pdf')
           
           
% system('pdfcrop obliqueSW.pdf obliqueSW.pdf');
          
          

• 1.
• 2.
• 3.
• 4.
• 5.
• 6.
• 7.
• 8.
• 9.
• 10.
• 11.
• 12.
• 13.
• 14.
• 15.
• 16.
• 17.
• 18.
• 19.
• 20.
• 21.
• 22.
• 23.
• 24.
• 25.
• 26.
• 27.
• 28.

3 仿真结果

4 参考文献

[1]王宇, 王世芬. 马赫数对后掠激波和湍流边界层干扰特性的影响[J]. 空气动力学学报, 1994, 012(003):313-319.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。