CN114257293B - 基于时空割集的ngso星座对geo卫星通信地球站干扰建模方法 - Google Patents

Abstract

本发明公开了一种基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法，包括以下步骤：计算NGSO星座卫星状态分布周期；分析NGSO星座对GEO卫星通信地球站可见干扰卫星数量及各干扰卫星数量时的干扰持续时间；构建NGSO卫星对GEO卫星通信地球站接收干扰信号数学模型；计算NGSO星座对GEO卫星通信地球站接收干扰信号；计算不同NGSO星座对GEO卫星通信地球站可见干扰卫星数量时GEO卫星通信地球站接收到干扰信号的均值；构建时空域分割的NGSO星座对GEO卫星地球站干扰模型。本发明在充分考虑GEO卫星通信地球站地理空间分布和NGSO星座状态分布时间周期等因素的基础上，能够全面、精确地分析和计算GEO卫星通信地球站在全球不同位置、不同时间受NGSO通信星座干扰的情况。

Description

基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模 方法

技术领域

本发明涉及卫星通信技术领域，尤其涉及一种基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法。

背景技术

近年来，卫星通信技术蓬勃发展，卫星通信正在由单一轨道GEO通信系统向高中低轨结合的多轨道卫星通信系统发展、由单颗卫星中继通信向星地一体化组网通信发展，同时随着卫星制造技术的发展以及卫星发射成本的降低，越来越多的卫星通信系统被提出并开始建设。然而，卫星通信的频谱资源是非常有限的。面对卫星通信快速发展所导致的频谱资源枯竭，不同卫星通信系统之间进行同频共用是目前解决上述问题的常用手段，而这势必将带来卫星通信系统间的同频干扰问题。因此，研究非静止轨道星座(简称NGSO星座)对现有GEO卫星通信系统的干扰显得尤为重要。然而，NGSO星座对GEO卫星通信地球站干扰分析存在以下难点，一是地面上不同位置的地球站受到NGSO星座干扰的情况不同，难以对不同GEO卫星通信地球站的受干扰情况进行统一建模；二是NSGO卫星相对于同一GEO卫星通信地球站的位置变化迅速，导致对同一GEO卫星通信地球站干扰的NSGO卫星数量、链路动态变化，难以对同一GEO卫星通信地球站在不同时刻的受干扰情况进行精确建模。因此，针对NGSO星座对GEO卫星通信地球站干扰场景，如何对全球不同位置的GEO卫星通信地球站在不同时段受整个NGSO星座干扰情况进行统一、准确地分析，成为了卫星通信系统间干扰分析、干扰监测中亟待解决的问题。

发明内容

针对NGSO星座对GEO卫星通信系统的干扰分析复杂，干扰情况随NSGO卫星相对于GEO卫星通信地球站的位置快速变化的问题，公开了一种基于时空割集的NGSO星座对GEO卫星地球站干扰建模方法，定义并计算了NGSO星座干扰分析周期，能够对NGSO卫星星座对GEO卫星通信地球站的干扰信号进行精确建模，同时考虑了GEO卫星通信地球站位置对干扰模型的影响，能够对整个NGSO星座对全球不同位置GEO卫星通信地球站进行干扰情况分析。

本发明公开了一种基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法，其步骤包括：

S1，计算非静止轨道星座,即NGSO星座，的卫星状态分布周期；

对于地球表面上任一定点，定义同一NGSO星座在运行过程中出现前后两次的卫星分布完全相同情况的时间差为NGSO星座卫星状态分布周期，记作T_s；对于一个确定的NGSO星座，根据星座的结构和参数，通过计算得到T_s；

当NGSO星座采用倾斜圆轨道星座中的Walker Delta星座时，根据描述其星座结构的Walker代码N/P/F，其中N表示星座的卫星总数，P表示星座的轨道面数量，F表示相位因子，计算NGSO星座卫星状态分布周期，其计算公式为：

其中[，]表示求两者的最小公倍数，T_E表示地球自转周期，T表示NGSO星座的一个轨道周期；

