CN109547075B - 一种多波束卫星通信系统用户分组的方法及其装置 - Google Patents

一种多波束卫星通信系统用户分组的方法及其装置 Download PDF

Info

Publication number
CN109547075B
CN109547075B CN201811569312.9A CN201811569312A CN109547075B CN 109547075 B CN109547075 B CN 109547075B CN 201811569312 A CN201811569312 A CN 201811569312A CN 109547075 B CN109547075 B CN 109547075B
Authority
CN
China
Prior art keywords
user
group
users
grouping
executing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811569312.9A
Other languages
English (en)
Other versions
CN109547075A (zh
Inventor
戚晨皓
陈华健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201811569312.9A priority Critical patent/CN109547075B/zh
Publication of CN109547075A publication Critical patent/CN109547075A/zh
Application granted granted Critical
Publication of CN109547075B publication Critical patent/CN109547075B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Radio Relay Systems (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种多波束卫星通信系统用户分组的方法及其装置，属于无线通信领域。其方法包括初始化系统参数，初始化用户小组，并为每个小组分配初始用户，按照设定的用户分组选择标准为用户分配小组，输出用户小组集合，设计预编码矩阵。本发明针对卫星服务的地面用户进行分组，确定了卫星每个时隙服务的用户，能够大大提升了系统和速率，并且可以服务更多的用户。

Description

一种多波束卫星通信系统用户分组的方法及其装置
技术领域
本发明涉及一种多波束卫星通信系统用户分组的方法及其装置，属于无线通信技术领域。
背景技术
随着智能设备的普及和地面移动通信速率的不断提升，卫星通信作为地面移动通信的延伸和补充，亟需提高数据传输速率以迎合更多的应用场景。多输入多输出(Multiple-Input Multiple-Output,MIMO)技术在地面的成功应用，促使研究人员积极探索MIMO技术在卫星通信中的适用性。多波束卫星通信系统是MIMO在卫星系统应用中的最成功的方案之一，目前商用的卫星系统均采用多波束卫星的方案。多波束卫星通信通过卫星的多波束天线形成几十甚至几百个能量高的窄点波束，覆盖地面的一定区域，实现卫星与该区域用户间的广播或多播通信服务。
为了支持太比特通信，多波束卫星通信系统需要使用尽可能大的用户带宽，这使得全频复用成为必然，但相邻波束间的干扰严重制约了系统性能。因此，需要采用在信号传输之前进行预编码。此外，由于地面用户数目远远大于卫星天线的数目，卫星在每个时隙只能选择一部分用户传输数据，用户分组的技术是十分必要的，它直接影响力整个卫星系统的传输速率表现。
发明内容
本发明的目的在于克服上述多波束卫星通信系统的不足，提供一种在信号传输之前进行预编码的多波束卫星通信系统用户分组的方法及其装置。
本发明的目的是这样实现的
本发明一种多波束卫星通信系统用户分组的方法及其装置，其根据卫星服务用户场景的不同，对卫星服务的地面用户进行分组，分为两个方案：方案一只选择大部分用户被服务的情况，以实现小组最大化和速率，考虑的是最大化原则；方案二选择所有用户都被服务的情况，采用公平性原则。
本发明提供的方案一：一种多波束卫星通信系统用户分组的方法，其包括以下步骤：
步骤1、初始化系统参数，所述系统参数包括用户集合、小组总数和小组循环次数；
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤9；
步骤3、选择用户集合中用户信道向量的模值最大的用户，初始化用户小组，并更新用户集合，所述小组循环次数加1，计算当前小组和速率；
步骤4、判断小组用户个数是否小于等于波束个数，如果满足条件，执行步骤5，否则执行步骤2；
步骤5、根据用户分组选择标准，为上述完成初始化的小组集合分配合适的用户；
步骤6、将所选合适的用户加入当前小组构成临时小组，计算临时小组和速率；
步骤7、判断临时小组和速率是否大于当前小组和速率，如果满足条件，则执行步骤8，否则执行步骤4；
步骤8、更新当前小组集合和当前小组和速率，小组循环次数加1，执行步骤2；
步骤9、输出最后的用户分组集合；
步骤10、输出每个时隙的预编码矩阵。
本发明步骤5中，所述用户分组选择标准包括采用最大化和速率方法，所述合适的用户：
Figure BDA0001915067960000021
其中，
Figure BDA0001915067960000022
表示
Figure BDA0001915067960000023
小组中用户信道向量组成的信道矩阵。
本发明步骤5中，所述用户分组选择标准包括采用最大化信道容量方法，所述合适的用户：
Figure BDA0001915067960000024
其中，初始化矩阵A＝γI，用户m对应的信道向量hm，每一次迭代更新为
Figure BDA0001915067960000025
本发明提供的方案二：一种多波束卫星通信系统用户分组的方法，其包括以下步骤：
步骤1、初始化系统参数，所述系统参数包括用户集合、小组总数和小组循环次数；
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤4；
步骤3、选择用户集合中用户信道向量的模值最大的用户，初始化用户小组，并更新用户集合，所述小组循环次数加1，执行步骤2；
步骤4、判断用户集合是否为空集，如果满足条件，执行步骤9，否则执行步骤5；
步骤5、选择用户集合中的用户和相对应的信道向量，初始化小组循环次数；
步骤6、判断小组循环次数是否小于等于小组总数，如果满足条件，执行步骤7，否则执行步骤8；
步骤7、计算当前小组加入新用户增益，小组循环次数加1，执行步骤6；
步骤8、更新增益最大的小组，更新用户集合，执行步骤4；
步骤9、输出最后的用户分组集合；
步骤10、输出每个时隙的预编码矩阵。
本发明步骤7中，所述当前小组加入新用户后所得增益，包括第t小组加入新用户所得的增益V(t)，即
Figure BDA0001915067960000031
其中，
Figure BDA0001915067960000032
为第t用户小组为
Figure BDA0001915067960000035
对应信道矩阵。
本发明步骤7中，步骤7中，所述当前小组新用户后所得增益为第t小组加入新用户所得的增益V(t)，即
Figure BDA0001915067960000033
其中，初始化矩阵A(t)＝γI，hm为用户m对应的信道向量。
将用户m加入新用户所得增益最大的小组
Figure BDA0001915067960000034
每一次迭代更新矩阵为
Figure BDA0001915067960000041
其中，增益最大的小组序号为t*＝max1≤t≤TV(t)。
本发明一种多波束卫星通信系统用户分组的装置，其特征在于，其包括：
初始化模块，用于初始化用户集合、小组总数和小组循环次数；
小组初始化模块，用于选择用户集合中模值最大的用户，初始化小组集合；
小组分配模块，用于按照设定的用户分组选择标准，为完成初始化的小组集合分配合适的用户；
小组输出模块，用于输出最后的用户分组；
预编码设计模块，用于设计每个时隙的预编码矩阵。
其中，所述小组分配模块包括：
用户选择模块，用于选择用户集合中的用户和相对应的信道向量；
计算模块，用于计算用户加入小组后带来的增益；
选择更新模块，用于选择合适的用户，并且更新小组集合和用户集合。
进一步地，所述用户分组选择标准包括选择大部分用户被服务，实现小组最大化和速率。
进一步地，所述用户分组选择标准包括选择所有用户都被服务，实现小组最大化增益。
有益效果
本发明充分利用信道信息，针对卫星服务的地面用户场景的不同进行分组，提出了最大化系统和速率的方案和服务所有用户的方案，以确定卫星每个时隙服务用户，相比于现有的用户分组方案，大大提升了系统和速率，并且可以服务更多的用户。
附图说明
图1是本发明所适用的多波束卫星通信系统的结构示意图；
图2是本发明多波束卫星通信系统的用户分组方法的装置的结构示意图；
图3A是本发明多波束卫星通信系统的用户分组方法的实施例一的流程图；
图3B是本发明多波束卫星通信系统的用户分组方法的实施例二的流程图；
图4A和图4B是本发明与现有方案在不同发射功率下系统和速率的变化对比图；
图5A和图5B是本发明与现有方案在不同时隙数目下系统和速率的变化对比图；
图6是本发明与现有方案在不同时隙数目下服务用户总数的对比变化图。
具体实施方式
下面结合实施例和说明书附图对本发明做详细的说明，但不限于此。
本发明适用于多波束卫星通信场景，其系统结构如图1所示。假设一个多波束卫星通信系统包含一颗宽带卫星和K个地面用户。卫星采用单馈源单波束架构，反射面将N个馈源信号转换N个发送信号，每个信号在地面覆盖一片区域形成一个波束。我们使用全频复用技术，以提升频谱效率。
假设K个用户均匀分布在所有波束中，每个波束中都会有多个用户，因此有N＜＜K。为了能够同时服务多个用户，多波束卫星同时采用时分复用接入(Time DivisionMultiple Access，TDMA)和空分复用接入(Space Division Multiple Access,SDMA)。根据TDMA，卫星通信以帧单位为多个用户服务，每个帧包含T个时隙。在第t(t＝1,2…,T)时隙，有Lt个用户，每个波束中最多选择1个用户以保证所选的Lt个用户能被卫星通过SDMA的方式服务。假设由K个用户构成的索引集为
Figure BDA0001915067960000051
即
Figure BDA0001915067960000052
第t个时隙所选的Lt用户构成的用户集合为St，因此有
Figure BDA0001915067960000053
为了保证每个用户的公平性，每个用户在每一帧最多被服务一次，因此我们设置：
Figure BDA0001915067960000054
多波束卫星信道
Figure BDA0001915067960000055
可以建模成：
H＝ΦC (2)
其中，矩阵
Figure BDA0001915067960000056
代表用户和卫星之间不同传播路径引起的相位差异；实矩阵
Figure BDA0001915067960000057
代表多波束天线模式。由于卫星天线馈源间距与传播路径相比很小，所以在直视传播环境(line-of-sight，LOS)中，用户与卫星不同馈源之间的相位通常认为是相同的。因此，矩阵Φ为对角阵，对角元素为[Φ]ii＝ejφi,i＝1,…,K，φi为(0,2π]内均匀分布的变量。矩阵C包含接收天线增益、卫星辐射模式、路径损耗和噪声功率的影响。C矩阵第(k,n)个元素为
Figure BDA0001915067960000061
其中，GR为用户接收天线增益，Gk,n为第n个馈源和第k个用户之间的增益。dk表示第k个用户和卫星之间的距离，λ为波长，k为玻尔兹曼常数，TR为晴空环境下接收噪声温度，BW为载波带宽。
假设第t个时隙选出的Lt个用户的索引为
Figure BDA0001915067960000062
m为用户集合
Figure BDA0001915067960000063
中的用户，那么有St＝{m1,m2,…,mt}。St中Lt个用户接收信号表示为
Figure BDA0001915067960000064
那么第t个时隙的传输模型可以表示为
Figure BDA0001915067960000065
其中
Figure BDA0001915067960000066
表示为多波束卫星与Lt个用户之间的信道，
Figure BDA0001915067960000067
表示多波束卫星与ml用户之间的信道，l＝1,…,Lt。
Figure BDA0001915067960000068
表示第t个时隙的传输信号。
Figure BDA0001915067960000069
表示均值为0噪声方差为σ2独立同分布加性高斯白噪声同分布的加性高斯白噪声(Additive White Gaussian Noise，AWGN)，即
Figure BDA00019150679600000610
为了消除由信道矩阵引起的波束间干扰，引入预编码矩阵进行预处理来对抗信道失真。假设
x＝Ws (5)
其中
Figure BDA00019150679600000611
表示预编码矩阵，
Figure BDA00019150679600000612
表示发送信号，
Figure BDA00019150679600000616
表示给用户ml发送的信号，l＝1,…,Lt。假设其相互独立且模值为1，即
Figure BDA00019150679600000613
将式(5)带入式(4)，得
Figure BDA00019150679600000614
用户ml的接收信号可以表示为
Figure BDA00019150679600000615
用户ml的信干噪比(Signal to Interference plus Noise Ratio,SINR)为
Figure BDA0001915067960000071
用户ml的可达速率表示为
Figure BDA0001915067960000072
预编码采用ZF预编码，可以表示为
W＝αF (10)
其中F表示未归一化的预编码矩阵，α是归一化因子，分别表示为
Figure BDA0001915067960000073
ρ表示发射功率，σ2表示噪声方差，那么用户ml的SINR表示为
Figure BDA0001915067960000074
那么用户ml的可达速率可以表示为
Figure BDA0001915067960000075
那么对于含有K个用户的小组，总和速率为
Figure BDA0001915067960000076
本发明一种多波束卫星通信系统用户分组的方法的装置，如图2所示，其包括：初始化模块、小组初始化模块、小组分配模块、小组输出模块和预编码设计模块。具体地，初始化模块用于初始化参数的设置，初始化参数包括初始化用户集合、设置小组总数和初始化小组循环次数。小组初始化模块用于选择用户集合中模值最大的用户，初始化设定的小组集合。小组分配模块用于按照设定的用户分组选择标准，为完成初始化的小组集合分配合适的用户，其包括用户选择模块、计算模块和选择更新模块，用于为用户按照所设计的准则分配小组。其中，用户选择模块用于选择用户集合中的用户及其对应的信道向量；计算模块用于计算用户加入小组后带来的增益；选择更新模块，用于选择合适的用户，并且更新小组集合和用户集合。小组输出模块用于输出最后的用户分组，并进一步通过预编码设计模块设计预编码矩阵，通过预处理来对抗多波束卫星与用户之间的信道失真。
本发明一种多波束卫星通信系统用户分组的方法，其根据卫星服务用户场景的不同，对卫星服务的地面用户进行分组，一般分为两组：一组只选择大部分用户被服务的情况，以实现小组最大化和速率，考虑的是最大化原则；另一组选择所有用户都被服务的情况，采用公平性原则。
针对上述方法，本发明提出了两种方案：最大化系统和速率的方案和服务所有用户的方案。具体如下。
实例一：最大化系统和速率方案，其多波束卫星通信系统在最大化和速率场景下的用户分组方法的流程图，见图3A所示。
设计原则：在指定时隙总数的条件下，最大化小组的系统和速率。
随着小组中用户数目的增多，小组的系统和速率呈先升后降的趋势。因此每一个小组，都有一个最优的用户数目，这个最优用户数目与系统波束的个数相关。为了最大化系统和速率，使每个小组选择最优的用户数目，由此提出最大化和速率用户分组算法(Sum-rate Maximization User Grouping，SMUG)和低复杂度最大化和速率用户分组算法(LowComplexity Sum-rate Maximization User Grouping，LCSMUG)，其中，低复杂度最大化和速率用户分组算法采用最大化信道容量方法。
方法1：SMUG算法
步骤1、初始化用户集合
Figure BDA0001915067960000081
设置小组总数为T，初始化小组循环次数t＝0。
步骤2、判断小组循环次数t是否小于等于小组总数T，若满足条件，执行步骤3，否则执行步骤9。
步骤3、循环计算用户集合K中所有用户信道向量的模值，并选择模值最大的用户u*，即
Figure BDA0001915067960000082
初始化小组用户
Figure BDA0001915067960000091
将u*加入小组用户
Figure BDA0001915067960000092
即
Figure BDA0001915067960000093
从用户集合
Figure BDA0001915067960000094
中除去用户u*，即
Figure BDA0001915067960000095
小组循环次数t加1。通过式(14)计算当前小组的和速率
Figure BDA0001915067960000096
步骤4、判断小组
Figure BDA0001915067960000097
用户个数是否小于等于波束个数N，即
Figure BDA0001915067960000098
如果满足条件，执行步骤5，否则执行步骤2。
步骤5、根据当前用户小组
Figure BDA0001915067960000099
创建临时用户小组
Figure BDA00019150679600000910
m为用户集合
Figure BDA00019150679600000911
中的用户，选择合适的用户m*，使系统最大化和速率，即选择用户分组标准为
Figure BDA00019150679600000912
其中，
Figure BDA00019150679600000913
表示临时用户小组
Figure BDA00019150679600000914
中用户信道向量组成的信道矩阵。
根据式(14)，在为当前用户小组
Figure BDA00019150679600000915
分配新的用户时，
Figure BDA00019150679600000916
是固定的，因此最大化
Figure BDA00019150679600000917
相当于最大化
Figure BDA00019150679600000918
等价于最小化
Figure BDA00019150679600000919
因此在选择合适用户的时候，可以将用户集合
Figure BDA00019150679600000920
中的用户m临时加入到当前用户
Figure BDA00019150679600000921
构成临时用户小组
Figure BDA00019150679600000922
然后计算临时用户小组
Figure BDA00019150679600000923
的
Figure BDA00019150679600000924
选择最小的值对应的用户m*，对应于式(16)。
步骤6、将合适的用户m*加入当前用户小组
Figure BDA00019150679600000925
构成临时用户小组
Figure BDA00019150679600000926
即
Figure BDA00019150679600000927
根据式(14)计算临时小组和速率
Figure BDA00019150679600000928
步骤7、判断临时小组和速率
Figure BDA00019150679600000929
是否大于当前小组和速率
Figure BDA00019150679600000930
如果满足条件，则执行步骤8，否则执行步骤2。
步骤8、将当前小组更新为临时小组，即
Figure BDA00019150679600000931
当前小组和速率更新为临时小组和速率，即
Figure BDA00019150679600000932
将用户m*从用户集合
Figure BDA00019150679600000933
中除去，即
Figure BDA00019150679600000934
执行步骤4。
步骤9、输出
Figure BDA00019150679600000935
即所求得的用户分组集合。
步骤10、根据式(10)和式(11)输出每个时隙的预编码矩阵。
方法2：LCSMUG算法
SMUG算法在选择用户的过程中需要多次计算高维矩阵的逆，因此复杂度很高。我们针对这一缺陷做出改进，从信道容量的角度出发，对于
Figure BDA0001915067960000101
的信道传输模型，其信道容量可以定义为
Figure BDA0001915067960000102
其中γ＝ρ/σ2，表示发射功率与噪声方差的比值。在当前用户小组
Figure BDA00019150679600001018
的基础上，现在想要选择一个新的用户m，定义形成的新的用户小组为
Figure BDA0001915067960000103
那么式(17)可以展开分解为
Figure BDA0001915067960000104
其中，IN表示N维单位矩阵。
步骤1、初始化用户集合
Figure BDA0001915067960000105
设置小组总数为T，初始化小组循环次数t＝0。
步骤2、判断小组循环次数t是否小于等于小组总数T，若满足条件，执行步骤3，否则执行步骤4。
步骤3、初始化矩阵A＝γI。循环计算用户集合
Figure BDA00019150679600001019
中所有用户信道向量的模值，并选择模值最大的用户u*，见式(15)。初始化用户小组
Figure BDA0001915067960000106
将u*加入用户小组
Figure BDA0001915067960000107
即
Figure BDA0001915067960000108
从用户集合
Figure BDA0001915067960000109
中除去用户u*，即
Figure BDA00019150679600001010
小组循环次数t加1。通过式(14)计算当前小组的和速率
Figure BDA00019150679600001011
通过式(19)更新矩阵A。
Figure BDA00019150679600001012
步骤4、判断用户小组
Figure BDA00019150679600001013
是否小于等于波束个数N，即
Figure BDA00019150679600001014
如果满足条件，执行步骤5，否则执行步骤2。
步骤5、选择用户集合
Figure BDA00019150679600001015
中的用户m，以及对应的信道向量hm；根据当前用户小组
Figure BDA00019150679600001016
计算加入新用户所得的增益值，其中，
Figure BDA00019150679600001017
为第t(t＝1,2…,T)个时隙所选的Lt个用户构成的用户集合，T为卫星通信单位帧所包含的时隙个数。
根据每个用户加入当前小组所得的增益值选择合适的用户m*，用户分组选择标准为最大化信道容量方法，即
Figure BDA0001915067960000111
根据式(18)，最大化
Figure BDA0001915067960000112
等效于最大化：
Figure BDA0001915067960000113
进一步等效于最大化：
Figure BDA0001915067960000114
由于
Figure BDA0001915067960000115
可以预先计算，因此定义
Figure BDA0001915067960000116
那么在选择合适的用户时可以计算用户集合
Figure BDA0001915067960000117
中的用户m的
Figure BDA0001915067960000118
选择最大值对应的用户m*，对应于式(20)。
步骤6、将合适的用户m*加入当前用户小组
Figure BDA0001915067960000119
构成临时用户小组
Figure BDA00019150679600001110
即
Figure BDA00019150679600001111
通过式(14)计算临时小组和速率
Figure BDA00019150679600001112
通过式(19)更新矩阵A。
步骤7、判断临时小组和速率
Figure BDA00019150679600001113
是否大于当前小组和速率
Figure BDA00019150679600001114
如果满足条件，则执行步骤8，否则执行步骤2。
步骤8、将当前用户小组更新为临时用户小组，即
Figure BDA00019150679600001115
当前小组和速率更新为临时小组和速率，即
Figure BDA00019150679600001116
将用户m*从用户集合
Figure BDA00019150679600001117
中除去，即
Figure BDA00019150679600001118
执行步骤4。
步骤9、输出
Figure BDA00019150679600001119
即所求得的用户分组集合。
步骤10、根据式(10)和式(11)输出每个时隙的预编码矩阵。
实例二：服务所有用户方案，其多波束卫星通信系统在服务所有用户场景下的用户分组方法的流程图，见图3B所示。
设计标准：在指定时隙总数下，服务地面所有用户。
由于实例一考虑最大化原则，为了小组最大化和速率，只能选择大部分用户。而卫星信道改变比较缓慢，即剩下的部分用户可能永远不会被服务。针对这一缺陷，实例二采用公平性原则，考虑所有用户都被服务的情况。由此提出服务所有用户的最大化和速率用户分组算法(SMUG with full selection of users，SMUG-FS)和服务所有用户的低复杂度最大化和速率用户分组算法(LCSMUG with full selection of users，LCSMUG-FS)。
方法3：SMUG-FS算法
步骤1、初始化用户集合
Figure BDA00019150679600001120
设置小组总数为T，初始化小组循环次数t＝0。
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤4。
步骤3、计算用户集合
Figure BDA0001915067960000121
中所有用户信道向量的模值
Figure BDA0001915067960000122
并选择模值最大的用户u*，见式(15)。将u*加入用户小组
Figure BDA0001915067960000123
从用户集合
Figure BDA0001915067960000124
中除去用户u*，即
Figure BDA0001915067960000125
小组循环次数t加1，执行步骤2。
步骤4、判断用户集合
Figure BDA0001915067960000126
是否为空集，如果满足条件，执行步骤9，否则执行步骤5；
步骤5、选择用户集合
Figure BDA0001915067960000127
中的用户m，对应的信道向量为hm，初始化小组循环次数t为0；
步骤6、判断小组循环次数t是否小于等于小组总数T，如果满足条件，执行步骤7，否则执行步骤8；
步骤7、第t用户小组为
Figure BDA0001915067960000128
对应信道矩阵
Figure BDA0001915067960000129
计算第t小组加入新用户所得的增益V(t)，即
Figure BDA00019150679600001210
小组循环次数加1，执行步骤6。
步骤8、选择加入新用户所得增益最大的用户小组
Figure BDA00019150679600001218
即t*＝max1≤t≤TV(t)，
将用户m加入临时用户小组
Figure BDA00019150679600001211
从用户集合
Figure BDA00019150679600001212
中除去用户m，即
Figure BDA00019150679600001213
执行步骤4。
步骤9、输出
Figure BDA00019150679600001214
即所求得的用户分组集合。
步骤10、根据式(10)和式(11)输出每个时隙的预编码矩阵。
方法4：LCSMUG-FS算法
步骤1、初始化用户集合
Figure BDA00019150679600001215
设置小组总数为T，初始化小组循环次数t＝0。
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤4。
步骤3、初始化矩阵
Figure BDA00019150679600001216
计算用户集合
Figure BDA00019150679600001217
中所有用户信道向量的模值
Figure BDA0001915067960000131
并选择模值最大的用户u*，见式(15)。将u*加入用户小组
Figure BDA0001915067960000132
Figure BDA0001915067960000133
通过式(19)更新矩阵A(t)，从用户集合
Figure BDA0001915067960000134
中除去用户u*，即
Figure BDA0001915067960000135
小组循环次数t加1，执行步骤2。
步骤4、判断用户集合是否为空集，如果满足条件，执行步骤9，否则执行步骤5；
步骤5、选择用户集合
Figure BDA0001915067960000136
中的用户m，用户m对应的信道向量为hm，初始化小组循环次数t为0；
步骤6、判断小组循环次数t是否小于等于小组总数T，如果满足条件，执行步骤7，否则执行步骤8；
步骤7、计算第t小组加入新用户所得的增益V(t)，
Figure BDA0001915067960000137
小组循环次数t加1，执行步骤6。
步骤8、选择加入新用户所得增益最大的小组
Figure BDA0001915067960000138
即t*＝max1≤t≤TV(t)，将用户m加入新用户所得增益最大的小组
Figure BDA0001915067960000139
通过式(19)更新矩阵A(t*)，再选择合适的小组t*每一次迭代更新矩阵为
Figure BDA00019150679600001310
从用户集合
Figure BDA00019150679600001311
中除去用户m，即
Figure BDA00019150679600001312
执行步骤4。
步骤9、输出
Figure BDA00019150679600001313
即所求得的用户分组集合。
步骤10、根据式(10)和式(11)输出每个时隙的预编码矩阵。
为了验证本发明方法的效果，将本发明方法与现有的多天线下行正交分组算法(multiple antenna downlink orthogonal clustering，MADOC)方案进行仿真对比验证实验：假设多波束卫星在同步轨道运行，通过反射面将16个馈源产生16个波束服务地面用户，地面上320个用户均匀分布在16个波束中。卫星载波频率为20GHz(Ka波段)，用户链路带宽B为500MHz，用户天线增益为41.7dBi，地面用户接收品质因数G/T为17.68dB/K，玻尔兹曼常数κ＝1.38×10-23J/K。定义小组平均和速率即系统平均和速率为
Figure BDA0001915067960000141
图4A和图4B比较了本发明用户分组方法和现有用户分组方法MADOC算法在不同发射功率下小组平均和速率对比曲线图。设置发送功率为12dBW到24dBW，时隙数目为20。从图4A中可以看出，在发送功率为12dBW到24dBW，SUMG算法性能要优于MADOC算法；从图4B中可以看出，在发送功率为12dBW到24dBW，LCSMUG算法性能要优于MADOC算法；从图4A中可以看出，当发射功率大于20dBW时，SUMG-FS算法的性能优于MADOC算法；从图4B中可以看出，当发射功率大于19dBW时，LCSMUG-FS算法的性能也优于MADOC算法。
图5A和图5B比较了本发明用户分组方法和现有用户分组方法MADOC算法在不同时隙数量下小组平均和速率对比曲线图。设置放射功率为16dBW。从图5A和图5B中可以看出，对于SUMG、LCSMUG和MADOC算法，随着时隙数目的增加，小组平均和速率均呈下降趋势；从图5A和图5B中可以看出，而对于SUMG-FS和LCSMUG-FS算法，小组平均和速率呈先上升后下降的趋势。从图中可以看出，SUMG和LCSMUG算法性能要优于MADOC算法，当时隙数目大于22时，SUMG-FS和LCSMUG-FS算法的性能要优于MADOC算法，当时隙数目大于26时，SUMG-FS和LCSMUG-FS算法的性能要优于SUMG和LCSMUG算法。
图6比较了本发明用户分组方法和现有用户分组方法MADOC算法在不同时隙数量下服务用户总数对比曲线图。由于SUMG-FS和LCSMUG-FS算法考虑服务所有的用户，因此没有参与仿真。可以看出，随着时隙数目的增加，SUMG、LCSMUG和MADOC算法服务用户总数都呈上升趋势，SUMG和LCSMUG算法相比于MADOC算法可以服务更多的用户，此外在时隙数目小于24时，SMUG算法可以比LCSMUG算法服务更多的用户。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步地详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种多波束卫星通信系统用户分组的方法，其包括以下步骤：
步骤1、初始化系统参数，所述系统参数包括用户集合、小组总数和小组循环次数；
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤9；
步骤3、选择用户集合中用户信道向量的模值最大的用户，初始化用户小组，并更新用户集合，所述小组循环次数加1，计算当前小组和速率；
步骤4、判断当前小组用户数量是否小于等于波束个数，如果满足条件，执行步骤5，否则执行步骤2；
步骤5、根据用户分组选择标准，为上述完成初始化的小组集合分配合适的用户；
步骤6、将所选合适的用户加入当前小组构成临时小组，计算临时小组和速率；
步骤7、判断临时小组和速率是否大于当前小组和速率，如果满足条件，则执行步骤8，否则执行步骤2；
步骤8、更新当前小组集合和当前小组和速率，执行步骤4；
步骤9、输出最后的用户分组集合；
步骤10、输出每个时隙的预编码矩阵。
2.根据权利要求1所述的方法，其特征在于，步骤5中，所述用户分组选择标准包括采用最大化和速率方法，所述合适的用户：
Figure FDA0003133988070000011
其中，
Figure FDA0003133988070000012
表示临时小组个数
Figure FDA0003133988070000013
中用户信道向量组成的信道矩阵。
3.根据权利要求1所述的方法，其特征在于，步骤5中，所述用户分组选择标准包括采用最大化信道容量方法，所述合适的用户：
Figure FDA0003133988070000014
其中，初始化矩阵A＝γI，hm为用户m对应的信道向量hm，在选择合适的用户m*后，每一次迭代更新为：
Figure FDA0003133988070000021
4.一种多波束卫星通信系统用户分组的方法，其包括以下步骤：
步骤1、初始化系统参数，所述系统参数包括用户集合、小组总数和小组循环次数；
步骤2、判断小组循环次数是否小于等于小组总数，若满足条件，执行步骤3，否则执行步骤4；
步骤3、选择用户集合中用户信道向量的模值最大的用户，初始化用户小组，并更新用户集合，所述小组循环次数加1，执行步骤2；
步骤4、判断用户集合是否为空集，如果满足条件，执行步骤9，否则执行步骤5；
步骤5、选择用户集合中的用户和相对应的信道向量，初始化小组循环次数；
步骤6、判断小组循环次数是否小于等于小组总数，如果满足条件，执行步骤7，否则执行步骤8；
步骤7、计算当前小组加入新用户后所得增益，小组循环次数加1，执行步骤6；
步骤8、更新增益最大的小组，更新用户集合，执行步骤4；
步骤9、输出最后的用户分组集合；
步骤10、输出每个时隙的预编码矩阵和小组和速率。
5.根据权利要求4所述的方法，其特征在于，步骤7中，所述当前小组加入新用户后所得增益为第t小组加入新用户所得的增益V(t)，即
Figure FDA0003133988070000022
其中，
Figure FDA0003133988070000023
为第t用户小组的小组用户个数
Figure FDA0003133988070000024
对应信道矩阵。
6.根据权利要求4所述的方法，其特征在于，步骤7中，所述当前小组加入新用户后所得增益为第t小组加入新用户所得的增益V(t)，即
Figure FDA0003133988070000025
其中，初始化矩阵A(t)＝γI，hm为用户m对应的信道向量；
将用户m加入新用户所得增益最大的小组
Figure FDA0003133988070000026
每一次迭代更新矩阵为
Figure FDA0003133988070000031
其中，增益最大的小组序号为t*＝max1≤t≤TV(t)。
7.一种多波束卫星通信系统用户分组的方法的装置，其特征在于，其包括：
初始化模块，用于初始化用户集合、小组总数和小组循环次数；
小组初始化模块，用于选择用户集合中模值最大的用户，初始化小组集合；
小组分配模块，根据卫星服务的地面用户场景的不同进行分组，提出最大化系统和速率方案和服务所有用户方案，根据卫星服务用户场景的不同，对卫星服务的地面用户进行分组，一般分为两组：一组只选择大部分用户被服务的情况，以实现小组最大化和速率，考虑的是最大化原则；另一组选择所有用户都被服务的情况，采用公平性原则；为完成初始化的小组集合分配合适的用户；
小组输出模块，用于输出最后的用户分组；
预编码设计模块，用于设计每个时隙的预编码矩阵。
8.根据权利要求7所述的装置，其特征在于，所述小组分配模块包括：
用户选择模块，用于选择用户集合中的用户和相对应的信道向量；
计算模块，用于计算用户加入小组后带来的增益；
选择更新模块，用于选择合适的用户，并且更新小组集合和用户集合。
9.根据权利要求7所述的装置，其特征在于，所述用户分组选择标准包括选择大部分用户被服务，实现小组最大化和速率。
10.根据权利要求7所述的装置，其特征在于，所述用户分组选择标准包括选择所有用户都被服务，实现小组最大化增益。
CN201811569312.9A 2018-12-21 2018-12-21 一种多波束卫星通信系统用户分组的方法及其装置 Active CN109547075B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811569312.9A CN109547075B (zh) 2018-12-21 2018-12-21 一种多波束卫星通信系统用户分组的方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811569312.9A CN109547075B (zh) 2018-12-21 2018-12-21 一种多波束卫星通信系统用户分组的方法及其装置

Publications (2)

Publication Number Publication Date
CN109547075A CN109547075A (zh) 2019-03-29
CN109547075B true CN109547075B (zh) 2021-09-07

Family

ID=65856072

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811569312.9A Active CN109547075B (zh) 2018-12-21 2018-12-21 一种多波束卫星通信系统用户分组的方法及其装置

Country Status (1)

Country Link
CN (1) CN109547075B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110518961B (zh) 2019-08-29 2020-12-11 东南大学 大规模mimo卫星移动通信方法及系统
CN112751594B (zh) * 2020-12-17 2021-12-14 天地信息网络研究院(安徽)有限公司 基于和速率最大化的卫星通信用户分组方法
CN113177561B (zh) * 2021-04-27 2022-06-24 深圳市慧鲤科技有限公司 多用户集合方法及装置、电子设备和存储介质
CN113595612B (zh) * 2021-06-18 2022-11-15 中国人民解放军63921部队 适于卫星移动通信系统的双星用户调度方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634244A (zh) * 2013-12-03 2014-03-12 东南大学 密集分布式无线网络动态导频复用信道估计与导频分配方法
CN103796319A (zh) * 2014-01-16 2014-05-14 北京大学 一种多波束卫星移动通信下行链路频率复用方法
CN104320174A (zh) * 2014-09-16 2015-01-28 南昌航空大学 一种基于部分信道信息的卫星多波束协作传输方法
CN105007129A (zh) * 2015-06-16 2015-10-28 东南大学 多波束卫星移动通信系统容量估计及用户调度
CN107864506A (zh) * 2017-10-25 2018-03-30 南京理工大学 基于萤火虫算法的Ka多波束卫星通信系统资源分配方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8693970B2 (en) * 2009-04-13 2014-04-08 Viasat, Inc. Multi-beam active phased array architecture with independant polarization control

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103634244A (zh) * 2013-12-03 2014-03-12 东南大学 密集分布式无线网络动态导频复用信道估计与导频分配方法
CN103796319A (zh) * 2014-01-16 2014-05-14 北京大学 一种多波束卫星移动通信下行链路频率复用方法
CN104320174A (zh) * 2014-09-16 2015-01-28 南昌航空大学 一种基于部分信道信息的卫星多波束协作传输方法
CN105007129A (zh) * 2015-06-16 2015-10-28 东南大学 多波束卫星移动通信系统容量估计及用户调度
CN107864506A (zh) * 2017-10-25 2018-03-30 南京理工大学 基于萤火虫算法的Ka多波束卫星通信系统资源分配方法

Also Published As

Publication number Publication date
CN109547075A (zh) 2019-03-29

Similar Documents

Publication Publication Date Title
CN109547075B (zh) 一种多波束卫星通信系统用户分组的方法及其装置
Zheng et al. Generic optimization of linear precoding in multibeam satellite systems
EP2474104B1 (en) Tomlinson harashima precoding with additional receiver processing in a multi-user multiple-input multiple-output wireless transmission system
Guidotti et al. Geographical scheduling for multicast precoding in multi-beam satellite systems
Li et al. Optimizing channel-statistics-based analog beamforming for millimeter-wave multi-user massive MIMO downlink
US20230261710A1 (en) Pragmatic solutions for massive mimo for broadband telecommunication systems
Wang et al. NOMA-enabled multi-beam satellite systems: Joint optimization to overcome offered-requested data mismatches
CN113938183B (zh) 多波束卫星系统下基于非正交多址的通信资源分配方法
Guidotti et al. Design trade-off analysis of precoding multi-beam satellite communication systems
Guidotti et al. User clustering for multicast precoding in multi-beam satellite systems
Christopoulos et al. Frame based precoding in satellite communications: A multicast approach
Riviello et al. Joint graph-based user scheduling and beamforming in LEO-MIMO satellite communication systems
Boussemart et al. On the achievable rates for the return-link of multi-beam satellite systems using successive interference cancellation
Ahmad et al. Graph-based user scheduling algorithms for LEO-MIMO non-terrestrial networks
Christopoulos et al. Multibeam joint precoding: Frame-based design
Shi et al. Beamforming with superposition coding in multiple antenna satellite communications
Arvola et al. Two-layer precoding for dimensionality reduction in massive MIMO
US10498410B2 (en) Methods and apparatuses for transmit weights
Vázouez et al. Closed-form multicast precoding for satellite flexible payloads
Vidal et al. Joint power, frequency and precoding optimisation in a satellite sdma communication system
CN113258985A (zh) 单站多星mimo上注系统能效优化方法
Wang et al. Beamforming and power allocation in noma-based multibeam satellite systems with outage constraint
Kaulich et al. Efficient rate splitting method using successive beamforming techniques
Ahmad et al. Improved Graph-Based User Scheduling For Sum-Rate Maximization in LEO-NTN Systems
Nguyen et al. Resource allocation for multibeam MISO satellite systems: Sum rate versus proportional fair optimization

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant