卫星星历是描述卫星运动轨道的信息。也可以说卫星星历就是一组对应某一时刻的轨道参数及其变率。有了卫星星历就可以计算出任意时刻的卫星位置及其速度。GPS卫星星历分为预报星历和后处理星历。预报星历又称广播星历。

GPS广播星历参数共有16个,其中包括1个参考时刻,6个对应参考时刻的开普勒轨道参数和9个反映摄动力影响的参数。这些参数通过GPS卫星发射的含有轨道信息的导航电文传递给用户。

1.星历参考时刻 :

2.长半轴平方根:

3.偏心率:

4.参考历元下平近点角:

5.近地点角距:

6.轨道倾角:

7.本周初始历元的升交点赤经:

8.平运动差(由精密星历计算得到的卫星平均角速度与按给定参数计算所得的平均角速度之差):n

9.轨道倾角变化率(弧度/秒):

10.升交点赤经变化率(弧度/秒):

11.纬度幅角的余弦调和项改正的振幅(弧度)

12.纬度幅角的正弦调和项改正的振幅(弧度):

13.轨道半径的余弦调和项改正的振幅(m):

14.轨道半径的正弦调和项改正的振幅(m):

15.轨道倾角的余弦调和项改正的振幅(弧度):

16.轨道倾角的正弦调和项改正的振幅(弧度):

以下为matlab程序:

%the homework 1 of chapater 1
%student :Taylen
%time :2018/9/25
%程序功能:根据所提供的星历参数,计算此卫星在信号发射时刻t( GPS时间)239050.7223s
%时的时空位置%GPS广播星历参数共有16个,其中包括1个参考时刻,6个对应参考时刻的开普勒轨道参数
%和9个反映摄动力影响的参数。这些参数通过GPS卫星发射的含有轨道信息的导航电文传递给用户。
%时间参数
pra1 = 244800;           %te---星历参考时刻,即星历表参考历元(s)%开普勒六参数
pra2 = -1.064739758;     %M0---按参考历元te计算的平近点角(弧度)
pra3 = 0.005912038265;   %e---轨道偏心率
pra4 = 5153.65531;       %sqrt(a)长半轴平方根
pra5 = -1.717457876;     %w0---近地点角距
pra6 = 0.9848407943;     %i0---按参考历元计算的轨道倾角(弧度)
pra7 = 1.038062244;      %Ω0---本周初始历元的升交点赤经(弧度)%轨道摄动九参数
pra8  = 4.249105564e-9;      % Δ?---平运动差(弧度)
pra9  = 7.422851197e-51;     % dot_i---轨道倾角变化率
pra10 = -8.151768125e-9;     % dot_Omega---升交点赤经变化率
pra11 = 3.054738045e-7;      % Cuc---升交点角距的调和改正项振幅
pra12 = 2.237036824e-6;      % Cus---升交点角距的调和改正项振幅
pra13 = 350.53125;           % Crc---卫星地心距的调和改正项振幅
pra14 = 2.53125;             % Crs---卫星地心距的调和改正项振幅
pra15 = -8.381903172e-8;     % Cic---轨道倾角的调和改正项振幅
pra16 = 8.940696716e-8;      % Cis---轨道倾角的调和改正项振幅t=239050.7223;  %此卫星在信号发射时刻t ( GPS时间)miu = 3.986005e14;      % 地心引力常数
we = 7.2921151467e-5;   % 地球自转角速度
F = -4.442807633e-10;   % Constant, [sec/(meter)^(1/2)]% 计算归化时间
te = pra1;
tk = t - te;
% 计算卫星的平均角速度
a = pra4.^2;
n0 = sqrt(miu/a.^3);
n = n0 + pra8;
% 计算平近点角
Mk = pra2 + n * tk;
% 计算偏近点角Ek
%迭代计算:相邻两次计算差之绝对值值<1e-15时结束迭代计算,Ek的迭代初始值为0
Ek0=0;Ek=Mk;
while abs(Ek0-Ek)>1e-15Ek0 = Ek;Ek=Mk + pra3 * sin(Ek0);
end
% 计算卫星钟差相对论校正值 Compute relativistic correction term% 计算真近点角:取值在(-pi,pi]
v1=sqrt(1-pra3^2)*sin(Ek);
v2=cos(Ek)-pra3;
vk=atan(v1/v2);
if v1>0 & v2<0
Vk = pi + vk;%第二象限
elseif v1<0 & v2<0
Vk = vk - pi;%第三象限
end
% 计算升交点角距
Faik = Vk + pra7;
% 计算升交点角距改正值
Sigmauk=pra11*cos(2*Faik)+pra12*sin(2*Faik);
% 计算卫星地心向径改正值
Sigmark=pra13*cos(2*Faik)+pra14*sin(2*Faik);
% 计算卫星轨道倾角改正值
Sigmaik=pra15*cos(2*Faik)+pra16*sin(2*Faik);
% 修正升交点角距
uk=Faik+Sigmauk;
% 修正卫星地心相径
rk=a*(1-pra3*cos(Ek))+Sigmark;
% 修正卫星轨道倾角
ik=pra6+Sigmaik+tk*pra9;
% 计算卫星轨道平面直角坐标
xk=rk*cos(uk);
yk=rk*sin(uk);
% 计算升交点赤经
Omegak=pra5+(pra10-we)*tk-we*te;
% 将卫星坐标由轨道平面转换到ECEF坐标系
Xk=xk*cos(Omegak)-yk*cos(ik)*sin(Omegak);
Yk=xk*sin(Omegak)+yk*cos(ik)*cos(Omegak);
Zk=yk*sin(ik);

