CN115622594A - 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 - Google Patents
一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 Download PDFInfo
- Publication number
- CN115622594A CN115622594A CN202211232668.XA CN202211232668A CN115622594A CN 115622594 A CN115622594 A CN 115622594A CN 202211232668 A CN202211232668 A CN 202211232668A CN 115622594 A CN115622594 A CN 115622594A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- satellite
- sensor
- edge
- axis
- 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.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 title claims abstract description 16
- 239000013598 vector Substances 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 39
- 230000009466 transformation Effects 0.000 claims description 16
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 240000001549 Ipomoea eriocarpa Species 0.000 description 1
- 235000005146 Ipomoea eriocarpa Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0408—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas using two or more beams, i.e. beam diversity
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/18502—Airborne stations
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18515—Transmission equipment in satellites or space-based relays
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Radio Relay Systems (AREA)
Abstract
本发明公开了一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,涉及卫星通信仿真技术领域。本发明中从传感器的视角出发,根据卫星传感器姿态和卫星姿态准确的计算出卫星通信点波束的对地覆盖边缘,本发明能够根据传感器的指向建立波束边缘单位观测矢量,通过一系列的坐标系转换后联合直线方程和地球椭球体方程快速计算得到卫星通信点波束覆盖边缘的空间直角坐标,本发明的计算方法可接入实时的卫星遥测数据,根据卫星的天线指向方位角和俯仰角以及卫星姿态数据直接计算出卫星点波束边缘与地球交点的经纬度。
Description
技术领域
本发明涉及卫星通信仿真技术领域,尤其涉及一种基于传感器视角的卫星通信点波束覆盖边缘计算方法。
背景技术
卫星通信点波束具有覆盖范围更小、功率更强、抗干扰能力强等特点,可以实现较高的数据传输速率,其被广泛应用于通信卫星的通信波束设计。卫星通信点波束对地覆盖区域是卫星传感器形成的锥状波束与地球表面相交的区域,其覆盖位置和大小与卫星传感器波束半角、传感器姿态和卫星姿态相关。
在现有的点波束覆盖区域计算方法文献(李德治,窦朝晖,吕波等.卫星点波束覆盖区域算法研究[J].载人航天,2009,15(4))中是已知波束指向中心点求解波束覆盖边缘点,该计算方法没有考虑传感器姿态和卫星姿态的影响。在文献(薛永娇,王德民,齐春东等.载荷侧摆条件下卫星点波束覆盖区域算法研究[J].电子技术应用,2015,41(3))中建立以星下点为原点的站心坐标系,利用站心坐标系中波束的旋转角确定唯一的圆锥母线作为波束中心线,进而得到卫星波束中心与地球的交点,最终确定卫星覆盖范围,该计算方法可计算得到波束中心点和覆盖边缘点,同样没有考虑卫星传感器姿态和卫星姿态的影响。在专利(中国人民解放军海军工程大学.卫星通信点波束覆盖区域计算方法:中国,104348539[P].2015-02-11)中采用向量积公式求解覆盖边缘点,该方法计算时需按经纬度将地球划分为若干个点,需遍历卫星可视范围内的所有点进行计算,步进越小计算量越大,并且没有考虑卫星传感器姿态和卫星姿态的影响。
鉴于以上情况,本专利提出一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,将传感器坐标系中的波束覆盖边缘观测矢量转换为地心地固坐标系中的矢量进行计算,考虑了卫星传感器姿态和卫星姿态的影响,能够实时准确的计算出点波束的覆盖边缘区域,在卫星点波束覆盖特性分析中具有较高的使用价值。
发明内容
本发明的目的在于:为了解决上述问题,而提出的一种基于传感器视角的卫星通信点波束覆盖边缘计算方法。
为了实现上述目的,本发明采用了如下技术方案:
一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,包括以下步骤:
S1:建立传感器波束边缘单位观测矢量;以波束指向建立坐标系,坐标原点定义为传感器中心,X轴为传感器波束指向,Y轴垂直于X轴,向上为正,Z轴与X、Y轴所成平面符合右手定则,传感器中心原点为O1,M为波束边缘任意一点,对矢量O1M进行归一化得到波束边缘单位观测矢量Vunit;
S2:坐标系转换;进行波束边缘点坐标计算之前需将波束边缘单位观测矢量Vunit转换至地心地固坐标系,在同一个坐标系下进行波束边缘点坐标计算,坐标系包括波束坐标系、传感器坐标系、卫星本体坐标系、卫星轨道坐标系、地心惯性坐标系和地心地固坐标系;
S3:波束边缘点坐标计算;经过坐标转换得到地心地固坐标系下的波束边缘单位观测矢量:
Vecr=ReiRioRobRbsRsdVunit
波束边缘单位观测矢量经过卫星所在的位置,且波束边缘单位观测矢量Vecr和卫星位置矢量Secr为已知量:
Vecr=[v1,v2,v3]T
Secr=[s1,s2,s3]T
与波束边缘单位观测矢量共线的直线上一点的坐标为(X,Y,Z),与波束边缘单位观测矢量共线的直线方程为:
(X,Y,Z)为地固坐标系下地球上任一点的坐标,则地球椭球体方程式为:
其中a为地球的长半轴,b为地球的短半轴,联合上述直线方程和地球椭球体方程式可计算出与波束边缘单位观测矢量共线的直线与地球表面交点的空间直角坐标(X,Y,Z),将单位观测矢量Vunit绕X轴按步进旋转N次求解出卫星通信点波束覆盖边缘的空间直角坐标,最后将空间直角坐标(X,Y,Z)转换成大地坐标系下的经纬度和高程以便于绘图显示。
优选地,卫星本体坐标系至卫星轨道坐标系的转换矩阵Rob由当前的卫星姿态得到,卫星滚动角roll、俯仰角pitch,偏航角yaw构成欧拉旋转矩阵为Rz(yaw)Rx(roll)Ry(pitch)。
优选地,卫星轨道坐标系到地心惯性坐标系的转换矩阵Rio由当前卫星星历计算得到,卫星星历包括卫星位置矢量和速度矢量。
优选地,地心惯性坐标系到地心地固坐标系的转换矩阵Rei由极移矩阵EP、地球旋转矩阵ER、章动矩阵NR和岁差矩阵PR共同计算得到。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1、本发明能够根据传感器的指向建立波束边缘单位观测矢量,通过一系列的坐标系转换后联合直线方程和地球椭球体方程快速计算得到卫星通信点波束覆盖边缘的空间直角坐标。
2、本发明考虑了卫星传感器姿态和卫星姿态的影响,计算精度高,在卫星点波束指向调整规划和卫星点波束实时覆盖区域展示方面具有良好的应用前景。
3、本发明的计算方法可接入实时的卫星遥测数据,根据卫星的天线指向方位角和俯仰角以及卫星姿态数据直接计算出卫星点波束边缘与地球交点的经纬度。
附图说明
图1示出了根据本发明实施例提供的基于传感器视角的卫星通信点波束覆盖边缘示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1,本发明提供一种技术方案:
一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,包括以下步骤:
S1:建立传感器波束边缘单位观测矢量;以波束指向建立坐标系,坐标原点定义为传感器中心,X轴为传感器波束指向,Y轴垂直于X轴,向上为正,Z轴与X、Y轴所成平面符合右手定则,传感器中心原点为O1,M为波束边缘任意一点,对矢量O1M进行归一化得到波束边缘单位观测矢量Vunit,V为传感器视角的波束边缘观测矢量V=[v1,v2,v3]T,设|O1C|=d,圆锥截面圆半径|CM|=r,则O1M=[d,r,0]T,观测矢量V表达式如下:
归一化:
S2:坐标系转换;进行波束边缘点坐标计算之前需将波束边缘单位观测矢量Vunit转换至地心地固坐标系,在同一个坐标系下进行波束边缘点坐标计算,坐标系包括波束坐标系、传感器坐标系、卫星本体坐标系、卫星轨道坐标系、地心惯性坐标系和地心地固坐标系;
S3:波束边缘点坐标计算;经过坐标转换得到地心地固坐标系下的波束边缘单位观测矢量:
Vecr=ReiRioRobRbsRsdVunit
波束边缘单位观测矢量经过卫星所在的位置,且波束边缘单位观测矢量Vecr和卫星位置矢量Secr为已知量:
Vecr=[v1,v2,v3]T
Secr=[s1,s2,s3]T
与波束边缘单位观测矢量共线的直线上一点的坐标为(X,Y,Z),与波束边缘单位观测矢量共线的直线方程为:
(X,Y,Z)为地固坐标系下地球上任一点的坐标,则地球椭球体方程式为:
其中a为地球的长半轴,b为地球的短半轴,联合上述直线方程和地球椭球体方程式可计算出与波束边缘单位观测矢量共线的直线与地球表面交点的空间直角坐标(X,Y,Z),取其中离卫星较近的交点,将单位观测矢量Vunit绕X轴按步进旋转N次求解出卫星通信点波束覆盖边缘的空间直角坐标,最后将空间直角坐标(X,Y,Z)转换成大地坐标系下的经纬度和高程L、B、H。
其中a为地球椭球的长半径,e为椭球第一偏心率。
具体的,如图1所示,卫星本体坐标系至卫星轨道坐标系的转换矩阵Rob由当前的卫星姿态得到,卫星滚动角roll、俯仰角pitch,偏航角yaw构成欧拉旋转矩阵为Rz(yaw)Rx(roll)Ry(pitch)。
Rob=Rz(yaw)Rx(roll)Ry(pitch)
具体的,如图1所示,卫星轨道坐标系到地心惯性坐标系的转换矩阵Rio由当前卫星星历计算得到,卫星星历包括卫星位置矢量和速度矢量,卫星位置矢量和速度矢量由轨道六根数推算得到。
具体的,如图1所示,地心惯性坐标系到地心地固坐标系的转换矩阵Rei由极移矩阵EP、地球旋转矩阵ER、章动矩阵NR和岁差矩阵PR共同计算得到,其中地球旋转矩阵ER由格林尼治恒星时角GAST计算得出。
Rei=EP·ER·NR·PR
本发明中从传感器的视角出发,根据卫星传感器姿态和卫星姿态准确的计算出卫星通信点波束的对地覆盖边缘,本发明能够根据传感器的指向建立波束边缘单位观测矢量,通过一系列的坐标系转换后联合直线方程和地球椭球体方程快速计算得到卫星通信点波束覆盖边缘的空间直角坐标,本发明的计算方法可接入实时的卫星遥测数据,根据卫星的天线指向方位角和俯仰角以及卫星姿态数据直接计算出卫星点波束边缘与地球交点的经纬度。
实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (6)
1.一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,其特征在于,包括以下步骤:
S1:建立传感器波束边缘单位观测矢量;以波束指向建立坐标系,坐标原点定义为传感器中心,X轴为传感器波束指向,Y轴垂直于X轴,向上为正,Z轴与X、Y轴所成平面符合右手定则,传感器中心原点为O1,M为波束边缘任意一点,对矢量O1M进行归一化得到波束边缘单位观测矢量Vunit;
S2:坐标系转换;进行波束边缘点坐标计算之前需将波束边缘单位观测矢量Vunit转换至地心地固坐标系,在同一个坐标系下进行波束边缘点坐标计算,坐标系包括波束坐标系、传感器坐标系、卫星本体坐标系、卫星轨道坐标系、地心惯性坐标系和地心地固坐标系;
S3:波束边缘点坐标计算;经过坐标转换得到地心地固坐标系下的波束边缘单位观测矢量:
Vecr=ReiRioRobRbsRsdVunit
波束边缘单位观测矢量经过卫星所在的位置,且波束边缘单位观测矢量Vecr和卫星位置矢量Secr为已知量:
Vecr=[v1,v2,v3]T
Secr=[s1,s2,s3]T
与波束边缘单位观测矢量共线的直线上一点的坐标为(X,Y,Z),与波束边缘单位观测矢量共线的直线方程为:
(X,Y,Z)为地固坐标系下地球上任一点的坐标,则地球椭球体方程式为:
4.根据权利要求1所述的一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,其特征在于,卫星本体坐标系至卫星轨道坐标系的转换矩阵Rob由当前的卫星姿态得到,卫星滚动角roll、俯仰角pitch,偏航角yaw构成欧拉旋转矩阵为Rz(yaw)Rx(roll)Ry(pitch)。
5.根据权利要求1所述的一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,其特征在于,卫星轨道坐标系到地心惯性坐标系的转换矩阵Rio由当前卫星星历计算得到,卫星星历包括卫星位置矢量和速度矢量。
6.根据权利要求1所述的一种基于传感器视角的卫星通信点波束覆盖边缘计算方法,其特征在于,地心惯性坐标系到地心地固坐标系的转换矩阵Rei由极移矩阵EP、地球旋转矩阵ER、章动矩阵NR和岁差矩阵PR共同计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211232668.XA CN115622594B (zh) | 2022-10-10 | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211232668.XA CN115622594B (zh) | 2022-10-10 | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115622594A true CN115622594A (zh) | 2023-01-17 |
CN115622594B CN115622594B (zh) | 2024-07-02 |
Family
ID=
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116975504A (zh) * | 2023-09-22 | 2023-10-31 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819019A (zh) * | 2012-07-20 | 2012-12-12 | 西安空间无线电技术研究所 | 一种卫星波束与地球交点坐标的确定方法 |
CN104348539A (zh) * | 2014-07-23 | 2015-02-11 | 中国人民解放军海军工程大学 | 卫星通信点波束覆盖区域计算方法 |
CN106197425A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于卫星姿态角的地面目标点位置的计算方法 |
US20170288769A1 (en) * | 2010-05-02 | 2017-10-05 | Viasat, Inc. | Flexible beamforming for satellite communications |
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170288769A1 (en) * | 2010-05-02 | 2017-10-05 | Viasat, Inc. | Flexible beamforming for satellite communications |
CN102819019A (zh) * | 2012-07-20 | 2012-12-12 | 西安空间无线电技术研究所 | 一种卫星波束与地球交点坐标的确定方法 |
CN104348539A (zh) * | 2014-07-23 | 2015-02-11 | 中国人民解放军海军工程大学 | 卫星通信点波束覆盖区域计算方法 |
CN106197425A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于卫星姿态角的地面目标点位置的计算方法 |
Non-Patent Citations (1)
Title |
---|
张大伟 等: "大椭圆轨道卫星姿态基准坐标系与运动学研究", 上海航天, no. 02, 25 April 2017 (2017-04-25) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116975504A (zh) * | 2023-09-22 | 2023-10-31 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
CN116975504B (zh) * | 2023-09-22 | 2023-12-15 | 中科星图测控技术股份有限公司 | 一种用于卫星侦察覆盖区域目标的快速计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108508918B (zh) | 一种静轨遥感卫星数传天线高精度实时对地指向控制方法 | |
CN108427427B (zh) | 一种航天器对地表定向目标姿态角计算方法 | |
CN111142099B (zh) | 解决球面相控阵天线跟踪过顶盲捕目标的方法 | |
EP1772742B1 (en) | Correction of the distance between phase centres of two directional antenneas of a navigational satellite | |
CN109781059B (zh) | 星载点波束天线对地指向精度评估系统 | |
CN107450582A (zh) | 一种基于星上实时规划的相控阵数传引导控制方法 | |
CN112130590A (zh) | 一种基于瞬时惯性系下速度补偿的星载天线对地指向确定方法 | |
CN103675760A (zh) | 一种星载地球同步轨道合成孔径雷达姿态导引方法 | |
CN114142234A (zh) | 相控阵天线与卫星的对星方法、系统、电子设备及介质 | |
CN103697893B (zh) | 利用大气偏振光的三维定姿方法 | |
CN112255606A (zh) | 基于单反射面天线的Geo-SAR卫星正侧视成像姿态角计算方法 | |
CN111924133A (zh) | 适应空中信号高精度三维定位的编队构形设计方法及系统 | |
CN111121769A (zh) | 机械扫描便携式合作目标对空监视设备及自主校北方法 | |
CN113386979B (zh) | 一种自适应规避太阳的数传姿态规划方法 | |
EP2080033A1 (de) | Navigationseinrichtung und verfahren zum ermitteln von orientierungen | |
CN106250684B (zh) | 基于地固系数据的卫星过境时间快速计算方法 | |
CN115622594A (zh) | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 | |
CN115622594B (zh) | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 | |
US6283415B1 (en) | Simplified yaw steering method for satellite antenna beam control | |
CN115632699A (zh) | 一种目标天线的指向控制系统、方法、设备及存储介质 | |
CN112607056B (zh) | 雷达卫星对目标观测自主启动触发方法和系统 | |
CN111983612B (zh) | 一种sar滑动聚束模式方位去斜方法 | |
CN115604851B (zh) | 一种基于固定截距的大规模卫星星座频率干扰规避方法 | |
CN118138107A (zh) | 一种异轨星间天线链路指向算法的构建方法 | |
CN111427382B (zh) | 一种复杂月面地形下基准变化的中继指向规划方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |