高低轨同频共存卫星通信系统的共线干扰估计方法

技术领域

本发明属于卫星通信领域，尤其涉及一种高低轨同频共存卫星通信系统的共线(In-Line)干扰估计方法。

背景技术

随着卫星通信技术的不断发展，LEO低轨宽带卫星通信星座成为全球各国家和机构的研究热点。由于空间频率资源受限，LEO星座与高通量GEO卫星之间不可避免的存在多系统同频工作的情况。对于高低轨同频共存卫星通信系统来说，在频谱共享过程中，部分卫星或终端相距较近的情况下，会出现通信链路重合的情况，此时会产生共线(In-Line)干扰，对正常的卫星通信造成极大的影响。因此，如何针对不同的高低轨同频共存卫星通信系统应用场景，估计LEO(低轨道)卫星通信系统和GEO(高轨道)卫星通信系统之间的共线干扰，成为了当前构建大规模宽带LEO星座系统过程中亟待解决的问题。

目前，有关同频干扰估计的研究主要集中在单一卫星通信系统中多用户接入所引起的同频干扰，尚无针对高低轨同频共存卫星通信系统共线干扰的研究。公开号为CN110572192的中国专利中，提出了一种低轨卫星的频率规避方法和装置，提出了LEO与GEO卫星的共线干扰规避策略，并没有研究针对LEO与GEO系统In-Line条件下的干扰分析方法。

发明内容

针对高低轨同频共存卫星通信系统应用场景下如何估计LEO卫星通信系统和GEO卫星通信系统之间的In-Line干扰问题，本发明公开了一种高低轨同频共存卫星通信系统的共线干扰估计方法，包括如下步骤：

S1，根据高低轨卫星通信系统的星历数据、地面站经纬度，构建两者之间的干扰分析场景；

S2，根据星历数据，采样不同时刻低轨卫星的运动轨迹，筛选出高低轨卫星通信系统间的共线干扰链路；

S3，根据筛选出的共线干扰链路，区分干扰链路类型；

S4，配置地面站、LEO卫星和GEO卫星的通信终端参数、卫星星历及地球站位置参数；

S5，分别估计接收端收到的来自信号源的载波功率C、噪声功率N和来自干扰源的载波功率I；

所述共线(In-Line)干扰，为部分卫星或终端相距较近时出现的干扰链路与通信链路重合所导致的干扰，即GEO卫星地面站与LEO卫星地面站相距较近或同址时，当LEO卫星运行到地面站至GEO链路附近时，导致地面站至GEO链路与地面站至LEO链路重合，产生共线干扰。

步骤S1中所述的高低轨卫星通信系统间的干扰分析场景，高轨卫星通信系统为高轨卫星星座与其地面站构成的卫星通信系统，低轨卫星通信系统为低轨卫星星座与其地面站构成的卫星通信系统；两者之间的干扰分析场景构建，是基于两者的星历参数和地面站位置，完成卫星星历数据与地球站经纬度的转换，将卫星星历数据转换为经度、纬度、高度信息，完成轨道坐标系与J2000地球惯性坐标系之间的转换，将卫星星历所依托的轨道坐标系转换至J2000地球惯性坐标系，再完成高低轨卫星通信系统的干扰分析场景构建。

所述的将卫星星历数据转换为经度、纬度、高度信息，具体包括：

轨道坐标系六根数换算，用于将轨道六根数换算为轨道坐标系中的坐标值。建立轨道坐标系，该坐标系是空间直角坐标系，该坐标系的坐标原点位于地心，X轴和Y轴位于轨道平面上，Z轴和轨道平面的法线矢量

所述的完成轨道坐标系与J2000地球惯性坐标系之间的转换，将卫星星历所依托的轨道坐标系转换至J2000地球惯性坐标系，具体是完成用J2000地球惯性坐标系所表示的坐标点，与经纬度和大地高坐标系所表示点之间的换算；轨道坐标系坐标轴经三次旋转即可和大地坐标系的三个坐标轴

