CN102590831B - 载体旋转条件下判断卫星可见性的方法和装置 - Google Patents

载体旋转条件下判断卫星可见性的方法和装置 Download PDF

Info

Publication number
CN102590831B
CN102590831B CN2012100317718A CN201210031771A CN102590831B CN 102590831 B CN102590831 B CN 102590831B CN 2012100317718 A CN2012100317718 A CN 2012100317718A CN 201210031771 A CN201210031771 A CN 201210031771A CN 102590831 B CN102590831 B CN 102590831B
Authority
CN
China
Prior art keywords
satellite
carrier
coordinate system
angle
carrier coordinate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2012100317718A
Other languages
English (en)
Other versions
CN102590831A (zh
Inventor
刘学维
郑瑞锋
童长海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huali Zhixin (Chengdu) integrated circuit Co., Ltd
Original Assignee
Beijing HWA Create Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing HWA Create Co Ltd filed Critical Beijing HWA Create Co Ltd
Priority to CN2012100317718A priority Critical patent/CN102590831B/zh
Publication of CN102590831A publication Critical patent/CN102590831A/zh
Application granted granted Critical
Publication of CN102590831B publication Critical patent/CN102590831B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种在载体旋转条件下判断卫星可见性的方法和装置,该方法主要包括:建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。本发明实施例通过建立载体坐标系,并将卫星在地平坐标系中的位置转换到载体坐标系,并计算卫星仰角,从而可以在载体旋转的情况下根据卫星仰角有效地判断卫星可见性。

Description

载体旋转条件下判断卫星可见性的方法和装置
技术领域
本发明涉及通信技术领域,尤其涉及一种载体旋转条件下判断卫星可见性的方法和装置。
背景技术
继GPS(Global Positioning System,全球定位系统)后,各国致力于研制自己的卫星导航系统,我国在北斗一代卫星导航系统的基础上,开创了北斗二号卫星导航系统的新局面,并取得了阶段性成果。然而,北斗卫星导航系统的发展却面临着众多挑战。为了开展北斗制导技术研究,旋转体制的设备在高速飞行的同时还进行高速旋转运动,导致旋转体制的设备的可见星数目、接收卫星导航信号的幅度、多普勒、载波相位快速变化,因此,对接收机的捕获、跟踪性能提出了更高的要求。
现有技术中还没有在载体旋转的情况下,有效地判断卫星可见的方法。
发明内容
本发明的实施例提供了一种载体旋转条件下判断卫星可见性的方法和装置,以实现在载体旋转的情况下有效地判断卫星可见性。
一种在载体旋转条件下判断卫星可见性的方法,包括:
建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;
根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。
一在载体旋转条件下判断卫星可见性的装置,包括:
卫星位置获取模块,用于建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;
卫星可见性判断模块,用于根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。
由上述本发明的实施例提供的技术方案可以看出,本发明实施例通过建立载体坐标系,并将卫星在地平坐标系中的位置转换到载体坐标系,并计算卫星仰角,从而可以在载体旋转的情况下根据卫星仰角有效地判断卫星可见性。
附图说明
图1为本发明实施例一提供的一种在载体旋转条件下判断卫星可见性的方法的处理流程图;
图2为本发明实施例二提供的一种在载体旋转条件下判断卫星可见性的装置的具体结构图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图,以对本发明实施例进行清楚、完整地描述人员。
实施例一
本发明实施例通过建立载体坐标系,并将卫星在地平坐标系中的位置转换到载体坐标系,从而计算卫星仰角,根据卫星仰角判断卫星可见性。
该实施例提供的在载体旋转条件下判断卫星可见性的方法的处理流程如图1所示,包括如下的处理步骤:
步骤11、建立载体坐标系。
空间中的载体的高速旋转运动会造成载体对接收天线的遮挡,造成卫星可见性变化,在载体旋转的同时,载体的接收信号的多普勒和信号功率也随之发生变化。
本发明实施例建立的载体坐标系的原点O为旋转载体的天线相位中心,Y轴为载体轴线方向,Z轴为载体天线所在方向,X轴与Y、Z轴构成右手直角坐标系。载体航向角、俯仰角、横滚角分别为载体坐标系三个坐标轴与地平坐标系三轴的夹角。理论上认为天线辐射的信号以天线相位中心为圆心,向外辐射。
步骤12、计算载体旋转角度。
设载体沿轴线方向顺时针转动,自旋角速度为ω,角速度一阶变率为ω′,角速度二阶变率为ω″,仿真时间间隔为Δt,则载体旋转角度,即横滚角α计算公式如下:
α = ωΔt + 1 2 ω ′ Δ t 2 + 1 6 ω ′ ′ Δ t 3
步骤13、基于载体旋转角度计算载体坐标系下的卫星位置。
设地平坐标系下卫星位置为[e n u]T,上述上标T代表转置,载体坐标系原点在地平坐标系下的坐标为[eo no uo]T,则基于上述载体的横滚角α,卫星在上述载体坐标系下的位置[x y z]T的计算公式如下:
x y z = R [ e n u - e o n o u o ] = cos α 0 - sin α 0 1 0 sin α 0 cos α [ e n u - e o n o u o ]
其中,R为地平坐标系到载体坐标系的转换矩阵,是由地平坐标系和载体坐标系以及他们之间的旋转角度α得到的,将地平坐标系绕Y轴旋转α得到载体坐标系,其转换矩阵即为R。
步骤14:基于载体坐标系下的卫星位置计算卫星仰角,根据卫星仰角和设定的卫星屏蔽角判断卫星可见性。
卫星在载体坐标系下的仰角计算公式如下。
θ = arctan z x 2 + y 2
设卫星屏蔽角为
Figure BDA0000135409530000042
则当
Figure BDA0000135409530000043
时,判定卫星可见,当时卫,则判定卫星不可见。
实施例二
该实施例提供了一种在载体旋转条件下判断卫星可见性的装置,其具体结构如图2所示,包括如下的模块:
卫星位置获取模块21,用于建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;
卫星可见性判断模块22,用于根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。
具体的,所述的卫星位置获取模块21可以包括:
载体旋转角度计算模块211,用于建立载体坐标系,所述载体坐标系的原点O为旋转载体的天线相位中心,Y轴为载体轴线方向,Z轴为载体天线所在方向,X轴与Y、Z轴构成右手直角坐标系,所述载体坐标系与地平坐标系之间的转动角度即为载体航向角、俯仰角和横滚角;
设载体沿轴线方向顺时针转动,自旋角速度为ω,所述ω的一阶变率为ω′,所述ω的二阶变率为ω″,仿真时间间隔为Δt,所述载体的旋转角度,即横滚角α计算公式如下:
α = ωΔt + 1 2 ω ′ Δ t 2 + 1 6 ω ′ ′ Δ t 3 ;
卫星位置计算模块212,用于设地平坐标系下卫星位置为[e n u]T,载体坐标系原点在地平坐标系下的坐标为[eo no uo]T,则基于所述载体的横滚角α,卫星在所述载体坐标系下的位置[x y z]T的计算公式如下:
x y z = R [ e n u - e o n o u o ] = cos α 0 - sin α 0 1 0 sin α 0 cos α [ e n u - e o n o u o ]
其中,所述R为地平坐标系到载体坐标系的转换矩阵。
具体的,所述的卫星可见性判断模块22,还用于设卫星在所述载体坐标系下的位置为[x y z]T,所述卫星在所述载体坐标系下的仰角计算公式如下:
θ = arctan z x 2 + y 2
设卫星屏蔽角为
Figure BDA0000135409530000053
则当
Figure BDA0000135409530000054
时,则判定所述卫星可见,当
Figure BDA0000135409530000055
时,则判定所述卫星不可见。
应用本发明实施例的装置判断卫星可见性的方法的具体处理过程与前述方法实施例类似,此处不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
综上所述,本发明实施例通过建立载体坐标系,并将卫星在地平坐标系中的位置转换到载体坐标系,并计算卫星仰角,从而可以在载体旋转的情况下根据卫星仰角有效地判断卫星可见性。
本发明实施例提出了一种载体旋转条件下卫星可见性的判断方法,有效的解决了由于飞机或导弹在姿态变化时,信号被自身遮挡引起的卫星可见性变化问题。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (5)

1.一种在载体旋转条件下判断卫星可见性的方法,其特征在于,包括:
建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;
其中,所述的建立载体坐标系,计算载体旋转角度,具体包括:建立的载体坐标系的原点O为旋转载体的天线相位中心,Y轴为载体轴线方向,Z轴为载体天线所在方向,X轴与Y、Z轴构成右手直角坐标系,所述载体坐标系与地平坐标系之间的转动角度即为载体航向角、俯仰角和横滚角;设载体沿轴线方向顺时针转动,自旋角速度为ω,所述ω的一阶变率为ω′,所述ω的二阶变率为ω″,仿真时间间隔为Δt,所述载体的旋转角度,即横滚角α计算公式如下: α = ωΔt + 1 2 ω ′ Δt 2 + 1 6 ω ′ ′ Δ t 3 ;
根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。
2.根据权利要求1所述的在载体旋转条件下判断卫星可见性的方法,其特征在于,所述的根据所述载体旋转角度计算所述载体坐标系下的卫星位置,包括:
设地平坐标系下卫星位置为[e n u]T,载体坐标系原点在地平坐标系下的坐标为[eo no uo]T,则基于所述载体的横滚角α,卫星在所述载体坐标系下的位置[x y z]T的计算公式如下:
x y z = R [ e n u - e o n o u o ] = cos α 0 - sin 0 1 0 sin α 0 cos α [ e n u - e o n o u o ]
其中,所述R为地平坐标系到载体坐标系的转换矩阵。
3.根据权利要求1或2所述的在载体旋转条件下判断卫星可见性的方法,其特征在于,所述的根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性,包括:
设卫星在所述载体坐标系下的位置为[x y z]T,所述卫星在所述载体坐标系下的仰角计算公式如下:
θ = arctan z x 2 + y 2
设卫星屏蔽角为
Figure FDA0000363335110000022
,则当
Figure FDA0000363335110000023
时,则判定所述卫星可见,当
Figure FDA0000363335110000024
时,则判定所述卫星不可见。
4.一在载体旋转条件下判断卫星可见性的装置,其特征在于,包括:
卫星位置获取模块,用于建立载体坐标系,计算载体旋转角度,根据所述载体旋转角度计算所述载体坐标系下的卫星位置;
卫星可见性判断模块,用于根据所述载体坐标系下的卫星位置计算卫星仰角,将所述卫星仰角和设定的卫星屏蔽角进行比较,根据比较结果判断卫星可见性。
其中,所述的卫星位置获取模块包括:
载体旋转角度计算模块,用于建立载体坐标系,所述载体坐标系的原点O为旋转载体的天线相位中心,Y轴为载体轴线方向,Z轴为载体天线所在方向,X轴与Y、Z轴构成右手直角坐标系,所述载体坐标系与地平坐标系之间的转动角度即为载体航向角、俯仰角和横滚角;
设载体沿轴线方向顺时针转动,自旋角速度为ω,所述ω的一阶变率为ω′,所述ω的二阶变率为ω″,仿真时间间隔为Δt,所述载体的旋转角度,即横滚角α计算公式如下:
α = ωΔt + 1 2 ω ′ Δt 2 + 1 6 ω ′ ′ Δ t 3 ;
卫星位置计算模块,用于设地平坐标系下卫星位置为[e n u]T,载体坐标系原点在地平坐标系下的坐标为[eo no uo]T,则基于所述载体的横滚角α,卫星在所述载体坐标系下的位置[xyz]T的计算公式如下:
x y z = R [ e n u - e o n o u o ] = cos α 0 - sin 0 1 0 sin α 0 cos α [ e n u - e o n o u o ]
其中,所述R为地平坐标系到载体坐标系的转换矩阵。
5.根据权利要求4所述的在载体旋转条件下判断卫星可见性的方法,其特征在于,所述的卫星可见性判断模块,还用于设卫星在所述载体坐标系下的位置为[x y z]T,所述卫星在所述载体坐标系下的仰角计算公式如下:
θ = arctan z x 2 + y 2
设卫星屏蔽角为
Figure FDA0000363335110000033
,则当时,则判定所述卫星可见,当
Figure FDA0000363335110000035
时,则判定所述卫星不可见。
CN2012100317718A 2012-02-13 2012-02-13 载体旋转条件下判断卫星可见性的方法和装置 Active CN102590831B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100317718A CN102590831B (zh) 2012-02-13 2012-02-13 载体旋转条件下判断卫星可见性的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100317718A CN102590831B (zh) 2012-02-13 2012-02-13 载体旋转条件下判断卫星可见性的方法和装置

Publications (2)

Publication Number Publication Date
CN102590831A CN102590831A (zh) 2012-07-18
CN102590831B true CN102590831B (zh) 2013-11-13

Family

ID=46479742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100317718A Active CN102590831B (zh) 2012-02-13 2012-02-13 载体旋转条件下判断卫星可见性的方法和装置

Country Status (1)

Country Link
CN (1) CN102590831B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104180810B (zh) * 2014-08-26 2016-08-31 上海微小卫星工程中心 一种卫星测控出入境判定方法及装置
CN106324646B (zh) * 2016-08-23 2018-11-16 西安电子科技大学 基于卫星导航的导弹测姿方法
CN106446446B (zh) * 2016-10-14 2020-02-18 上海微小卫星工程中心 星上高精度自主出入境方法
CN106970398B (zh) * 2017-03-27 2020-05-19 中国电建集团西北勘测设计研究院有限公司 顾及卫星遮挡条件的卫星可见性分析及星历预报方法
CN107255823A (zh) * 2017-06-20 2017-10-17 北京遥测技术研究所 一种载体多天线可用导航卫星信号的产生方法
CN107677266B (zh) * 2017-09-03 2023-06-20 陈应天 基于自旋-仰角跟踪理论的星光导航系统及其解算方法
CN109813306A (zh) * 2019-01-17 2019-05-28 李子月 一种无人车规划轨迹卫星定位数据可信度计算方法
CN111830536A (zh) * 2020-08-10 2020-10-27 南京林业大学 一种结合地形影响的卫星可见性判断方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261317A (zh) * 2008-04-25 2008-09-10 浙江大学 高动态多模式卫星导航信号源模拟方法及其装置
CN101493511A (zh) * 2008-01-23 2009-07-29 中国科学院微电子研究所 一种全球定位系统卫星搜索调度方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6556832B1 (en) * 2000-02-04 2003-04-29 Qualcomm Incorporated Method and apparatus for evaluation of position location performance

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101493511A (zh) * 2008-01-23 2009-07-29 中国科学院微电子研究所 一种全球定位系统卫星搜索调度方法
CN101261317A (zh) * 2008-04-25 2008-09-10 浙江大学 高动态多模式卫星导航信号源模拟方法及其装置

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
刘旭东 等.旋转多天线条件下导航解算方法.《中国惯性技术学报》.2010,第18卷(第4期),466-470,475.
刘旭东 等.载体旋转条件下GPS中频信号生成方法.《飞行器测控学报》.2009,第28卷(第6期),91页左栏第1行至92页左栏第18行.
旋转多天线条件下导航解算方法;刘旭东 等;《中国惯性技术学报》;20100831;第18卷(第4期);466-470,475 *
旋转条件GPS接收信号频率和相位变化分析;申强 等;《北京理工大学学报》;20090131;第29卷(第1期);35-37 *
申强 等.旋转条件GPS接收信号频率和相位变化分析.《北京理工大学学报》.2009,第29卷(第1期),35-37.
载体旋转条件下GPS中频信号生成方法;刘旭东 等;《飞行器测控学报》;20091231;第28卷(第6期);91页左栏第1行至92页左栏第18行 *

Also Published As

Publication number Publication date
CN102590831A (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102590831B (zh) 载体旋转条件下判断卫星可见性的方法和装置
US10355351B2 (en) Antenna array pointing direction estimation and control
CN103557876B (zh) 一种用于天线跟踪稳定平台的捷联惯导初始对准方法
CN103335649B (zh) 一种惯性导航系统极区导航参数解算方法
CN103972654B (zh) 直升机旋翼遮挡下的动中通天线卫星对星跟踪装置
CN106410410B (zh) 一种具有物理水平平台的vsat天线系统卫星捕获跟踪方法
CN102736091B (zh) 星球表面广域探测的卫星导航方法及系统
CN103900608B (zh) 一种基于四元数ckf的低精度惯导初始对准方法
CN102707080B (zh) 一种星敏感器模拟捷联惯导陀螺的方法
CN102157790A (zh) 一种用于动中通的天线跟踪系统
CN107607947A (zh) 基于卡尔曼滤波的星载雷达成像参数在线估计方法
CN103744052B (zh) 一种应用于空中目标定位的双星测时差测向方法及装置
CN104374388A (zh) 一种基于偏振光传感器的航姿测定方法
CN101629826A (zh) 基于单轴旋转的光纤陀螺捷联惯性导航系统粗对准方法
CN103675760B (zh) 一种星载地球同步轨道合成孔径雷达姿态导引方法
CN102117952A (zh) 四轴结构天线稳定框架及控制方法
CN104049241B (zh) 目标位置坐标未知的双基地合成孔径雷达的空间同步方法
CN103743413A (zh) 倾斜状态下调制寻北仪安装误差在线估计与寻北误差补偿方法
CN105116430A (zh) 用于动中通的伪航向的基于卡尔曼滤波的海泊态搜星方法
CN109765530A (zh) 一种运动平台雷达波束解耦方法
CN111121769B (zh) 机械扫描便携式合作目标对空监视设备及自主校北方法
CN110187400B (zh) 基于航向跟踪的海空重力扰动水平分量测量误差调制方法
CN206202675U (zh) 无人飞行器
CN105115519A (zh) 应用于动中通系统的惯导系统初始对准方法
CN108981703A (zh) 火星着陆器联合位置估计方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20200716

Address after: 610000 Sichuan city of Chengdu province Tianfu Zheng Xing Shun Shing Street No. 172

Patentee after: Huali Zhixin (Chengdu) integrated circuit Co., Ltd

Address before: 100094 No. 8, No. 18, West Wang Road, Haidian District, Beijing

Patentee before: HWA CREATE Corp.,Ltd.

TR01 Transfer of patent right