- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
运用matlab计算gps卫星的坐标位置
clear
clc
formatlong
tp=input(tp=);
toc=input(toc=);
a0=input(a0=);
a1=input(a1=);
a2=input(a2=);
toe=input(toe=);
M0=input(M0=);
a=input(长半径a=);
deltan=input(卫星平均角速度之差deltan=);
e=input(e=);
w=input(w=);
Cuc=input(Cuc=);
Cus=input(Cus=);
Cic=input(Cic=);
Cis=input(Cis=);
Crc=input(Crc=);
Crs=input(Crs=);
i0=input(i0=);
I=input(轨道倾角变化率I=);
OM0=input(OM0=);
OM=input(升交点赤径变化率OM=);
tt=a0+a1*(tp-toc)+a2*(tp-toc);
t=tp-tt;
tk=t-toe;
u=3.986005e14;
n0=(sqrt(u))/(a*a*a);
n=n0+deltan;
Mk=M0+n*tk;
{
n=MK;
ek0=0;
ek1=n+e*sin(ek0);
ek2=n+e*sin(ek1);
EK=ek2;
}
Dk=1;
Ek=0;
n1=0;
whileabs(Ek-Dk)0.0000000001
n1=n1+1;
Ek=Dk;
Dk=Mk+e*sin(Ek);
end
Ek=Dk;
Vk=atan((sqrt(1-e*e)*sin(Ek)/(cos(Ek)-e));
ifsin(Ek)0cos(Ek)-e0
Vk=pi-Vk;
elseifsin(Ek)0cos(Ek)-e0
Vk=pi+Vk;
elseifsin(Ek)0cos(Ek)-e0
Vk=2*pi-Vk;
end
Faik=Vk+w;
SigmaU=Cuc*cos(2*Faik)+Cus*sin(2*Faik);
SigmaR=Crc*cos(2*Faik)+Crs*sin(2*Faik);
SigmaI=Cic*cos(2*Faik)+Cis*sin(2*Faik);
Uk=Faik+SigmaU;
Rk=a*(1-e*cos(Ek))+SigmaR;
Ik=i0+SigmaI+I*tk;
X0=Rk*cos(Uk);
Y0=Rk*sin(Uk);
we=7-5;
OMK=OM0+(OM-we)*tk-we*toe;
Xk=X0*cos(OMK)-Y0*cos(Ik)*sin(OMK);
Yk=X0*sin(OMK)+Y0*cos(Ik)*cos(OMK);
Zk=Y0*sin(Ik);
disp([卫星钟差改正dt=,num2str(tt)])
disp([归化时刻tk=,num2str(tk)])
disp([平均运行角速度n=,num2str(n)])
disp([卫星平近点角Mk=,num2str(Mk)])
disp([偏近点角Ek=,num2str(Ek)])
disp([真近点角Vk=,num2str(Vk)])
disp([升交距角Faik=,num2str(Faik)])
disp([摄动改正项SigmaU=,num2str(SigmaU)])
disp([摄动改正项SigmaR=,num2str(SigmaR)])
disp([摄动改正项SigmaI=,num2str(SigmaI)])
disp(经过摄动改正项:)
disp([升交距角Uk=,num2str(Uk)])
disp([卫星矢径Rk=,num2str(Rk)])
disp([轨道倾角Ik=,num2str(Ik)])
disp(卫星在轨道平面坐标系的坐标)
disp([X0=,num2str(X0)])
disp([Y0=,num2str(Y0)])
disp([观测时刻升交点经度OMK=,num2str(OMK)])
disp(卫星在地心固定坐标系中的直角坐标)
disp([Xk=,num2str(Xk)])
disp([Yk=,num2str(Yk)])
disp([Zk=,num2str
您可能关注的文档
- 心理健康趣味活动策划(10篇).pdf
- 南开大学考试卷_大学语文.pdf
- 二类医疗器械产品生产实验室要求.pdf
- 项目服务投标文件方案(成本管理及配置方案).pdf
- 酒店餐饮督察部岗位职责.pdf
- 高中英语主谓一致用法详解文档.pdf
- 一道美丽的风景线作文600字范文精选.pdf
- 全册教学反思部编版九年级历史下册.pdf
- 22春东北大学《数控机床与编程》在线作业一答案参考9.pdf
- 虚拟现实技术在博物馆展览中的应用探究.pdf
- 第三单元解决问题的策略 三年级下册数学高频考点重难点讲义(苏教版)(1).docx
- 4.4 解比例 数学六年级下册同步培优讲义(苏教版).docx
- 第四单元《根据方向与距离确定物体位置》(教案)五年级下册数学青岛版.docx
- (奥数典型题)第三讲 分解质因数 五年级下册数学思维拓展提升讲义(人教版).docx
- 2.2:2、5、3的倍数-人教版五年级数学下册第二单元:因数和倍数.docx
- 第六单元正比例和反比例 六年级下册数学高频考点重难点讲义(苏教版).docx
- 第二单元除数是一位数的除法 三年级下册数学高频考点重难点讲义(人教版).docx
- 第二单元-认识三角形和四边形 四年级数学下册提升(北师大版).docx
- 第一单元观察物体(三)五年级下册数学高频考点重难点讲义(人教版).docx
- 第九单元 数学广角—推理 二年级数学下册重难点知识点(人教版).docx
最近下载
- 涂药法-中医护理技术操作规程.ppt VIP
- 弱势群体社会医学(医学课件).ppt VIP
- 小学语文四年级上册《语文园地六》作业设计.docx
- 2024-2030年中国鹿行业经营风险与发展前景趋势预测分析研究报告.docx
- GB_T27065-2015_合格评定产品、过程和服务认证机构要求.doc
- 个人能力就业分析报告.pdf VIP
- 成语故事《三顾茅庐》故事简介PPT.pptx VIP
- 现代遗传学(厦大)中国大学MOOC慕课 章节测验期中期末考试答案.docx
- 北京市西城区2022-2023学年六年级上学期英语期末试卷(音频暂未更新).docx VIP
- 在役聚乙烯PE地埋燃气管道失效可能性评分表、后果评分表.pdf VIP
文档评论(0)