matlab绘制星下点轨迹,卫星星下点轨迹Matlab仿真.docx

使用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;endlmd = deltalmd - we*t;% use symetry to generate the other datafor j = 1:6 lmd(j+7) = -lmd(7-j); fai(j+7) = fai(7-j);endfor j = 1:12 lmd(j+13) = lmd(13-j); fai(j+13) = -fai(13-j);end h = geoshow('landareas.shp', 'FaceColor', [1 1 1]);grid onhold onplot(lmd, fai); title(['GEO1ìμà￡?????i=', num2str(i)])2.回归轨道卫星，回归周期1天，倾角分别为60度，周期为4h。clc; clear;t = [0 1/3 1/2 2/3 4/5 1];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; % use symetry to generate the other datafor j = 1:5 lmd(j+6) = lmd(6) + ( lmd(6) - lmd(6-j) ); fai(j+6) = fai(6-j);endfor j = 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);endcnt = 1;for m = 1:5 for j = 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; % record when tranverse from east to west 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 load stillh = geoshow('landareas.shp', 'FaceColor', [1 1 1]);grid onhold onplot(lmd1(2:20), fai1(2:20), 'b--'); % earth stillplot(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));for m = 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');endtitle(['D???μ?1ì?￡￡oT=4h1ìμà￡?????i=', num2str(i)])地球不转时的星下点clc; clear;t = [0 1/3 1/2 2/3 4/5 1];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;% earth still% use symetry to generate the other dat