基于布谷鸟搜索及蚁群优化的波束选择方法及装置

1.本发明涉及移动通信技术技术领域，更具体的说是涉及一种基于布谷鸟搜索及蚁群优化的波束选择方法及装置。


背景技术：

2.目前，毫米波大规模mimo系统具有大带宽、高数据传输速率、强频谱扩展能力和抗干扰能力等优势，成为5g中的一项极具潜力的关键技术。然而传统的mimo系统需要为每根天线配备一条专用射频(radio frequency,rf)链，这导致毫米波大规模mimo系统具有较高的硬件成本和能源消耗。为了减小所需rf链的数量，提出了基于离散透镜阵列(discrete lens array,dla)的波束空间mimo系统，利用dla将传统空间信道转换成波束空间信道，而波束空间信道具有稀疏性，只有少量主导波束对系统性能有较大贡献，因此有效的波束选择对于减少rf链数量和提高系统性能至关重要。
3.但是，现有技术中通常采用蚁群优化算法实现毫米波大规模mimo系统的波束选择，首先为各用户分配一个备选波束集，然后在每次迭代中计算当前用户备选波束集中各波束被选择的概率，找出被选概率最大的波束分配给当前用户，然而aco算法中存在的正反馈机制会使次优解占据优势，从而易使得算法陷入局部最优，从而较难寻求到全局最优解。
4.因此，如何提供一种能够解决上述问题的波束选择方法是本领域技术人员亟需解决的问题。


技术实现要素：

