NGSO卫星地球探测系统对深空探测系统的干扰确定方法

1.NGSO卫星地球探测系统对深空探测系统的干扰确定方法，其特征在于步骤如下：
第一步，明确两个时变系统，即NGSO卫星地球探测系统和深空探测系统的轨迹参数NGSO卫星地球探测系统和深空探测系统的轨迹参数包括历元时刻、卫星轨道高度、轨道倾角、升交点赤经、近地点幅角、真近点角和椭圆轨道的偏心率；
第二步，统一两个时变系统的参照系
将深空探测系统的轨迹参数使用地心坐标系坐标来表示，得到两个时变系统的空间位置信息，即方位角、仰角和距离
第三步，明确两个时变系统的性能参数
所述性能参数包括卫星发射和接收系统特性、地球站发射和接收系统特性、载波功率特性、频率范围、地球站站址、业务链路特性及通信系统保护门限要求；
第四步，确定干扰分析关注的时间弧段
NGSO卫星地球探测系统对深空探测系统的干扰链路建立时间，与深空探测系统链路建立时间的重合部分，即为干扰分析过程中关注的时间弧段；
第五步，确定干扰分析总时间及采样频度
NGSO卫星地球探测系统和深空探测系统同时在轨的工作时刻为时间起点，NGSO卫星地球探测系统和深空探测系统的一方中卫星系统的寿命到期为时间终点，即干扰分析总时间；采样频度根据能够捕捉到的干扰发生数量决定；
第六步，确定干扰链路上的发射和接收增益及干扰链路距离长度
根据第五步确定的干扰分析总时间及采样频度，得到两个时变系统对应于干扰分析总时间内相应采样频度确定的每一时刻的空间位置信息，求得NGSO卫星地球探测系统对深空探测系统干扰链路上的发射和接收增益及干扰链路距离长度；
第七步，确定被干扰链路接收系统的功率谱密度
根据第六步计算得到的干扰链路上的发射和接收增益及干扰链路距离长度，通过传输方程计算相应采样频度确定的每一时刻的干扰信号到达被干扰链路接收系统，即卫星或地球站的功率谱密度；
PSD＝P+G(θ)+G0(φ)-FSL-10log10(BWIa) (1)
其中，
P——干扰信号的发射功率，dBW；
G(θ)——干扰信号在被干扰链路接收方向的天线增益，dBi；
G0(φ)——被干扰链路接收系统在干扰信号方向的接收增益，dBi；
FSL——干扰信号到达被干扰链路接收系统的自由空间传播路径损失，dB；
BWIa——干扰信号载波带宽，Hz；
第八步，确定同频干扰
根据被干扰链路接收系统的干扰保护标准，即可接受的干扰信号的功率谱密度的限值，将第七步计算得到的功率谱密度与限值要求做比较，若超过限值要求则说明该采样频度确定的时刻存在不可接受的干扰。
2.根据权利要求1所述的NGSO卫星地球探测系统对深空探测系统的干扰确定方法，其特征在于：所述第一步中，历元时刻的选取与所述两个时变系统的轨迹参数初始状态直接相关，对干扰技术结果有直接影响，通常选取NGSO卫星地球探测系统进入轨道的时间起点。

