function?ydot=bili（ty）
%?y（1）=theta??y（2）=V??y（3）=x??y（4）=y?y（5）=xm??y（6）=ym??y（7）=m
T=288.15-0.0065*y（4）;%溫度與高度的關(guān)系，用于算當(dāng)?shù)匾羲?br/>gamma=1.4;
r=287;
a=sqrt（T*r*gamma）;%當(dāng)?shù)匾羲?br/>Re=6370000;%地球半徑
g=9.8*（Re/（Re+y（4）））^2;%當(dāng)?shù)刂亓铀俣?br/>rho=1.225*exp（-0.00015*y（4））;%當(dāng)?shù)乜諝饷芏?br/>s=1.94;???????%機(jī)翼面積
M=y（2）/a;%馬赫數(shù)
Vm=-300;%目標(biāo)速度
K=3;%比例系數(shù)
e=180/pi;
q=atan（（y（6）-y（4））/（y（5）-y（3）））;
eta=q*e-y（1）;
R=（y（6）-y（4））/sin（q）;
%導(dǎo)彈推力，質(zhì)量與時(shí)間的關(guān)系
t1=2.79;t2=t1+40.82;
dotm=219.45*（t<=t1）+7.42*（t>t1&tp=40159*g*（t<=t1）+5135.94*g*（t>t1&t%升阻力系數(shù)的擬合
Mi=[1.52.02.53.03.54.04.5];?%給出的馬赫數(shù)據(jù)
Clai=[0.03370.03240.0300.0240.0250.02330.0244];?%給出的數(shù)據(jù)?
Cla=interp1（MiClaiM‘spline‘）;%線性插值求升力線斜率
alpha=（K*y（7）*y（2）*（y（2）*sin（eta/e）-Vm*sin（q））/R+y（7）*g*cos（y（1）/e））/（0.5*rho*y（2）^2*s*Cla+p/e）;%求攻角
alphai=0:2:8;????????????%給出的數(shù)據(jù)
Mi2=[1.52.02.53.03.54.0];???%給出的數(shù)據(jù)
Cxi=[0.0430.0460.0540.0670.084;0.0370.03950.0470.05850.0745;0.03170.03350.0400.05050.0645;...
????0.0270.0290.0340.04350.056;0.0230.0250.02950.0370.0485;0.02050.02250.0260.03250.042];
Cx=interp2（alphaiMi2CxialphaM‘spline‘）;?%線性插值求阻力系數(shù)
%微分方程組
%?y（1）=theta??y（2）=V??y（3）=x??y（4）=y?y（5）=xm??y（6）=ym??y（7）=m
ydot=[（p*alpha/57.3+0.5*rho*y（2）^2*s*Cla*alpha-y（7）*g*cos（y（1）/e））/（y（7）*y（2））;...
??????p*cos（alpha/e）/y（7）-0.5*rho*y（2）*y（2）*Cx*s/y（7）-g*sin（y（1）/e）;...
??????y（2）*cos（y（1）/e）;...
??????y（2）*sin（y（1）/e）;...
??????Vm;...
??????0;...
??????-dotm];

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1604??2018-11-30?21:01??比例導(dǎo)引法1\bili.m

?????文件????????181??2018-11-30?21:02??比例導(dǎo)引法1\bilievents.m

?????文件????????438??2018-11-30?21:02??比例導(dǎo)引法1\zhu.m

?????文件???????2426??2018-12-02?22:07??比例導(dǎo)引法1\三個(gè)全部集合在這個(gè)文件運(yùn)行時(shí)請點(diǎn)擊zhu其他為函數(shù)文件.m

?????目錄??????????0??2018-12-02?22:08??比例導(dǎo)引法1

-----------?---------??----------?-----??----

?????????????????4649????????????????????5