CN113433571B - 一种下行链路最差干扰快速估算方法 - Google Patents

一种下行链路最差干扰快速估算方法 Download PDF

Info

Publication number
CN113433571B
CN113433571B CN202110555538.9A CN202110555538A CN113433571B CN 113433571 B CN113433571 B CN 113433571B CN 202110555538 A CN202110555538 A CN 202110555538A CN 113433571 B CN113433571 B CN 113433571B
Authority
CN
China
Prior art keywords
satellite
gso
ngso
coordinate system
ngso satellite
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
CN202110555538.9A
Other languages
English (en)
Other versions
CN113433571A (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.)
CASIC Space Engineering Development Co Ltd
Original Assignee
CASIC Space Engineering Development 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 CASIC Space Engineering Development Co Ltd filed Critical CASIC Space Engineering Development Co Ltd
Priority to CN202110555538.9A priority Critical patent/CN113433571B/zh
Publication of CN113433571A publication Critical patent/CN113433571A/zh
Application granted granted Critical
Publication of CN113433571B publication Critical patent/CN113433571B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • 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/18519—Operations control, administration or maintenance
    • 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/19—Earth-synchronous 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/195—Non-synchronous stations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Signal Processing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种下行链路最差干扰快速估算方法，包括：S101：初始化NGSO卫星的纬度幅角；S102：计算NGSO卫星在J2000坐标系的三维位置矢量、赤经和J2000坐标系到轨道坐标系的坐标转换矩阵；S103：初始化GSO卫星的赤经；S104：计算GSO卫星在J2000坐标系的位置矢量；S105：计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量；S106：判断GSO卫星到NGSO卫星的延长线是否与地球相交；S107：计算与地球交点在J2000坐标系的矢量表达式；S108：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益；S109：计算与地球交点处GSO地球站的EPFD；S110：存储EPFD、对应NGSO卫星纬度幅角、GSO卫星赤经和地球表面交点的坐标；S111：更新GSO的赤经；S112：更新NGSO的纬度幅角。

Description

一种下行链路最差干扰快速估算方法
技术领域
本发明涉及下行链路最差干扰快速估算领域。更具体地，涉及一种下行链路最差干扰快速估算方法。
背景技术
随时随地的互联网接入是当前信息时代最重要的特征，近年来，国内外多家企业都提出了非对地静止卫星轨道(NGSO)卫星系统，计划为全球提供随时随地的互联网接入服务。NGSO卫星系统的主要用频包括Ku频段、Ka频段等，根据《无线电规则》，NGSO卫星系统的空间电台的发射，在对地静止卫星轨道(GSO)可视的地球表面任何点上产生的等效功率通量密度(EPFD)在给定的百分比时间内均不得超过给定的限值，并且给出了EPFD限值对应的地球站基准天线和基准带宽。
对于参考频段为17.8-18.6GHz、地球站基准天线参考天线直径为1m、参考辐射模式为ITU-R S.1428-1建议书、基准带宽为40kHz，NGSO卫星系统发射的EPFD不得超过的限值如表1所示，可以看出表1的第二列给出了几个“限值点”，第三列给出了每个“限值点”对应的不能超出的时间百分比，某NGSO卫星系统发射的EPFD超出相应值的时间百分比与表1的限值点关系的示意图如图1所示。第二列最后一行的EPFD限值为“-164”，对应的不超出时间百分比为“100％”，即限值“-164”是完全不能超过的。对于其他不同的频段、地球站参考天线直径或参考辐射模式，NGSO卫星系统发射的EPFD不得超过的限值的格式与表1的限值格式类似，即都存在一个完全不能超过的限值。
由于EPFD限值是以时间百分比的形式给出，所以在进行限值检查计算时，通常采用“时间切片”的方法完成，即设置一定的仿真总时长和仿真步长，分别计算不同仿真时刻的EPFD，并统计EPFD超出相应值的时间百分比。为了获得可信度高的仿真结果，仿真总时长应尽量长、仿真步长应尽量短，这带来的问题是仿真计算耗时长，无法快速确定系统间的干扰是否满足要求，从而影响系统方案的设计进度。
发明内容
本发明的一个目的在于提供一种下行链路最差干扰快速估算方法，用于快速确定NGSO卫星系统对GSO卫星系统的下行链路最差干扰，从而可以快速判断NGSO卫星系统对GSO卫星系统的下行链路干扰是否满足限值要求。
为达到上述目的，本发明采用下述技术方案：
一种下行链路最差干扰快速估算方法，包括以下步骤：
S101：初始化NGSO卫星的纬度幅角；
S102：计算NGSO卫星在J2000坐标系的三维位置矢量、赤经和J2000坐标系到轨道坐标系的坐标转换矩阵；
S103：初始化GSO卫星的赤经；
S104：计算GSO卫星在J2000坐标系的位置矢量；
S105：计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量；
S106：判断GSO卫星到NGSO卫星的延长线是否与地球表面相交，如果相交，转到S107，不相交转到S111；
S107：计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式；
S108：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益；
S109：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD；
S110：存储EPFD，以及对应NGSO卫星纬度幅角、GSO卫星赤经以及GSO卫星到NGSO卫星的延长线与地球表面的交点的坐标；
S111：更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算；
S112：更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算。
在一个具体示例中，所述初始化NGSO卫星的纬度幅角包括：
S1011：将NGSO卫星的纬度幅角取为u＝u0，其中，u0为NGSO卫星的初始纬度幅角u0＝0°。
在一个具体示例中，所述计算NGSO卫星在J2000坐标系的位置矢量、赤经以及J2000坐标系到轨道坐标系的坐标转换矩阵包括：
S1021：根据纬度幅角u和近地点幅角ω，计算真近点角：
θ＝u-ω
S1022：根据所述真近点角计算NGSO卫星轨道半径为：
Figure BDA0003077035430000031
其中，a是NGSO卫星半长轴，e是NGSO卫星偏心率；
S1023：NGSO卫星在J2000坐标系的位置矢量为：
Figure BDA0003077035430000032
其中，
Figure BDA0003077035430000033
为单位向量，计算公式为：
Figure BDA0003077035430000041
其中，Ω是NGSO卫星升交点赤经，i是NGSO卫星轨道倾角；
S1024：将NGSO卫星在J2000坐标系的位置矢量的三维分量记为xN、yN和zN，NGSO卫星在J2000坐标系的赤经为：
Figure BDA0003077035430000042
S1025：J2000坐标系到轨道坐标系的坐标转换矩阵为：
Figure BDA0003077035430000043
其中，Mx()、My()和Mz()分别为绕X、Y和Z轴的坐标转换矩阵计算公式为：
Figure BDA0003077035430000044
Figure BDA0003077035430000045
Figure BDA0003077035430000046
在一个具体示例中，所述初始化GSO卫星的赤经包括：
S1031：初始化GSO卫星的赤经：
LG＝LN-ΔLG
其中，ΔLG为GSO卫星赤经区间宽度的一半。
在一个具体示例中，所述计算GSO卫星在J2000坐标系的位置矢量包括：
S1041：GSO卫星在J2000坐标系的位置矢量rG为：
Figure BDA0003077035430000051
其中，Re是地球平均半径，BG为GSO卫星的赤纬，HG为GSO卫星的高度。
在一个具体示例中，所述计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量包括：
S1051：GSO卫星到NGSO卫星的延长线在J2000坐标系的单位矢量为：
Figure BDA0003077035430000052
S1052：GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量为：
Figure BDA0003077035430000053
在一个具体示例中，所述判断GSO卫星到NGSO卫星的延长线是否与地球表面相交包括：
S1061：NGSO卫星到地球的切线与轨道坐标系+Z轴方向的夹角为：
Figure BDA0003077035430000054
S1062：GSO卫星到NGSO卫星的延长线与轨道坐标系+Z轴方向的夹角为：
Figure BDA0003077035430000055
其中，zo是GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量do第三维分量；
S1063：当
Figure BDA0003077035430000056
时，则延长线与地球表面相交；否则，延长线不与地球表面相交。
在一个具体示例中，所述计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式包括：
S1071：GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式为：
rd＝rN+kdi
其中，系数k计算公式为：
Figure BDA0003077035430000061
在一个具体示例中，所述计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD包括：
S1091：NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地面的交点处GSO地球站的EPFD为：
EPFD＝Pt+Gt-10log10(4πd2)+Gr-Grmax
其中，Pt是NGSO卫星空间电台的发射功率，Gt为NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益，Gr是GSO地球站的接收增益，Grmax是GSO地球站的最大接收增益，在GSO地球站处于GSO卫星与NGSO卫星延长线方向的情况下，Gr＝Grmax，d是NGSO卫星空间电台与GSO卫星到NGSO卫星的延长线与地球表面的交点的距离，即d＝|rd-rN|。
在一个具体示例中，更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算包括：
S1111：更新GSO卫星的赤经：
首先计算LG1：
LG1＝LG+Ls
更新LG：
LG＝LG1
其中，Ls为GSO卫星赤经的更新步长；
S1112：判断LG≤LN+ΔLG，如果符合，转到S104，否则转到S112。
在一个具体示例中，所述更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算包括：
S1121：更新NGSO卫星的纬度幅角：
首先计算u1：
u1＝u+us
更新u：
u＝u1
其中，us为NGSO卫星纬度幅角的更新步长；
S1122：判断u≤360°，若是，转到S102，否则结束。
本发明的有益效果如下：
本发明的一种下行链路最差干扰快速估算方法，可以用于快速确定NGSO卫星系统对GSO卫星系统的下行链路最差干扰，从而可以快速判断NGSO卫星系统对GSO卫星系统的下行链路干扰是否满足限值要求。
本发明的一种下行链路最差干扰快速估算方法获得的最差场景的卫星和地球站位置，可以作为“时间切片”仿真分析的初值，从而可以为“时间切片”仿真分析的仿真初值选取提供依据。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为某NGSO卫星系统发射的EPFD超出相应值的时间百分比与限值点关系的示意图；
图2为最差干扰场景对应GSO卫星、NGSO卫星以及地球站的空间几何关系示意图；
图3为本发明实施例一种下行链路最差干扰快速估算方法的流程图；
图4为最差场景(即NGSO卫星处于GSO卫星与GSO地球站连线上)时，不采取干扰规避措施的情况下EPFD随纬度幅角的变化；
图5为最差场景(即NGSO卫星处于GSO卫星与GSO地球站连线上)时采取干扰规避措施的情况下EPFD随纬度幅角的变化。
具体实施方式
为使本发明的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
NGSO卫星的运动轨迹可以由六个自由变量定义，分别为半长轴a、偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω和纬度幅角u，其中前五个变量是不变量或慢变量，最后一个变量纬度幅角是快变量，从而NGSO卫星位置的变化主要取决于纬度幅角；同时，GSO卫星的位置可以采用定点经度描述。
基于此，在寻找最差干扰时，将NGSO卫星的纬度幅角和GSO卫星经度作为变量分别确定NGSO卫星的位置和GSO卫星的位置，根据最差场景分析，GSO地球站位于GSO卫星与NGSO卫星的延长线与地球表面交点，从而确定了GSO地球站的位置。在确定GSO地球站位置后，从而可以确定NGSO卫星与GSO地球站的距离、NGSO卫星发射天线在NGSO卫星到GSO地球站的方向的发射增益，进而可以计算NGSO卫星在地球表面NGSO地球站处产生的EPFD。同时，采用最差干扰快速估算方法获得的NGSO卫星位置、GSO卫星位置以及GSO地球站位置，可以为作为进一步的“时间切片”仿真分析的初始值，从而可以为“时间切片”仿真分析的初值选取提供依据。
本发明一种下行链路最差干扰快速估算方法，涉及非对地静止卫星轨道卫星系统NGSO对GSO卫星系统的下行链路最差干扰快速估算方法，该方法首先根据EPFD的计算公式的特征，分析确定了最差干扰场景；然后建立了最差干扰场景下，NGSO纬度幅角、GSO卫星经度取不同值时EPFD的计算方法，进而确定EPFD最大值。
根据EPFD限值格式，每种场景都存在一个“100％”不能超过的限值，即完全不能超过的值。如图1和表1所示，当NGSO卫星系统发射的EPFD最大值超过时间百分比为“100％”对应的限值时，NGSO卫星系统发射的EPFD不满足要求，从而不需要进行进一步判断是否满足其他限值，从而可以实现快速初步判断，因此下行链路最差干扰的快速计算具有重要意义。
Figure BDA0003077035430000091
表1
根据EPFD的计算公式，EPFD的大小主要由发射功率、发射天线增益、接收天线增益和传播距离决定。其中，发射功率为卫星系统的设计参数，通常为固定值或变化范围较小，并且由传播距离导致的EPFD改变量变化、发射偏离角导致的发射天线增益变化相对于接收偏离角导致的接收增益变化都较小，如图2所示，因此，下行链路的最差干扰通常出现在NGSO卫星处于GSO卫星与GSO地球站连线上。
本发明的一个实施例提供一种下行链路最差干扰快速估算方法，如图3所示，包括以下步骤：
S101：初始化NGSO卫星的纬度幅角；
S102：计算NGSO卫星在J2000坐标系的三维位置矢量、赤经和J2000坐标系到轨道坐标系的坐标转换矩阵；
S103：初始化GSO卫星的赤经；
S104：计算GSO卫星在J2000坐标系的位置矢量；
S105：计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量；
S106：判断GSO卫星到NGSO卫星的延长线是否与地球表面相交，如果相交，转到S107，不相交转到S111；
S107：计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式；
S108：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益；
S109：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD；
S110：存储EPFD，以及对应NGSO卫星纬度幅角、GSO卫星赤经以及GSO卫星到NGSO卫星的延长线与地球表面的交点的坐标；
S111：更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算；
S112：更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算。
在一个具体示例中，所述初始化NGSO卫星的纬度幅角包括：
S1011：将NGSO卫星的纬度幅角取为u＝u0，其中，u0为NGSO卫星的初始纬度幅角u0＝0°。
在一个具体示例中，所述计算NGSO卫星在J2000坐标系的位置矢量、赤经以及J2000坐标系到轨道坐标系的坐标转换矩阵包括：
S1021：根据纬度幅角u和近地点幅角ω，计算真近点角：
θ＝u-ω
S1022：根据所述真近点角计算NGSO卫星轨道半径为：
Figure BDA0003077035430000111
其中，a是NGSO卫星半长轴，e是NGSO卫星偏心率；
S1023：NGSO卫星在J2000坐标系的位置矢量为：
Figure BDA0003077035430000112
其中，
Figure BDA0003077035430000113
为单位向量，计算公式为：
Figure BDA0003077035430000114
其中，Ω是NGSO卫星升交点赤经，i是NGSO卫星轨道倾角；
S1024：将NGSO卫星在J2000坐标系的位置矢量的三维分量记为xN、yN和zN，NGSO卫星在J2000坐标系的赤经为：
Figure BDA0003077035430000121
S1025：J2000坐标系到轨道坐标系的坐标转换矩阵为：
Figure BDA0003077035430000122
其中，Mx()、My()和Mz()分别为绕X、Y和Z轴的坐标转换矩阵计算公式为：
Figure BDA0003077035430000123
Figure BDA0003077035430000124
Figure BDA0003077035430000125
在一个具体示例中，所述初始化GSO卫星的赤经包括：
S1031：初始化GSO卫星的赤经：
LG＝LN-ΔLG
其中，ΔLG为GSO卫星赤经区间宽度的一半。
在一个具体示例中，所述计算GSO卫星在J2000坐标系的位置矢量包括：
S1041：GSO卫星在J2000坐标系的位置矢量rG为：
Figure BDA0003077035430000126
其中，Re是地球平均半径，BG为GSO卫星的赤纬，HG为GSO卫星的高度。
在一个具体示例中，所述计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量包括：
S1051：GSO卫星到NGSO卫星的延长线在J2000坐标系的单位矢量为：
Figure BDA0003077035430000131
S1052：GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量为：
Figure BDA0003077035430000132
在一个具体示例中，所述判断GSO卫星到NGSO卫星的延长线是否与地球相交包括：
S1061：NGSO卫星到地球的切线与轨道坐标系+Z轴方向的夹角为：
Figure BDA0003077035430000133
S1062：GSO卫星到NGSO卫星的延长线与轨道坐标系+Z轴方向的夹角为：
Figure BDA0003077035430000134
其中，zo是GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量do第三维分量即在Z轴上的分量；
S1063：当
Figure BDA0003077035430000135
时，则延长线与地球表面相交；否则，延长线不与地球表面相交。
在一个具体示例中，所述计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式包括：
S1071：GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式为：
rd＝rN+kdi
其中，系数k计算公式为：
Figure BDA0003077035430000141
在一个具体实施例中，所述计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益包括：
S1081：NGSO卫星空间电台发射天线的天线方向图函数通常可以在轨道坐标系中给出，根据S1052给出的GSO卫星到NGSO卫星的延长线在轨道坐标系得单位矢量do，可以确定NGSO卫星空间电台发射天线在GSO卫星到NGSO卫星的延长线方向的发射增益Gt。
在一个具体示例中，所述计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD包括：
S1091：NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD为：
EPFD＝Pt+Gt-10log10(4πd2)+Gr-Grmax
其中，Pt是NGSO卫星空间电台的发射功率，通常为一给定的固定值，Gt为NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益，Gr是GSO地球站的接收增益，Grmax是GSO地球站的最大接收增益，由天线的特性决定，在GSO地球站处于GSO卫星与NGSO卫星延长线方向的情况下，Gr＝Grmax，d是NGSO卫星空间电台与GSO卫星到NGSO卫星的延长线与地球表面的交点距离，即d＝|rd-rN|。
在一个具体示例中，所述更新GSO卫星的赤经经判断后转到S104，直至完成所有赤经的计算包括：
S1111：更新GSO卫星的赤经：
首先计算LG1：
LG1＝LG+Ls
更新LG：
LG＝LG1
其中，Ls为GSO卫星赤经的更新步长；
S1112：判断LG≤LN+ΔLG，如果符合，转到S104，否则转到S112。
在一个具体示例中，所述更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算包括：
S1121：更新NGSO卫星的纬度幅角：
首先计算u1：
u1＝u+us
更新u：
u＝u1
其中，us为NGSO卫星纬度幅角的更新步长；
S1122：判断u≤360°，若是，转到S102，否则结束。
其中，所述GSO卫星赤经的更新步长和NGSO卫星纬度幅角的更新步长为仿真中的自定义取值，例如取1秒左右。
本发明提供一个具体实施例为NGSO卫星系统的轨道高度为1200km、轨道倾角为87°、轨道偏心率为0；NGSO卫星波束在轨道坐标系的X轴方向的张角为±24°。
在不采取干扰规避措施的情况下，采用上述步骤计算NGSO卫星处于不同纬度幅角、GSO卫星处于不同经度时，某NGSO卫星系统发射的EPFD。GSO卫星经度与NGSO卫星经度差为0°，NGSO卫星发射的EPFD随纬度幅角的变化，以及GSO卫星经度取不同值时NGSO卫星发射的EPFD最大值随纬度幅角的变化如图4所示，从图中可以看出：单颗NGSO卫星发射的最大EPFD为-137.5dB(W/m2)/40kHz，超过了-164dB(W/m2)/40kHz，无需进行采用时间切片的方法进行进一步的仿真，即可确定NGSO卫星系统发射的EPFD超出了EPFD限值。
在采取干扰规避措施的情况下，采用上述步骤计算NGSO卫星处于不同纬度幅角、GSO卫星处于不同经度时，某NGSO卫星系统发射的EPFD。GSO卫星经度与NGSO卫星经度差为0°，NGSO卫星发射的EPFD随纬度幅角的变化，以及GSO卫星经度取不同值时NGSO卫星发射的EPFD最大值随纬度幅角的变化如图5所示，从图中可以看出：单颗NGSO卫星发射的最大EPFD为-165.84dB(W/m2)/40kHz，没有超过-164dB(W/m2)/40kHz，从而可以初步判断采用干扰规避措施的情况下，NGSO卫星系统发射的EPFD满足EPFD限值，从而可以快速判断干扰规避措施的有效性。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

Claims (11)

1.一种下行链路最差干扰快速估算方法，其特征在于，包括以下步骤：
S101：初始化NGSO卫星的纬度幅角；
S102：计算NGSO卫星在J2000坐标系的三维位置矢量、赤经和J2000坐标系到轨道坐标系的坐标转换矩阵；
S103：初始化GSO卫星的赤经；
S104：计算GSO卫星在J2000坐标系的位置矢量；
S105：计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量；
S106：判断GSO卫星到NGSO卫星的延长线是否与地球表面相交，如果相交，转到S107，不相交转到S111；
S107：计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式；
S108：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益；
S109：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD；
S110：存储EPFD，以及对应NGSO卫星纬度幅角、GSO卫星赤经以及GSO卫星到NGSO卫星的延长线与地球表面的交点的坐标；
S111：更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算；
S112：更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算；
S113：完成所述计算，获得所有NGSO卫星纬度幅角、所有GSO卫星赤经对应的所有EPFD值，根据EPFD限值格式，每种场景都存在一个“100％”不能超过的限值，即完全不能超过的值，进一步计算所有EPFD值的最大值，即为最差干扰值。
2.根据权利要求1所述的方法，其特征在于，所述初始化NGSO卫星的纬度幅角包括：
S1011：将NGSO卫星的纬度幅角取为u＝u0，其中，u0为NGSO卫星的初始纬度幅角u0＝0°。
3.根据权利要求2所述的方法，其特征在于，所述计算NGSO卫星在J2000坐标系的位置矢量、赤经以及J2000坐标系到轨道坐标系的坐标转换矩阵包括：
S1021：根据纬度幅角u和近地点幅角ω，计算真近点角：
θ＝u-ω
S1022：根据所述真近点角计算NGSO卫星轨道半径为：
Figure FDA0003690097700000021
其中，a是NGSO卫星半长轴，e是NGSO卫星偏心率；
S1023：NGSO卫星在J2000坐标系的位置矢量为：
Figure FDA0003690097700000022
其中，
Figure FDA0003690097700000023
为单位向量，计算公式为：
Figure FDA0003690097700000024
其中，Ω是NGSO卫星升交点赤经，i是NGSO卫星轨道倾角；
S1024：将NGSO卫星在J2000坐标系的位置矢量的三维分量记为xN、yN和zN，NGSO卫星在J2000坐标系的赤经为：
Figure FDA0003690097700000031
S1025：J2000坐标系到轨道坐标系的坐标转换矩阵为：
Figure FDA0003690097700000032
其中，Mx()、My()和Mz()分别为绕X、Y和Z轴的坐标转换矩阵计算公式为：
Figure FDA0003690097700000033
Figure FDA0003690097700000034
Figure FDA0003690097700000035
4.根据权利要求3所述的方法，其特征在于，所述初始化GSO卫星的赤经包括：
S1031：初始化GSO卫星的赤经：
LG＝LN-ΔLG
其中，ΔLG为GSO卫星赤经区间宽度的一半。
5.根据权利要求4所述的方法，其特征在于，所述计算GSO卫星在J2000坐标系的位置矢量包括：
S1041：GSO卫星在J2000坐标系的位置矢量rG为：
Figure FDA0003690097700000041
其中，Re是地球平均半径，BG为GSO卫星的赤纬，HG为GSO卫星的高度。
6.根据权利要求5所述的方法，其特征在于，所述计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量包括：
S1051：GSO卫星到NGSO卫星的延长线在J2000坐标系的单位矢量为：
Figure FDA0003690097700000042
S1052：GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量为：
Figure FDA0003690097700000043
7.根据权利要求6所述的方法，其特征在于，所述判断GSO卫星到NGSO卫星的延长线是否与地球表面相交包括：
S1061：NGSO卫星到地球的切线与轨道坐标系+Z轴方向的夹角为：
Figure FDA0003690097700000044
S1062：GSO卫星到NGSO卫星的延长线与轨道坐标系+Z轴方向的夹角为：
Figure FDA0003690097700000045
其中，zo是GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量do第三维分量；
S1063：当
Figure FDA0003690097700000051
时，则延长线与地球表面相交；否则，延长线不与地球表面相交。
8.根据权利要求7所述的方法，其特征在于，所述计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式包括：
S1071：GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式为：
rd＝rN+kdi
其中，系数k计算公式为：
Figure FDA0003690097700000052
9.根据权利要求8所述的方法，其特征在于，所述计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD包括：
S1091：NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD为：
EPFD＝Pt+Gt-10log10(4πd2)+Gr-Grmax
其中，Pt是NGSO卫星空间电台的发射功率，Gt为NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益，Gr是GSO地球站的接收增益，Grmax是GSO地球站的最大接收增益，在GSO地球站处于GSO卫星与NGSO卫星延长线方向的情况下，Gr＝Grmax，d是NGSO卫星空间电台与GSO卫星到NGSO卫星的延长线与地球表面的交点的距离，即d＝|rd-rN|。
10.根据权利要求9所述的方法，其特征在于，更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算包括：
S1111：更新GSO卫星的赤经：
首先计算LG1：
LG1＝LG+Ls
更新LG：
LG＝LG1
其中，Ls为GSO卫星赤经的更新步长；
S1112：判断LG≤LN+ΔLG，如果符合，转到S104，否则转到S112。
11.根据权利要求10所述的方法，其特征在于，所述更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算包括：
S1121：更新NGSO卫星的纬度幅角：
首先计算u1：
u1＝u+us
更新u：
u＝u1
其中，us为NGSO卫星纬度幅角的更新步长；
S1122：判断u≤360°，若是，转到S102，否则结束。
CN202110555538.9A 2021-05-21 2021-05-21 一种下行链路最差干扰快速估算方法 Active CN113433571B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110555538.9A CN113433571B (zh) 2021-05-21 2021-05-21 一种下行链路最差干扰快速估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110555538.9A CN113433571B (zh) 2021-05-21 2021-05-21 一种下行链路最差干扰快速估算方法