NGSO卫星地球探测系统对深空探测系统的干扰确定方法\n技术领域\n[0001] 本发明涉及一种NGSO(非地球静止轨道)卫星地球探测系统对深空探测系统使用相同频率的干扰确定方法，适用于气象、海洋、资源、环境等NGSO卫星地球探测系统与火星探测、金星探测等深空探测系统之间的干扰确定。\n背景技术\n[0002] 根据无线电频率使用的国际规则要求，深空探测指对距离地球表面200万公里外的天体目标进行探测，由于其距离远、信号传播路径长，导致其信号非常微弱，对外来干扰极其敏感。同时，随着当今卫星地球探测及近地空间探测等业务的蓬勃发展，各种相关无线电业务的卫星系统不断增多，对相应卫星系统在深空探测专用频段内产生干扰的可能性和干扰程度进行评估就变得日益重要。\n[0003] 干扰判断准则和干扰确定方法是进行上述评估的必要技术基础。\n[0004] 在干扰判断准则方面，国际电联已提出了对深空探测业务系统在相应无线电频段内的保护标准，即\n[0005] 1)对深空探测地球站的保护标准：\n[0006] -222dB(W/Hz)，(适用频段：2GHz附近)；\n[0007] -221dB(W/Hz)，(适用频段：8GHz附近)；\n[0008] -220dB(W/Hz)，(适用频段：13GHz附近)；\n[0009] -217dB(W/Hz)，(适用频段：32GHz附近)；\n[0010] 2)对深空探测器的保护标准：\n[0011] -193dB(W/20Hz)，(适用频段：2GHz附近)；\n[0012] -190dB(W/20Hz)，(适用频段：7GHz附近)；\n[0013] -186dB(W/20Hz)，(适用频段：17GHz附近)；\n[0014] -183dB(W/20Hz)，(适用频段：34GHz附近)；\n[0015] 至于干扰确定方法方面，目前国际上的研究主要集中在GSO卫星系统之间的干扰评估方法，并形成了协调弧、ΔT/T、C/I等较为成熟的干扰确定方法，但这类方法不适用于时变系统之间的干扰分析。\n[0016] 由于发射深空探测系统目前仍局限于极少数宇航业务领域较发达国家，在我国该业务如火星探测等项目也仅处于刚刚起步阶段，各项研究主要集中在关键技术突破、系统任务实现和实验验证等方面，而对来自其他NGSO卫星系统的同频干扰问题的研究则十分薄弱，尤其在干扰确定方法方面，目前国际上并没有统一而权威的方法，基本处于空白状态。本发明正是基于上述背景，为了解决如何确定NGSO卫星地球探测系统对深空探测系统在相应频段内的干扰提出的方法。\n发明内容\n[0017] 本发明的技术解决问题：克服现有技术的不足，提供一种NGSO卫星地球探测系统对深空探测系统的干扰确定方法，能够实现对NGSO卫星地球探测系统对深空探测系统在相应频段内可能产生的干扰准确的判断，从而对提高深空探测系统的抗干扰能力，并且工程实现容易。\n[0018] 本发明技术解决方案：NGSO卫星地球探测系统对深空探测系统的同频干扰确定方法，其特点在于步骤如下：\n[0019] 第一步，明确两个时变系统，即NGSO卫星地球探测系统和深空探测系统的轨迹参数\n[0020] NGSO卫星地球探测系统和深空探测系统的轨迹参数包括历元时刻、卫星轨道高度、轨道倾角、升交点赤经、近地点幅角、真近点角和椭圆轨道的偏心率；其中，历元时刻的选取与两个时变系统的轨迹参数初始状态直接相关，对干扰技术结果有直接影响，通常选取卫星系统进入轨道的时间起点；\n[0021] 第二步，统一两个时变系统的参照系\n[0022] 将深空探测器的轨迹参数使用地心坐标系坐标来表示，得到两个时变系统的空间位置信息，即方位角、仰角和距离\n[0023] 第三步，明确两个时变系统的性能参数\n[0024] 所述性能参数包括卫星发射和接收系统特性、地球站发射和接收系统特性、载波功率特性、频率范围、地球站站址、业务链路特性及通信系统保护门限要求；\n[0025] 第四步，确定干扰分析关注的时间弧段\n[0026] NGSO卫星地球探测系统对深空探测系统的干扰链路建立时间，与深空探测系统被干扰链路建立时间的重合部分，即为干扰分析过程中关注的时间弧段；\n[0027] 第五步，确定干扰分析总时间及采样频度\n[0028] NGSO卫星地球探测系统和深空探测系统同时在轨的工作时刻为时间起点，NGSO卫星地球探测系统和深空探测系统的一方中卫星系统的寿命到期为时间终点，即干扰分析总时间；采样频度根据能够捕捉到的干扰发生数量决定；\n[0029] 第六步，确定干扰链路上的发射和接收增益及干扰链路距离长度[0030] 根据第五步确定的干扰分析总时间及采样频度，得到两个时变系统对应于干扰分析总时间内相应采样频度确定的每一时刻的空间位置信息，求得NGSO卫星地球探测系统对深空探测系统干扰链路上的发射和接收增益及干扰链路距离长度；\n[0031] 第七步，确定被干扰链路接收系统的功率谱密度\n[0032] 根据第六步计算得到的干扰链路上的发射和接收增益及干扰链路距离长度，通过传输方程计算相应采样频度确定的每一时刻的干扰信号到达被干扰链路接收系统，即卫星或地球站的功率谱密度；\n[0033] 第八步，确定同频干扰\n[0034] 根据被干扰链路接收系统的干扰保护标准，即可接受的干扰信号的功率谱密度的限值，将第七步计算得到的功率谱密度与限值要求做比较，若超过限值要求则说明该采样频度确定的时刻存在不可接受的干扰。\n[0035] 所述第一步中，历元时刻的选取与所述两个时变系统的轨道参数初始状态直接相关，对干扰技术结果有直接影响，通常选取卫星系统进入轨道的时间起点。\n[0036] 本发明与现有技术相比的优点在于：通过本发明上述干扰确定方法，可以实现对NGSO卫星地球探测系统对深空探测系统在相应频段内可能产生的干扰进行准确的判断。例如，对干扰发生的频度、干扰持续时间、干扰强度及对整个深空探测系统在一定保护标准之下的干扰发生的时间百分比等得到数字化的评估结果，对NGSO卫星地球探测系统的干扰抑制要求提供指导性参考，从而对提高深空探测系统的抗干扰能力；本发明填补了国际上关于NGSO卫星地球探测系统与深空探测系统之间干扰确定方法的空白，并且工程实现容易。\n附图说明\n[0037] 图1为NGSO卫星地球探测系统与深空探测系统之间干扰确定方法的流程图；\n[0038] 图2为NGSO卫星地球探测系统对火星探测地球站接收系统的干扰模型示意图；\n[0039] 图3为火星探测系统地球站接收天线旁瓣模式示意图。\n具体实施方式\n[0040] 如图1所示，本发明具体实现步骤如下：\n[0041] 第一步，明确两个时变系统的轨迹参数。NGSO卫星地球探测系统和深空探测系统的轨迹参数包括历元时刻、卫星轨道高度、轨道倾角、升交点赤经、近地点幅角、真近点角和椭圆轨道的偏心率。其中，历元时刻的选取与两个时变系统的轨迹初始状态直接相关，对干扰技术结果有直接影响，通常选取卫星系统进入轨道的时间起点。历元时刻的选取与所述两个时变系统的轨道参数初始状态直接相关，对干扰技术结果有直接影响，通常选取卫星系统进入轨道的时间起点。\n[0042] 第二步，统一两个时变系统的参照系\n[0043] 将深空探测器的轨迹参数使用地心坐标系坐标来表示，得到两个时变系统的空间位置信息，即方位角、仰角和距离。\n[0044] 第三步，明确两个时变系统的性能参数\n[0045] 所述性能参数包括卫星发射和接收系统特性、地球站发射和接收系统特性、载波功率特性、频率范围、地球站站址、业务链路特性及通信系统保护门限要求。\n[0046] 第四步，确定干扰分析关注的时间弧段\n[0047] NGSO卫星地球探测系统对深空探测系统的干扰链路建立时间，与深空探测系统被干扰链路建立时间的重合部分，即为干扰分析过程中关注的时间弧段。\n[0048] 第五步，确定干扰分析总时间及采样频度\n[0049] NGSO卫星地球探测系统和深空探测系统同时在轨的工作时刻为时间起点，NGSO卫星地球探测系统和深空探测系统的一方中卫星系统的寿命到期为时间终点，即干扰分析总时间；采样频度根据能够捕捉到的干扰发生数量决定。\n[0050] 第六步，确定干扰链路上的发射和接收增益及干扰链路距离长度[0051] 根据第五步确定的干扰分析总时间及采样频度，得到两个时变系统对应于干扰分析总时间内相应采样频度确定的每一时刻的空间位置信息，求得NGSO卫星地球探测系统对深空探测系统干扰链路上的发射和接收增益及干扰链路距离长度。\n[0052] 第七步，确定被干扰链路接收系统的功率谱密度\n[0053] 根据第六步计算得到的干扰链路上的发射和接收增益及干扰链路距离长度，通过传输方程计算相应采样频度确定的每一时刻下的干扰信号到达被干扰链路接收系统，即卫星或地球站的功率谱密度。\n[0054] 式一：PSD＝P+G(θ)+G0(φ)-FSL-10log10(BWIa)\n[0055] 其中，\n[0056] P——干扰信号的发射功率，dBW；\n[0057] G(θ)——干扰信号在被干扰链路接收方向的天线增益，dBi；\n[0058] G0(φ)——被干扰链路接收系统在干扰信号方向的接收增益，dBi；\n[0059] FSL——干扰信号到达被干扰链路接收系统的自由空间传播路径损失，dB；\n[0060] BWIa——干扰信号载波带宽，Hz。\n[0061] 第八步，确定同频干扰\n[0062] 根据被干扰链路接收系统的干扰保护标准，即可接受的干扰信号的功率谱密度的限值，将第七步计算得到的功率谱密度与限值要求做比较，若超过限值要求则说明该采样频度确定的时刻存在不可接受的干扰。\n[0063] 下面以NGSO卫星地球探测系统和火星探测系统下行业务(空对地方向)在\n8400-8450GHz频段的干扰确定方法为例，对本发明做进一步说明。\n[0064] 实施例1\n[0065] NGSO卫星地球探测系统对火星探测地球站接收系统的干扰模型示意图如图2所示。其中，卫星地球探测系统地球站接收来自卫星地球探测系统的传输信号，深空探测系统地球站接收来自深空探测系统的传输信号。由于卫星地球探测系统发射天线在深空探测系统地球站方向也有一定的增益，从而导致干扰信号被深空探测系统地球站接收。使用本发明方法确定这一干扰。\n[0066] 1.明确两个卫星系统的轨迹参数\n[0067] NGSO卫星地球探测系统和火星探测系统轨迹参数如表1所示。\n[0068] 表1 NGSO卫星地球探测系统与火星探测系统轨迹参数\n[0069] \n[0070] 2 统一两个时变系统的参照系\n[0071] 将NGSO卫星地球探测系统和火星探测系统的轨迹参数使用地心坐标系表示，如表2所示。\n[0072] 表2 统一NGSO卫星地球探测系统与火星探测系统的参照系\n[0073] \n[0074] 3 明确双方卫星系统性能参数\n[0075] NGSO卫星地球探测系统和火星探测系统参数如表3所示。\n[0076] 表3 NGSO卫星地球探测系统与火星探测系统参数\n[0077] \n[0078] \n[0079] 4.确定干扰分析关注的时间弧段\n[0080] 火星探测系统地球站起始跟踪仰角：7°-88°，即认为仰角为7°-88°时地球站都可以跟踪捕获到火星探测器，进而建立被干扰链路；NGSO卫星地球探测系统地球站起始跟踪仰角：5°，即认为仰角大于或等于5°时地球站都可以跟踪捕获到NGSO卫星，进而建立干扰链路。NGSO卫星地球探测系统对深空探测系统的干扰链路建立时间，与火星探测系统被干扰链路建立时间的重合部分，即为干扰分析过程中关注的时间弧段。\n[0081] 5.确定干扰分析总时间及采样频度\n[0082] 干扰分析时间：2012.10.8 05:00:00-2013.10.8 05:00:00；\n[0083] 干扰计算参数采样频度(仿真步长)：5秒。\n[0084] 6.求得每个采样时刻对应的干扰卫星所处的空间位置，干扰卫星在被干扰地球站方向的天线增益、被干扰地球站在干扰卫星方向的接收增益，以及干扰链路距离长度。\n[0085] 7.计算干扰信号到达被干扰链路地球站接收系统的功率谱密度。\n[0086] 干扰信号的功率谱密度(PSD，单位dB(W/Hz))按照下式计算：\n[0087] PSD＝P+G(θ)+G0(φ)-FSL-10log10(BWIa)\n[0088] 式中，\n[0089] P——干扰卫星的发射功率，dBW；\n[0090] G(θ)——干扰卫星在被干扰地球站方向的天线增益，dBi；\n[0091] G0(φ)——被干扰地球站在干扰卫星方向的接收增益，dBi；\n[0092] FSL——干扰信号到达被干扰地球站的自由空间传播路径损失，dB；\n[0093] BWIa——干扰信号载波带宽，Hz。\n[0094] 8.将第7步计算所得的每个采样时刻的干扰信号功率谱密度与深空探测器保护标准(7GHz附近，干扰信号功率谱密度不得超过-190dB(W/20Hz))进行比较，可以对干扰情况进行确定。\n[0095] 通过计算机仿真实现上述干扰确定方法实施案例，可以看到NGSO卫星地球探测系统与火星探测系统可以兼容工作，评估结果与实际状态相符。因此，本发明干扰确定方法可以广泛应用于NGSO卫星地球探测系统与深空探测系统之间的干扰确定，从而提高深空探测系统的抗干扰能力。