S2，在一个NGSO星座卫星状态分布周期的时间内，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量及各干扰卫星数量的干扰持续时间；

设定GEO卫星通信地球站在大地坐标系中坐标为(l,w,h),其中l表示该地球站的大地经度，w表示该地球站的大地纬度，h表示该地球站的大地高度；

建立地心直角坐标系，原点O设在地球的质量中心，x轴与首子午面与赤道面的交线重合，向东为正，z轴与地球旋转轴重合，向北为正，y轴与xoz平面垂直，三个坐标轴构成右手系；

将GEO卫星通信地球站的大地坐标系坐标转换为地心直角坐标系坐标(x₀,y₀,z₀)，其转换公式为

其中N₀为卯酉圈的半径，

a为地球椭球长半轴，b为地球椭球短半轴，e表示地球椭球曲率；

GEO卫星通信地球站天线方向始终指向GEO卫星，GEO卫星在地心直角坐标系的坐标为(x_S,y_S,0)，将GEO卫星通信地球站的可见范围看作以GEO卫星通信地球站为顶点、以GEO卫星通信地球站到GEO卫星连线为轴、张角为θ的一个圆锥体，其中

此圆锥体的锥面S在地心直角坐标系中的表达式为

其中x_t、y_t、z_t为待计算的中间参数，(x,y,z)为锥面S上任意一点在地心直角坐标系下的坐标；

NGSO星座轨道所在球面在地心直角坐标系中的表达式U为

U:x²+y²+z²＝h_s ²，

其中h_s表示NGSO星座轨道半径，(x,y,z)为NGSO星座轨道所在球面上任意一点在地心直角坐标系下的坐标；

将锥面S的表达式与NGSO星座轨道所在球面的表达式U联立，即得到GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围边界曲线表达式L，记GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围表达式为M；

在一个NGSO星座卫星状态分布周期T_s的时间内，根据NGSO星座轨道参数，将NGSO星座的NGSO卫星k在时间t，在地心直角坐标系的坐标表示为(x_l(t),y_l(t),z_l(t))_k，其中k表示NGSO卫星编号且k∈[1,2,…,N]，t表示时间且0＜t≤T_s；

当NGSO卫星k的坐标(x_l(t),y_l(t),z_l(t))_k∈M时，判定NGSO卫星k为NGSO星座对GEO卫星通信地球站的可见干扰卫星；

根据NGSO星座轨道参数与可见范围M，计算得到位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行期间的最大可见干扰卫星数量，记为n(l,w)；根据NGSO星座在一个NGSO星座卫星状态分布周期T_s的时间内的运行时刻t，得到NGSO星座对GEO卫星通信地球站的在时刻t的可见干扰卫星数量n＝N(l,w,t)，N(l,w,t)为位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行时刻t的可见干扰卫星数量，且N(l,w,t)∈[1，2,…,n(l,w)]；

根据不同时刻的NGSO星座对GEO卫星通信地球站的可见干扰卫星的数量，将一个NGSO星座卫星状态分布周期T_s分割为n(l,w)个种类的时间片段，每一个种类的时间片段构成一个时间片段集合，则可见干扰卫星数量为n时的时间片段集合T_n表达式为：

其中i＝1,2,…,I_nmax，n∈[1，2,…,n(l,w)]，I_nmax表示可见干扰卫星数量为n时的时间片段总数量，(t_2i-1,t_2i)_n表示可见干扰卫星数量为n时的第i个时间片段，t_2i-1表示n个可见干扰卫星的第i次出现时刻，t_2i表示n个可见干扰卫星的第i次出现后的消失时刻；

因此，NGSO星座对GEO卫星通信地球站的各可见干扰卫星数量为n时的持续时间T_total ⁿ表示为：

从而得到各干扰卫星数量的干扰持续时间。

S3，计算GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号；

将运动到GEO卫星通信地球站可见范围内的一个NGSO星座卫星所发射的干扰信号I_s(t)表示为：

其中A_s表示干扰信号幅度，f为干扰信号频率，

