rightf.m

这是比例导引的MATLAB程序,其中龙格库塔积分的程序很简单且通用性强

%右函数定义的，用于单步龙格库塔法计算。
%如果需要使用全局变量，要使用global定义。
% 这种格式的，也完全可以被ode45函数调用。

function dy = rightF(t,y) % right functions
global XNC VM;

dy = zeros(5,1);    % a column vector

dy(1) = XNC/VM;
dy(2) = -XNC*sin(y(1)); %过载垂直于theta（弹道倾角）
dy(3) =1*XNC*cos(y(1));

dy(4) =y(2);
dy(5) =y(3);

% 下面的，为方便调用程序，把stepRK.m的源代码放在这

%函数说明，t表示自变量（一般为时间），y表示变量的初值（是一个行向量），h表示步长
% %运行的结果，也是一个行向量
% 
% function dy = stepRK(t,y,h)
% K1=rightF(t,y);
% K2=rightF(t+h/2,y'+h*K1/2);
% K3=rightF(t+h/2,y'+h*K2/2);
% K4=rightF(t+h,y'+h*K3);
% dy=y'+(K1+2*K2+2*K3+K4)*h/6;
% dy=dy';
%