[0041] 图1是本发明方法一种多果蝇群协同优化小波常模盲均衡方法原理图。
[0042] 图2是本发明方法仿真结果：(a)三种不同方法的收敛曲线，（b）WTCMA的输出星座图，(c)SFOA-WTCMA的输出星座图，(d)本发明方法MFOA-WTCMA的输出星座图。具体实施方式[0043] 下面结合附图对本发明的技术方案做进一步的详细说明：
[0044] 如图1所示，一种多果蝇群协同优化小波常模盲均衡方法原理，MFOA（Multi-fruit FliesOptimization Algorithm）为多果蝇群协同优化方法，图中除去MFOA部分，其余部分就是小波常模盲均衡方法WTCMA。本发明首先给出单果蝇群优化方法SFOA并分析其性能；其次，将协同进化引入SFOA中，提出多果蝇群协同优化方法MFOA；再次，用MFOA对波常模盲均衡方法WTCMA进行优化。
[0045] 小波常模盲均衡方法WTCMA
[0046] 如图1所示，a(k)是发射信号，h(k)是信道脉冲响应，n(k)是信道噪声，y(k)是正交小波变换器的输入信号；R(k)是正交小波变换器的输出信号，f(k)是小波常模盲均衡方法WTCMA的权向量，z(k)是小波常模盲均衡器的输出信号。
[0047] 正交小波变换器的输入信号为
[0048] y(k)=a(k)h(k)+n(k) (1)
[0049] 正交小波变换器的输出信号为
[0050] R(k)=Qy(k) (2)
[0051] 式中，Q为正交小波变换矩阵。
[0052] 小波常模盲均衡器的输出信号为
[0053] z(k)=fT(k)R(k) (3)
[0054] 式中，上标T表示转置。
[0055] z(k)经误差生成器得到的误差信号e(k)为
[0056]
[0057] 式中， 为发射信号a(k)的统计模值，是一个正常数。
[0058] WTCMA的权向量更新公式为
[0059]
[0060] 式中， 为发射信号a(k)的统计模值，是一个正常数，E表示数学期望；*表示取共轭，μ为步长，
diag表示对角矩阵； (k)表示对小波变换系数rI,K(k)的平
均功率估计， 表示尺度变换系数sI+1,K(k)的平均功率估计，其迭代公式为
[0061]
[0062] 式中，β为平滑因子，且0<β<1，一般取略小于1的数，rI,K(k)是尺度参数为I，平移参数为K的小波变换系数，I、K均取正整数；sI+1,K(k)为小波分解的最大尺度为I+1，平移参数为K的尺度变换系数；以上式子构成了小波常模盲均衡方法WTCMA。
[0063] 单果蝇群优化方法SFOA是模拟果蝇群觅食行为的过程，利用果蝇群中果蝇个体协作机制和信息共享机制不断迭代搜索更新味道浓度和位置，寻求待求解问题的全局最优解。协同进化源于自然界中种群间的协作进化机制，它是利用种群间所具有的自动调节和自动适应原理来相互驱使、影响和制约各个种群的进化；不仅可提高种群各自的搜索性能还能使全局性能得到极大的优化，突破了单种群体搜索最优解的局限；它将一个种群分为多个子种群或者直接采用多个种群，每个种群独立进化的同时相互共享和交互各种群之间的信息，协同进化整个系统直至搜索到全局最优解。
[0064] 针对SFOA搜索复杂大规模区域时局部搜索能力不足、易出现早熟现象的缺点以及协同进化的优点，将协同进化和正反反馈与实时信息共享机制与SFOA相结合，提出了多果蝇群协同优化方法MFOA。MFOA是一种基于SFOA的多果蝇群并行协同进化、具有正反反馈机制的改进果蝇群优化方法，它采用多个果蝇群并行拓扑结构，各果蝇群不仅利用从外界获得的食物味道浓度信息和位置信息指导自身的搜索完成各自的独立进化，同时把本果蝇群体搜索获得的经验与其他果蝇群分享完成系统的协同进化并通过评价各个果蝇群的最优解来指导各果蝇的进化；MFOA全面提高了局部的和全局的搜索能力、收敛效率和求解精度。
[0065] MFOA的操作过程为：（1）采用多果蝇群并行的拓扑结构进行搜索，将N个果蝇群作为N个独立进化的种群，由N个果蝇群同时搜索整个解空间；(2)实行评价机制，找出各果蝇群全局最优解；以食物味道浓度为衡量优劣的标准，N个果蝇群先各自搜索到各自种群的最优解，然后评价N个最优解的味道浓度，从N个果蝇群中找出全局最优解；(3)采用正反反馈与实时信息共享机制选取N个果蝇群中的最优解来指导各果蝇群的进化。N个果蝇群相互跟踪全局最优解，评价各果蝇群的全局最优解的优劣，选取其中的最优解指导下一步搜索，来实现信息的多向传递和实时共享。
[0066] 本发明利用MFOA优化小波常模盲均衡方法WTCMA，其原理如图1所示。图1所示的一种多果蝇群协同优化小波常模盲均衡方法MFOA-WTCMA，包括如下步骤：
[0067] 步骤一、将发射信号a(k)经过信道脉冲响应h(k)后加入信道噪声n(k)，得到正交小波变换器的输入信号y(k)：y(k)=a(k)h(k)+n(k)，其中，k为时间序列且为正整数；
[0068] 步骤二、将步骤一所述的正交小波变换器的输入信号y(k)经正交小波变换器进行正交小波变换后，得正交小波变换器的输出信号R(k)：R(k)=Qy(k)，其中，Q为正交小波变换矩阵；
[0069] 步骤三、将步骤二所述的正交小波变换器的输出信号R(k)经过小波常模盲均衡方H法WTCMA处理后得到小波常模盲均衡器的输出信号z(k)：z(k)=f(k)R(k)，其中f(k)为小波常模盲均衡方法WTCMA的权向量，上标H表示共轭转置；
[0070] 步骤四、将步骤二所述的正交小波变换器输出信号R(k)作为多果蝇群协同优化方法的输入信号，以获得WTCMA的权向量f(k)的初始优化权向量fopt(0)，该多果蝇群协同优化方法以食物味道浓度为性能评价标准，将最低食物味道浓度称为最优食物味道浓度，将与最优食物味道浓度对应的位置向量作为最优解；N个果蝇群先各自搜索到各自种群的最优解，然后比较N个最优解的食物味道浓度，从N个果蝇群中找出全局最优解；用该全局最优解来指导各果蝇群的搜索过程，使N个果蝇群跟踪全局最优解并进行下一次搜索，再次获取N个果蝇群的当前全局最优食物味道浓度；当当前全局最优食物味道浓度不变时，则与之对应的位置向量就是N个果蝇群协同搜索得到的当前全局最优位置向量，将此全局最优位置向量作为小波常模盲均衡方法WTCMA的初始优化权向量；
[0071] 步骤五、在步骤四获得初始优化权向量fopt(0)后，对小波常模盲均衡方法WTCMA的权向量f(k)进行更新。
[0072] 所述步骤四中获得WTCMA的权向量f(k)的初始优化权向量fopt(0)具体步骤如下：
[0073] (401)果蝇群初始化：随机产生N个果蝇群，每个果蝇群包括M只果蝇个体；N个果蝇群的初始位置向量c=(c1,c2...,cd,...,cD)，其中cd为N个果蝇群的第d维初始位置；第n个果蝇群的位置向量Cn=(Cn1,Cn2...,Cnd,...,CnD)，其中Cnd为第n个果蝇群第d维位置；第n个果蝇群中第m只果蝇的位置向量Cnm=(Cnm1,Cnm2,...,Cnmd,...,CnmD)，其中Cnmd为第n个果蝇群中第m只果蝇的第d维位置；n=1,2,...,N，m=1,2,...,M，d=1,2,...,D，N、M、D为正整数，D为位置向量的维数；第n个果蝇群第m只果蝇的位置向量处的食物味道浓度记为S(Cnm)；
[0074] (402)确定食物味道浓度：第n个果蝇群中第m只果蝇的位置向量Cnm对应于小波常模盲均衡方法WTCMA的一个初始权向量f(0)，将WTCMA的代价函数作为第n个果蝇群第m只果蝇的位置向量处的食物味道浓度S(Cnm)，即有
[0075] S(Cnm)=JWTCMA(f(0))；
[0076] 式中，JWTCMA(f(0))是由WTCMA的初始权向量f(0)确定的代价函数；
[0077] (403)在N个果蝇群中，获取每个果蝇群的位置向量：每个果蝇群的位置向量均按第n个果蝇群的位置向量更新公式同时进行更新，即
[0078] Cn=c+Vn；
[0079] 式中，Vn=(Vn1,Vn2,...,Vnd,...,VnD)为第n个果蝇群的步进向量，Vnd=unrand(0,1)为区间[0,un]内的随机数，un为第n个果蝇群的步进因子；
[0080] (404)在N个果蝇群中，获取每个果蝇群中第m只果蝇的位置向量：每个果蝇群中第m只果蝇的位置向量，均按第n个果蝇群中第m只果蝇的位置向量更新公式同时进行更新，即[0081] Cnm=Cn+Unm；
[0082] 式中，Unm=(Unm1,Unm2,...,Unmd,...,UnmD)为第n个果蝇群中第m只果蝇的步进向量，Unmd=unmrand(0,1)为区间[0,unm]内的随机数，unm为第n个果蝇群中第m只果蝇的步进因子；
[0083] (405)采用比较法从N个果蝇群中获取搜索到最低食物味道浓度的果蝇群及其对应的位置向量，获取过程如下：
[0084] 首先获得第n个果蝇群中第m只果蝇个体的食物味道浓度Snm，即
[0085] Snm=S(Cnm)；
[0086] 比较第n个果蝇群中M个果蝇个体的食物味道浓度大小，将其中的最低食物味道浓度记为Sn-opt，且
[0087]
[0088] 式中， 表示从M个值中选取最小值；第n个果蝇群中最低食物味道浓度Sn-opt称为该果蝇群的最优食物味道浓度，与该果蝇群的最优食物味道浓度对应的果蝇个体称为该果蝇群的最优个体；该最优个体的位置向量Cn-opt称为该果蝇群的最优位置向量；
[0089] 再次，从N个果蝇群的最优食物味道浓度中获取全局最低食物味道浓度Sopt，即[0090]
[0091] 该全局最低食物味道浓度称为N个果蝇群的全局最优食物味道浓度，与N个果蝇群的全局最优食物味道浓度对应的位置向量Copt称为N个果蝇群的全局最优位置向量，并作为下一次N个果蝇群协同搜索时，N个果蝇群的初始位置向量，即
[0092] c=Copt；
[0093] (406)重复(403)到(405)的步骤，进行下一次的N个果蝇群协同搜索，再次获取N个果蝇群的当前全局最优食物味道浓度；当当前全局最优食物味道浓度不变时，则与之对应的位置向量就是N个果蝇群协同搜索得到的当前全局最优位置向量，将此当前全局最优位置向量作为小波常模盲均衡方法WTCMA的初始优化权向量fopt(0)。
[0094] 所述步骤五中对小波常模盲均衡方法WTCMA的权向量f(k)进行更新方法如下：即[0095]
[0096] 式中， 为发射信号a(k)的统计模值， 是一个正常数，E表示数学期望；*表示取共轭，μ为步长，
diag表示对角矩阵； (k)表示对小波变换系数rI,K(k)的
平均功率估计， 表示尺度变换系数sI+1,K(k)的平均功率估计，其迭代公式为
[0097]
[0098] 式中，β为平滑因子，且0<β<1，rI,K(k)是尺度参数为I，平移参数为K的小波变换系数，I、K均取正整数；sI+1,K(k)为小波分解的最大尺度为I+1，平移参数为K的尺度变换系数。
[0099] 实施实例：
[0100] 为了检验本发明方法MFOA-WTCMA的性能，以最小均方误差曲线(MSE)作为衡量性能指标，以WTCMA、SFOA-WTCMA为比较对象进行仿真实验。
[0101] 发射信号为16QAM，均衡器权长为16，信噪比为25dB，采用最小相位水声信道h=[0.9656,-0.0906,0.0578,0.2368]；果蝇种群数为2，果蝇种群规模100，果蝇初始化位置[-0.1,0.1]，果蝇种群迭代步进值[-0.05,0.05]，最大进化代数为200；SFOA-WTCMA的步长μSFOA-WTCMA=0.0035，MFOA-WTCMA的步长μMFOA-WTCMA=0.0035；采用DB2小波进行分解，分解层数为2层；初始功率设置为8，遗忘因子β=0.99，600次蒙特卡诺仿真，结果如图2所示。
[0102] 图2中的(a)为最小均方误差曲线（MSE）收敛曲线，表明本发明方法MFOA-WTCMA的稳态误差比WTCMA减小约3.5dB、比SFOA-WTCMA减小约1dB。图2中的(b)、(c)、(d)表明，本发明方法MFOA-WTCMA的输出星座图比WTCMA、SFOA-WTCMA相比，更为紧凑、清晰。因此，本发明方法MFOA-WTCMA的性能优于SFOA-WTCMA，本发明方法具有收敛速度最快、均方误差最小和全局性能最优的特点，在通信技术领域有很强的实用价值。