表示干扰信号的初相位，I_s(t)表示t时刻的干扰信号；

根据GEO卫星通信地球站天线接收增益、NGSO星座卫星天线发射增益和所发射的干扰信号的传输损耗，得到GEO卫星通信地球站接收到的NGSO卫星k的干扰信号幅度函数A_rk(t)，A_rk(t)为随时间t变化的函数，其表达式为：

其中A_rk(t)和A_s取dB为单位，D为NGSO星座卫星天线直径，单位为米，α(t)表示NGSO星座卫星相对于GEO卫星通信地球站的在时间t的仰角，θ为GEO卫星通信地球站天线仰角，

为NGSO卫星天线的3dB波瓣宽度的一半，d_k为NGSO星座卫星k到GEO卫星通信地球站的通信链路长度，单位为km，ω_s为NGSO星座卫星的移动角速度，单位为度每秒，α₀表示NGSO星座卫星相对于GEO卫星通信地球站的初始仰角；

将GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号表示为

其中A_rk(t)表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号幅度随时间变化的函数，

表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号初相位；

S4，构建NGSO星座对GEO卫星通信地球站接收的干扰信号分析模型；

当NGSO星座对GEO卫星通信地球站的可见干扰卫星数量n＝N(l,w,t)时，此GEO卫星通信地球站接收到的所有可见干扰卫星造成的叠加信号R(l,w,t)表示为

