abc.rar_matlab变轨_卫星轨道计算_变轨_地月轨道_地球月球matlab

close all;clear all; ap=0; %相机视角采用自动设置 load topo %装载图象 topo 和相应的图象色图 topomap1 figure('colormap',topomap1,'Color',[.0 .0 .0]); %设置图形窗的色图和背景色 [x,y,z] = sphere(50); %绘制50格网状球面 azzy.DataAspectRatio = [1 1 1];azzy.PlotBoxAspectRatioMode = 'auto'; fa = axes('Visible','off', azzy); %设置轴的数据宽高比和坐标框三度比 szzy1.AmbientStrength = 0.1;szzy1.DiffuseStrength = 1; szzy1.SpecularColorReflectance = .9; szzy1.SpecularExponent = 30; szzy1.SpecularStrength = 1; s=surface(x,y,z,szzy1,'FaceLighting','phong','FaceColor','texture',... 'EdgeColor','none','Cdata',topo,'Parent',fa); %设置面对象并且赋给S light( 'position' ,[-1 0 1], 'color' ,[0.5 1 0.5]); light( 'position' ,[-1.5 0.5 -0.5], 'color' ,[.6 .2 .2 ]); light( 'Position' ,[1.5 1.5 -1]); light( 'Position' ,[0 -1.5 0], 'color' ,[0.6 0.6 1]); lighting gouraud set(gca,'CameraViewAngleMode','manual') ang=get(gca,'CameraViewAngle'); %获取相机视角 ang=.1*ang; %改变相机视角（跟上句都可不用） hold on; h=[1.46 1.5 1.54 1.57 1.60]; r0=1.45; for l=1:5 v0=h(l)*(2*pi/r0); %角动量守恒 r=[r0,0,0];v=[0,-v0,0]; %初始状态 state=[r(1),r(2),r(3),v(1),v(2),v(3)]; GM=4*pi^2; tau=-.002; nStep=[2290 2940 4015 5380 7805]; time=0; for i=1:nStep(l) state=rk6(state,time,tau,'gravrk',GM); r=[state(1) state(2) state(3)]; v=[state(4) state(5) state(6)]; time=time+tau; Sx(i)=r(1); Sy(i)=r(2); Sz(i)=r(3); timeS(i)=time; end plot3(Sx,Sy,Sz,'c-'),axis([-15,1,-10,10,-10,10]);axis off; t=line('Color',[1 0.4 0.4],'Marker','.','MarkerSize',20);%给卫星的点 n=length(Sx);i=1;j=1; view(166,21) while 1 set(gca,'CameraViewAngle',ang) drawnow rotate(s,[0 0 1],.3); set(t,'xdata',Sx(i),'ydata',Sy(i),'zdata',Sz(i)); pause(0.0001); i=i+10; if i>n break end end end