基于Matlab模拟斜激波的 Theta-Beta-MachNo 关系附完整代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

基于Matlab模拟斜激波的 Theta-Beta-MachNo 关系附完整代码

⛄ 完整代码

close all; clear all; clc

g = 1.4;                                 % Specific heat ratio for air

beta = 0:(pi/180):(pi/2);                % Range for shock wave angle

m = 0;    

                                         % theta (=flow deflection angle)

for M1 = 1:0.01:5                        % Upstream Mach Number

    m = m+1;

    %================================

    %     theta-beta-M relation

    %================================

    Nr = ((M1^2)*((sin(beta)).^2))-1;    

    Dr = ((g+(cos(2*beta)))*M1^2)+2;

    theta = atan(2*cot(beta).*Nr./Dr);

    %================================

    %      max. theta for a M1

    %================================

    a(m) = max(theta);                   % max theta for the Mach No.

    b(m) = beta(find(theta==a(m)));      % find the beta for max. theta

    plot(theta,beta,'-b')

    hold on

end

plot(a,b,'-r','Linewidth',1.5)

xlabel('\theta')

ylabel('\beta')

axis([0 42*pi/180 0 pi/2])

⛄ 运行结果

⛄ 参考文献

[1]王晨星. 利用多介质PPM方法数值模拟斜激波与物质交界面的相互作用[D]. 中国工程物理研究院, 2004​

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

❤️ 关注我领取海量matlab电子书和数学建模资料