基于滑模控制的含有攻击时间和攻击角度约束的制导律代码

基于滑模控制（sliding model control）的含有攻击时间和攻击角度约束的制导律仿真代码，含有弹道图、攻击角度图、攻角图等等，仿真环境：MATLAB

function [ty] = RK4（ftbegantendy1h）
n=（tend-tbegan）/h; 
y（:1）=y1;
t（1）=tbegan;
global V;       %导弹速度
global qd;      %λd 撞击角度
global td;      %撞击时间
global flag;
flag=0;
i=1;
while（i<=n&&t（i）    K1=h*f（t（i）y（1:4i））;
    K2=h*f（t（i）+h/2y（1:4i）+K1/2）;
    K3=h*f（t（i）+h/2y（1:4i）+K2/2）;
    K4=h*f（t（i）+hy（1:4i）+K3）;
    t（i+1）=t（i）+h;
    y（1:4i+1）=y（1:4i）+（K1+2*K2+2*K3+K4）/6;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    y（5i+1）=td-t（i）-y（1i）;
    if（y（5i+1）<0.0000001）
       flag=1;
    end
    r=y（1i）*V;
    lamuda=y（2i）+qd;
    y（6i+1）=-r*cos（lamuda）;
    y（7i+1）=-r*sin（lamuda）;
    y（8i+1）=lamuda;
    y（9i+1）=V*（（y（3i+1）-y（3i））/h-V*sin（y（3i））/r）;
    
    i=i+1;   
end

end

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

     文件        798  2018-12-02 11:03  code\RK4.m

     文件        674  2018-12-02 11:00  code\solvefunc.m

     文件       1262  2018-12-02 11:03  code\solver.m

     文件      51128  2018-12-02 10:37  code\td=55s，λd改变，a变化曲线.fig

     文件      49072  2018-12-02 10:37  code\td=55s，λd改变，λ变化曲线.fig

     文件      50347  2018-12-02 10:37  code\td=55s，λd改变，σ变化曲线.fig

     文件      58204  2018-12-02 10:29  code\td=55s，λd改变，弹道图.fig

     文件      55334  2018-12-02 10:56  code\λd=-120°，td改变，a变化曲线.fig

     文件      55292  2018-12-02 10:53  code\λd=-120°，td改变，λ变化曲线.fig

     文件      54307  2018-12-02 10:54  code\λd=-120°，td改变，σ变化曲线.fig

     文件      64458  2018-12-02 10:52  code\λd=-120°，td改变，弹道图.fig

     目录          0  2019-09-02 11:24  code

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

               440876                    12