衛(wèi)星坐標(biāo)計(jì)算實(shí)驗(yàn)驕陽(yáng)教學(xué)

1、西南交通大學(xué)衛(wèi)星坐標(biāo)計(jì)算實(shí)驗(yàn)實(shí)驗(yàn)報(bào)告專(zhuān) 業(yè)： 班 級(jí)： 姓 名： 學(xué) 號(hào)： 成 績(jī)： 2017年3月20日習(xí)題與教育一、實(shí)驗(yàn)步驟：1、將星歷文件中需要的數(shù)據(jù)作為矩陣導(dǎo)入matlab中的變量。2、編寫(xiě)計(jì)算程序代碼1）計(jì)算平均角速度：2）改正平角速度n：3）計(jì)算平近點(diǎn)角：4）通過(guò)迭代計(jì)算，計(jì)算偏近點(diǎn)角：5）計(jì)算真近點(diǎn)角：6）計(jì)算升交距角：7）計(jì)算衛(wèi)星軌道攝動(dòng)項(xiàng)改正數(shù)：8）計(jì)算改正后的真近點(diǎn)角：9）計(jì)算改正后的向徑：10）計(jì)算改正后的傾角：11）計(jì)算軌道平面內(nèi)的坐標(biāo)：12）計(jì)算觀測(cè)瞬間升交點(diǎn)的經(jīng)度：13）計(jì)算旋轉(zhuǎn)矩陣： 14）衛(wèi)星坐標(biāo)：15）計(jì)算鐘差：二、實(shí)現(xiàn)代碼：function pos=spos

2、(br,dt)gm=3.986005e+14;we=7.2921151467e-5;%a0 a1 a2a0=br(1,2);a1=br(1,3);a2=br(1,4);% iode crs dn m0crs=br(2,2); dn=br(2,3); m0=br(2,4);% cuc e cus sqacuc=br(3,1);e=br(3,2); cus=br(3,3);sqa=br(3,4);% toe cic om0 cistoe =br(4,1);cic =br(4,2); om0 =br(4,3);cis =br(4,4);% i0 crc w domi0 =br(5,1);crc =b

3、r(5,2); w =br(5,3);dom =br(5,4);% di weekdi=br(6,1);n0=sqrt(gm)/sqa3;n=n0+dn;mk=m0+n*dt;ek=mk;te=inf;while abs(te-ek)1e-14 te=ek; ek=mk+e*sin(ek);endxx=cos(ek)-e;yy=sqrt(1-e*e)*sin(ek);v=atan2(yy,xx);u=v+w;du=cuc*cos(2*u)+cus*sin(2*u);dr=crc*cos(2*u)+crs*sin(2*u);di=cic*cos(2*u)+cis*sin(2*u);uk=u+du;rk=sqa*sqa*(1-e*cos(ek)+dr;ik=i0+di+di*dt;x=rk*cos(uk);rk*sin(uk);0;rx=1 0 0;0 cos(ik) -sin(ik);0 sin(ik) cos(ik);lt=om0+(dom-we)*(toe+dt)-dom*toe;rz=cos(lt) -sin(lt) 0;sin(lt) cos(lt) 0;0 0 1;%計(jì)算鐘差ddt=a0+a1*(