马赫数和速度相互转换函数（matlab插值计算）

function [Mach] = v2Mach(v,h)
% 速度与马赫数转换函数，插值计算
% 输入:速度v 单位 m/s ，高度h 单位m  输出：Mach
H=[0 100 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 ...
    8500 9000 9500 10000 11000 12000 13000 14000 15000 16000 17000 18000 19000 20000 25000 30000 ...
    35000 40000 45000 50000 55000 60000 65000 70000 75000 80000]';
a=[340.3 340.0 338.4 336.4 334.5 332.5 330.6 328.6 326.6 324.6 322.6 320.5 318.5 316.5 314.4 312.3 310.2 308.1 ...
    306.0 303.8 301.7 299.5 295.1*ones(1,10) 298.4 301.7 308.3 317.2 325.3 329.8 326.7 320.6 310.1 297.1 283.6 269.4]';
Vsound=interp1(H,a,h,'linear');
Mach=v/Vsound;
end


function [V] = Mach2V(Ma,h)
% 速度与马赫数转换函数，插值计算
% 输入:Mach ，高度h 单位m  输出：V 速度 单位 m/s
H=[0 100 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 ...
    8500 9000 9500 10000 11000 12000 13000 14000 15000 16000 17000 18000 19000 20000 25000 30000 ...
    35000 40000 45000 50000 55000 60000 65000 70000 75000 80000]';
a=[340.3 340.0 338.4 336.4 334.5 332.5 330.6 328.6 326.6 324.6 322.6 320.5 318.5 316.5 314.4 312.3 310.2 308.1 ...
    306.0 303.8 301.7 299.5 295.1*ones(1,10) 298.4 301.7 308.3 317.2 325.3 329.8 326.7 320.6 310.1 297.1 283.6 269.4]';
Vsound=interp1(H,a,h,'linear');
V=Ma*Vsound;
end