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

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


背景技术：

2.随着非静止轨道（ngso）卫星日益增多，合适轨位、频谱资源有限，波束密集、频谱重叠、位置动态，由此产生的ngso星座间同频干扰愈趋严重，波束间同频干扰链路多，干扰发生概率高，干扰情形变化快。高强度、集总性、时变性是当前ngso通信星座系统间同频干扰呈现的特点，此前的单链路、静态化干扰场景已经不再适用于当前复杂的干扰环境。
3.干噪比（i/n）、载干比（c/i）、等效功率通量密度（epfd）等是描述ngso通信星座间同频互扰的常用评价指标。现有技术中干噪比概率分布计算方法可分为外推法、数值法、基准卫星快照法。外推法通常采用外推卫星轨道位置、统计不同干扰值出现时间占比的方法获得，一般仿真时长为数月，随着卫星数量不断增多，干扰组合大幅增加，计算总量大，普通计算机难以有效仿真；数值法有两种，其中一种由fortes提出，通过推导卫星概率密度函数，在全球范围内积分计算干噪比概率分布，并采用自适应步长和高斯求积规则减少计算工作量，此方法在实施细节上尚未给出具体分析；另一种由niloofar okati等提出，利用随机几何研究大规模星座通信性能的思路，通过将大规模星座模拟为均匀分布（bpp模型），利用随机几何知识得出相关指标概率分布，数学解析表达式较为复杂，不易进行分析，仅应用于求解下行链路中断概率，尚未在干扰性能上进行应用；基准卫星快照法由清华大学靳瑾及林子翘等提出，主要通过确定基准卫星，划分区域，联合轨道参数以及卫星星座构型参数，得到卫星星座快照，从而计算星座间干噪比概率分布，该方法仍需要快照仿真，难以分析卫星数量、高度、地球站纬度对于干噪比的影响。
4.在巨型星座场景下，卫星数量庞大，沿用外推方法将导致计算量大幅上升，甚至出现普通pc机无法完成仿真，亟需大幅降低计算复杂度。一般的仿真方法不能给出星座高度、数量、地球站纬度与干噪比的关系，不能对星座设计给出参考的问题。针对上述存在的问题，亟需提供一种计算效率高，能够快速分析星座数量、高度、地球站纬度与干噪比之间关系，对星座设计给出参考的方法。


技术实现要素：

5.本发明提供了一种计算大规模非静止轨道星座下行链路干噪比分布的方法，通过用斐波那契网格法来表示干扰星座卫星的坐标以及接入卫星的概率位置坐标，能够推导出非静止星座系统间下行链路干噪比的初等数学表达式，最终得到了非静止轨道星座间的干噪比概率分布。该方法与外推方法和基准卫星快照法相比，计算结果准确度相似，计算效率大幅提升，并且能够快速的分析星座数量、高度、地球站纬度与干噪比之间的关系，可为星座间下行链路互扰分析及星座设计提供一定的参考。
6.具体技术方案如下：计算大规模非静止轨道星座下行链路干噪比分布的方法，包
括如下步骤：s01：建立卫星星座下行链路干扰模型，用于模拟干扰系统对受扰系统的干扰，所述受扰系统包括地球站和a星座，所述干扰系统包括b星座，所述b星座与a星座同频共存；s02：建立干扰评价模型，用于评价所述受扰系统受到的干扰水平；s03：利用斐波那契网格法建立卫星空间位置模型，用于求解卫星空间位置坐标，利用空域卫星出现概率获得等效卫星数；s04：建立卫星星座下行链路干噪比公式，用于计算地球站集总干噪比，根据所述干噪比的大小进行归类处理，获得干噪比分布。
7.优选的，所述卫星空间位置坐标包括b星座内干扰卫星坐标：；式中，为b星座干扰卫星高度最小值，r为地球半径，为b星座轨道高度，为b星座卫星个数。
8.优选的，所述b星座干扰卫星高度最小值的计算公式为：；式中，为b星座干扰卫星对地球站最小仰角。
9.优选的，所述b星座干扰卫星包括b星座位于干扰区域内的卫星。
10.优选的，所述等效卫星数为b星座等效卫星数，计算公式为：；式中，为b星座等效卫星数，为b星座卫星个数，p为b星座卫星在空域内出现的概率，r为地球半径，为b星座轨道高度，为b星座干扰卫星高度最小值，为b星座卫星倾角，l为计算空域中心的纬度，为b星座卫星对地球站最小仰角。
11.优选的，所述地球站集总干噪比的计算公式为：；式中，为所述地球站集总干噪比，为b星座等效卫星数，为a星座与b星座之间的重叠带宽的发射功率，为第n颗干扰卫星的天线发射增益，为地球站天线接收增益，为第n颗干扰卫星与地球站连线偏离地球站天线主轴夹角，为第n颗干
扰卫星与地球站连线偏离卫星天线主轴的夹角，为第n颗干扰卫星到地球站的距离，k为玻耳兹曼常量，t为地球站接收端的系统噪声温度，b为卫星通信带宽，为b星座干扰卫星高度最小值，r为地球半径，为b星座轨道高度，为通信载波波长。
12.优选的，所述a星座与b星座之间的重叠带宽的发射功率计算公式为：；；式中，为b星座发射功率，f为重叠带宽因子，为a星座的通信带宽，为a星座的中心频率，为b星座的通信带宽，为b星座的中心频率。
13.优选的，包括计算不同的t值对应的值，根据大小进行归类处理。
14.优选的，所述归类处理包括以1db为区间粒度，统计干噪比值在区间［e-1，e）的概率，计算公式为：；式中，e为统计值。
15.优选的，所述干扰评价模型采用的评价指标为干噪比、载干比或等效功率通量密度中的一种。
16.本发明的有益效果如下：本发明提供的计算大规模非静止轨道星座下行链路干噪比分布的方法，计算效率高，在仿真条件相同的前提下本发明提供方法的仿真时间仅为30s，与轨道外推方法相比误差不大，轨道外推方法时间约为500s；本系统能够提供星座高度、数量、地球站纬度对干噪比概率分布影响的分析，可为星座间下行链路互扰分析及星座设计提供一定的参考。
附图说明
17.图1示出了本发明实施例中下行链路干扰模型ngso星座a与ngso星座b同频共存示意图；图2示出了本发明实施例中采用的斐波那契网格法示意图；图3示出了本发明实施例中ngso星座b中干扰卫星的位置坐标示意图；图4示出了本发明实施例中ngso星座b中第n颗卫星干扰接入卫星示意图；图5示出了轨道外推方法和本发明提供方法的测试结果示意图；图6示出了采用轨道外推方法a星座高度变化范围为800km~1400km对i/n值的影响示意图；图7示出了采用本发明提供方法a星座高度变化范围为800km~1400km对i/n值的影响示意图；图8示出了轨道外推方法和本发明提供方法在a星座高度变化范围为800km~
1400km的i/n值对比示意图；图9示出了采用轨道外推方法b星座高度变化范围为800km~2000km对i/n值的影响示意图；图10示出了采用本发明提供方法b星座高度变化范围为800km~2000km对i/n值的影响示意图；图11示出了轨道外推方法和本发明提供方法在b星座高度变化范围为800km~2000km的i/n值对比示意图；其中，图5、图8和图11中图例“＋”代表轨道外推方法结果，“o”代表本发明方法结果。
具体实施方式
18.下面对本技术的实施方式作进一步地详细描述，显然，所描述的实施例仅是本技术的一部分实施例，而不是所有实施例的穷举。需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。
19.说明书和权利要求书中的术语“第一”、“第二”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备，不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
20.应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
21.实施例
22.本发明旨在提供一种计算大规模非静止轨道星座下行链路干噪比分布的方法，以缓解利用外推方法仿真巨型星座的干噪比概率分布计算量大的技术问题、缓解快照分析方法无法实现分析星座参数对干噪比影响的问题。本发明对干噪比、载干比、等效功率通量密度epfd等ngso通信星座间同频互扰的常用评价指标均是适用的。本发明实施例中选用的指标是干噪比，通过计算星座卫星间的干噪比概率分布，评估星座间的干扰情况。主要通过建立下行链路干扰模型、干扰评价模型、卫星空间位置模型推导出卫星下行链路干噪比公式，最终通过计算得到干噪比值的概率分布。
23.一、建立卫星星座下行链路干扰模型
24.图1示出了本发明实施例中下行链路干扰模型当ngso星座a与ngso星座b同频共存示意图，图中包括某地区的地球站e；ngso星座a为e的通信星座，简称为a星座；ngso星座b简称为b星座，b星座与a星座同频共存。其中，受扰系统由a星座中一颗卫星和地球站e组成，该卫星称为接入卫星；干扰系统由b星座中多颗卫星组成，此类卫星称为干扰卫星。规定仰角
大于的空域中的b星座卫星对地球站e能造成干扰，相应的空域称为干扰区域，记为。优选的，b星座为大规模星座，a星座与b星座其运行轨道视为圆轨道，卫星在其轨道面上均匀分布，即每颗卫星与其他卫星之间的最短距离相等。地球半径为r，a星座卫星个数为，a星座卫星轨道高度为，b星座卫星个数为，b星座卫星轨道高度为，b星座卫星倾角为。
25.二、建立干扰评价模型
26.通用的卫星通信干扰评价指标一般有三种，分别是干噪比、载干比、等效功率通量密度epfd。优选的，本实施例利用干噪比来对干扰水平进行评价。单颗卫星通信系统的干扰评价指标，计算公式：；式中：为干扰卫星重叠带宽的发射功率，为卫星天线发射增益，为卫星与地球站连线偏离卫星天线主轴的夹角，为地球站天线接收增益，为卫星与地球站连线偏离地球站天线主轴的夹角，k为玻耳兹曼常量，t为地球站接收端的系统噪声温度，b为卫星通信带宽，d为干扰卫星发射机到地球站的距离，为通信载波波长。
27.根据itu-rs.1432-1建议书，当一个卫星通信系统与另一个卫星系统同频共存，其受到的干噪比值不应超过-12.2db，即系统噪声的6%。
28.三、建立卫星空间位置模型
29.从单颗卫星通信系统的干扰评价指标的计算公式可以看出，对地面固定点来说，给定发射功率、卫星通信带宽、卫星天线发射增益、接收增益等特性后，由、、d确定，即由卫星位置坐标唯一决定，为了给出卫星位置坐标，使用斐波那契网格法。
30.斐波那契网格法是用来解决在球面上均匀分布多个点的问题的一种数学方法。图2示出了本发明实施例中采用的斐波那契网格法示意图，以球心为原点建立空间直角坐标系，在单位球面上，均匀分布n个点，各点的坐标如下式所示：；式中，n为均匀分布的点，n为点的个数，。
31.优选的，卫星空间位置坐标的求解方法，以地心为原心，以地心与地球站连线为z
轴，建立空间直角坐标系，分别给出干扰卫星位置坐标以及接入卫星位置坐标。
32.优选的，干扰卫星位置坐标的求解方法，利用斐波那契网格法，b星座卫星坐标如下所示为：；；式中，为b星座卫星个数，r为地球半径，为b星座轨道高度。
33.由于只有区域内卫星对地球站造成干扰，图3示出了本发明实施例中ngso星座b中干扰卫星的位置坐标示意图，计算卫星坐标，考虑空间在xoz平面的投影，b星座干扰卫星高度最小值记为，若卫星的z轴坐标高于，则说明该卫星出现在区域内，根据空间几何学的知识，计算出，计算公式如下：；；令，得。
34.干扰区域内干扰卫星坐标表达式为：；式中，为b星座干扰卫星高度最小值，r为地球半径，为b星座轨道高度，为b星座卫星个数，为b星座干扰卫星对地球站最小仰角。
35.以上星座模型是均匀分布的，无论地球站位于何处，其干扰区域内的卫星数量都是相等的。而对于当前的星座而言，卫星并非均匀的分布在整个球面上，其分布与地球站的纬度、卫星轨道的倾角有关。对于固定轨道倾角高度的星座，为有效描述不同纬度下干扰卫星个数的不同，引入等效卫星数。
36.对某空域，星座的卫星在中出现概率为：；；
式中：a为空域的球面积；γ为卫星轨道的倾角；l为计算空域中心的纬度；、分别为空域的俯仰角差、方位角差对应的地心角差。
37.对于b星座，在干扰区域内出现的卫星个数为。假设卫星均匀分布，可得干扰卫星个数为，式中，s为的表面积。
38.为使两种方法在干扰区域内的卫星数相等，调整b星座卫星总数等于等效卫星数，即。
39.b星座等效卫星数为：；式中，为b星座等效卫星数，为b星座卫星个数，p为b星座卫星在空域内出现的概率，r为地球半径，为b星座轨道高度，为b星座干扰卫星高度最小值，为b星座卫星倾角，l为计算空域中心的纬度，为b星座卫星对地球站最小仰角。优选的，l为计算空域中心的纬度，l等于地球站纬度。
40.四、卫星星座下行链路干噪比公式
41.图4示出了本发明实施例中ngso星座b中第n颗卫星干扰接入卫星示意图，对于接入卫星每一个可能出现的位置坐标入卫星每一个可能出现的位置坐标，该位置的接入卫星也用表示，为t的函数，不同的t时刻值（即不同的接入卫星位置）对应不同的干噪比值。地球站的集总干噪比计算式：；式中，为a星座与b星座之间的重叠带宽的发射功率，为第n颗干扰卫星的天线发射增益，为地球站天线接收增益，为第n颗干扰卫星与地球站连线偏离地球站天线主轴夹角，为第n颗干扰卫星与地球站连线偏离卫星天线主轴的夹角，为第n颗干扰卫星到地球站的距离，k为玻耳兹曼常量，t为地球站接收端的系统噪声温度，b为卫星通信带宽，为b星座干扰卫星高度最小值记，r为地球半径，为b星座轨道高度，为通信载波波长。
42.a星座与b星座之间的重叠带宽的发射功率由下式计算：；
；式中，f为重叠带宽因子，为a星座与b星座之间的重叠带宽的发射功率，为a星座的通信带宽，为a星座的中心频率，为b星座的通信带宽，为b星座的中心频率，为b星座发射功率。
43.低地球轨道系统卫星天线，如没有实测的方向图数据，可根据itu-r s.1528建议书提供的方向图，得：；；式中，为天线直径，为通信载波波长，为天线的轴外角度，为为天线半功率波束宽度，为主瓣最大增益，为0dbi远端副瓣电平。
44.地球站，如果没有实测的方向图数据，可根据itu-rs.465给出的ngso地球站天线方向图得：；；式中，为天线等效半径，为通信载波波长，为天线的轴外角度。
45.假定干扰卫星发射天线的中心方向是指向地心的。干扰卫星的坐标为，地球站坐标（0，0，r），地心坐标（0，0，0），第n颗干扰卫星与地球站连线偏离卫星天线主轴的夹角，；地球站坐标（0，0，r），干扰卫星坐标为，接入卫星坐标为，可得第n颗干扰卫星与地球站连线偏离地球站天线主轴夹角：。
46.第n颗干扰卫星到地球站的距离为：
；即可获得干扰星座b对地球站在接入卫星时的值。
47.为t的函数，不同的t值（即不同的接入卫星位置）对应不同的干噪比值。通过计算不同的t值对应的值并进行归类处理，可最终得到干噪比值的概率分布。
48.设以1db为区间粒度，统计在区间［e-1，e）的概率，其中e为统计的值，则统计方法为：。
49.对所有干扰值与概率值重复此过程，最终得到的概率分布。在其他参数固定的前提下，该表达式与四个位置参数相关，分别是干扰卫星的个数、高度、接入卫星的高度以及地球站的纬度l。这样，通过控制变量法，即可得到干扰水平与卫星个数、高度、地球站纬度之间的关系。
50.五、系统测试结果评价
51.图5示出了采用本发明提供系统的测试结果示意图，从图中可以看出，本发明提供的非静止轨道星座下行链路干噪比分布计算系统，在仿真条件相同的前提下，与外推方法相比误差不大，计算效率高，步长为5s，模拟仿真10天的数据，轨道外推方法时间约为500秒，本方法仿真时间仅用时30秒。
52.图6示出了采用轨道外推方法a星座高度变化范围为800km~1400km对i/n值的影响示意图，图7示出了采用本发明提供方法星座a星座变化范围为800km~1400km对i/n值的影响示意图，图8示出了轨道外推方法和本发明提供方法在a星座高度变化范围为800km~1400km的i/n值对比示意图。从图6~图8可以看出，随着a星座卫星高度不断增加，概率期望值不断增大，说明b星座对a星座的干扰不断增强。
53.图9示出了采用轨道外推方法b星座高度变化范围为800km~2000km对i/n值的影响示意图,图10示出了采用本发明提供方法b星座高度变化范围为800km~2000km对i/n值的影响示意图,图11示出了轨道外推方法和本发明提供方法在b星座高度变化范围为800km~2000km的i/n值对比示意图。从图9~图11可以看出，随着b星座卫星高度不断增加，概率期望值先增大后减小，说明b星座对a星座的干扰先增强后减弱。
54.本发明实施例中提供的计算大规模非静止轨道星座下行链路干噪比分布的方法，通过用斐波那契网格法来表示干扰星座卫星的坐标以及接入卫星的概率位置坐标，推导出非静止星座系统间下行链路干噪比的初等数学表达式，得到非静止轨道星座间的干噪比概率分布。该方法与现有技术相比，计算结果准确度相似，计算效率大幅提升，并且能够快速的分析星座数量、高度、地球站纬度与干噪比值之间的关系，可为星座间下行链路互扰分析及星座设计提供一定的参考。
55.显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

