Matlab实现斜激波受偏向角、冲击角和马赫数的影响关系
Matlab实现斜激波受偏向角、冲击角和马赫数的影响关系
TT_Matlab
博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，完整matlab代码或者程序定制加qq1575304183。
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, $ heta$ (rad)’
,
’interpreter’
,
’latex’
)
ylabel(
"Shock wave angle, $eta$ (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
,
’$ heta= heta_{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代码问题可私信交流。
部分理论引用网络文献，若有侵权联系博主删除。