认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度方法

认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度方法
【专利摘要】认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度方法，涉及信息与通信技术领域，为了解决地面组件与卫星组件的同频干扰等问题。本发明所构建的系统模型由一个多波束卫星，多个地面辅助组件和地球站组成，主要是针对点波束内地球站对相邻点波束内的地面用户的上行链路造成干扰的场景分析。为了减少点波束边缘的CGC所受的干扰，本发明采用了禁区原则，利用改进的频分复用方案，对认知星地一体化系统中同频干扰及频谱有效性进行了理论分析与仿真验证，同时，本发明联合考虑频谱效率与干扰抑制两方面，提出了能使系统在可接收的干扰水平下实现最大的频谱效率的最优EA策略。
【专利说明】
认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区 宽度方法
技术领域
[0001] 本发明设及信息与通信技术领域，具体设及在认知星地一体化系统中对既能保证 频谱有效性，又能抑制组件间干扰的最优EA策略的技术。
【背景技术】
[0002] 太空域和信息域是国家经济发展和军事斗争两个重要战略制高点，目前，我国已 经有近百颗的空间飞行器在轨运行，覆盖了通信、导航、中继、深空探测、载人航天等所有领 域，同时配套使用的地面网系统也日益成熟，包括卫星运控站、测控站、上行站、监测站等。 卫星通信已经成为信息化时代不可缺少的通信方式。运种通信方式具有宽频带、大容量、广 覆盖的优势，但也存在阴影衰落严重等问题。
[0003] 针对卫星通信的缺点，地面无线通信网络由于系统容量大，通信质量好，频带利用 率高，手机发射功率低，而且频率规划灵活，可W覆盖城市及室内等等优点，可W作为卫星 通信的必要补充，但是地面通信在偏远农村地区覆盖不彻底，而且在城市环境中，会有严重 的噪声和干扰，卫星通信恰巧可W弥补地面通信的不足。随着卫星通信技术的日益成熟及 地面互联网技术的迅速发展，针对W上提到的卫星与地面通信的优缺点，需要综合考虑两 种通信方式，利用卫星连接远程地面互联网，建立空天地一体化的信息网络。
[0004] 国际电信联盟（ITU)提出了认知星地一体化的概念，认为它是能发挥重要作用的 下一代的网络。此外，在欧洲技术平台肥TW0RLD2020,卫星工作组颁布了关于在5G中卫星作 用的白皮书。卫星网络能够为人口低密度区域提供最佳且最综合的覆盖，而地面网络或者 是地面组件能够为高密度人口区域，如城市环境，提供最高的带宽和最低成本的覆盖，所 W，认知星地一体化网络的设计是有广阔前景的。
[0005] 认知星地一体化系统是由基于多波束空间卫星网络SBN(Satellite Based 化twork)和辅助地面组件CGC(Complementary Ground Components)构成。CGC并不等同于 移动服务运营商所使用的独立的地面基站，它们与卫星组件使用相同的频率资源，并受相 同的网络管理机制制约，与集成的用户设备(肥)通信。CGC也可W复用除所在点波束频率之 外的卫星频率，CGC可W增强卫星由于严重阴影衰落而难W覆盖的区域，特别是在高楼密集 的城市。
[0006] 在实际通信过程中，认知星地一体化网络将会带来很多优势，但也会面临很多挑 战，特别是干扰问题。无线通信的快速发展，增加了对频谱资源的需求，频谱稀缺的问题变 得日益严重。由于频谱资源的匿乏，在卫星和地面网络中利用频分复用技术，能够有效改善 频谱的有效性，扩大系统容量。然而，当卫星网络和辅助地面组件CGCW相同的频率同时发 送信息时，将会导致同频干扰的问题，如何在满足系统要求的干扰水平下，使频谱利用率达 到最优是一体化的卫星地面通信系统亟待解决的关键问题。
[0007] 近几年来，国内外有很多研究设及到认知星地一体化网络中干扰模型建立及频率 复用策略。Deslandes等人构造了一个由多波束卫星、地面基站W及混合的地面和卫星用户 终端组成的模型，主要分析了地面用户上行信号干扰卫星用户的上行信号的情形。Kang等 人构造了一个由多波束卫星、地面辅助组件(CGCs)和移动地球站组成的干扰模型主要针对 地面组件CGC对卫星上行链路的干扰进行分析。Deslandes、KangW及化e-sub化等人阐述 了禁区EZ的概念，它是位于蜂窝边缘的一个特定的区域，其中的地面基站不能使用同蜂窝 和邻近蜂窝的频率波段，有效的改善了系统的性能。
[0008] 但是在认知星地一体化系统中，存在很多同频干扰的场景，除上面提到的两种干 扰场景外，还有地球站对相邻波束内的地面用户的上行链路造成干扰、CGC的下行信号链路 干扰卫星下行信号等很多场景，而目前的大多数研究都集中对地面用户上行信号干扰卫星 用户上行信号的场景分析，而很少设及到其他干扰场景，同时在认知星地一体化系统中也 缺乏对所有干扰场景的概括分析，此外，对于上述研究所提到的禁区，目前的研究只是停留 在仿真验证EA可W减少系统所受的干扰水平，而且会降低系统的频谱利用率，对如何在权 衡系统可接受的干扰水平与频率利用率的前提下，在不同的频率复用下寻求最佳的EA取值 使系统的性能最优而缺乏研究。

