比例导引律Matlab程序,是否有关于比例导引方面的仿真程序？

本人想学习这方面的知识，上网看了一个类似代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%    PNGE.m

%    Proportional Navigation Guidance Emluator

%    导弹拦截的比例制导仿真程序

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [IT,MD,PT]=PNGE(MT,MTSD,ILP,ILD,IV,IAL,S)

%用比例导引律仿真拦截器拦截来袭导弹的时间步长模拟算法

%  ChengAihua,PLA Information Engineering University,ZhengZhou,China

%  Email:aihuacheng@gmail.com

%  All rights reserved

%比例导引律：拦截器速度向量的偏转角速度(法向过载)与目标视线的旋转角速度成正比

%输入参数列表

% MT    Missile Track                    来袭弹道导弹运动轨迹(N×3矩阵)

% MTSD  Missile Track's Sample Disdance  导弹轨迹采样时间间隔，兼仿真时间步长

% ILP   Interceptor's Launch Position    拦截器发射地点(1×3矩阵)

% ILD   Interceptor's Launch Direction   拦截器发射方向归一化向量(1×3矩阵)

% IV    Interceptor's Veloci