Publications (2)

Publication Number Publication Date
CN113433571A CN113433571A (zh) 2021-09-24
CN113433571B true CN113433571B (zh) 2022-09-23

Family

ID=77802743

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110555538.9A Active CN113433571B (zh) 2021-05-21 2021-05-21 一种下行链路最差干扰快速估算方法

Country Status (1)

Country Link
CN (1) CN113433571B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6678520B1 (en) * 1999-01-07 2004-01-13 Hughes Electronics Corporation Method and apparatus for providing wideband services using medium and low earth orbit satellites
US6511020B2 (en) * 2000-01-07 2003-01-28 The Boeing Company Method for limiting interference between satellite communications systems
CN110417460B (zh) * 2019-08-16 2021-03-30 国家无线电监测中心 一种非静止轨道卫星对静止轨道卫星干扰的分析方法
CN111541477B (zh) * 2019-11-25 2022-03-04 航天科工空间工程发展有限公司 低轨星座系统内部频率干扰抑制方法及装置
CN111934739B (zh) * 2020-07-02 2022-04-19 航天科工空间工程发展有限公司 一种降低ngso卫星俯仰角度的干扰规避方法
CN112054838B (zh) * 2020-09-11 2022-05-17 中国科学院微小卫星创新研究院 一种ngso卫星偏置方案的设计方法

