卫星坐标计算实验_第1页
卫星坐标计算实验_第2页
卫星坐标计算实验_第3页
卫星坐标计算实验_第4页
卫星坐标计算实验_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

西南交通大学卫星坐标计算实验实验报告专 业： 班 级： 姓 名： 学 号： 成 绩： 2017年3月20日一、实验步骤：1、将星历文件中需要的数据作为矩阵导入MATLAB中的变量。2、编写计算程序代码1）计算平均角速度：2）改正平角速度n：3）计算平近点角：4）通过迭代计算，计算偏近点角：5）计算真近点角：6）计算升交距角：7）计算卫星轨道摄动项改正数：8）计算改正后的真近点角：9）计算改正后的向径：10）计算改正后的倾角：11）计算轨道平面内的坐标：12）计算观测瞬间升交点的经度：13）计算旋转矩阵： 14）卫星坐标：15）计算钟差：二、实现代码：function pos=SPOS(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 =br(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;%计算钟差ddt=a0+
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论