Matlab根据广播星历表计算卫星坐标

clear,clc;%读取卫星星历文件
filename='20040130.04N.txt';
fid=fopen(filename);%按行读取星历文件的每一个元素，从PRN7号卫星开始
data=textscan(fid,'%f %f %f %f %f %f %f %f %f
%f\n','HeaderLines',11);%将读取的数组转换成数值型矩阵
B1=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B2=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B3=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B4=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B5=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B6=cell2mat(data);
data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
B7=cell2mat(data);
data=textscan(fid,'%f\n','HeaderLines',1);
B8=cell2mat(data);%定义常量GM,Ve
GM=3.986004e14;
Ve=7.2921151467e-5;%求出轨道长半轴a
sqrta=B3(4);
a=sqrta^2;
e=B3(2);%求出平均角速度n0
n0=sqrt(GM/(a^3));
XYZ=[];%计算从2004年1月30日16：00开始，每隔1分钟的卫星坐标
for i=0:1:20%计算当前时刻的儒略日JulianDay=fix(365.25*(2004-1))+fix(30.6001*(1+13))+30+(16+i/60+0/3600)/24+1720981.5;%将儒略日转换成GPS时t2=(JulianDay-2444244.5)/7-fix((JulianDay-2444244.5)/7);t=t2*604800;%计算卫星钟差改正后的当前时刻的GPS时t=t-(B1(8)+B1(9)*(t-B4(1))+B1(10)*(t-B4(1))^2);%计算参考时刻的儒略日JulianDayoe=fix(365.25*(2004-1))+fix(30.6001*(1+13))+30+(8+0/60+0/3600)/24+1720981.5;%将儒略日转换成GPS时t2oe=(JulianDayoe-2444244.5)/7-fix((JulianDayoe-2444244.5)/7);toe=t2oe*604800;%计算时间差tktk=t-toe;%计算改正平角速度n=n0+B2(3);%计算平近点角EkMk=B2(4)+n*tk;%迭代计算，求出偏近点角EkEk=Mk;Ek0=inf;whileabs(Ek0-Ek)>1e-100Ek0=Ek;Ek=Mk+e*sin(Ek0);end%计算真近点角Sinfk=(sqrt(1-e^2)*sin(Ek))/(1-e*cos(Ek));Cosfk=(cos(Ek)-e)/(1-e*cos(Ek));fk=atan2(Sinfk,Cosfk);%计算升交角距phik=fk+B5(3);%从卫星星历文件中读取各个方向的振幅Crs,Cuc,Cus,Cic,Cis,CrcCrs=B2(2);Cuc=B3(1);Cus=B3(3);Cic=B4(2);Cis=B4(4);Crc=B5(2);%计算卫星轨道摄动项改正数deltauk=Cus*sin(2*phik)+Cuc*cos(2*phik);deltark=Crs*sin(2*phik)+Crc*cos(2*phik);deltaik=Cis*sin(2*phik)+Cic*cos(2*phik);%计算改正后的升交角距ukuk=phik+deltauk;%计算改正后的向径rkrk=a*(1-e*cos(Ek))+deltark;%计算改正后的倾角ikik=B5(1)+deltaik+B6(1)*tk;%计算观测瞬间升交点的经度LkLk=B4(3)+(B5(4)-Ve)*tk-Ve*B4(1);%计算卫星在轨道平面内的坐标x,y,zxyz=[rk*cos(uk);rk*sin(uk);0];%计算旋转矩阵的值RzRx=[cos(Lk),-sin(Lk)*cos(ik),sin(Lk)*sin(ik);sin(Lk),cos(Lk)*cos(ik),-cos(Lk)*sin(ik);0,sin(ik),cos(ik)];%计算卫星在协议地球坐标系中的位置XYZ=vpa([XYZ;1e-3*(RzRx*xyz)'],10);
endXYZ

卫星星历文件

Matlab根据广播星历表计算卫星坐标相关推荐

用Java编程卫星坐标_2 分钟掌握卫星星座图
话不多说,助您 2 分钟掌握卫星星座图,知晓原理,并能绘制效果基础什么是卫星星座? 详细科普不在本文浪费时间了,完全不懂的同学瞄一下卫星星座-百度百科就可以了. 实现储备知识简单解释一下所需 ...
利用广播星历计算卫星在ECEF下的坐标
#define _CRT_SECURE_NO_WARNINGS 1 #include <stdio.h> #include <stdlib.h> #include <ma ...
GPS卫星坐标用计算机计算,GPS卫星位置的计算
[摘要]针对GPS卫星导航电文的特点和里边各种参数含义及用法,阐述了利用导航电文提供的数据计算卫星位置的方法,为进行基线解算时提供准确的卫星坐标数据. 0前言全球卫星定位系统(GPS)是一个实时的 ...
卫星过顶计算matlab,一种圆轨道低轨卫星过顶剩余可见时长预测方法
一种圆轨道低轨卫星过顶剩余可见时长预测方法 [技术领域] [0001] 本发明涉及卫星通信领域,特别是涉及一种圆轨道低轨卫星过顶剩余可见时长预测方法. [背景技术] [0002] 圆轨道低轨卫星是指 ...
基于MATLAB计算卫星位置
matlab卫星定位认识星历文件观测文件及位置计算星历文件观测文件 matlab面对对象卫星位置计算卫星高度角方位角计算(待更新...) 认识星历文件观测文件及位置计算初学者,欢迎指正批评. ...
卫星过顶计算matlab,基于SGP4模型的卫星轨道计算.docx
基于 SGP4 模型的卫星轨道计算刁宁辉 ① , 基于 SGP4 模型的卫星轨道计算刁宁辉 ① ,刘建强 ① ,孙从容 ① ,孟鹏 ② (① 国家卫星海洋应用中心 ,北京 ...
c语言编程gps卫星坐标计算公式,C语言计算GPS卫星位置演示教学.doc
C语言计算GPS卫星位置演示教学.doc C语言计算GPS卫星位置1 概述在用GPS信号进行导航定位以及制订观测计划时,都必须已知GPS卫星在空间的瞬间位置.卫星位置的计算是根据卫星电文所提供的轨道参 ...
基于卫星星历计算卫星在CGCS2000大地坐标系中的坐标
目录一.北斗系统概述 1.空间星座 2.坐标系统 3.时间系统二.实验目的三.实验内容四.实验过程五.实验结果一.北斗系统概述 1.空间星座北斗卫星导航系统简称北斗系统,英文缩写为 BD ...
matlab多项式计算题目,MATLAB数据分析与多项式计算-习题答案
<MATLAB数据分析与多项式计算-习题答案>由会员分享,可在线阅读,更多相关<MATLAB数据分析与多项式计算-习题答案(4页珍藏版)>请在人人文库网上搜索. 1.第6章 M ...
matlab 水文,MATLAB在水文水能计算中的应用
第 25卷第 1期 2006年第 1期红水河 HongShui River V01．25．No．1 No．1．20o6 MATLAB在水文水能计算中的应用黄飞仁 ,黄汉球2,曾英先 (1．广西水电 ...

Matlab根据广播星历表计算卫星坐标

Matlab根据广播星历表计算卫星坐标

卫星星历文件

Matlab根据广播星历表计算卫星坐标相关推荐

最新文章

热门文章