带落角约束的变结构末制导律（线性）源代码

最新推荐文章于 2017-03-07 19:34:29 发布

51星系

最新推荐文章于 2017-03-07 19:34:29 发布

阅读量1.7k

点赞数 2

CC 4.0 BY-SA版权

分类专栏：科研航空航天仿真算法文章标签： matlab 源代码

本文链接：https://blog.youkuaiyun.com/stk10/article/details/52814902

本文介绍了带落角约束的变结构末制导律设计，并提供了线性化的源代码实现。参考了周荻的著作，探讨了变结构导引律在实际应用中的关键步骤和数学模型。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

参考文献：周荻那本变结构导引律的设计

clear

%===================================================================================
%--------------------------------情景说明------------------------------------------
%目标正弦机动,At=Asin(*),速度大小不变
%导弹速度大小为常量
%==========================================================================
Vt=500;
Xt=10000;
Yt=8000;
ThetaT=pi/9;
Xt1=[]; %存放中间数据，用来分析结果，如作图
Yt1=[];
At1=[];

Vm=2000;
Xm=0;
Ym=0;
ThetaM=2*pi/9;
Xm1=[]; %存放中间数据，用来分析结果，如作图
Ym1=[];
Am1=[];
AngleV=[];
DotQ1=[];
Q1=[];

n=1;
Xt1(n)=Xt;
Yt1(n)=Yt;

Xm1(n)=Xm;
Ym1(n)=Ym;

Q=atan((Yt-Ym)/(Xt-Xm));
DotQ=((Vt*sin(ThetaT)-Vm*sin(ThetaM))*(Xt-Xm)-(Yt-Ym)*(Vt*cos(ThetaT)-Vm*cos(ThetaM)))/((Xt-Xm)^2+(Yt-Ym)^2);
DotR=((Xm-Xt)*(Vm*cos(ThetaM)-Vt*cos(ThetaT))+(Ym-Yt)*(Vm*sin(ThetaM)-Vt*sin(ThetaT)))/sqrt((Xm-Xt)^2+(Ym-Yt)^2);

t=0;
Dt=0.01;
k=10;
E1=40;
E2=0.01;
c=2;
Qd=-pi/2;
At=20;
while(DotR<=0)
        DotQ1(n)=DotQ;
        Q1(n)=Q;
        
        %if(abs(t-3)<=1e