其中：

所述的完成J2000地球惯性坐标系与地球站位置所依托的地理经纬度坐标系间之间的转换，其具体过程为，

对于将地理经纬度坐标系转换为J2000地球惯性坐标系，其公式为：

其中，

对于将J2000地球惯性坐标系转换为地理经纬度坐标系，其公式为：

步骤S2所述的根据星历数据，采样不同时刻低轨卫星的运动轨迹，筛选出高低轨卫星卫星通信系统间的共线干扰链路，其具体包括：

S21，根据地面站位置及GEO卫星星历参数，设定地面站至GEO链路参数；

S22，根据地面站至GEO链路的链路类型及其链路的波束范围，估计该链路在LEO轨道所在高度上的投影区域，其过程具体为：

S221，基于J2000地球惯性坐标系，假设GEO卫星坐标为(x

S222，假设LEO卫星高度为H

S23，估计当前时刻的LEO卫星位置，判断是否有LEO卫星出现在步骤S222所述的覆盖范围内；

S24，若LEO卫星出现在所述覆盖范围内，则估计共线干扰，并估计下一时刻的LEO卫星位置，重新进行步骤S23的判断；

S25，若LEO卫星未出现在所述覆盖范围内，则估计下一时刻的LEO卫星位置，重新进行步骤S23的判断。

步骤S3所述的干扰链路类型，包括地面站至LEO卫星上行链路、地面站至GEO卫星上行链路、LEO卫星至地面站下行链路和GEO至地面站下行链路。

对于地面站至LEO卫星上行链路，信号源为地面站通信终端，接收端为LEO卫星通信终端，干扰源为地面站；

对于地面站至GEO卫星上行链路，信号源为地面站通信终端，接收端为GEO卫星通信终端，干扰源为地面站；

对于LEO卫星至地面站下行链路，信号源为LEO卫星通信终端，接收端为地面站通信终端，干扰源为GEO卫星；

对于GEO卫星至地面站下行链路，信号源为GEO卫星通信终端，接收端为地面站通信终端，干扰源为LEO卫星。

步骤S4所述的地面站、LEO卫星和GEO卫星的通信终端参数和卫星星历及地球站位置参数具体为，通信终端分为上行链路通信终端和下行链路通信终端，下行链路通信终端参数包括其工作频率(GHz)、卫星天线口径(m)、卫星发射功率(W)、天线效率(％)、用户天线口径(m)、用户热噪声温度(K)、信号带宽(MHz)，上行链路通信终端参数包括其工作频率(GHz)、卫星天线口径(m)、卫星热噪声温度(K)、天线效率(％)、用户发射功率(W)、用户天线口径(m)、信号带宽(MHz)；卫星星历及地球站位置参数包括半长轴、偏心率、倾角、升交点赤经、近地点幅角、真近点角，地球站的经度、纬度、高度。

步骤S5所述的分别估计接收端收到的来自信号源的载波功率C、噪声功率N

载波功率C的估计方法为

噪声功率N

来自干扰源的载波功率I的估计方法为，

针对地面站至LEO卫星上行链路和地面站至GEO卫星上行链路，其发射天线增益G

其中，μ＝2.07123sin(θ)/sin(θ

在不同干扰链路类型下，对应的信号源载波功率、噪声功率和干扰源载波功率的估计的具体公式如下，其中，脚标GS、LEO、GEO分别代表地面站、LEO卫星和GEO卫星：

对于地面站至LEO卫星上行链路，其信号源载波功率C

其中，P

对于地面站至GEO卫星上行链路，其信号源载波功率C

其中，G

对于LEO卫星至地面站下行链路，其信号源载波功率C

其中，P

对于GEO卫星至地面站下行链路，其信号源载波功率C

其中，G

进一步地，所述的在共线干扰条件下接收端的载噪比C/(N

本方法的有益效果为：

一是对于GEO与NGSO系统的共线干扰问题，针对性的提出了一种干扰估计方法，该方法区分了不同链路类型，并在不同干扰链路类型下分别给出对应的信号源载波功率、噪声功率和干扰源载波功率等参数的具体计算方法，提高了共线干扰估计的准确度和针对性。

二是本发明所提的计算共线干扰方法，主要应用于高低轨卫星间的干扰分析，计算卫星系统间链路共线时的干扰情况，即为最差干扰情况，准确量化了不同时刻卫星系统间的干扰程度，可以应用于S、Ku、Ka等多种频段卫星通信系统设计和空间频率资源的申报协调。

附图说明

图1为本发明提供的一种高低轨同频共存卫星通信系统的共线干扰估计方法的流程示意图；

图1中，模块101对应步骤S1，模块102对应步骤S2，模块103对应步骤S3，模块104对应步骤S4，模块105和106对应步骤S5；

图2为构建的干扰分析场景示意图；

图3为地面站至LEO卫星上行链路干扰地面站至GEO卫星上行链路示意图；

图4为筛选高低轨卫星通信系统间的共线干扰链路示意图；

图5为筛选高低轨卫星通信系统间的共线干扰链路流程图。

具体实施方式

为了使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明公开了一种高低轨同频共存卫星通信系统的共线干扰估计方法，包括如下步骤：

S1，根据高低轨卫星通信系统的星历数据、地面站经纬度，构建两者之间的干扰分析场景；

S2，根据星历数据，采样不同时刻低轨卫星的运动轨迹，筛选出高低轨卫星通信系统间的共线干扰链路；

S3，根据筛选出的共线干扰链路，区分干扰链路类型；

S4，配置地面站、LEO卫星和GEO卫星的通信终端参数、卫星星历及地球站位置参数；

S5，分别估计接收端收到的来自信号源的载波功率C、噪声功率N和来自干扰源的载波功率I。

图1为本发明提供的一种高低轨同频共存卫星通信系统的共线干扰估计方法的流程示意图。

所述共线(In-Line)干扰，为部分卫星或终端相距较近时出现的干扰链路与通信链路重合所导致的干扰，即GEO卫星地面站与LEO卫星地面站相距较近或同址时，当LEO卫星运行到地面站至GEO链路附近时，导致地面站至GEO链路与地面站至LEO链路重合，产生共线干扰。

步骤S1中所述的高低轨卫星通信系统间的干扰分析场景，高轨卫星通信系统为高轨卫星星座与其地面站构成的卫星通信系统，低轨卫星通信系统为低轨卫星星座与其地面站构成的卫星通信系统；两者之间的干扰分析场景构建，是基于两者的星历参数和地面站位置，完成卫星星历数据与地球站经纬度的转换，将卫星星历数据转换为经度、纬度、高度信息，完成轨道坐标系与J2000地球惯性坐标系之间的转换，将卫星星历所依托的轨道坐标系转换至J2000地球惯性坐标系，再完成高低轨卫星通信系统的干扰分析场景构建。图2为构建的干扰分析场景示意图。

将卫星星历数据转换为经度、纬度、高度信息，具体包括：

轨道坐标系六根数换算，用于将轨道六根数换算为轨道坐标系中的坐标值。轨道六根数具体包括轨道半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角和真近角。建立轨道坐标系，该坐标系是空间直角坐标系，该坐标系的坐标原点位于地心，X轴和Y轴位于轨道平面上，Z轴和轨道平面的法线矢量

所述的完成轨道坐标系与J2000地球惯性坐标系之间的转换，将卫星星历所依托的轨道坐标系转换至J2000地球惯性坐标系，具体是完成用J2000地球惯性坐标系所表示的坐标点，与经纬度和大地高坐标系所表示点之间的换算；轨道坐标系坐标轴经三次旋转即可和大地坐标系的三个坐标轴

其中：

所述的完成J2000地球惯性坐标系与地球站位置所依托的地理经纬度坐标系间之间的转换，其具体过程为，

对于将地理经纬度坐标系转换为J2000地球惯性坐标系，其公式为：

其中，

对于将J2000地球惯性坐标系转换为地理经纬度坐标系，其公式为：

步骤S2所述的根据星历数据，采样不同时刻低轨卫星的运动轨迹，筛选出高低轨卫星卫星通信系统间的共线干扰链路，筛选流程如图5所示，其具体包括：

S21，根据地面站位置及GEO卫星星历参数，设定地面站至GEO链路参数；

S22，根据地面站至GEO链路的链路类型及其链路的波束范围，估计该链路在LEO轨道所在高度上的投影区域，其过程具体为：

S221，基于J2000地球惯性坐标系，假设GEO卫星坐标为(x

S222，假设LEO卫星高度为H

S23，估计当前时刻的LEO卫星位置，判断是否有LEO卫星出现在步骤S222所述的覆盖范围内；

S24，若LEO卫星出现在所述覆盖范围内，则估计共线干扰，并估计下一时刻的LEO卫星位置，重新进行步骤S23的判断；

S25，若LEO卫星未出现在所述覆盖范围内，则估计下一时刻的LEO卫星位置，重新进行步骤S23的判断。

步骤S3所述的干扰链路类型，包括地面站至LEO卫星上行链路、地面站至GEO卫星上行链路、LEO卫星至地面站下行链路和GEO至地面站下行链路。

对于地面站至LEO卫星上行链路，信号源为地面站通信终端，接收端为LEO卫星通信终端，干扰源为地面站；

对于地面站至GEO卫星上行链路，信号源为地面站通信终端，接收端为GEO卫星通信终端，干扰源为地面站；

对于LEO卫星至地面站下行链路，信号源为LEO卫星通信终端，接收端为地面站通信终端，干扰源为GEO卫星；

对于GEO卫星至地面站下行链路，信号源为GEO卫星通信终端，接收端为地面站通信终端，干扰源为LEO卫星。

图3为地面站至LEO卫星上行链路干扰地面站至GEO卫星上行链路示意图。图4为筛选高低轨卫星通信系统间的共线干扰链路示意图；图5为筛选高低轨卫星通信系统间的共线干扰链路流程图。

步骤S4所述的地面站、LEO卫星和GEO卫星的通信终端参数和卫星星历及地球站位置参数具体为，通信终端分为上行链路通信终端和下行链路通信终端，下行链路通信终端参数包括其工作频率(GHz)、卫星天线口径(m)、卫星发射功率(W)、天线效率(％)、用户天线口径(m)、用户热噪声温度(K)、信号带宽(MHz)，上行链路通信终端参数包括其工作频率(GHz)、卫星天线口径(m)、卫星热噪声温度(K)、天线效率(％)、用户发射功率(W)、用户天线口径(m)、信号带宽(MHz)；卫星星历及地球站位置参数包括半长轴、偏心率、倾角、升交点赤经、近地点幅角、真近点角，地球站的经度、纬度、高度。

步骤S5所述的分别估计接收端收到的来自信号源的载波功率C、噪声功率N

载波功率C的估计方法为

噪声功率N

接收端收到的来自干扰源的载波功率I的估计方法与估计来自信号源的载波功率C的方法一致，来自干扰源的载波功率I的估计方法为，

针对地面站至LEO卫星上行链路和地面站至GEO卫星上行链路，其发射天线增益G

其中，μ＝2.07123sin(θ)/sin(θ

在不同干扰链路类型下，对应的信号源载波功率、噪声功率和干扰源载波功率的估计的具体公式如下，其中，脚标GS、LEO、GEO分别代表地面站、LEO卫星和GEO卫星：

对于地面站至LEO卫星上行链路，其信号源载波功率C

其中，P

对于地面站至GEO卫星上行链路，其信号源载波功率C

其中，G

对于LEO卫星至地面站下行链路，其信号源载波功率C

其中，P

对于GEO卫星至地面站下行链路，其信号源载波功率C

其中，G

进一步地，所述的在共线干扰条件下接收端的载噪比C/(N