卫星导航定位技术二:由星历参数求解卫星时空位置相关推荐

  1. 全球卫星导航定位技术的现状

    全球卫星导航定位技术的现状与发展 [摘 要]随着卫星导航定位系统的发展壮大,全球定位系统进去了与定位.导航.授时.移动通信和因特网等多源信息载体融合的新时期,一直以来我国对全球卫星导航定位系统的探索使 ...

  2. 青岛科技大学|物联网工程|物联网定位技术(第二讲)|15:00

    目录 物联网定位技术(第二讲) 1. 卫星的轨道高度与覆盖区域有何关系,试画图给予说明覆盖区地心角与覆盖面积的关系 2. 试给出实际的卫星地面覆盖区和用户空间可视区所对应的半地心角的公式并请给予解释 ...

  3. 室内定位技术方案---Wifi、RFID、bluetooth、Zigbee

    室内定位技术方案 一.室内定位系统原理 国内常规定位技术有Wifi定位.RFID定位.蓝牙定位.Zigbee定位等几种方式,这几种定位方式有什么不同,各有什么利弊,下面我们将详细分析这几种定位原理.优 ...

  4. UWB定位,新一代的精确定位技术

    超宽带(UWB)无线定位技术由于功耗低.抗多径效果好.安全性高.系统复杂度低,尤其是能提供非常精确的定位精度等优点,而成为未来无线定位技术的热点和首选.文章将UWB定位技术与其他无线定位技术就定位环境 ...

  5. 【转载】导航定位技术方案

    早在15世纪,当人类开始探索海洋的时候,定位技术也随之催生.当时的定位方法十分粗糙,就是是运用航海图和星象图以确定自己的位置. 随着社会的进步和科技的发展,定位技术在技术手段.定位精度.可用性等方面均 ...

  6. 10种室内定位技术原理深度解析

    在高度城市化的今天,室内空间越来越庞大复杂.人类战胜了大自然,却在自己构筑的钢筋水泥中迷了路. 东晋高僧法显西行求法乘船回国时写道:大海弥漫无边,不识东西,唯望日.月.星宿而进.大意就是说,大海辽阔无 ...

  7. 灵思科电子科技—室内定位技术有哪些_七大室内定位技术详解

    1.蓝牙室内定位技术 蓝牙室内技术是利用在室内安装的若干个蓝牙局域网接入点,把网络维持成基于多用户的基础网络连接模式,并保证蓝牙局域网接入点始终是这个微微网(piconet)的主设备,然后通过测量信号 ...

  8. 老司机揭秘手机定位技术,实在是精彩!

    本文来源于鲜枣课程 今天我们所处的移动互联网时代,手机成了每个人的生活标配. 这些手机里,安装了形形色色的APP,提供了各种服务,彻底改变了我们的生活. 这些服务里面,就包括我们今天的主角--定位. ...

  9. 智慧楼宇篇 3 —— 室内定位技术(三) - UWB室内定位技术

    UWB室内定位技术 1  引言 本文探讨室内定位技术中的一种:UWB室内定位技术,并在定位技术系列最后对各种定位技术进行总结,敬请关注微信公众号"智物客"后续文章. 2  概述 U ...

最新文章

  1. 【高并发】32位多核CPU并发读写long型数据为何会出现诡异问题?看完这篇我懂了!...
  2. 四、垃圾收集之垃圾收集算法
  3. 寄存器和存储器的区别_寄存器、累加器、暂存器都是什么?它们有什么区别?...
  4. FD.io/VPP — VPP Agent — Telemetry Plugin
  5. Centos7 安装 nginx 服务器的两种方式
  6. HttpClient接口Post方式
  7. sql安装目录下log文件夹_Linux安装Hive数据仓库工具
  8. 移民申请--申请澳大利亚技术移民的职业评估有技巧(2)
  9. html用jq设置动态效果,jQuery实现基本动画效果的方法详解
  10. Redis - Redis command timed out nested exception is io.lettuce.core.RedisCommandTimeoutException
  11. python中forward的作用_基于numpy的前馈神经网络(feedforward neural network)
  12. python dataframe索引转成列_Pandas之DataFrame对象的列和索引之间的转化
  13. Windows Embedded Webcast 2008年1月预告
  14. js中的call和apply方法的区别
  15. ASP.NET vNext MVC 6 电商网站开发实战
  16. [系统安全] 十五.Chrome密码保存功能渗透解析、Chrome蓝屏漏洞及音乐软件漏洞复现
  17. 设置电脑的保护色(绿豆沙色)
  18. Hubble数据库x某股份制商业银行:一款混布数据库的“交易+分析”提速之路
  19. PRD文档范例,千万收藏的产品经理写作手册
  20. 学驾照,科目一计分题

热门文章

  1. 鼠标控制两台台式机_台式机和控制台上的真棒实时GI
  2. 如何复制word中带修订标记的段落
  3. python是开源的跨平台编程语言_Python是一种跨平台、开源、免费的动态编程语言。...
  4. 具有不安全、不正确或缺少SameSite属性的Cookie
  5. 软件测试教程 自动化测试appium篇
  6. java基础巩固-宇宙第一AiYWM:为了维持生计,架构知识+分布式微服务+高并发高可用高性能知识序幕就此拉开(一:总览篇)~整起
  7. Vue3项目引入 vue-quill 编辑器组件并封装使用
  8. QBC和QBE查询方法
  9. 什么是日志服务(事件日志服务)
  10. mysql 查重名现象_mysql查询重复数据【转】