技术领域:
[0001] 本
发明涉及
数据采集和
信号处理技术领域,具体讲是一种基于混合猴群优化的小波加权多模盲均衡方法。背景技术:
[0002] 目前,
水声通信是人们普遍接受的水下通信方式,水声信道中的多径传播、高背景噪声等因素会使信号在传输过程中产生严重的码间干扰(Inter-symbol Interference,ISI),通信
质量无法得到保证。为解决这一问题,各种均衡技术应运而生。加权多模方法(WMMA)利用判决符号的指数幂构成加权项,可以自适应调
制模值,能有效降低模型误差,但在获取方法非凸性代价函数的全局最优解时,仍沿用了梯度下降的思想,无法克服方法易陷入局部极值的问题,难以进一步提高均衡效果,将
小波变换(WT)引入WMMA,可降低
信噪比,但只能一定程度上提高均衡效果。
[0003] 猴群优化方法(MA)具有良好的全局随机搜索能
力,但局部搜索能力不够强,也容易陷入局部最优解,结合具有优秀搜索路径的列维飞行模式(LF)和具有极强局部搜索能力的模拟
退火优化方法(SA),可提高全局和局部意义下的搜索能力。发明内容:
[0004] 本发明要解决的技术问题是,提供一种可以大大降低均衡器的稳态均方误差,加快收敛速度的基于混合猴群优化的小波加权多模盲均衡方法。
[0005] 本发明的技术解决方案是,提供一种基于混合猴群优化的小波加权多模盲均衡方法,包括基于小波变换的加权多模盲均衡方法以及混合猴群优化方法,所述方法结合在一起后的具体步骤如下:
[0006] 步骤①初始化阶段:先设置基于混合猴群优化的小波加权多模盲均衡方法中的相关参数,接着随机产生初始种群,生成信号y(k),然后确定混合猴群优化方法的适应度函数,将y(k)作为混合猴群优化方法的
输入信号,最后初始化猴群全局最优
位置向量;
[0007] 步骤②混合猴群优化方法的
迭代寻优阶段:通过混合猴群优化方法获得猴群的全局最优位置向量,并将其作为基于小波变换的加权多模盲均衡方法的初始权向量;
[0008] 步骤③信号均衡输出阶段:对上述步骤①中的输入信号y(k)通过小波加权多模盲均衡方法均衡输出。
[0009] 优选地,根据本发明所述的一种基于混合猴群优化的小波加权多模盲均衡方法,其中,初始化阶段的具体步骤如下:
[0010] 步骤a设置基于混合猴群优化的小波加权多模盲均衡方法中所有相关参数;
[0011] 步骤b生成信号y(k):信号经信道传输后将发生一定变化。
[0012] y(k)=cT(k)a(k)+b(k) (1)
[0013] y(k)=cT(k)a(k)+b(k)
[0014] 式(1)中,a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列。
[0015] 步骤c随机产生初始种群:在n维空间创建一个规模为m的猴群X=[X1,X2,…,Xm],则第i只猴子的位置可以用一个n维的向量Xi=(xi1,xi2,…,xin),i=1,2,…,m来表示,每只人工猴的位置向量与基于小波变换的加权多模盲均衡方法中的权向量设置为相同形式。猴群初始位置的分配采用的是随机形式,过程如下:
[0016] xij=xmin,j+(xmax,j-xmin,j)rand (2)
[0017] 式(2)中,i=1,2,…,m,j=1,2,…,n;xij为第i只猴子在第j维的实际位置,xmin,j和xmax,j分别表示搜索空间第j维的下界和上界,rand产生一个在区间[0,1]上的实数;
[0018] 步骤d确定适应度函数:将混合猴群优化方法中适应度函数的倒数对应于小波加权多模盲均衡方法(WT-WMMA)的代价函数,两者关系如下:
[0019]
[0020] 式(2)中,L(Xi)为基于小波变换的加权多模盲均衡方法的代价函数,f(Xi)为混合猴群优化方法的适应度函数。基于小波变换的加权多模盲均衡方法的最佳权向量是要使其代价函数取得最小值,而混合猴群优化方法中最终取得的全局最优解是适应度函数的最大值,故利用混合猴群优化方法求解均衡器的代价函数时将两者作以上对应;
[0021] 步骤e初始化猴群全局最优位置向量:将y(k)作为混合猴群优化方法的输入信号,根据适应度函数,分别计算初始状态中每只猴子的位置向量对应的适应度函数值,比较结果,将该群体中的最大适应度函数值定义为全局最优适应度函数值f(X*)的初始值,该值对应的位置向量定义为全局最优位置向量X*的初始值。
[0022] 优选地,根据本发明所述的一种基于混合猴群优化的小波加权多模盲均衡方法,其中,混合猴群优化方法的迭代寻优阶段的具体步骤如下(以第i只猴子为例):
[0023] 步骤a爬过程,具体如下:
[0024] 1)在第t次的爬行中,随机生成一个向量ΔXi=(Δxi1,Δxi2,…,Δxin),i=1,2,…,m,分量Δxij以相同的概率0.5取爬步长λi(λi>0)或-λi,爬步长将采用列维飞行模式确定:
[0025]
[0026] Levy~ε=t-1-β(0<β≤2) (5)
[0027] 式(4)中,α为步长控制因子(α>0,具体数值可根据实际情况设置); 为点对点乘法;L(β)为随机搜索路径,它移动的步长服从列维分布。
[0028] 2)计算伪梯度
[0029]
[0030] 式(6)中,i=1,2,…,m,j=1,2,…,n,向量f'i(xi)=(f'i1(xi),f'i2(xi),…f'in(xi))称为适应度函数在点Xi处的伪梯度;
[0031] 3)设向量Y=(y1,y2,…yn),向量中各分量为
[0032] yj=xij+λi·sign(f'ij(xi)) (7)
[0033] 式(7)中,i=1,2,…,m,j=1,2,…,n。
[0034] 4)若Y满足搜索空间范围[a,b]内,则更新位置,Xi←Y,否则,保持Xi不变。
[0035] 5)返回爬过程的步骤1),继续进行迭代,直到达到设置的最大爬次数或者前后两次迭代过程中适应度函数值无明显变化,爬过程结束,转入望-跳过程;
[0036] 步骤b望-跳过程,具体步骤如下:
[0037] 1)经过一段时间的攀爬,每只人工猴在其
视野范围(xij-γ,xij+γ),i=1,2,…,m,j=1,2,…,n内进行多次眺望以探索周围是否有更好的位置,每次眺望的结果用Y=(y1,y2,…yn)表示,若Y在搜索空间范围[a,b]内并有f(Y)>f(Xi),则Xi←Y;否则,保持Xi不变;
[0038] 2)重复望-跳过程步骤1)直到有满足条件的Y产生时或者达到设定次数终止;
[0039] 3)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
[0040] 步骤c模拟退火操作,具体步骤如下:
[0041] 1)在第i只猴子当前位置一定范围内进行随机扰动,生成一个满足约束条件的位置向量Y=(y1,y2,…yn);
[0042] 2)若f(Y)>f(Xi),则将其作为新解,更新人工猴位置,否则根据Metropolis准则,以状态转移概率P决定是否接受新解;
[0043]
[0044] 式(8)中:fk为第k次迭代所得的最优解的适应度函数值,k=0,1,…,K;K为
马尔科夫链的长度;fk+1为第(k+1)次迭代所得的最优解的适应度函数值;P(Tk+1)为
温度为Tk+1时的接收概率;
[0045] 3)根据约束条件判断是否结束,若结束,转入步骤4);否则转入步骤1);
[0046] 4)若未达到冷却状态,则按式(9)进行降温处理,
[0047] Ti+1=Ti·α (9)
[0048] 式(9)中:为温度冷却系数.若达到冷却状态,转入步骤5);否则转入步骤1);
[0049] 5)计算所有猴子在模拟退火操作中所有位置向量对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
[0050] 6)若达到设定次数,进入翻过程;否则,返回爬过程,进行迭代寻优;
[0051] 步骤d翻过程的主要目的是迫使猴群开拓新的搜索区域。经过以上迭代,每只猴子在翻区间[c,d]内沿当前位置指向
支点的方向(当θ>0时)或反方向(当θ<0时)翻跳以更新当前的位置,具体步骤如下:
[0052] 1)在翻区间[c,d]内随机产生一个实数θ,令
[0053] yj=xij+θ(pj-xij) (10)
[0054] 式(10)中,i=1,2,…,m,j=1,2,…,n,支点 设Y=(y1,y2,…yn),若Y在搜索空间范围[a,b]内并有f(Y)>f(Xi),则Xi←Y;否则,保持Xi不变;
[0055] 2)重复翻过程步骤1)直到有满足条件的Y产生或者达到最大设置次数时终止;
[0056] 3)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f* *(X),用该值对应的位置向量更新X;
[0057] 3)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
[0058] 4)检验是否满足结束条件,若满足则转入步骤e;否则,转到爬过程;
[0059] 步骤e输出全局最优位置向量 该向量即为小波加权多模盲均衡方法中权向量F(k)的初始值,迭代寻优阶段结束。
[0060] 优选地,根据本发明所述的一种基于混合猴群优化的小波加权多模盲均衡方法,其中,信号均衡输出阶段的具体步骤如下:
[0061] 步骤a将信号y(k)分为
实部和
虚部,并分别对其进行小波变换。
[0062] 根据小波变换理论,若设Q为小波变换矩阵,则经过小波变换后均衡器的输入信号R(k)为:
[0063] y(k)=yRe(k)+j·yIm(k) (11)
[0064] R(k)=RRe(k)+i·RIm(k)=yRe(k)Q+i·yIm(k)Q (12)
[0065]
[0066]
[0067]
[0068]
[0069] 式中,下标Re和Im分别代表参数的实部和虚部,uj,m(k)、sJ,m(k)分别表示相应的小波和尺度变换系数,均衡器长度为L,J为小波分解的最大尺度,L=2J,k=0,1,…,L-1,kj=L/2j-1(j=1,2,…,J)为在尺度j下小波函数的最大平移;φJ,m(k)为小波函数; 为尺度函数。
[0070] 步骤b信号均衡输出:
[0072]
[0073] WT-WMMA误差函数为
[0074]
[0075]
[0076] WT-WMMA代价函数为
[0077]
[0078]
[0079]
[0080] 式中,加权因子λRe,λIm∈[0,2],对方形QAM信号
星座图,选择λRe=λIm,λRe和λIm分别是信号实部和虚部的加权因子,仿真实验中λRe和λIm的值是根据仿真实例中经验得到;和 分别是判决信号 的实部与虚部, 与z(k)近似相等。
[0081] 权向量F(k)的迭代
[0082] F(k)=FRe(k)+jFIm(k) (21a)
[0083]
[0084]
[0085]
[0086]
[0087] 式(21b)-(21c)中:μ为权向量的迭代步长,R*(k)为R(k)的共轭,表示对尺度变换系数mI,n
(k)、小波变换系数ri,n(k)的
平均功率估计,其迭代公式为:
[0088]
[0089]
[0090] 式(22)和(23)中,η(0<η<1)为平滑因子。
[0091] 本发明的有益效果是:加权多模方法(WMMA)中的加权项在均衡器系数的迭代更新过程中,能够自适应地调整模值,同等条件下均衡效果更好,结合小波变换WT降
低信号自相关性后,均衡效果进一步提高。猴群优化方法(MA)采用列维飞行模式(LF)确定爬步长,并与模拟退火优化方法(SA)有机结合后,提高了方法的全局寻优能力,提高了寻优结果的
精度,用该方法来搜寻WT-WMMA中代价函数的全局最优解,能有效逃离局部极值,均衡器性能再一次得到提高,收敛速度加快的同时稳态均方误差也得到了有效的控制。
附图说明:
[0092] 图1是本发明的原理图;
[0094] 图3是使用本发明方法后的收敛曲线图;
[0095] 图4是使用本发明方法后的WT-MMA输出星座图;
[0096] 图5是使用本发明方法后的WT-WMMA输出星座图;
[0097] 图6是使用本发明方法后的MA-WT-WMMA输出星座图;
[0098] 图7是使用本发明方法后的LFSA-MA-WT-WMMA输出星座图。具体实施方式:
[0099] 下面结合附图和具体实施方式对本发明一种基于混合猴群优化的小波加权多模盲均衡方法作进一步说明:
[0100] 本发明的原理
框图如图1所示,图1中,a(k)是信源发射复信号;c(k)是长度为L的信道脉冲响应向量,c(k)=[c(k),L,c(k-L+1)]T;b(k)是通信过程中的加性高斯白噪声;y(k)是长度为L的经信道传输后的复信号序列,y(k)=[y(k+L),L,y(k-L)]T,R(k)是y(k)经过小波变换后的复信号;F(k)为均衡器的权向量,F(k)=[F0(k),L,FL(k)]T(上标T表示转置);z(k)为均衡器的输出复信号;下标Re和Im分别代表参数的实部和虚部。
[0101] 如图2所示,本发明一种基于混合猴群优化的小波加权多模盲均衡方法的具体步骤如下:
[0102] 一、初始化阶段
[0103] 具体步骤如下:
[0104] 步骤a设置基于混合猴群优化的小波加权多模盲均衡方法中所有相关参数;
[0105] 步骤b生成信号y(k):信号经信道传输后将发生一定变化。
[0106] y(k)=cT(k)a(k)+b(k) (1)
[0107] 式(1)中,a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列。
[0108] 步骤c随机产生初始种群:在n维空间创建一个规模为m的猴群X=[X1,X2,…,Xm],则第i只猴子的位置可以用一个n维的向量Xi=(xi1,xi2,…,xin),i=1,2,…,m来表示,每只人工猴的位置向量与基于小波变换的加权多模盲均衡方法中的权向量设置为相同形式。猴群初始位置的分配采用的是随机形式,过程如下:
[0109] xij=xmin,j+(xmax,j-xmin,j)rand (2)
[0110] 式(2)中,i=1,2,…,m;j=1,2,…,n,xij为第i只猴子在第j维的实际位置,xmin,j和xmax,j分别表示搜索空间第j维的下界和上界,rand产生一个在区间[0,1]上的实数;
[0111] 步骤d确定适应度函数:将混合猴群优化方法中适应度函数的倒数对应于小波加权多模盲均衡方法(WT-WMMA)的代价函数,两者关系如下:
[0112]
[0113] 式(3)中,L(Xi)为基于小波变换的加权多模盲均衡方法的代价函数,f(Xi)为混合猴群优化方法的适应度函数。基于小波变换的加权多模盲均衡方法的最佳权向量是要使其代价函数取得最小值,而混合猴群优化方法中最终取得的全局最优解是适应度函数的最大值,故利用混合猴群优化方法求解均衡器的代价函数时将两者作以上对应;
[0114] 步骤e初始化全局最优位置向量:在本方法中,每只猴子的位置向量对应盲均衡器的一个权向量。将y(k)作为混合猴群优化方法的输入信号,根据适应度函数,分别计算每只猴子的初始位置向量对应的适应度函数值,比较结果,将该群体中的最大适应度函数值定义为全局最优适应度值f(X*)的初始值,该值对应的位置向量定义为全局最优位置向量X*的初始值。
[0115] 二、混合猴群优化方法的迭代寻优阶段
[0116] 具体步骤如下(以第i只猴子为例):
[0117] 步骤a爬过程。猴群优化方法的爬过程主要是利用伪梯度思想,计算适应度函数在当前位置相邻两点的函数值,继而进行比较移动的过程,具体如下:
[0118] 1)在第t次的爬行中,随机生成一个向量ΔXi=(Δxi1,Δxi2,…,Δxin),i=1,2,…,m,分量Δxij以相同的概率0.5取爬步长λi(λi>0)或-λi,爬步长将采用列维飞行模式确定:
[0119]
[0120] Levy~ε=t-1-β(0<β≤2) (5)
[0121] 式(4)中,α为步长控制因子(α>0,具体数值可根据实际情况设置);为点对点乘法;L(β)为随机搜索路径,它移动的步长服从列维分布。
[0122] 2)计算伪梯度
[0123]
[0124] 式(6)中,i=1,2,…,m,j=1,2,…,n,向量f'i(xi)=(f'i1(xi),f'i2(xi),…f'in(xi))称为适应度函数在点Xi处的伪梯度。
[0125] 3)设向量Y=(y1,y2,…yn),向量中各分量为
[0126] yj=xij+λ·sign(f'ij(xi)) (7)
[0127] 式(7)中,i=1,2,…,m,j=1,2,…,n。
[0128] 4)若Y满足搜索空间范围[a,b]内,则更新位置,Xi←Y,否则,保持Xi不变。
[0129] 5)返回爬过程的步骤1),继续进行迭代,直到达到设置的最大爬次数或者前后两次迭代过程中适应度函数值无明显变化,爬过程结束,转入望-跳过程;
[0130] 步骤b望-跳过程,每只猴子通过爬过程达到各自山峰的最高位置,在视野长度为γ的范围内,通过向四周多次眺望,寻找周围是否存在比当前山峰更好的位置,如果有,则跳离当前位置到更好的位置上去。具体步骤如下:
[0131] 1)在视野范围(xij-γ,xij+γ),i=1,2,…,m,j=1,2,…,n内多次眺望,每次眺望结果用随机产生Y=(y1,y2,…yn)来表示,若Y在搜索空间范围[a,b]内并有f(Y)>f(Xi),则Xi←Y;否则,保持Xi不变。
[0132] 2)重复望-跳过程步骤1)直到有满足条件的Y产生时或者达到设定次数终止。
[0133] 3)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*。
[0134] 步骤c模拟退火操作,具体步骤如下:
[0135] 1)在第i只猴子当前位置一定范围内进行随机扰动,生成一个满足约束条件的位置向量Y=(y1,y2,…yn),
[0136] 2)若f(Y)>f(Xi),则将其作为新解,更新人工猴位置,否则根据Metropolis准则,以状态转移概率P决定是否接受新解。
[0137]
[0138] 式(8)中:fk为第k次迭代所得的最优解的适应度函数值,k=0,1,…,K;K为马尔科夫链的长度;fk+1为第(k+1)次迭代所得的最优解的适应度函数值;P(Tk+1)为温度为Tk+1时的接收概率;
[0139] 3)根据约束条件判断是否结束,若结束,转入步骤4);否则转入步骤1);
[0140] 4)若未达到冷却状态,则按式(9)进行降温处理,
[0141]
[0142] 式(9)中:为温度冷却系数.若达到冷却状态,转入步骤5);否则转入步骤1);
[0143] 5)计算所有猴子在模拟退火操作中所有位置向量对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
[0144] 6)若达到设定次数,进入翻过程;否则,返回爬过程,进行迭代寻优;
[0145] 步骤d翻过程的主要目的是迫使猴群开拓新的搜索区域。经过以上迭代,每只猴子在翻区间[c,d]内沿当前位置指向支点的方向(当θ>0时)或反方向(当θ<0时)翻跳以更新当前的位置,具体步骤如下:
[0146] 1)在翻区间[c,d]内随机产生一个实数θ,令
[0147] yj=xij+θ(pj-xij) (10)
[0148] 式(10)中,i=1,2,…,m,j=1,2,…,n,支点 设Y=(y1,y2,…yn),若Y在搜索空间范围[a,b]内并有f(Y)>f(Xi),则Xi←Y;否则,保持Xi不变。
[0149] 2)重复翻过程步骤1)直到有满足条件的Y产生或者达到最大设置次数时终止。
[0150] 3)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
[0151] 4)检验是否满足结束条件,若满足则转入步骤e;否则,转到爬过程;
[0152] 步骤e输出全局最优位置向量 该向量即为小波加权多模盲均衡方法中权向量F(k)的初始值,迭代寻优阶段结束。
[0153] 三、信号均衡输出阶段
[0154] 具体步骤如下:
[0155] 步骤a将信号y(k)分为实部和虚部,并分别对其进行小波变换。
[0156] 根据小波变换理论,若设Q为小波变换矩阵,则经过小波变换后均衡器的输入信号R(k)为:
[0157] y(k)=yRe(k)+j·yIm(k) (11)
[0158] R(k)=RRe(k)+i·RIm(k)=yRe(k)Q+i·yIm(k)Q (12)
[0159]
[0160]
[0161]
[0162]
[0163] 式中,下标Re和Im分别代表参数的实部和虚部,uRej,m(k)、uImj,m(k)、sReJ,m(k)、sImJ,m(k)分别为相应的小波和尺度变换的实部和虚部;均衡器长度为L,J为小波分解的最大尺度,L=2J,k=0,1,…,L-1,kj=L/2j-1(j=1,2,…,J)为在尺度j下小波函数的最大平移;φJ,m(k)为小波函数; 为尺度函数。
[0164] 步骤b信号均衡输出:
[0165] 均衡器的输出信号为
[0166]
[0167] WT-WMMA误差函数为
[0168]
[0169]
[0170] WT-WMMA代价函数为
[0171]
[0172]
[0173]
[0174] 式中,加权因子λRe,λIm∈[0,2],对方形QAM信号星座图,选择λRe=λIm,λRe和λIm分别是信号实部和虚部的加权因子,仿真实验中λRe和λIm的值是根据仿真实例中经验得到;和 分别是判决信号 的实部与虚部, 与z(k)近似相等。
[0175] 权向量F(k)的迭代
[0176] F(k)=FRe(k)+jFIm(k) (21a)
[0177]
[0178]
[0179]
[0180]
[0181] 式(21b)-(21c)中:μ为权向量的迭代步长,R*(k)为R(k)的共轭,dRej,m、dImj,m、vReJ,m和vImJ,m分别为均衡器权系数的实部和虚部,表示对尺度变换系数mI,n
(k)、小波变换系数ri,n(k)的平均功率估计,其迭代公式为:
[0182]
[0183]
[0184] 式(22)和(23)中,η(0<η<1)为平滑因子。
[0185] 四、仿真试验
[0186] 为了验证LFSA-MA-WT-WMMA的有效性,以MMA、WT-WMMA、MA-WT-WMMA为比较对象,进行仿真实验。仿真中,水声信道采用h=[0.9656,-0.0906,0.0578,0.2368],发射信号为256QAM信号,信号
采样点均为10000点,盲均衡器的权长均为16,信噪比均25dB;对信道的输入信号采用DB2小波进行分解,分解层次是2层,功率初始值设置为4,遗忘因子β=0.99;加权因λ=1.3,其步长μMMA=0.0000075,μWT-WMMA=0.000034,μMA-WT-WMMA=0.000045,μLFSA-MA-WT-WMMA=0.0000051,500次蒙特卡诺仿真结果如图所示。从图3、图4、图5、图6和图7中明显看出,LFSA-MA-WT-WMMA的稳态均方误差(MSE)比MA-WT-WMMA小近2dB,比WT-WMMA小近
3dB,比MMA小近4dB,MSE得到了有效控制;LFSA-MA-WT-WMMA的收敛速度比MMA、WT-WMMA、MA-WT-WMMA明显加快,并LFSA-MA-WT-WMMA的输出星座图最为清晰、紧凑,恢复出的传输信号最准确。
[0187] 以上所述的实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种
变形和改进,均应落入本发明
权利要求书确定的保护范围内。