关于导引法

1）数值积分项为 R,q,xm,ym,x,y,ny,Eta(比例导引法涉及，并需要给出初值)

2）关于matlab的ode45使用方法：

[t,y,Tend,Yend,evennum] = ode45(@dif_missle_trace1,tspan,y0,op);        变步长，无法设置步长项

输入：

tspan:[t0,tf]初试时间及结束时间。

y0:初始值

dif_missle_trace1：数值积分函数。 dy=dif_missle_trace1(t,y)   y与dy的顺序及项数需要相等。

op:终止条件 op = odeset('Events',@eventfun); 其中eventfun函数如下

function [value, isterminal,direction] = eventfun(t,y)
value = y(2)-1e-5;%精度，判断停止，距离近乎为0时停运。
isterminal = 1;%为1时表示该时间会触发停止条件（value=0）
direction = 0;%判断方向：1为在零点递增，-1为递减，0任何情况都会触发。
end

输出：

t:d.   输出的时间矩阵

y: 对应时间矩阵的值矩阵

Tend:触发终止条件时的时间

Yend:终止时间的y值

evennum:触发时指针值 (好像是值isterminal)

3）swicth i

      case 1:

     case 2:

otherwise:

       end

关于现代制导规律：线性最优、自适应制导、微分对策制导  。建立原则：使导弹在接近目标的过程中付出能量或其它性能指标最小。