优秀参赛集题许fj matlab程序.pdfVIP

附录：

问题一：

(1)main

loadsatinfo.txt

satellite_06.initial=satinfo(7,:);

satellite_09.initial=satinfo(10,:);

[satellite_06.a,satellite_06.b,satellite_06.c,satellite_06.x,satellite_06.y,satellite_06.z]=satellite_RK4(satellite_06.in

itial,0.01);

[satellite_09.a,satellite_09.b,satellite_09.c,satellite_09.x,satellite_09.y,satellite_09.z]=satellite_RK4(satellite_09.in

itial,0.01);

(2)satellite_RK4

function[a,b,]=satellite_RK4(satellite,h)

GM=3.986005*10^14;

n=(250-0)/h+1;

x=zeros(1,n);

y=zeros(1,n);

z=zeros(1,n);

a=zeros(1,n);

b=zeros(1,n);

c=zeros(1,n);

x(1)=satellite(1,1);

y(1)=satellite(1,2);

z(1)=satellite(1,3);

a(1)=satellite(1,4);

b(1)=satellite(1,5);

c(1)=satellite(1,6);

fori=2:1:n

K1=(-1)*GM/(x(i-1)^2+y(i-1)^2+z(i-1)^2)^(3/2)*x(i-1);

L1=(-1)*GM/(x(i-1)^2+y(i-1)^2+z(i-1)^2)^(3/2)*y(i-1);

M1=(-1)*GM/(x(i-1)^2+y(i-1)^2+z(i-1)^2)^(3/2)*z(i-1);

N1=a(i-1);

P1=b(i-1);

Q1=c(i-1);

K2=(-1)*GM/((x(i-1)+h*N1/2)^2+(y(i-1)+h*P1/2)^2+(z(i-1)+h*Q1/2)^2)^(3/2)*(x(i-1)+h*N1/2);

L2=(-1)*GM/((x(i-1)+h*N1/2)^2+(y(i-1)+h*P1/2)^2+(z(i-1)+h*Q1/2)^2)^(3/2)*(y(i-1)+h*P1/2);

M2=(-1)*GM/((x(i-1)+h*N1/2)^2+(y(i-1)+h*P1/2)^2+(z(i-1)+h*Q1/2)^2)^(3/2)*(z(i-1)+h*Q1/2);

N2=a(i-1)+h*K1/2;

P2=b(i-1)+h*L1/2;

Q2=c(i-1)+h*M1/2;

K3=(-1)*GM/((x(i-1)+h*N2/2)^2+(y(i-1)+h*P2/2)^2+(z(i-1)+h*Q2/2)^2)^(3/2)*(x(i-1)+h*N2/2);

L3=(-1)*GM/((x(i-1)+h*N2/2)^2+(y(i-1)+h*P2/2)^2+(z(i-1)+h*Q2/2)^2)^(3/2)*(y(i-1)+h*P2/2);

M3=(-1)*GM/((x(i-1)+h*N2/2)^2+(y(i-1)+h*P2/2)^2+(z(i-1)+h*Q2/2)^2)^(3/2)*(z(i-1)+h*Q2/2);

N3=a(i-1)+h*K2/2;

P3=b(i-1)+h*L2/2;

Q3=c(i-1)+h*M2/2;

K4=(-1)*GM/((x(i-1)+h*N3)^2+(y(i-1)+h*P3)^2+(z(i-1)+h*Q3)^2)^(3/2)*(x(i-1)+h*N3);

L4=(-1)*GM/((x(i-1)+h*N3)^2+(y(i-1)+h*P3)^2+(z(i-1)+h*Q3)^2)^(3/2)*(y(i-1)+h*P3);