這幾天在做飛彈質量因數大作業和火箭主動段飛行大作業時,在計算飛彈公升力與阻力時,都遇到了計算不同飛彈高度時對應的空氣密度。
通常這種問題可直接查閱標準大氣表可以得到準確的值,但是在使用matlab進行計算時,利用查表法顯然效率太低,這是就需要乙個海拔高度與密度對應的函式來直接進行計算,但是,在剛開始做時,我在網上並沒有找到這樣的式子。
還是回去翻了翻之前的《空氣動力學》課本,找到了對應關係,記錄一下
1 rho0=1.225;t0=288.15;2if 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下就行