其中i表示GEO卫星通信地球站上空的第i颗可见干扰卫星，k(i)表示该第i颗可见干扰卫星在NGSO星座中的卫星编号，n(t)_k(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号时的通信链路噪声，其为高斯白噪声，A_rk(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的幅度，

表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的初相位；

S5，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量在不同取值时，GEO卫星通信地球站接收的干扰信号的均值；

对于GEO卫星通信地球站，将其在一个NGSO星座卫星状态分布周期T_s的时间内接收到的干扰信号均值表示为一维向量

其中，u_p表示可见干扰卫星数量为p时，GEO卫星通信地球站接收到的干扰信号的均值，p＝1,2,...,n(l,w)，其计算公式为

其中，I_pmax表示可见干扰卫星数量为p时的时间片段总数量。

S6，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型；

将GEO卫星通信地球站以其在大地坐标系中经纬度坐标，进行空间上的分割，对每个GEO卫星通信地球站得到一个独立的空间间隔，将一个NGSO星座卫星状态分布周期T_s以不同时间内的GEO卫星通信地球站的可见干扰卫星数量进行时间上的分割，得到若干个时间片段，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型,其表达式为：

其中，u(l,w,t)表示NGSO星座对地球上经纬度为(l,w)的GEO卫星通信地球站在一个NGSO星座卫星状态分布周期内时间片段(t_2i-1,t_2i)上的干扰信号均值。

本发明具有以下优点：

(1)本发明定义了NGSO星座卫星状态分布周期，将NGSO星座对GEO卫星通信地球站干扰分析的限定在一个卫星状态分布周期内，明确了分析需要进行的时长，既能完整地分析整个NGSO星座对GEO卫星通信地球站的干扰情况，又能避免时长外的重复分析；

(2)本发明将GEO卫星通信地球站的可见范围和NGSO星座轨道建模到地心直角坐标系中，通过立体几何表示出可见干扰卫星范围，能够形象、快捷地判断出可见干扰卫星数量；

(3)本发明将GEO卫星通信地球站以其在大地坐标系中的经纬度进行空间上的分割，将一个卫星状态分布周期时间以GEO卫星通信地球站的可见干扰卫星数量进行时间上的分割，能够分析GEO卫星通信地球站在全球不同位置、不同时段的干扰情况。

附图说明

图1为本发明中基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法实现流程图。

具体实施方式

为了更好的了解本发明内容，这里给出一个实施例。

下面结合附图，对本发明进行详细描述。

图1为本发明中基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法实现流程图。

如图1所示，本发明描述了一种基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法，具体步骤如下：

S1，计算非静止轨道星座,即NGSO星座，的卫星状态分布周期；

对于地球表面上任一定点，定义同一NGSO星座在运行过程中出现前后两次的卫星分布完全相同情况的时间差为NGSO星座卫星状态分布周期，记作T_s；对于一个确定的NGSO星座，根据星座的结构和参数，通过计算或者仿真系统仿真可以得到T_s；

当NGSO星座采用倾斜圆轨道星座中的Walker Delta星座时，根据描述其星座结构的Walker代码N/P/F，其中N表示星座的卫星总数，P表示星座的轨道面数量，F表示相位因子，计算NGSO星座卫星状态分布周期，其计算公式为：

其中[，]表示求两者的最小公倍数，T_E表示地球自转周期，T表示NGSO星座的一个轨道周期；

S2，在一个NGSO星座卫星状态分布周期的时间内，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量及各干扰卫星数量的干扰持续时间；

设定GEO卫星通信地球站在大地坐标系中坐标为(l,w,h),其中l表示该地球站的大地经度，w表示该地球站的大地纬度，h表示该地球站的大地高度；

建立地心直角坐标系，原点O设在地球的质量中心，x轴与首子午面与赤道面的交线重合，向东为正，z轴与地球旋转轴重合，向北为正，y轴与xoz平面垂直，三个坐标轴构成右手系；

将GEO卫星通信地球站的大地坐标系坐标转换为地心直角坐标系坐标(x₀,y₀,z₀)，其转换公式为

其中N₀为卯酉圈的半径，

a为地球椭球长半轴，b为地球椭球短半轴，e表示地球椭球曲率；

GEO卫星通信地球站天线方向始终指向GEO卫星，GEO卫星在地心直角坐标系的坐标为(x_S,y_S,0)，将GEO卫星通信地球站的可见范围看作以GEO卫星通信地球站为顶点、以GEO卫星通信地球站到GEO卫星连线为轴、张角为θ的一个圆锥体，其中

此圆锥体的锥面S在地心直角坐标系中的表达式为

其中x_t、y_t、z_t为待计算的中间参数，(x,y,z)为锥面S上任意一点在地心直角坐标系下的坐标；

NGSO星座轨道所在球面在地心直角坐标系中的表达式U为

U:x²+y²+z²＝h_s ²，

其中h_s表示NGSO星座轨道半径，(x,y,z)为NGSO星座轨道所在球面上任意一点在地心直角坐标系下的坐标；

将锥面S的表达式与NGSO星座轨道所在球面的表达式U联立，即得到GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围边界曲线表达式L，记GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围表达式为M；

在一个NGSO星座卫星状态分布周期T_s的时间内，根据NGSO星座轨道参数，将NGSO星座的NGSO卫星k在时间t，在地心直角坐标系的坐标表示为(x_l(t),y_l(t),z_l(t))_k，其中k表示NGSO卫星编号且k∈[1,2,…,N]，t表示时间且0＜t≤T_s；

当NGSO卫星k的坐标(x_l(t),y_l(t),z_l(t))_k∈M时，判定NGSO卫星k为NGSO星座对GEO卫星通信地球站的可见干扰卫星；

根据NGSO星座轨道参数与可见范围M，计算得到位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行期间的最大可见干扰卫星数量，记为n(l,w)；根据NGSO星座在一个NGSO星座卫星状态分布周期T_s的时间内的运行时刻t，得到NGSO星座对GEO卫星通信地球站的在时刻t的可见干扰卫星数量n＝N(l,w,t)，N(l,w,t)为位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行时刻t的可见干扰卫星数量，且N(l,w,t)∈[1，2,…,n(l,w)]；

根据不同时刻的NGSO星座对GEO卫星通信地球站的可见干扰卫星的数量，将一个NGSO星座卫星状态分布周期T_s分割为n(l,w)个种类的时间片段，每一个种类的时间片段构成一个时间片段集合，则可见干扰卫星数量为n时的时间片段集合T_n表达式为：

其中i＝1,2,…,I_nmax，n∈[1，2,…,n(l,w)]，I_nmax表示可见干扰卫星数量为n时的时间片段总数量，(t_2i-1,t_2i)_n表示可见干扰卫星数量为n时的第i个时间片段，t_2i-1表示n个可见干扰卫星的第i次出现时刻，t_2i表示n个可见干扰卫星的第i次出现后的消失时刻；

因此，NGSO星座对GEO卫星通信地球站的各可见干扰卫星数量为n时的持续时间T_total ⁿ表示为：

S3，计算GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号；

将运动到GEO卫星通信地球站可见范围内的一个NGSO星座卫星所发射的干扰信号I_s(t)表示为：

其中A_s表示干扰信号幅度，f为干扰信号频率，

表示干扰信号的初相位，I_s(t)表示t时刻的干扰信号；

根据GEO卫星通信地球站天线接收增益、NGSO星座卫星天线发射增益和所发射的干扰信号的传输损耗，得到GEO卫星通信地球站接收到的NGSO卫星k的干扰信号幅度函数A_rk(t)，A_rk(t)为随时间t变化的函数，其表达式为：

其中A_rk(t)和A_s取dB为单位，D为NGSO星座卫星天线直径，单位为米，α(t)表示NGSO星座卫星相对于GEO卫星通信地球站的在时间t的仰角，θ为GEO卫星通信地球站天线仰角，

为NGSO卫星天线的3dB波瓣宽度的一半，d_k为NGSO星座卫星k到GEO卫星通信地球站的通信链路长度，单位为km，ω_s为NGSO星座卫星的移动角速度，单位为度每秒，α₀表示NGSO星座卫星相对于GEO卫星通信地球站的初始仰角；

将GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号表示为

其中A_rk(t)表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号幅度随时间变化的函数，

表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号初相位；

S4，构建NGSO星座对GEO卫星通信地球站接收的干扰信号分析模型；

当NGSO星座对GEO卫星通信地球站的可见干扰卫星数量n＝N(l,w,t)时，此GEO卫星通信地球站接收到的所有可见干扰卫星造成的叠加信号R(l,w,t)表示为

其中i表示GEO卫星通信地球站上空的第i颗可见干扰卫星，k(i)表示该第i颗可见干扰卫星在NGSO星座中的卫星编号，n(t)_k(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号时的通信链路噪声，其为高斯白噪声，A_rk(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的幅度，

表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的初相位；

S5，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量在不同取值时，GEO卫星通信地球站接收的干扰信号的均值；

对于GEO卫星通信地球站，将其在一个NGSO星座卫星状态分布周期T_s的时间内接收到的干扰信号均值表示为一维向量

其中，u_p表示可见干扰卫星数量为p时，GEO卫星通信地球站接收到的干扰信号的均值，p＝1,2,...,n(l,w)，其计算公式为

其中，I_pmax表示可见干扰卫星数量为p时的时间片段总数量。

S6，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型；

将GEO卫星通信地球站以其在大地坐标系中经纬度坐标，进行空间上的分割，对每个GEO卫星通信地球站得到一个独立的空间间隔，将一个NGSO星座卫星状态分布周期T_s以不同时间内的GEO卫星通信地球站的可见干扰卫星数量进行时间上的分割，得到若干个时间片段，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型,其表达式为：

其中，u(l,w,t)表示NGSO星座对地球上经纬度为(l,w)的GEO卫星通信地球站在一个NGSO星座卫星状态分布周期内时间片段(t_2i-1,t_2i)上的干扰信号均值。

以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (2)

1.一种基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法，其特征在于，其具体步骤包括：

S1，计算非静止轨道星座,即NGSO星座，的卫星状态分布周期；

对于地球表面上任一定点，定义同一NGSO星座在运行过程中出现前后两次的卫星分布完全相同情况的时间差为NGSO星座卫星状态分布周期，记作T_s；对于一个确定的NGSO星座，根据星座的结构和参数，通过计算得到T_s；

S2，在一个NGSO星座卫星状态分布周期的时间内，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量及各干扰卫星数量的干扰持续时间；

S3，计算GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号；

S4，构建NGSO星座对GEO卫星通信地球站接收的干扰信号分析模型；

S5，计算NGSO星座对GEO卫星通信地球站的可见干扰卫星数量在不同取值时，GEO卫星通信地球站接收的干扰信号的均值；

S6，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型；

所述的步骤S2，其具体包括：

设定GEO卫星通信地球站在大地坐标系中坐标为(l,w,h),其中l表示该地球站的大地经度，w表示该地球站的大地纬度，h表示该地球站的大地高度；

建立地心直角坐标系，原点O设在地球的质量中心，x轴与首子午面与赤道面的交线重合，向东为正，z轴与地球旋转轴重合，向北为正，y轴与xoz平面垂直，三个坐标轴构成右手系；

将GEO卫星通信地球站的大地坐标系坐标转换为地心直角坐标系坐标(x₀,y₀,z₀)，其转换公式为

其中N₀为卯酉圈的半径，

a为地球椭球长半轴，b为地球椭球短半轴，e表示地球椭球曲率；

GEO卫星通信地球站天线方向始终指向GEO卫星，GEO卫星在地心直角坐标系的坐标为(x_S,y_S,0)，将GEO卫星通信地球站的可见范围看作以GEO卫星通信地球站为顶点、以GEO卫星通信地球站到GEO卫星连线为轴、张角为θ的一个圆锥体，其中

此圆锥体的锥面S在地心直角坐标系中的表达式为

S：

其中x_t、y_t、z_t为待计算的中间参数，(x,y,z)为锥面S上任意一点在地心直角坐标系下的坐标；

NGSO星座轨道所在球面在地心直角坐标系中的表达式U为

U:x²+y²+z²＝h_s ²，

其中h_s表示NGSO星座轨道半径，(x,y,z)为NGSO星座轨道所在球面上任意一点在地心直角坐标系下的坐标；

将锥面S的表达式与NGSO星座轨道所在球面的表达式U联立，即得到GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围边界曲线表达式L，记GEO卫星通信地球站在NGSO星座轨道所在球面上的可见范围表达式为M；

在一个NGSO星座卫星状态分布周期T_s的时间内，根据NGSO星座轨道参数，将NGSO星座的NGSO卫星k在时间t，在地心直角坐标系的坐标表示为(x_l(t),y_l(t),z_l(t))_k，其中k表示NGSO卫星编号且k∈[1,2,…,N]，t表示时间且0＜t≤T_s；

当NGSO卫星k的坐标(x_l(t),y_l(t),z_l(t))_k∈M时，判定NGSO卫星k为NGSO星座对GEO卫星通信地球站的可见干扰卫星；

根据NGSO星座轨道参数与可见范围M，计算得到位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行期间的最大可见干扰卫星数量，记为n(l,w)；根据NGSO星座在一个NGSO星座卫星状态分布周期T_s的时间内的运行时刻t，得到NGSO星座对GEO卫星通信地球站的在时刻t的可见干扰卫星数量n＝N(l,w,t)，N(l,w,t)为位于经度l、纬度w的GEO卫星通信地球站在NGSO星座运行时刻t的可见干扰卫星数量，且N(l,w,t)∈[1，2,…,n(l,w)]；

根据不同时刻的NGSO星座对GEO卫星通信地球站的可见干扰卫星的数量，将一个NGSO星座卫星状态分布周期T_s分割为n(l,w)个种类的时间片段，每一个种类的时间片段构成一个时间片段集合，则可见干扰卫星数量为n时的时间片段集合T_n表达式为：

其中i＝1,2,…,I_{n max}，n∈[1，2,…,n(l,w)]，I_{n max}表示可见干扰卫星数量为n时的时间片段总数量，(t_2i-1,t_2i)_n表示可见干扰卫星数量为n时的第i个时间片段，t_2i-1表示n个可见干扰卫星的第i次出现时刻，t_2i表示n个可见干扰卫星的第i次出现后的消失时刻；

因此，NGSO星座对GEO卫星通信地球站的各可见干扰卫星数量为n时的持续时间T_total ⁿ表示为：

从而得到各干扰卫星数量的干扰持续时间；

所述的步骤S3，其具体包括：

将运动到GEO卫星通信地球站可见范围内的一个NGSO星座卫星所发射的干扰信号I_s(t)表示为：

其中A_s表示干扰信号幅度，f为干扰信号频率，

表示干扰信号的初相位，I_s(t)表示t时刻的干扰信号；

根据GEO卫星通信地球站天线接收增益、NGSO星座卫星天线发射增益和所发射的干扰信号的传输损耗，得到GEO卫星通信地球站接收到的NGSO卫星k的干扰信号幅度函数A_rk(t)，A_rk(t)为随时间t变化的函数，其表达式为：

其中A_rk(t)和A_s取dB为单位，D为NGSO星座卫星天线直径，单位为米，α(t)表示NGSO星座卫星相对于GEO卫星通信地球站的在时间t的仰角，θ为GEO卫星通信地球站天线仰角，

为NGSO卫星天线的3dB波瓣宽度的一半，d_k为NGSO星座卫星k到GEO卫星通信地球站的通信链路长度，单位为km，ω_s为NGSO星座卫星的移动角速度，单位为度每秒，α₀表示NGSO星座卫星相对于GEO卫星通信地球站的初始仰角；

将GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号表示为

其中A_rk(t)表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号幅度随时间变化的函数，

表示GEO卫星通信地球站接收到的NGSO卫星k发射的干扰信号初相位；

所述的步骤S4，其具体包括：

当NGSO星座对GEO卫星通信地球站的可见干扰卫星数量n＝N(l,w,t)时，此GEO卫星通信地球站接收到的所有可见干扰卫星造成的叠加信号R(l,w,t)表示为

其中i表示GEO卫星通信地球站上空的第i颗可见干扰卫星，k(i)表示该第i颗可见干扰卫星在NGSO星座中的卫星编号，n(t)_k(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号时的通信链路噪声，其为高斯白噪声，A_rk(i)表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的幅度，

表示GEO卫星通信地球站接收第i颗可见干扰卫星发射的干扰信号的初相位；

所述的步骤S5，对于GEO卫星通信地球站，将其在一个NGSO星座卫星状态分布周期T_s的时间内接收到的干扰信号均值表示为一维向量

其中，u_p表示可见干扰卫星数量为p时，GEO卫星通信地球站接收到的干扰信号的均值，p＝1,2,...,n(l,w)，其计算公式为

其中，I_pmax表示可见干扰卫星数量为p时的时间片段总数量；

所述的步骤S6,将GEO卫星通信地球站以其在大地坐标系中经纬度坐标，进行空间上的分割，对每个GEO卫星通信地球站得到一个独立的空间间隔，将一个NGSO星座卫星状态分布周期T_s以不同时间内的GEO卫星通信地球站的可见干扰卫星数量进行时间上的分割，得到若干个时间片段，构建时空域分割的NGSO星座对GEO卫星通信地球站干扰模型,其表达式为：

其中，u(l,w,t)表示NGSO星座对地球上经纬度为(l,w)的GEO卫星通信地球站在一个NGSO星座卫星状态分布周期内时间片段(t_2i-1,t_2i)上的干扰信号均值。

2.如权利要求1所述的基于时空割集的NGSO星座对GEO卫星通信地球站干扰建模方法，其特征在于，所述的对于一个确定的NGSO星座，根据星座的结构和参数，通过计算得到T_s，其具体包括：

当NGSO星座采用倾斜圆轨道星座中的Walker Delta星座时，根据描述其星座结构的Walker代码N/P/F，其中N表示星座的卫星总数，P表示星座的轨道面数量，F表示相位因子，计算NGSO星座卫星状态分布周期，其计算公式为：

其中[，]表示求两者的最小公倍数，T_E表示地球自转周期，T表示NGSO星座的一个轨道周期。

Images