认识星历文件观测文件及位置计算

初学者，欢迎指正批评。qigezhubajie@163.com

星历文件

星历文件表头部分不做介绍。
第一张图片为导航数据文件中数据记录部分说明，第二章为数据记录部分。
第二张图片中的数据部分与第一张图片中的参数一一对应。
第二张图片中加灰部分为数据有效范围。

观测文件

观测文件表头部分不做介绍。

matlab面对对象卫星位置计算

时间转化模块不做介绍，自行查阅文件。
一定要自己敲一遍，理解的更深刻。
classdef satlite < handle
properties
PRN; af2; M0; sqrtA; delta_n; e ;
omega; Cuc; Cus; Crc; Crs; i0;
i_dot; Cic; Cis; OMEGA; OMEGA_dot;
toe ; af0; af1; toc;
X ; Y ; Z;
miu = 3.986005e14; % 地心引力常数
we = 7.2921151467e-5; % 地球自转角速度
tg%观测时间
end
methods
function self = satlite(PRN, af2, M0, sqrtA, delta_n, e, omega, …
Cuc, Cus, Crc, Crs, i0, i_dot, Cic, Cis, OMEGA, OMEGA_dot, …
toe, af0, af1, toc,tg)
self.PRN = PRN;self.af2 = af2;self.M0 = M0;
self.sqrtA = sqrtA;self.delta_n = delta_n;self.e = e;
self.omega = omega;self.Cuc = Cuc;self.Cus = Cus;
self.Crc = Crc;self.Crs = Crs;self.i0 = i0;
self.i_dot = i_dot;self.Cic = Cic;self.Cis = Cis;
self.OMEGA = OMEGA;self.OMEGA_dot = OMEGA_dot;
self.toe = toe;self.af0 = af0;
self.af1 = af1;self.toc = toc;
self.tg = tg;
end
function position(self)
%计算平均角速度
a = self.sqrtA^2;
n0 = sqrt(self.miu)/(self.sqrtA^3);
n = n0 + self.delta_n;
%计算归化时间
tb = self.af0 + self.af1*(self.tg - self.toc) + self.af2*(self.tg - self.toc)^2;
ts = self.tg - tb;
tk = ts - self.toe;
% 计算平近点角
Mk = self.M0 + n * tk;
% 计算偏近点角Ek
%迭代计算：相邻两次计算差之绝对值值<1e-15时结束迭代计算,Ek的迭代初始值为0
Mk = rem(Mk+2pi,2pi);
Ek0=0;Ek=Mk;
while abs(Ek0-Ek)>1e-15
Ek0 = Ek;
Ek=Mk + self.e * sin(Ek0);
end
% 计算真近点角:取值在（-pi,pi]
v1=sqrt(1-self.e^2)sin(Ek);
v2=cos(Ek)-self.e;
vk=atan2(v1,v2);
% 计算升交点角距
Faik = vk + self.omega;
% 计算升交点角距改正值摄动改正项计算
Sigmauk=self.Cuccos(2Faik)+self.Cussin(2Faik);
% 计算卫星地心向径改正值
Sigmark=self.Crccos(2Faik)+self.Crssin(2Faik);
% 计算卫星轨道倾角改正值
Sigmaik=self.Ciccos(2Faik)+self.Cissin(2Faik);
% 修正升交点角距
uk=Faik+Sigmauk;
% 修正卫星地心相径
rk=a(1-self.ecos(Ek))+Sigmark;
% 修正卫星轨道倾角
ik=self.i0+Sigmaik+tkself.i_dot;

        % 计算卫星轨道平面直角坐标xk=rk*cos(uk);yk=rk*sin(uk);% 计算升交点赤经Omegak=self.OMEGA + (self.OMEGA_dot-self.we)*tk-self.we*self.toe;% 将卫星坐标由轨道平面转换到ECEF坐标系self.X=xk*cos(Omegak)-yk*cos(ik)*sin(Omegak);self.Y=xk*sin(Omegak)+yk*cos(ik)*cos(Omegak);self.Z=yk*sin(ik);endend

end

卫星高度角方位角计算（待更新。。。）

基于MATLAB计算卫星位置相关推荐

matlab计算空间坐标,通过matlab计算卫星位置
卫星星历是描述卫星运动轨道的信息.也可以说卫星星历就是一组对应某一时刻的轨道参数及其变率.有了卫星星历就可以计算出任意时刻的卫星位置及其速度.GPS卫星星历分为预报星历和后处理星历.预报星历又称广播星 ...
python读取导航电文并计算卫星位置
python简单计算卫星位置前言一.思路总的可分为两个部分:获取导航参数和计算卫星位置. ①获取导航参数:首先讲导航星历中的数据切片,存入csv文件中,再读取csv文件的数据并赋值给各参数 ②计 ...
如何用Matlab建立信用风险模型,基于Matlab计算的KMV模型在商业银行信用风险管理中的实践应用.pdf...
科技论坛 2014.6 基于Matlab计算的KMV模型在商业银行信用风险管理中的应用李园 (天津大学管理与经济学部,天津,300072) 摘要:社会经济的不断发展,金融行业也处于不断发展之中, ...
用matlab计算卫星的瞬时位置,用MATLAB计算GPS卫星位置.PDF
用MATLAB计算GPS卫星位置高新技术用MATLAB计算GPS卫星位置罗利娟杨乐 (西安翻译学院陕西西安 710061) 摘要:本文主要介绍了GPs测量数据的常用格式RINEX标准 ...
用c++根据轨道六根数计算卫星位置
轨道六根数是描述卫星轨道的一组参数,包括: 轨道长半径(a):卫星轨道的半径,表示卫星到地球中心的平均距离. 轨道偏心率(e):卫星轨道的偏心率,表示轨道的椭圆程度. 轨道倾角(i):卫星轨道与地球赤 ...
用matlab计算卫星的瞬时位置,用MATLAB计算GPS卫星位置
GPS定位的基本原理简单来说就是在多时候需要把来自不同型号的接收机的RINEX标准文件里不是单一的一个文WGS-84空间直角坐标系中,确定未知点与数据放在一块进行处理,这就需要数据件,而是包括如下几种 ...
基于matlab计算aom衍射效率,傅里叶望远镜大气湍流模拟实验
引用本文曾志红, 罗秀娟, 王保峰, 夏爱利, 程志远, 司庆丹. 傅里叶望远镜大气湍流模拟实验.光子学报, 2014,43(6): 601002 ZENG Zhi-hong, LUO Xiu-ju ...
基于matlab计算aom衍射效率,实验七激光位相调制光谱特性检测-华东师范大学.ppt...
<实验七激光位相调制光谱特性检测-华东师范大学.ppt>由会员分享,可在线阅读,更多相关<实验七激光位相调制光谱特性检测-华东师范大学.ppt(9页珍藏版)>请在人人文库网上搜 ...
采用morison方程基于matlab计算大直径波浪力,用Morison方程计算分析悬浮隧道所受波浪力初探...
用Morison方程计算分析悬浮隧道所受波浪力初探用Morison方程计算分析悬浮隧道所受波浪力初探第16卷第3期石家庄铁道学院第卷第期后豕饮迥r子阮子很 2003年9月JOURNALOFSHIJ ...

基于MATLAB计算卫星位置

matlab卫星定位