在进行飞行力学建模时飞机的气动导数一般由拟合公式或插值表给出,拟合公式使用起来比较简单,插值表需要用插值方法计算,气动导数通常与两个变量有关,是二元插值问题,下面是Matlab计算气动导数一例:
function [ Cndr_0 ] = Cndr( alpha_0,beta_0 )
% 输入单位为度
% 计算指定攻角和侧滑角下的Cndr
alpha=[-10 -5 0 5 10 15 20 25 30 35 40 45]';
beta=[-30 -20 -10 0 10 20 30]';
Cndr1=[-0.018 -0.028 -0.037 -0.048 -0.043 -0.052 -0.062
-0.052 -0.051 -0.041 -0.045 -0.044 -0.034 -0.034
-0.052 -0.043 -0.038 -0.045 -0.041 -0.036 -0.027
-0.052 -0.046 -0.04 -0.045 -0.041 -0.036 -0.028
-0.054 -0.045 -0.04 -0.044 -0.04 -0.035 -0.027
-0.049 -0.049 -0.038 -0.045 -0.038 -0.028 -0.027
-0.059 -0.057 -0.037 -0.047 -0.034 -0.024 -0.023
-0.051 -0.052 -0.03 -0.048 -0.035 -0.023 -0.023
-0.03 -0.03 -0.027 -0.049 -0.035 -0.02 -0.019
-0.037 -0.033 -0.024 -0.045 -0.029 -0.016 -0.009
-0.026 -0.03 -0.019 -0.033 -0.022 -0.01 -0.025
-0.013 -0.008 -0.013 -0.016 -0.009 -0.014 -0.01]';
Cndr_0=interp2(alpha,beta,Cndr1,alpha_0,beta_0);
end
注:插值最大的问题是要注意插值范围,外插不可靠,有些角度只给出正的部分,负值部分对称过去即可(如β)。