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

(3)
[dYdt.m]
function Yd=DYdt(t,Y)
% t
% Y
global G ME %
xy=Y(1:2);Vxy=Y(3:4); %
r=sqrt(sum(xy.^2));
Yd=[Vxy;-G*ME*xy/r^3]; %
(4)
global G ME %
G=6.672e-11;ME=5.97e24;vy0=4000;x0=-4.2e7;t0=0;tf=60*60*24*9;
tspan=[t0,tf]; %
Y0=[x0;0;0;vy0]; %
[t,YY]=ode45('DYDt',tspan,Y0);%
X=YY(:,1); %
Y=YY(:,2); %
plot(X,Y,'b','Linewidth',2); hold on
axis('image') %
[XE,YE,ZE] = sphere(10); %
RE=0.64e7; %
XE=RE*XE;YE=RE*YE;ZE=0*ZE; %
mesh(XE,YE,ZE),hold off %
8
6
4
2
0
-2
-4
-6
-8
x 10 7
0 5 10 15 20
x 10 7
图 4.14-2
【例 4.14.2.1-2】上例中,程序间的参数(如 G 和 ME)传送,是依靠全局变量形式实现的。
一般说来,编写程序时,应尽量少用全局变量,以免引起混乱。本例演示参数如何在指令间
直接传送。
(1)
[DYDt2.m]
function Yd=DYDt2(t,Y,flag,G,ME)
% flag
%
switch flag
case '' %
X=Y(1:2);V=Y(3:4);r=sqrt(sum(X.^2));Yd=[V;-G*ME*X/r^3];
35