Also Published As

Publication number Publication date
CN113433571A (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN113131989B (zh) 一种ngso星座系统频谱共享仿真时间参数设计方法
CN105430664B (zh) 一种基于分类拟合预测传播路损的方法和装置
CN109587690A (zh) 无人机辅助认知卫星地面通信网络多点协作传输方法
CN108833041B (zh) 一种基于椭圆轨道的多波束低轨卫星信道仿真方法
CN112039564B (zh) 一种卫星大规模mimo广覆盖预编码传输方法
CN102111873B (zh) 选取可见基站的方法及装置、终端定位方法及装置
CN108123772A (zh) 一种基于梯度投影的无人机时频资源分配方法
CN111147127A (zh) 一种低轨卫星星座全球单重覆盖及频率复用算法
CN114039652B (zh) 基于建筑物几何分析的毫米波抗阻挡多无人机部署方法
CN116582202A (zh) 计算大规模非静止轨道星座下行链路干噪比分布的方法
CN113433571B (zh) 一种下行链路最差干扰快速估算方法
CN107395299B (zh) 一种450MHz频段卫星气象业务与铁路列调系统的干扰分析方法
Last et al. Effect of skywave interference on coverage of radiobeacon DGPS stations
CN111614418A (zh) 无人机数据链电磁环境效应预测方法及设备
CN115827185B (zh) 6g空中基站结合北斗空中避障的方法、存储介质及设备
CN115426027A (zh) 一种基于干扰分布的受扰系统地球站位置计算方法
US20110143657A1 (en) Method of establishing communication link between a mobile earth station and a satellite of mss and apparatus therefor
Hu et al. Analysis of marine wireless communication channel under high sea conditions
CN115173974A (zh) 基于空间资源数据库与轮廓线插值的下行干扰计算方法
Kerczewski et al. Parameter impact on sharing studies between UAS CNPC satellite transmitters and terrestrial systems
CN113660673B (zh) 一种基于面元法的超短波通信地图构建方法
CN111787508A (zh) 基于功率控制的雷达通信一体化uav网络效用优化方法
Wang et al. Performance analysis of spectrum sensing based on distributed satellite clusters under perturbation
US20040056800A1 (en) Base station, base station module and method for estimating parameters of uplink signals
CN118449581A (zh) 基于卫星干扰隔离角优化的ngso星座间干扰规避方法

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
GR01 Patent grant