winged-cone的六自由度部分开环模型

winged-cone的六自由度部分开环模型仿真，里面有一些物理参数可供参考

function out = cxfiv（in）
inputs    =   in（1:4）;
RE        =   inputs（1）;
LE        =   inputs（2）;
RUD       =   inputs（3）;
PLA       =   inputs（4）;

states    =   in（5:16）;
x         =   states（1）;
y         =   states（2）;
h         =   states（3）;
V         =   states（4）;
chi       =   states（5）;
gamma     =   states（6）;
alpha     =   states（7）;
beta      =   states（8）;
mu        =   states（9）;
p         =   states（10）;
q         =   states（11）;
r         =   states（12）;
%******************************************************
%飞行器相关的参数
%******************************************************
b         =   18.288*3.2808;
Jyy       =   9036000/14.594*3.2808^2;
Jzz       =   9036000/14.594*3.2808^2;
Jxx       =   915300/14.594*3.2808^2;
Mass      =   136000/14.594;
S         =   334.73*3.2808^2;
c         =   80;
g         =   9.8*3.2808;
%******************************************************
%纵向模型中的相关参数
%******************************************************
M         =   V/（（8.99e-9）*h^2-（9.16e-4）*h+996）;
rho       =   0.00238*exp（-h/24000）;
qb        =   rho*V^2/2;

ALPHA     =   alpha*180/pi;

beta      =   0;

%x_cg      =   100000;
%x_cg      =   4.4668*3.2808;
%******************************************************
%气动系数、气动力和气动力矩
%******************************************************
Clb       =   -1.402E-01+3.326E-02*M-7.590E-04*ALPHA+8.596E-06*（ALPHA*M）-3.794E-03*M^2+2.354E-06*ALPHA^2-1.044E-08*（ALPHA*M）^2+2.219E-04*M^3-8.964E-18*ALPHA^3-6.462E-06*M^4+3.803E-19*ALPHA^4+7.419E-08*M^5-3.353E-21*ALPHA^5;
ClRE      =   3.570E-04-9.569E-05*ALPHA-3.598E-05*M+1.170E-04*RE+2.794E-08*（ALPHA*M）*RE+4.950E-06*ALPHA^2+1.411E-06*M^2-1.160E-06*RE^2-4.641E-11*（（ALPHA*M）*RE）^2;
ClLE      =   -（3.570E-04-9.569E-05*ALPHA-3.598E-05*M+1.170E-04*LE+2.794E-08*（ALPHA*M）*LE+4.950E-06*ALPHA^2+1.411E-06*M^2-1.160E-06*LE^2-4.641E-11*（（ALPHA*M）*LE）^2）;
ClRUD     =   -5.0103E-19+6.2723E-20*ALPHA+2.3418E-20*M+0.00011441*RUD-2.6824E-06*（ALPHA*RUD）-3.4201E-21*（ALPHA*M）-3.5496E-06*（M*RUD）+5.5547E-08*（ALPHA*M）*RUD;
Clr       =   3.82E-01-1.06E-01*M+1.94E-03*ALPHA-8.15E-05*（ALPHA*M）+1.45E-02*M^2-9.76E-06*ALPHA^2+4.49E-08*（ALPHA*M）^2+1.02E-03*M^3-2.70E-07*ALPHA^3+3.56E-05*M^4+3.19E-08*ALPHA^4-4.81E-07*M^5-1.06E-09*ALPHA^5;
Clp       =   -2.99E-01+7.47E-02*M+1.38E-03*ALPHA-8.78E-05*（ALPHA*M）-9.13E-03*M^2-2.04E-04*ALPHA^2-1.52E-07*（ALPHA*M）^2+5.73E-04*M^3-3.86E-05*ALPHA^3-1.79E-05*M^4+4.21E-06*ALPHA^4+2.20E-07*M^5-1.15E-07*ALPHA^5;
Cl        =   Clb+ClRE+ClLE+ClRUD+Clr*r*b/2/V+Clp*p*b/2/V;
l_bar     =   qb*b*S*Cl;%滚转力矩系数

Cma       =   -2.192E-02+7.739E-03*M-2.260E-03*ALPHA+1.808E-04*（ALPHA*M）-8.849E-04*M^2+2.616E-04*ALPHA^2-2.880E-07*（ALPHA*M）^2+4.617E-05*M^3-7.887E-05*ALPHA^3-1.143E-06*M^4+8.288E-06*ALPHA^4+1.082E-08*M^5-2.789E-07*ALPHA^5;
CmRE      =   -5.67E-05-6.59E-05*ALPHA-1.51E-06*M+2.89E-04*RE+4.48E-06*（ALPHA*RE）-4.46E-06*（ALPHA*M）-5.87E-06*（M*RE）+9.72E-08*（ALPHA*M）*RE;
CmLE      =   -5.67E-05-6.59E-05*ALPHA-1.5

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8705  2011-07-01 20:39  开环特性\cxfiv.m

     文件      39015  2011-07-01 20:43  开环特性\cxo.mdl

     文件      40691  2011-07-01 20:41  开环特性\cxo1.mdl

     文件       9880  2011-07-01 20:41  开环特性\cxseven.m

     目录          0  2011-08-23 17:29  开环特性

----------- ---------  ---------- -----  ----

                98291                    5