這幾天在做導彈質量因數大作業和火箭主動段飛行大作業時，在計算導彈升力與阻力時，都遇到了計算不同導彈高度時對應的空氣密度。
通常這種問題可直接查閱標准大氣表可以得到准確的值，但是在使用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下就行