【发明内容】

[0009] 本发明是为了解决目前认知星地一体化系统中在频谱效率及干扰估计方面存在 的W下问题：
[0010] (1)、干扰场景单一，缺乏对其他干扰场景的建模研究；
[0011] (2)、缺乏在多种频率复用因子下的干扰估计与频谱有效性的研究；
[0012] (3)、没有寻找到最佳禁区EA的取值使系统性能最优；
[0013] 从而提出了一种认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽 度方法。
[0014] 认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度方法，
[0015] 设一个点波束内只有一个地球站，并且分布于点波束的中屯、，并能够与卫星通信， 同时设在点波束内的CGC是随机分布的，建立系统模型；
[0016] 它包括W下步骤：
[0017] 步骤A、估计组件间的干扰水平的步骤，具体为：
[0018] 步骤Al、设输入输入CGC数量M，点波束半径R，EA宽度值W，地球站的发射功率Pt,Ci， 地球站发射天线增益Gt, Gi，CGC的接收天线增益Gr ;
[0019]步骤A2、根据系统模型计算中屯、波束内M个CGC的位置矩阵Z;
[0020] 步骤A3、根据系统模型计算相邻地球站的位置矩阵ES_pos;
[0021] 步骤A4、根据EA宽度值W值，将CGC划分为屯决，得到位置矩阵Ri和元素个数m，i = 0,1,...,6;
[0022] 步骤A5、分割卫星频段成n部分{fi}，i = l，2, . . .，n，设外层循环迭代数的初始值i =0;
[0023] 步骤A6、等分Ri成n-1部分，分别分配频段。，j = l，2,. . .，n，设内层循环迭代数的 初始值P = I;
[0024] 当频率复用因子n = 7且P如0时，执行W下步骤：
[00巧]步骤A7、对Ri内每个坐标Z (P )，定位唯一同频地球站，计算距离d和路径损耗Lc;
[00%]步骤A8、根据公式：
[OO^ ] I j = Pt, ci+Gt, ci+G 广 Lci (dBW)
[002引计算每个CGC受到的干扰功率，存入矩阵10,令p = p+l;
[0029] 式中：j是CGC在中屯、点波束内的标号，i是频段的标号，l<i<6;Pt,Cl是指在相邻 的工作频率是fi的点波束内地球站的发射功率;Gt,Cl是指在相邻的工作频率是fi的点波束 内地球站指向第n个CGC的发射天线增益瓜是指CGC的接收天线增益;Lci是指从中屯、波束内 工作在fi的第j个CGC到相邻点波束同频地球站的路径损耗；
[0030] 步骤A9、当p〉m时，退出内层循环，令i = i+l;
[0031] 当i< 6时执行W下步骤：
[0032] 步骤AlO、等分Ri成n-2部分，分别分配频段扣j = 1，2，. .，5，设内层循环迭代数P = 1;
[0033] 当p<ni时，返回执行步骤A7:
[0034] 当i〉6时，执行W下步骤：
[00对步骤All、退出外层循环，得至IjM个CGC所受到的干扰功率{lm}，m=0，l，. . .，6;
[0036] 步骤A12、绘制复用因子为n时的干扰水平CDF曲线；
[0037] 步骤B、估计组件间的干扰水平的步骤，具体为：
[0038] 步骤B1、输入M，R，复用因子n，令i = l，W=0;
[0039] 当￥<則寸，执行W下步骤：
[0040] 步骤B2、根据W取值，把中屯、波束内的CGC分割成屯块，位置矩阵分别是Rm，m = 0， 1，...，6,得到区域Ri内CGC的数量m;
[0041 ]步骤B3、把工作频段划分成n段，计算子频带带宽Bsub, Nin, Nout;
[0042] 步骤B4、根据公式：
[0043]
[0044]
[0045] 式中，Nin和Nout分别指中屯、波束在EA内部和外部的CGC数量，n频率复用因子，Bsub是 指子频带的带宽，Ball是指系统使用的全部的频带带宽；
[0046] 步骤B5、循环停止后，绘制复用因子为n时的标准有效带宽Bavaii随W的变化曲线；
[0047] 步骤C、寻求最佳禁区化A)宽度的实现步骤：
[004引步骤Cl、输入M，R，链路有效性La，复用因子n;
[00例步骤C2、令i = l，W=0,计算p = [La -M+1]，并向负无穷方向取整；
[(K)加 ]当W < R时，执行W下步骤：
[0051 ]步骤C3、计算复用因子n，EA宽度值是W的所有CGC的干扰功率，存储于矩阵I;
[0052]步骤C4、对矩阵I中的元素由小到大排序，存储到矩阵I' ；
[00对步骤C5、取矩阵I'的第P个元素，存储于矩阵laccept(i)，令i = i + l，W = W+1500;
[0054] 步骤C6、循环停止后，得到在不同的La与n下Iaccept随W的变化曲线；取曲线转折点 处对应的W值，即：为最佳的禁区EA宽度。
[0055] 当复用因子n为3或4时，步骤A7和A8分别由步骤A71至步骤A73替换：
[0056] 步骤A71、对Ri内坐标为Z (P)的CGC，定位同频地球站个，计算距离d和路径损耗Lc，m =1;
[0化7] 当m < q时，循环执行步骤A72:
[0化引步骤A72、根据公式：
[0化9]
[0060] 计算坐标为Z(P)的CGC受到的第m个地球站的干扰功率，存入矩阵Km)，m=m+l;
[0061] 式中：q是指与中屯、波束第j个CGC同频的地球站的个数;下标Cip是指工作在fi频段 的第P个相邻点波束；
[0062] 当m〉q时，执行W下步骤：
[0063] 步骤A73、对矩阵I按行求和，得到坐标为Z(P)的CGCs受到的合计干扰功率，存入矩 阵 Ii，p = p+1;
[0064] 执行完W上S步后，执行步骤A9。
[0065] 本发明具有W下特点和显著进步：
[0066] 1、本发明概括说明了在认知星地一体化系统中可能存在的干扰场景，并对目前较 少研究的一种场景进行建模分析；
[0067] 2、相比于传统的频率复用方案，本发明所提出的改进的频率复用方案更能有效的 减弱CGC所受的干扰水平；
[0068] 3、针对目前大多数研究只在某一复用因子下进行，本发明研究了复用因子分别是 3，4，7时的干扰水平与频谱有效性随EA取值的变化情况，提供了更多的参考；
[0069] 4、本发明在兼顾系统频谱效率与可接受的干扰水平的前提下，在不同的复用因子 场景下，相应找到了使系统性能最优的EA宽度值W。
【附图说明】
[0070] 图1是认知星地一体化系统中组件间的干扰场景示意图；
[0071] 图2是卫星组件与CGC之间改进的频分复用方案的原理示意图；
[0072] 图3是中屯、波束内的CGC受到合计干扰的场景示意图；
[0073] 图4是计算中屯、波束内所有CGC所受的干扰水平的流程示意图；
[0074] 图5是在不同W下估计系统的频谱有效性的流程示意图；
[0075] 图6是仿真可接受干扰水平随W的变化情况的流程示意图；
[0076] 图7是复用因子是3时CGC所受到的干扰仿真示意图；
[0077] 图8是复用因子是4时CGC所受到的干扰仿真示意图；
[0078] 图9是复用因子是7时CGC所受到的干扰仿真示意图；
[0079] 图10是EA的宽度W对EA内CGC数量的影响仿真示意图；
[0080] 图11是S种频率复用因子下，系统标准有效带宽随W的变化情况仿真示意图；
[0081] 图12至图14是不同的复用因子及链路有效性下，可接受干扰水平随W的变化情况 仿真不意图；
【具体实施方式】
[0082]
【具体实施方式】、认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度 方法，认知星地一体化系统中可能存在的组件间的干扰场景如图I所示。
[0083] 图1中的颜色标识说明了CGC可W复用卫星频率。与地球站的发射功率相比，本发 明忽略了卫星用户的存在。场景①描述的是点波束内的地球站对相邻波束内的地面用户的 上行链路造成干扰，运会导致CGC接收的信号质量下降。由于点波束的覆盖范围远大于CGC 的覆盖范围，在一个点波束内会存在数W万计的同频CGC，位于点波束边缘的CGC将会遭受 更严重的干扰问题，因此同频的地球站对CGC的干扰不容忽视。场景①是本发明主要考虑的 情况。
[0084] 场景②描述的是CGC的下行信号链路干扰了卫星到地球站的下行信号。场景③描 述的是地面用户对相邻点波束内的地球站的上行信号造成干扰。在运种场景下，可W区分 两种干扰情况:单个实体的干扰问题（即当一个地面用户与卫星之间存在LOS且处于点波束 的边缘区域时，对地球站进行干扰)和合计干扰问题（即在点波束内所有同频的地面用户对 地球站的干扰）。场景④描述的是卫星的下行信号干扰了 CGC的下行信号。应当注意到在地 面蜂窝系统中有一个大的动态范围，因此通过稍微减小CGC的蜂窝半径，可W减轻干扰强 度。场景⑤和⑥描述的是地面蜂窝系统中的同频干扰问题。通过使用合适的频率复用策略 和适当调整CGC的发射功率，可W减轻地面系统的同频干扰问题。
[0085] 针对图1中干扰场景①中点波束边缘的CGC将会受到相邻点波束同频的ES的严重 干扰问题，对其进行建模，采取了禁区化A)的原则，在运一区域的CGC不可W使用所在点波 束与相邻点波束的频率。图2W频分复用因子是7为例，说明了卫星组件与CGC之间改进的频 分复用方案。图中点波束内的白色区域表示禁区，W是指EA的宽度，0 < W < R。
[0086] 本发明假设一个点波束内只有一个地球站，并且分布于点波束的中屯、与卫星通 信，同时假设在点波束内的CGC是随机分布的。建立了系统模型后，本发明在频率复用因子 分别是3，4，7时，进行组件间的干扰估计、频谱有效性分析W及寻找使系统性能最优的最佳 EA宽度W。具体实现步骤分S部分叙述如下：
[0087] (1)估计组件间的干扰水平由W下步骤实现：
[0088] 步骤一、输入CGC数量M，点波束半径R，EA宽度值W，地球站的发射功率Pt,Ci，地球站 发射天线增益Gt, Gi，CGC的接收天线增益Gr ;
[0089] 步骤二、从模型中计算中屯、波束内M个CGC的位置矩阵Z;
[0090] 步骤=、从模型中计算相邻地球站的位置矩阵ES_pos;
[0091] 步骤四、根据W值，划分CGC成屯决，得到位置矩阵Ri和元素个数m，i = 0，l，...，6;
[0092] 步骤五、分割卫星频段成n部分{fi}，i = l，2, . . .，n，迭代数数i = 0;
[0093] 步骤六、等分Ri成n-1部分，分别分配频段扣j = l，2,. . .，n，迭代数p = l;
[0094] 当频率复用因子11 = 7且9<11〇时，执行W下步骤：
[00M] 步骤屯、对Ri内每个Z(P)，定位唯一同频地球站，计算距离d和路径损耗Lc;
[0096] 步骤八、根据公式（1)计算每个CGC受到的干扰功率，存入矩阵Io，p = p+l;
[0097] 步骤九、当p〉m时，退出内层循环，i = i+l;
[009引当i < 6时执行W下步骤：
[0099] 步骤十、等分Ri成n-2部分，分别分配频段fj，j = l，2,. .，5,迭代数P = I;
[0100] 当P如财返回执行步骤屯：
[0101] 步骤十一、当i〉6时，退出外层循环，得到M个CGC所受到的干扰功率{lm}，m = 0， 1,---,6;
[0102] 步骤十二、画出复用因子是n时的干扰水平CDF曲线；
[0103] 当复用因子n是3或4时，不能再用公式（1)计算，因为对于每个CGC受到的不是单一 地球站的干扰，而是多个地球站的合计干扰，应用公式（3)来计算。上述步骤屯和八应更改 为：
[0104] 步骤一、对Ri内坐标为Z (P)的CGC，定位同频地球站q个，计算距离d和路径损耗Lc，m =1;
[01化]当m< q时，循环执行步骤二：
[0106] 步骤二、根据公式(3)计算坐标为Z(P)的CGC受到的第m个地球站的干扰功率，存入 矩阵I(m)，m=m+l;
[0107] 当m〉q时，执行W下步骤：
[0108] 步骤S、对矩阵I按行求和，得到坐标为Z(P)的CGCs受到的合计干扰功率，存入矩 阵 Ii，p = p+1;
[0109]执行完W上二步后，继续执行步骤九。
[0110] (2)估计系统的频谱有效性的实现步骤：
[0111 ] 步骤一、输入M，R，复用因子n，令i = l，W=0;
[0112]当W <則寸执行W下步骤：
[011引步骤二、根据W取值，把中屯、波束内的CGC分割成屯决，位置矩阵分别是Rm, m = 0， 1，...，6,得到区域Ri内CGC的数量m;
[0114] 步骤S、把工作频段划分成n段，计算子频带带宽Bsub，Nin，Nout;
[0115] 步骤四、根据公式（5)计算Bavaii，i = i+l，W=W+1500;
[0116] 步骤五、循环停止后，画出复用因子是n时的标准有效带宽Bavail随W的变化曲线。
[0117] (3)寻求最佳EA宽度的实现步骤：
[011引步骤一、输入M，R，链路有效性La，复用因子n;
[0119] 步骤二、令i = i，w=〇,计算p = [La ? M+1](取整）；
[0120] 当W< R时执行W下步骤：
[0121] 步骤S、计算复用因子n，EA宽度值是W的所有CGC的干扰功率，存储于矩阵I;
[0122] 步骤四、对矩阵I中的元素由小到大排序，存储到矩阵I'；
[0123] 步骤五、取矩阵I'的第P个元素，存储于矩阵Iaeeept( i )，i = i + 1，W = W+1500 ;
[0124] 步骤六、循环停止后，得到在不同的La与n下Iaccept随W的变化曲线。取曲线转折点 处对应的W值，即为最佳的EA宽度。
[0125] W上是该发明具体的实现步骤，下面从理论角度进行阐述分析。
[0126] 在频率复用因子为7的场景中，本发明W工作频段是f7的点波束作为中屯、波束，来 分析其内的CGC受到相邻点波束内的ES的干扰情况。如果令中屯、波束内工作在fi的第j个 CGC所受到的干扰为U那么。可W由下式计算：
[0127] Ij = Pt,ci+Gt'ci+G^Lci(dBW) (1)
[012引其中，j是CGC在中屯、点波束内的标号，i是频段的标号，1 < i <6。口*,Ci是指在相邻 的工作频率是fi的点波束内地球站的发射功率;Gt,Cl是指在相邻的工作频率是fi的点波束 内地球站指向第n个CGC的发射天线增益瓜是指CGC的接收天线增益;Lci是指从中屯、波束内 工作在fi的第j个CGC到相邻点波束同频地球站的路径损耗。
[0129] 在本发明中，为了简化仿真程序，假设地球站的仰角是高于48%此时用于干扰计 算的地球站朝向地平线方向的离轴天线增益Gt,Cl是一个常数，为-10地i JTU-R管理规范也 对地面站的EIRP进行了限定，W保护其他的地面站，在本发明仿真中研究最糟糕的情况，即 使用极限值40地W。
[0130] 在估计中，使用的传输模型是自由空间损耗，Lci可由下式计算得到：
[0131] Lci = 32.45+201g(dci)+201g(fi)(dB) (2)
[01创在上式中，dci是指中屯、波束内工作在fi的第j个CG巧师邻点波束同频地球站的距 离。du单位是km, f i单位是MHz。
[0133] 当频率复用因子n是3或4时，中屯、波束内的第j个CGC的干扰不能按照公式（1)计 算。因为此种情况下，第j个CGC不是受单个地球站的干扰，而是多个地球站的合计干扰。如 图3所示:假设中屯、波束的CGC工作在f2频段，则此时会同时受到相邻立个点波束C2、C4和Cs 的同频地球站的干扰。此时CGC所受到干扰水平是=个干扰值之和。
[0134] 由此，得到此时干扰水平的计算公式：
(3>
[0135]
[0136] 其中，q是指与中屯、波束第j个CGC同频的地球站的个数。下标Cip是指工作在fi频段 的第P个相邻点波束。
[0137] 假设中屯、波束内CGC的总数是M，为了得到M个CGC的干扰水平分布规律，W及得到 复用因子取值与干扰水平的关系，定义CGC的累积概率分布函数：
[01 測
（4)
[0139] 其中，e(t)是阶跃函数。
[0140] 为了在不同的复用因子下估计系统的频谱有效性，通过定义标准有效带宽来衡量 系统的频谱有效性：
(5)
[0141]
[01创上式中，Nin和Nout分别指中屯、波束在EA内部和外部的CGC数量，n频率复用因子，Bsub 是指子频带的带宽，Ball是指传统模式下中屯、波束内每个CGC允许使用的频带资源的总和。
[0143] 在联合考虑干扰水平与频谱有效性的同时，为了得到最佳的EA值，需要考虑系统 可接受干扰水平与EA大小的关系，定义链路有效性La:
[0144] (6)
[0145] 其中=Iaccept是指系统可接受的干扰水平，。是指中屯、波束内第j个CGC所受到的干 扰，N是中屯、波束内CGC的总数。
[0146] 本发明在Matlab平台上进行了仿真验证，仿真参数为：
[0147]
[0148] 具体的仿真结果如下：
[0149] (I)S种频率复用因子下，EA的宽度W对CGC所受干扰水平分布的影响，参见图7至 9。
[0150] 从仿真结果图可W看出，随着EA的增大，干扰水平的CDF曲线越快接近1，对中屯、波 束内CGC的干扰水平在不断减小。运是因为，在EA内的CGC不会接受来自最近的相邻波束的 ES的干扰。比如W频率复用因子是7时为例，在传统模式下，即W=O时，CGC接收的干扰值不 大于-IHdBW的概率是72.7%，而在改进模式中，当胖=20缸时，其概率是78%，当胖=40缸 时，其概率是82.5 %，当W=60km时，其概率值是85 %。
[0151] 比较图7至图9,也可W发现随着频率复用因子的增加，干扰水平在不断减小。Ww =20km为例，当干扰水平的CDF取0.別寸，在复用因子是3时，对应的干扰水平是-346.3dB，当 复用因子是4时，其值是-231.6地，当复用因子是7时，其值是-113.8地。
[0152] (2化A的宽度W对EA内CGC数量的影响，参见图10。
[0153] (3)S种频率复用因子下，系统标准有效带宽随W的变化情况，参见图11。
[0154] 从仿真结果图可W看出，系统标准有效带宽随W的变化趋势与CGC数量随EA的变化 趋势恰好相反，运是因为EA中的CGC不允许使用所在点波束和相邻点波束的频段，能够使用 除所在点波束之外的频率的CGC就越少，由图10曲线可W看出，提高频率复用因子，可W增 大频率利用率，同时W越大，系统的标准有效带宽越小，即频率利用率越小，而由图8可知EA 增大可W减小CGC所受到的干扰。因此在决定EA宽度W取值时，应充分权衡系统可接受的干 扰水平与频率利用率两个因素。
[0155] (4)不同的复用因子及链路有效性下，可接受干扰水平随W的变化情况，参见图12 至图14。
[0156] 由仿真结果图可W看出，随着W增大，系统可接受的干扰水平在不断减小，曲线前 半部分下降趋势明显，而后半部分则趋于平缓。运说明，可W通过选择一个W值，使系统可接 受的干扰水平维持稳定，W如果取值过大，将会造成系统频谱利用率下降（图12已经证明）， 应该联合考虑系统可接受的干扰水平与频谱有效性，可W得出在曲线由急剧下降到平缓的 转折点处所对应的W值就是最优的EA宽度。W图14为例具体说明。当要求系统的链路有效性 是99 %时，可取W = 37km，当链路有效性是95 %、90 %和85 %时，EA的宽度分别可设置成W = 45km、58km和68km。
[0157] 经上述仿真试验验证，本发明具有W下特点和显著进步：
[0158] 1、本发明概括说明了在认知星地一体化系统中可能存在的干扰场景，并对目前较 少研究的一种场景进行建模分析；
[0159] 2、相比于传统的频率复用方案，本发明所提出的改进的频率复用方案更能有效的 减弱CGC所受的干扰水平；
[0160] 3、针对目前大多数研究只在某一复用因子下进行，本发明研究了复用因子分别是 3，4，7时的干扰水平与频谱有效性随EA取值的变化情况，提供了更多的参考；
[0161] 4、本发明在兼顾系统频谱效率与可接受的干扰水平的前提下，在不同的复用因子 场景下，相应找到了使系统性能最优的EA宽度值W。
【主权项】
1.认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区宽度方法，其特征是： 设一个点波束内只有一个地球站，并且分布于点波束的中心，并能够与卫星通信，同时 设在点波束内的辅助地面组件(CGC)是随机分布的，建立系统模型； 该方法包括以下步骤： 步骤A、估计组件间的干扰水平的步骤，具体为： 步骤A1、设输入输入CGC数量Μ，点波束半径R，EA宽度值W，地球站的发射功率Pt, Ci，地球 站发射天线增益Gt, Ci，CGC的接收天线增益Gr; 步骤A2、根据系统模型计算中心波束内Μ个CGC的位置矩阵z; 步骤A3、根据系统模型计算相邻地球站的位置矩阵ES_pos; 步骤A4、根据EA宽度值W值，将CGC划分为七块，得到位置矩阵Ri和元素个数ru，i = 0， 1，· · ·，6; 步骤A5、分割卫星频段成η部分{h}，i = l，2,. . .，n，设外层循环迭代数的初始值i = 0; 步骤A6、等分1^成11-1部分，分别分配频段f」，j = l，2, . . .，n，设内层循环迭代数的初始 值P = 1; 当频率复用因子η = 7且p < no时，执行以下步骤： 步骤A7、对Ri内每个坐标z (p)，定位唯一同频地球站，计算距离d和路径损耗L。； 步骤A8、根据公式： Ij = Pt,ci+Gt>ci+Gr-Lci(dBff) 计算每个CGC受到的干扰功率，存入矩阵Ιο，令p = p+l; 式中：j是CGC在中心点波束内的标号，i是频段的标号，1 < i < 6;Pt,Ci是指在相邻的工作 频率是fi的点波束内地球站的发射功率;Gt, Ci是指在相邻的工作频率是fi的点波束内地球 站指向第η个CGC的发射天线增益;Gr是指CGC的接收天线增益;Ui是指从中心波束内工作在 h的第j个CGC到相邻点波束同频地球站的路径损耗； 步骤A9、当p>m时，退出内层循环，令i = i+l; 当6时执行以下步骤： 步骤A10、等分心成11-2部分，分别分配频段f」，j = 1，2，. .，5，设内层循环迭代数p = 1; 当p < m时，返回执行步骤A7: 当i >6时，执行以下步骤： 步骤All、退出外层循环，得到Μ个CGC所受到的干扰功率{Im}，m=0，l，...，6; 步骤A12、绘制复用因子为η时的干扰水平CDF曲线； 步骤B、估计组件间的干扰水平的步骤，具体为： 步骤B1、输入M，R，复用因子n，令i = l，W=0; 当W5R时，执行以下步骤： 步骤B2、根据W取值，把中心波束内的CGC分割成七块，位置矩阵分别是Rm，m = 0，l，...， 6，得到区域Ri内CGC的数量m; 步骤B3、把工作频段划分成η段，计算子频带带宽Bsub，Nin，N?t; 步骤B4、根据公式：计算 Bavaii，令i = i+l，W=W+1500; 式中，NiJPNcmt分别指中心波束在EA内部和外部的CGC数量，η频率复用因子，Bsub是指子 频带的带宽，Ball是指传统模式下中心波束内每个CGC允许使用的频带资源的总和； 步骤B5、循环停止后，绘制复用因子为η时的标准有效带宽Bavall随W的变化曲线； 步骤C、寻求最佳禁区(EA)宽度的实现步骤： 步骤C1、输入Μ，R，链路有效性La，复用因子η; 步骤〇2、令1 = 1，￥=0，计算口=[1^*1+1]，并向负无穷方向取整； 当W5R时，执行以下步骤： 步骤C3、计算复用因子η，ΕΑ宽度值是W的所有CGC的干扰功率，存储于矩阵I; 步骤C4、对矩阵I中的元素由小到大排序，存储到矩阵Γ ; 步骤C5、取矩阵Γ的第ρ个元素，存储于矩阵Iac^pt(i)，令i = i+1，W=W+1500; 步骤C6、循环停止后，得到在不同的Lgn下13。。_随1的变化曲线;取曲线转折点处对应 的W值，即：为最佳的禁区EA宽度。2. 根据权利要求1所述的认知星地一体化系统中联合频谱效率及干扰抑制的最优禁区 宽度方法，其特征在于当复用因子η为3或4时，步骤A7和A8分别由步骤A71至步骤A73替换： 步骤Α71、对Ri内坐标为z (ρ)的CGC，定位同频地球站的数量q，计算距离d和路径损耗L。， m=l; 当q时，循环执行步骤A72: 步骤A72、根据公式：计算坐标为z (ρ)的CGC受到的第m个地球站的干扰功率，存入矩阵I(m)，m=m+1; 式中：q是指与中心波束第j个CGC同频的地球站的个数;下标C1P是指工作在6频段的第 P个相邻点波束； 当m>q时，执行以下步骤： 步骤A73、对矩阵I按行求和，得到坐标为z(p)的多个地面辅助组件(CGCs)受到的合计 干扰功率，存入矩阵Ii，P = P+l; 执tx完以上二步后，执彳丁步骤A9。3. 根据权利要求1所述的根据权利要求1所述的认知星地一体化系统中联合频谱效率 及干扰抑制的最优禁区宽度方法，其特征在于该方法的应用场景为：点波束内的地球站对 相邻波束内的地面用户的上行链路造成干扰，导致CGC接收的信号质量下降的场合。
【文档编号】H04Q11/02GK105827301SQ201610231526
【公开日】2016年8月3日
【申请日】2016年4月14日
【发明人】贾敏, 蔺萍, 刘晓锋, 郭庆, 顾学迈
【申请人】哈尔滨工业大学