内容简介作者简介 - 科学与工程计算系

otherwise
error(['Unknown flag ''' flag '''.']);
end
(2)
[t,YY]=ode45('DYDt2',tspan,Y0,[],G,ME);
4.14.2.2 解算指令较复杂格式的使用示例
【例 4.14.2.2-1】带事件设置的 ODE 文件及主程序编写演示。本例将以较高精度计算卫星经
过近地点和远地点的时间,并在图上标志。
(1)
[DYDt3.m]
function varargout=DYDt3(t,Y,flag,G,ME,tspan,Y0)
%
%
%
%
%
%
%
switch flag
case '' %
varargout{1} = f(t,Y,G,ME); %
case 'init' %
[varargout{1:3}] = fi(tspan,Y0); %
case 'events' %
[varargout{1:3}] = fev(t,Y,Y0);
otherwise
error(['Unknown flag ''' flag '''.']);
end
% ------------------------------------------------------------------
function Yd = f(t,Y,G,ME)
%
X=Y(1:2);V=Y(3:4);r=sqrt(sum(X.^2));Yd=[V; -G*ME*X/r^3];
% ------------------------------------------------------------------
function [ts,y0,options] = fi(tspan,Y0)
%
ts=tspan;y0 = Y0;
options = odeset('Events','on','Reltol',1e-5,'Abstol',1e-4);
%
% ------------------------------------------------------------------
function [value,isterminal,direction] = fev(t,Y,Y0)
%
dDSQdt = 2 * ((Y(1:2)-Y0(1:2))' * Y(3:4));
%dDSQdt
value = [dDSQdt; dDSQdt]; %
direction = [1; -1]; %
isterminal = [1; 0]; %
36