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 onendplot(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');    

3 仿真结果

4 参考文献

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

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

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