卫星星下点轨迹Matlab仿真.docx

卫星星下点轨迹Matlab仿真.docx

《卫星星下点轨迹Matlab仿真.docx》由会员分享，可在线阅读，更多相关《卫星星下点轨迹Matlab仿真.docx（8页珍藏版）》请在冰豆网上搜索。

卫星星下点轨迹Matlab仿真

利用Matlab绘制卫星星下点轨迹

1.地球静止轨道卫星，倾角别离为0，30，90度。

clc;clear;

t=0:

1:

6;

we=360/24;

u=we*t;

i=30;

fai=asind（sind（i）*sind（u））;

deltalmd=atand（cosd（i）*tand（u））;

if（i==90）

deltalmd（end）=90;

end

lmd=deltalmd-we*t;

%usesymetrytogeneratetheotherdata

forj=1:

6

lmd（j+7）=-lmd（7-j）;

fai（j+7）=fai（7-j）;

end

forj=1:

12

lmd（j+13）=lmd（13-j）;

fai（j+13）=-fai（13-j）;

end

h=geoshow（'','FaceColor',[111]）;

gridon

holdon

plot（lmd,fai）;title（['GEO¹ìµÀ£¬Çã½Çi=',num2str（i）]）

2.回归轨道卫星，回归周期1天，倾角别离为60度，周期为4h。

clc;clear;

t=[01/31/22/34/51];

we=360/24;

w=180/2;

u=w*t;

i=60;

fai=asind（sind（i）*sind（u））;

deltalmd=atand（cosd（i）*tand（u））;

lmd=deltalmd-we*t;

%usesymetrytogeneratetheotherdata

forj=1:

5

lmd（j+6）=lmd（6）+（lmd（6）-lmd（6-j））;

fai（j+6）=fai（6-j）;

end

forj=1:

10

if（lmd（11）+（lmd（11）-lmd（11-j）））>180

lmd（j+11）=-180+rem（lmd（11）+（lmd（11）-lmd（11-j））,180）;

else

lmd（j+11）=lmd（11）+（lmd（11）-lmd（11-j））;

end

fai（j+11）=-fai（11-j）;

end

cnt=1;

form=1:

5

forj=1:

21

if（lmd（j+21*（m-1））+60）>180

lmd（j+21*m）=-180+rem（lmd（j+21*（m-1））+60,180）;

record（m,cnt）=j;%recordwhentranversefromeasttowest

cnt=cnt+1;

else

lmd（j+21*m）=lmd（j+21*（m-1））+60;

end

fai（j+21*m）=fai（j+21*（m-1））;

end

cnt=1;

end

loadstill

h=geoshow（'','FaceColor',[111]）;

gridon

holdon

plot（lmd1（2:

20）,fai1（2:

20）,'b--'）;%earthstill

plot（lmd（1:

6）,fai（1:

6）,'bo'）;

plot（lmd（21*6）,fai（21*6）,'bo'）;

plot（lmd（1:

13）,fai（1:

13））;plot（lmd（14:

21）,fai（14:

21））;

form=1:

5

plot（lmd（21*m+1:

record（m,1）+21*m-1）,fai（21*m+1:

record（m,1）+21*m-1））;plot（lmd（record（m,1）+21*m:

21*（m+1））,fai（record（m,1）+21*m:

21*（m+1）））;

plot（lmd（21*m）,fai（21*m）,'bo'）;

end

title（['ÐÇϵã¹ì¼££ºT=4h¹ìµÀ£¬Çã½Çi=',num2str（i）]）

地球不转时的星下点

clc;clear;

t=[01/31/22/34/51];

we=360/24;

w=180/2;

u=w*t;

i=60;

fai=asind（sind（i）*sind（u））;

deltalmd=atand（cosd（i）*tand（u））;

lmd=deltalmd;%earthstill

%usesymetrytogeneratetheotherdata

forj=1:

5

lmd（j+6）=lmd（6）+（lmd（6）-lmd（6-j））;

fai（j+6）=fai（6-j）;

end

forj=1:

10

if（lmd（11）+（lmd（11）-lmd（11-j）））>180

lmd（j+11）=-180+rem（lmd（11）+（lmd（11）-lmd（11-j））,180）;

else

lmd（j+11）=lmd（11）+（lmd（11）-lmd（11-j））;

end

fai（j+11）=-fai（11-j）;

end

forj=1:

21

if（lmd（j）+180）>180

lmd（j）=-180+rem（lmd（j）+180,180）;

else

lmd（j）=lmd（j）+180;

end

fai（j）=fai（j）;

end

lmd（11）=0;

lmd1=lmd;

fai1=fai;

savestilllmd1fai1

h=geoshow（'','FaceColor',[111]）;

gridon

holdon

plot（lmd（2:

20）,fai（2:

20））;

%plot（lmd（1:

13）,fai（1:

13））;plot（lmd（14:

21）,fai（14:

21））;

%form=1:

5

%plot（lmd（21*m+1:

record（m,1）+21*m-1）,fai（21*m+1:

record（m,1）+21*m-1））;plot（lmd（record（m,1）+21*m:

21*（m+1））,fai（record（m,1）+21*m:

21*（m+1）））;

%end

title（['T=4h¹ìµÀ£¬Çã½Çi=',num2str（i）]）