这几天在做导弹质量因数大作业和火箭主动段飞行大作业时，在计算导弹升力与阻力时，都遇到了计算不同导弹高度时对应的空气密度。
通常这种问题可直接查阅标准大气表可以得到准确的值，但是在使用matlab进行计算时，利用查表法显然效率太低，这是就需要一个海拔高度与密度对应的函数来直接进行计算，但是，在刚开始做时，我在网上并没有找到这样的式子。
还是回去翻了翻之前的《空气动力学》课本，找到了对应关系，记录一下
1 rho0=1.225;T0=288.15; 2 if h<=11000 3 T=T0-0.0065*h; 4 Rho=rho0*(T/T0)^4.25588; 5 elseif h>11000&&h<=20000 6 T=216.65; 7 Rho=0.36392*exp((-h+11000)/6341.62); 8 else 9 T=216.65+0.001*(h-20000); 10 Rho=0.088035*(T/216.65)^-35.1632; 11 end
其中，
高度0-11000,对流层
11000-20000m 平流层较低位置
20000-32000m 平流层较高位置
另外，也可直接使用查表法：查表数据 数据是.mat格式的，可用matlab直接打开，用的时候load下就行