卫星目标位置.rar_卫星位置仿真_坐标系_瞄准_瞄准点_轨迹坐标

%%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