技术特征：
1.一种计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，包括如下步骤：s01：建立卫星星座下行链路干扰模型，用于模拟干扰系统对受扰系统的干扰，所述受扰系统包括地球站和a星座，所述干扰系统包括b星座，所述b星座与a星座同频共存；s02：建立干扰评价模型，用于评价所述受扰系统受到的干扰水平；s03：利用斐波那契网格法建立卫星空间位置模型，用于求解卫星空间位置坐标，利用空域卫星出现概率获得等效卫星数；s04：建立卫星星座下行链路干噪比公式，用于计算地球站集总干噪比，根据所述干噪比的大小进行归类处理，获得干噪比分布。2.根据权利要求1所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述卫星空间位置坐标包括b星座内干扰卫星坐标：；式中，为b星座干扰卫星高度最小值，r为地球半径，为b星座轨道高度，为b星座卫星个数。3.根据权利要求2所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述b星座干扰卫星高度最小值的计算公式为：；式中，为b星座干扰卫星对地球站最小仰角。4.根据权利要求3所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述b星座干扰卫星包括b星座位于干扰区域内的卫星。5.根据权利要求1所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述等效卫星数为b星座等效卫星数，计算公式为：；式中，为b星座等效卫星数，为b星座卫星个数，p为b星座卫星在空域内出现的概率，r为地球半径，为b星座轨道高度，为b星座干扰卫星高度最小值，为b星座卫星倾角，l为计算空域中心的纬度，为b星座卫星对地球站最小仰角。6.根据权利要求1所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述地球站集总干噪比的计算公式为：
；式中，为所述地球站集总干噪比，为b星座等效卫星数，为a星座与b星座之间的重叠带宽的发射功率，为第n颗干扰卫星的天线发射增益，为地球站天线接收增益，为第n颗干扰卫星与地球站连线偏离地球站天线主轴夹角，为第n颗干扰卫星与地球站连线偏离卫星天线主轴的夹角，为第n颗干扰卫星到地球站的距离，k为玻耳兹曼常量，t为地球站接收端的系统噪声温度，b为卫星通信带宽，为b星座干扰卫星高度最小值，r为地球半径，为b星座轨道高度，为通信载波波长。7.根据权利要求6所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述a星座与b星座之间的重叠带宽的发射功率计算公式为：；；式中，为b星座发射功率，f为重叠带宽因子，为a星座的通信带宽，为a星座的中心频率，为b星座的通信带宽，为b星座的中心频率。8.根据权利要求6所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，包括计算不同的t值对应的值，根据大小进行归类处理。9.根据权利要求8所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述归类处理包括以1db为区间粒度，统计干噪比值在区间［e-1，e）的概率，计算公式为：；式中，e为统计值。10.根据权利要求1所述的计算大规模非静止轨道星座下行链路干噪比分布的方法，其特征在于，所述干扰评价模型采用的评价指标为干噪比、载干比或等效功率通量密度中的一种。

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


技术研发人员：张云帆 储飞黄 贾录良 万颖 苏英豪 曾庆尧
受保护的技术使用者：中国人民解放军战略支援部队航天工程大学
技术研发日：2023.07.13
技术公布日：2023/8/14