5.有鉴于此，本发明提供了一种基于布谷鸟搜索及蚁群优化的波束选择方法及装置，提出一种布谷鸟搜索(cuckoo search,cs)算法和aco算法相结合的波束选择方案csaco，首先基于全局搜索能力强的cs算法构建波束选择系统，然后引入局部搜索能力强的aco算法对cs算法的随机解进行局部寻优，从而加快cs算法的收敛速度，同时cs算法中解的随机性也可防止aco算法陷入局部最优，cs算法和aco算法联合优化可使可行解朝向全局最优的方向进化，最终达到接近最优的系统和速率性能。
6.为了实现上述目的，本发明采用如下技术方案：
7.一种基于布谷鸟搜索及蚁群优化的波束选择方法，应用于毫米波大规模mimo系统，所述毫米波大规模mimo系统包括多个用户，包括以下步骤：
8.s1：获取备选波束方向集合b＝{b1,b2,...,bk}∈r
k*c
，其中b表示所有用户的备选波束方向集合，bk表示用户k的备选波束方向集合，k表示总的用户数量，c表示各用户备选波束方向集合中波束的数量；
9.s2：引入位置索引集合xi＝{x
i,1
,x
i,2
,...,x
i,k
,...,x
i,k
}∈r
k*1
(i＝1,2,...,m)，其中x
i,k
表示个体i给用户k分配的波束为bk中的第x
i,k
个元素，m表示一个种群内的个体数量；
10.s3：对所述位置索引集合进行多次更新，通过获取x
i,k
与bk之间的对应关系，从所
述备选波束方向集合中寻找出为所述用户分配的波束方向；
11.s4：通过波束选择方法对个体进行迭代求解出全局最优解，完成波束选择。
12.优选的，所述s4具体包括：
13.s41：初始化种群个体，并通过levy飞行生成随机解；
14.s42：通过蚁群优化方法对所述随机解进行局部寻优，比较所有个体适应度，保留质量高的个体为第一个体；
15.s43：对所述第一个体进行舍弃并随机生成新个体；
16.s44：比较所述新个体以及所述第一个体的适应度，并保留质量高的个体为第二个体；
17.s45：对所述第二个体判断是否满足迭代停止条件，若满足则计算所述第二个体的适应度，并保留质量最高的个体，即为全局最优解。
18.优选的，所述s45：对第二个体判断是否满足迭代停止条件，若满足则计算第二个体的，所述s45中还包括：若不满足则返回至所述步骤s41。
19.优选的，所述步骤s41中，生成随机解的具体公式为：
[0020][0021]
式中，式中，l(β)是levy飞行的随机寻优路线，μ、v均是服从正态分布的随机数，β是参数偏度。
[0022]
进一步，本发明还提供一种根据上述任一项所述的一种基于布谷鸟搜索及蚁群优化的波束选择方法的装置，包括：
[0023]
获取模块，用于获取备选波束方向集合b＝{b1,b2,...,bk}∈r
k*c
，其中b表示所有用户的备选波束方向集合，bk表示用户k的备选波束方向集合，k表示总的用户数量，c表示各用户备选波束方向集合中波束的数量；以及引入位置索引集合xi＝{x
i,1
,x
i,2
,...,x
i,k
,...,x
i,k
}∈r
k*1
(i＝1,2,...,m)，其中x
i,k
表示个体i给用户k分配的波束为bk中的第x
i,k
个元素，m表示一个种群内的个体数量；
[0024]
更新模块，用于对所述位置索引集合进行多次更新，通过获取x
i,k
与bk之间的对应关系，从所述备选波束方向集合中寻找出为所述用户分配的波束方向；
[0025]
选择模块，用于通过波束选择方法对个体进行迭代求解出全局最优解，完成波束选择。
[0026]
经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种基于布谷鸟搜索及蚁群优化的波束选择方法及装置，将所有分配给用户的波束组合起来构成种群中的个体，而对所有个体进行优化可同时考虑到各用户间的干扰并直接以最大化系统和速率为优化目标进行寻优。cs算法依据levy飞行机制可搜索到较好的全局可行解，但其强随机性可能使cs算法在最优解附近震荡从而导致算法收敛速度慢；而aco算法依据其较强的局部寻优能力可搜索到cs算法随机解附近的最优解，从而加快算法的收敛速度，同时全局可行解的强随机性也可防止aco算法陷入局部最优。
附图说明
[0027]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现
有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
[0028]
图1为本发明提供的一种基于布谷鸟搜索及蚁群优化的波束选择方法的整体流程图；
[0029]
图2为本发明实施例提供的步骤s4的具体执行流程图；
[0030]
图3为本发明提供的一种基于布谷鸟搜索及蚁群优化的波束选择装置结构示意图；
[0031]
图4为本发明实施例提供的各算法在不同信噪比条件下的可达和速率性能对比曲线图。
具体实施方式
[0032]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0033]
参见附图1所示，本发明实施例公开了一种基于布谷鸟搜索及蚁群优化的波束选择方法，应用于毫米波大规模mimo系统，所述毫米波大规模mimo系统包括多个用户，包括以下步骤：
[0034]
s1：获取备选波束方向集合b＝{b1,b2,...,bk}∈r
k*c
，其中b表示所有用户的备选波束方向集合，bk表示用户k的备选波束方向集合，k表示总的用户数量，c表示各用户备选波束方向集合中波束的数量；
[0035]
s2：引入位置索引集合xi＝{x
i,1
,x
i,2
,...,x
i,k
,...,x
i,k
}∈r
k*1
(i＝1,2,...,m)，其中x
i,k
表示个体i给用户k分配的波束为bk中的第x
i,k
个元素，m表示一个种群内的个体数量；
[0036]
s3：对位置索引集合进行多次更新，通过获取x
i,k
与bk之间的对应关系，从备选波束方向集合中寻找出为用户分配的波束方向；
[0037]
s4：通过波束选择方法对个体进行迭代求解出全局最优解，完成波束选择。
[0038]
设本实施例提供的方法共有t
max
次迭代，在第t次迭代中，对于用户k有n个可选的波束方向，但由于信道具有稀疏性，各用户对应的主导波束数目远小于n，故本方法设定用户k在每次迭代时可选的波束数目为c＜＜n，将用户k的备选波束集表示为bk＝{i
k,1
,i
k,2
,...,i
k,c
,i
k,c+1
,...,i
k,c
}∈r
1*c
，其中i
k,c
为用户k的第c个备选波束，则所有用户的备选波束集b＝{b1,b2,...,bk}∈r
k*c
。本方法根据最大幅度准则为各用户预先求出bk，且不同备选波束对应的信道功率幅值关系为其中表示用户k的第c个波束所对应的信道功率幅值。
[0039]
为了保证每次迭代分配给用户k的波束均为bk中的元素，本方法引入一个新的位置索引集合，其中xi＝{x
i,1
,x
i,2
,...,x
i,k
,...,x
i,k
}∈r
k*1
(i＝1,2,...,m)表示一个种群内的个体数量(鸟巢数)，x
i,k
表示个体i给用户k分配的波束为bk中的第x
i,k
个元素。在本实施例提供的方法的迭代过程中，对集合xi不断进行更新，并且通过x
i,k
和bk的对应关系从中找
出个体i为用户k分配的波束方向g
i,k
，若gi＝{g
i,1
,g
i,2
,...,g
i,k
}∈r
k*1
使得的值越小，其中为波束空间信道矩阵，则认为个体i的质量越高，最终通过迭代求出质量最高的个体g。
[0040]
参见附图2所示，在一个具体的实施例中，s4具体包括：
[0041]
s41：初始化种群个体，并通过levy飞行生成随机解；
[0042]
s42：通过蚁群优化方法对随机解进行局部寻优，比较所有个体适应度，保留质量高的个体为第一个体；
[0043]
s43：对第一个体进行舍弃并随机生成新个体；
[0044]
s44：比较新个体以及第一个体的适应度，并保留质量高的个体为第二个体；
[0045]
s45：对第二个体判断是否满足迭代停止条件，若满足则计算第二个体的适应度，并保留质量最高的个体，即为全局最优解。
[0046]
在一个具体的实施例中，s45中还包括：若不满足则返回至步骤s41。
[0047]
在一个具体的实施例中，步骤s41中，生成随机解的具体公式为：
[0048][0049]
式中，l(β)是levy飞行的随机寻优路线，μ、v均是服从正态分布的随机数，β是参数偏度。
[0050]
具体的，对种群p＝{x1,x2,...,xm}∈r
k*m
进行更新，且在算法第一次迭代时给所有用户分配的波束均是使得其信道功率幅值最大所对应的波束方向，因此初始化种群p1为全1矩阵。levy飞行的随机游走策略满足其中表示个体i在第t代的取值，α指步长因子，表示点乘运算，l(β)是levy飞行的随机寻优路线，其满足μ、v均是服从正态分布的随机数，μ、v是参数偏度，一般取1.5。levy飞行利用强随机性在种群p1基础上进行随机游走并生成种群p2，然后把p2作为aco算法初始值，使用aco算法对p2进行局部寻优。为防止所有个体都进行局部寻优会使算法复杂度增加，故本方法仅从p2中抽取q个个体进行局部寻优并生成种群p3，然后分别对p2和p3中所有个体求适应度，并在p2中替换掉局部寻优后质量提升的个体构成种群p4。
[0051]
在现实环境中，布谷鸟蛋会有一定概率被宿主发现，故在算法中按照一定的概率p
α
舍弃p4中的某些个体并随机生成新个体构成种群p5，然后分别对p4和p5中的所有个体求适应度，并在p4中替换掉“舍弃再生成”操作后质量提升的个体构成种群p6，得到种群p6后判断是否满足停止迭代条件，若不满足，则令p1＝p6并重复上述步骤；若满足停止迭代条件，则找出p6中质量最高的个体，其中舍弃概率可以为0.25，“质量较高的个体”是能够使得d的值越小的个体。
[0052]
传统的停止迭代条件是判断是否达到最大迭代次数t
max
，而这会导致一些低效率的迭代计算。为了进一步提升所提算法的效率，本发明实施例中引入一个新的停止迭代条件，即经过迭代得到的最优解(质量最高的个体)连续e次都不变时就停止迭代，除此之外，为了避免算法一直不收敛，当达到最大迭代次数t
max
时也让其停止迭代。
[0053]
利用aco算法对levy飞行产生的随机解进行局部寻优，进一步提升levy飞行产生
的随机解的质量，较大程度避免cs算法在最优解附近发生震荡的问题，从而提高cs算法的收敛精度和收敛速度，同时随机性较强的levy飞行机制也可防止aco算法陷入局部最优。
[0054]
参见附图3所示，本发明实施例还提供一种根据上述实施例任一项的一种基于布谷鸟搜索及蚁群优化的波束选择方法的装置，包括：
[0055]
获取模块，用于获取备选波束方向集合b＝{b1,b2,...,bk}∈r
k*c
，其中b表示所有用户的备选波束方向集合，bk表示用户k的备选波束方向集合，k表示总的用户数量，c表示各用户备选波束方向集合中波束的数量；以及引入位置索引集合xi＝{x
i,1
,x
i,2
,...,x
i,k
,...,x
i,k
}∈r
k*1
(i＝1,2,...,m)，其中x
i,k
表示个体i给用户k分配的波束为bk中的第x
i,k
个元素，m表示一个种群内的个体数量；
[0056]
更新模块，用于对位置索引集合进行多次更新，通过获取x
i,k
与bk之间的对应关系，从备选波束方向集合中寻找出为用户分配的波束方向；
[0057]
选择模块，用于通过波束选择方法对个体进行迭代求解出全局最优解，完成波束选择。
[0058]
为评估本实施例所提供的方法的波束选择方法性能，同时对比分析csaco算法和现有的mm算法、ia算法、aco算法性能。
[0059]
仿真中，设置基站端天线数n＝256，用户数k＝32，rf链数n
rf
＝k，根据saleh valenzuela信道模型来生成空间信道数据，信道参数设置如下：(1)多径数l＝3；(2)a
k,l
～n(0,1)，其中a
k,l
代表第k个用户在第l条路径的复增益和延迟；(3)其它参数配置如表1所示。
[0060]
表1 csaco算法主要参数配置
[0061][0062]
各算法在不同信噪比条件下的可达和速率如图4所示，除全数字预编码外，在信噪比较小时，各算法的可达和速率相差不大，但在信噪比较大时，本文所提csacao算法可达和速率比现有算法有所提升，这是因为csaco算法结合了全局搜索和局部搜索的优点，能够使得可行解朝向全局最优的方向进化且不易陷入局部最优，从而相对其它算法能够搜索到质量较高的解。同时当用户数较少时，各用户之间的干扰也较小，故各算法为不同用户选择的波束方向基本一致，从而导致可达和速率相差较小，但在用户数增多时，用户间干扰随之增加，而csaco算法将所有用户的波束组合起来构成可进化的个体，使其在充分考虑用户间干扰条件下搜索到较好的全局可行解。
[0063]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0064]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。