%%draw satellite orbit&place
function drawsatellite(a,e,W,i,w,M0,angle,k)
% a卫星轨道的长半轴.
% e是椭圆的偏心率.
% W是轨道升交点赤经
% i是轨道倾角
% w是轨道近地点幅角
% M0是初始平近点角
% angle是天线视角（假定雷达天线瞄准线指向卫星运动的右侧）
% k是观测次数
GM=3.986005e14;
we=sqrt(GM/(a^3));
T=2*pi/we;
t1=0:0.1:T;
M1=M0*2*pi/360+we*t1; %平近点角
Dk1=1;
Ek1=0;
n1=0;
while abs(Ek1-Dk1)>0.0000000001
n1=n1+1;
Ek1=Dk1;
Dk1=M1+e*sin(Ek1);
end
E1=Dk1; %迭代求偏进点角
x1=a*(cos(E1)-e); %轨道坐标
y1=a*sqrt((1-e^2))*sin(E1);
z1=0*E1;
t2=linspace(0,T,k);
M2=M0*2*pi/360+we*t2; %平近点角
Dk2=1;
Ek2=0;
n2=0;
while abs(Ek2-Dk2)>0.0000000001
n2=n2+1;
Ek2=Dk2;
Dk2=M2+e*sin(Ek2);
end
E2=Dk2; %迭代求偏进点角
x2=a*(cos(E2)-e); %卫星坐标
y2=a*sqrt((1-e^2))*sin(E2);
z2=0*E2;
DtoR=2*pi/360;
W0=W*DtoR;%升交点的经度
i0=i*DtoR;%轨道的倾角
w0=w*DtoR;%近地点幅角
R123=orbit_earth(W0,i0,w0);
Ans1=R123*[x1;y1;z1]; %轨道坐标系-惯性坐标系
x10=Ans1(1,:);
y10=Ans1(2,:);
z10=Ans1(3,:);
plot3(x10,y10,z10,'r--','LineWidth',1);%绘制轨道
hold on;
Ans2=R123*[x2;y2;z2]; %轨道坐标系-惯性坐标系
x20=Ans2(1,:);
y20=Ans2(2,:);
z20=Ans2(3,:);
plot3(x20,y20,z20,'d','MarkerFaceColor','y');%绘制卫星坐标
axis equal;
axis off;
hold on;
Ans3=targetplace(a,e,W,i,w,M0,angle,k);
x30=Ans3(1,:);
y30=Ans3(2,:);
z30=Ans3(3,:);
plot3(x30,y30,z30,'o','MarkerFaceColor','c');%绘制卫星波束中心点（地球表面）
[X,Y,Z]=sphere(100);
RE=0.64e7;
X=RE*X;
Y=RE*Y;
Z=RE*Z;
hold on;
mesh(X,Y,Z,'EdgeColor','k'); %绘制地球
legend('轨道','卫星','目标','地球');
end
评论0