卫星星下点轨迹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）]）
升级会员 