[发明专利]计算大规模非静止轨道星座下行链路干噪比分布的方法

说明书

本发明公开了计算大规模非静止轨道星座下行链路干噪比分布的方法，包括：建立卫星星座下行链路干扰模型，用于模拟干扰系统对受扰系统的干扰，所述受扰系统包括地球站和A星座，所述干扰系统包括B星座，所述B星座与A星座同频共存；建立干扰评价模型，用于评价所述受扰系统受到的干扰水平；利用斐波那契网格法建立卫星空间位置模型，用于求解卫星空间位置坐标，利用空域卫星出现概率获得等效卫星数；建立卫星星座下行链路干噪比公式，用于计算地球站集总干噪比，根据所述干噪比的大小进行归类处理，获得干噪比分布。本发明方法能快速的分析星座数量、高度、地球站纬度与干噪比之间的关系，为星座间下行链路互扰分析及星座设计提供一定的参考。

技术领域

本发明涉及卫星通信技术领域，尤其涉及一种计算大规模非静止轨道星座下行链路干噪比分布的方法。

背景技术

随着非静止轨道（NGSO）卫星日益增多，合适轨位、频谱资源有限，波束密集、频谱重叠、位置动态，由此产生的NGSO星座间同频干扰愈趋严重，波束间同频干扰链路多，干扰发生概率高，干扰情形变化快。高强度、集总性、时变性是当前NGSO通信星座系统间同频干扰呈现的特点，此前的单链路、静态化干扰场景已经不再适用于当前复杂的干扰环境。

干噪比（I/N）、载干比（C/I）、等效功率通量密度（EPFD）等是描述NGSO通信星座间同频互扰的常用评价指标。现有技术中干噪比概率分布计算方法可分为外推法、数值法、基准卫星快照法。外推法通常采用外推卫星轨道位置、统计不同干扰值出现时间占比的方法获得，一般仿真时长为数月，随着卫星数量不断增多，干扰组合大幅增加，计算总量大，普通计算机难以有效仿真；数值法有两种，其中一种由Fortes提出，通过推导卫星概率密度函数，在全球范围内积分计算干噪比概率分布，并采用自适应步长和高斯求积规则减少计算工作量，此方法在实施细节上尚未给出具体分析；另一种由Niloofar Okati等提出，利用随机几何研究大规模星座通信性能的思路，通过将大规模星座模拟为均匀分布（BPP模型），利用随机几何知识得出相关指标概率分布，数学解析表达式较为复杂，不易进行分析，仅应用于求解下行链路中断概率，尚未在干扰性能上进行应用；基准卫星快照法由清华大学靳瑾及林子翘等提出，主要通过确定基准卫星，划分区域，联合轨道参数以及卫星星座构型参数，得到卫星星座快照，从而计算星座间干噪比概率分布，该方法仍需要快照仿真，难以分析卫星数量、高度、地球站纬度对于干噪比的影响。

在巨型星座场景下，卫星数量庞大，沿用外推方法将导致计算量大幅上升，甚至出现普通PC机无法完成仿真，亟需大幅降低计算复杂度。一般的仿真方法不能给出星座高度、数量、地球站纬度与干噪比的关系，不能对星座设计给出参考的问题。针对上述存在的问题，亟需提供一种计算效率高，能够快速分析星座数量、高度、地球站纬度与干噪比之间关系，对星座设计给出参考的方法。

发明内容

本发明提供了一种计算大规模非静止轨道星座下行链路干噪比分布的方法，通过用斐波那契网格法来表示干扰星座卫星的坐标以及接入卫星的概率位置坐标，能够推导出非静止星座系统间下行链路干噪比的初等数学表达式，最终得到了非静止轨道星座间的干噪比概率分布。该方法与外推方法和基准卫星快照法相比，计算结果准确度相似，计算效率大幅提升，并且能够快速的分析星座数量、高度、地球站纬度与干噪比之间的关系，可为星座间下行链路互扰分析及星座设计提供一定的参考。

具体技术方案如下：计算大规模非静止轨道星座下行链路干噪比分布的方法，包括如下步骤：

S01：建立卫星星座下行链路干扰模型，用于模拟干扰系统对受扰系统的干扰，所述受扰系统包括地球站和A星座，所述干扰系统包括B星座，所述B星座与A星座同频共存；

S02：建立干扰评价模型，用于评价所述受扰系统受到的干扰水平；

S03：利用斐波那契网格法建立卫星空间位置模型，用于求解卫星空间位置坐标，利用空域卫星出现概率获得等效卫星数；