发明内容
本发明目的是为了克服CMA的收敛速度慢、均方误差较大及对权向量的初始化的敏感性的缺陷,提出了一种混沌与最速下降法联合优化的正交小波常数模盲均衡方法。该方法通过对均衡器的接收信号进行归一化的正交小波变换,以降低输入信号的相关性,加快收敛速度;针对CMA的局部收敛问题,利用混沌优化的搜索过程按混沌运动自身的规律和特性进行、内在的随机性和遍历性进行高效全局寻优的特点,将其与最速下降法结合对均衡器的权向量进行优化,使优化后权向量的值位于最优点的邻域范围内,再切换至正交小波常数模盲均衡方法,利用该方法对权向量进行更新,最终使权向量收敛至全局最优解。
本发明为实现上述目的,采用如下技术方案:
本发明混沌与最速下降法联合优化的正交小波常数模盲均衡方法,包括如下步骤:
a.)将发射信号a(n)经过脉冲响应信道c(n)得到信道输出向量x(n),其中n为时间序列,下同;
b.)采用信道噪声w(n)和步骤a所述的信道输出向量x(n)得到正交小波变换器(WT)的输入信号y(n):y(n)=w(n)+x(n);
c.)将步骤b所述的正交小波变换器(WT)的输入信号y(n)经过正交小波变换后,则均衡器输入信号R(n):
R(n)=Qy(n) (1)
式中,Q为正交变换矩阵,均衡器输出为:
z(n)=fH(n)R(n) (2)
误差函数e(n)为
e(n)=|z(n)|2-Rcma (3)
式中,Rcma为小波CMA的模,且
Rcma=E{|a(n)|4}/E{|a(n)|2} (4)
代价函数J为
Jf(n)=E[e2(n)] (5)
由最速下降法,得到均衡器权向量的迭代公式为
式中,μ为步长因子,
为正交小波功率归一化矩阵。其中,diag[ ]表示对角矩阵,“*”表示共轭,
和
分别表示对小波系数r
l,k和尺度系数s
L,k的平均功率估计,r
l,k(n)表示小波空间l层分解的第k个信号,s
L,k(n)表示尺度空间中最大分解层数L时的第k个信号,可由下式递推得到
式中,β为平滑因子,且0<β<1。;
其特征在于:
将混沌优化算法与最速下降法结合,利用均衡器接收的一小段数据进行权向量初始化:
令α=[1,2,Λ m],采用16QAM时,m=16,则η=16QAM(α),即η表示对α正交幅度调制后的输出信号。根据式(6),令经过M1次最速下降法迭代优化后的权向量为f*1,均衡器的权长为L,n的取值范围为1ΛN,则
R(n)=Qy(n+L-1:-1:n) (8)
z(n)=fH(n)R(n) (9)
调制误差e(n)为
e(n)=min[|z(n)-η|2] (10)
式中,min( )表示函数最小化。
定义平均调制误差
式中,zk(n)表示第k次混沌优化时均衡器的输出信号z(n),AME(k)表示第k次混沌优化的平均调制误差。
为提高权向量的收敛精度,将权向量的实部和虚部分别作为优化变量进行优化,权向量的复数形式可表示为fi(n)=fi1(n)+jfi2(n),“j”为虚数虚单位,其中
fi1=ci+dixi1(n) (12)
fi2=ci+dixi2(n) (13)
式中,fi(n)表示权向量f中的第i个元素,ci,di为常数,xi1(n)和xi2(n)分别为第i个权向量的实部和虚部所对应的混沌变量。
利用式(12)、式(13)对式(9)中的权向量f的实部和虚部分别进行M2次迭代,每次迭代的过程中调制误差都随着权向量的不同而发生改变,除去e(n)中的最大值和最小值,式(11)修正为
将混沌优化过程中AME(k)的最小值赋给J*2,所对应的权向量为f*2,J*2为混沌优化过程中的代价函数值。优化后切换至正交小波常数模盲均衡方法的切换条件是:
J*2(i-1)-J*2(i)<ζ (15)
式中,ζ为一正数,J*2(i)表示在第i次的混合优化过程中,经过M2次混沌优化后得到的J*2,i的取值范围为1ΛM3。在优化过程中若满足式(15),切换到正交小波常数模盲均衡方法,若不满足此条件,则在进行M3次混合优化后,切换到正交小波常数模盲均衡方法。
传统的常数模方法(CMA)存在收敛速度慢、均方误差较大、不同位置的权向量初始化会收敛到不同的极小值点的缺点。针对其存在的缺点,发明了混沌与最速下降法联合优化的正交小波常数模盲均衡方法(Chaos optimization basedorthogonal Wavelet Transform Constant Modulus blind equalization Algorithm,CWTCMA)。该方法利用混沌优化算法全局寻优,使权向量收敛到全局最优解的邻域,然后切换到正交小波常数模盲均衡方法,进行局部寻优;对均衡器的输入信号进行正交小波变换,加快了权向量收敛速度、减小了均方误差。水声信道的仿真结果表明:与CMA,WT-CMA相比,本发明CWTCMA具有更好的均衡性能。
具体实施方式
由于混沌运动具有随机性、遍历性、规律性的特点,混沌搜索能在一定的范围内按其自身的规律不重复的遍历每一个状态。混沌优化算法就是根据其遍历性和规律性的特点,采用混沌变量在一定的范围内遍历每一个状态,能够使混沌变量的搜索跳出局部极值点。为避免CMA的权向量收敛到局部极小值点,本发明将权向量作为优化变量。
常见的混沌映射函数有Lorenz映射、Logistic映射、改进型的logistic映射、Henon映射等,其中改进的Logistic映射与其它映射相比,具有结构简单、计算量小、使用方便的特点,所以采用改进的logistic映射作为权向量混沌优化的迭代公式。
该映射的表达式为
x(n+1)=1-2x2(n) (1)
式中,x(n)表示混沌变量x在第n次迭代过程中的值,映射的相空间范围为(-1,1)。
由于CMA中权向量的取值范围与改进的logistic映射的遍历空间不同,式(2)将第i个混沌变量xi(n)映射到相应的第i个优化变量中,也就是优化变量fi(n),这样使混沌变量的取值范围“放大”到权向量的取值范围。
fi(n)=ci+dixi(n) (2)
式中,ci,di为常数。
为提高权向量的收敛精度,将权向量的实部和虚部分别作为优化变量进行优化,权向量的复数形式可表示为fi(n)=fi1(n)+jfi2(n),“j”为虚数虚单位,于是将式(2)改写成式(3)和式(4),
fi1(n)=ci+dixi1(n) (3)
fi2(n)=ci+dixi2(n) (4)
式中,xi1(n)和xi2(n)分别为第i个权向量的实部和虚部所对应的混沌变量。
混沌优化常数模盲均衡权向量
对CMA权向量进行混沌优化的基本步骤如下:
1:设置最大的迭代次数M2,并对权向量实部和虚部所对应的混沌变量xi1(n)、xi2(n)赋初值,令其为xi1(0)和xi2(0),其中i=1,ΛL,L为权向量的长度;
2:通过式(3)、式(4),将x
i1(0)和x
i2(0)映射到权向量的优化区间,得到权向量的实部f
i1(0)和虚部f
i2(0),令
所对应的CMA的代价函数J(0);
3:进行混沌搜索,得到f
i(n)和J(n),如果J(n)<J(0),则
J
*=J(n);
4:当n>M
2时,
保持不变,结束;否则令n=n+1,转到步骤3。
混沌与最速下降法联合优化常数模盲均衡权向量
利用混沌优化算法优化权向量理论上可以遍历所有的状态,但优化时间较长,由于CMA中权向量的迭代利用了最速下降法,而将混沌优化算法与最速下降法有机结合的混合算法具有全局收敛和快速收敛的特点。因此,可以将混沌优化算法与最速下降法结合,来优化均衡器权向量。首先通过混沌优化全局寻优,使均衡器权向量的值接近全局最优点,然后采用最速下降法在最优点的邻域范围内局部寻优。利用混合算法搜索,有利于权向量跳出局部最优点,接近全局最优点,并提高收敛的精度。
利用混合算法对盲均衡器的权向量f(n)进行优化,求解代价函数的最小值Jmin,步骤如下:
1:设最速下降法和混沌优化的最大迭代次数分别为M1、M2,混合搜索次数M3,令计数器k=0,初始权向量f(0);
2:以f(0)为初始点,进行M1次最速下降法搜索,得到优化后的均衡器权向量f*1(n)和代价函数值J*1;
3:以f*1(n)为初始点,进行M2次混沌优化搜索得到f*2(n)及J*2;
4:令k=k+1,如果k>M3,优化结束,否则转步骤5;
5:如果J*2<J*1,令f(n)=f*2(n),如果J*2≥J*1,令f(n)=f*1(n),转步骤1;
基于混沌优化算法和最速下降法的计算步骤,如图1所示。
混沌与最速下降法联合优化正交小波盲均衡权向量
为了克服CMA局部收敛的缺点,本发明采用混沌优化算法与最速下降法结合来对权向量进行初始化;采用正交小波函数对均衡器的输入信号进行变换,通过降低输入信号的相关性,来加快收敛速度。这样引入混沌与最速下降法联合优化算法的正交小波变换常数模盲均衡方法,称为混沌与最速下降法联合优化的正交小波常数模盲均衡方法(Chaos optimization based orthogonal Wavelet TransformConstant Modulus blind equalization Algorithm,CWTCMA)。其原理框图,如图2所示。
图2中,a(n)是输入发射序列,c(n)是长度为N
h的基带信道响应向量,w(n)为噪声向量,Q为正交小波变换矩阵,y(n)为信道的输出向量,R(n)为均衡器的输入信号,f(n)为均衡器的权向量,其初始值通过混沌与最速下降法的混合优化得到,z(n)为均衡器的输出,
为判决器的输出。过程如下:
1)正交小波变换盲均衡方法
图2中,正交小波变换后的输出为
R(n)=Qy(n) (5)
z(n)=fH(n)R(n) (6)
误差函数e(n)为
e(n)=|z(n)|2-Rcma (7)
式中,Rcma为小波CMA算法的模,且
Rcma=E{|a(n)|4}/E{|a(n)|2} (8)
代价函数J为
Jf(n)=E[e2(n)] (9)
由最速下降法,得到均衡器权向量的迭代公式为
式中,μ为步长因子,
为正交小波功率归一化矩阵。其中,diag[ ]表示对角矩阵,“*”表示共轭,
和
分别表示对小波系数r
l,k和尺度系数s
L,k的平均功率估计,r
l,k(n)表示小波空间l层分解的第k个信号,s
L,k(n)表示尺度空间中最大分解层数L时的第k个信号,可由下式递推得到
式中,β为平滑因子,且0<β<1,一般取略小于1的数。称式(5)~(11)为基于正交小波的常数模盲均衡方法(Wavelet Transform based CMA,WT-CMA)。
2)权向量的混沌与最速下降法联合优化
将混沌优化算法与最速下降法结合,利用均衡器接收的一小段数据进行权向量初始化。
令α=[1,2,Λm],采用16QAM时,m=16,则η=16QAM(α),即η表示对α正交幅度调制后的输出信号。根据式(10),令经过M1次最速下降法迭代优化后的权向量为f*1,均衡器的权长为L,n的取值范围为1ΛN,则
R(n)=Qy(n+L-1:-1:n) (12)
z(n)=fH(n)R(n) (13)
式中,min( )表示函数最小化。
调制误差e(n)为
e(n)=min[|z(n)-η|2] (14)
定义平均调制误差
式中,zk(n)表示第k次混沌优化时均衡器的输出信号z(n),AME(k)表示第k次混沌优化的平均调制误差。利用式(3)、式(4)对式(13)中的权向量f的实部和虚部分别进行M2次迭代,每次迭代的过程中调制误差都随着权向量的不同而发生改变,除去e(n)中的最大值和最小值,式(15)修正为
将混沌优化过程中AME(k)的最小值赋给J*2,所对应的权向量为f*2。优化后切换至正交小波盲均衡方法的切换条件是:
J*2(i-1)-J*2(i)<ζ (17)
式中,ζ为一正数,J*2(i)表示在第i次的混合优化过程中,经过M2次混沌优化后得到的J*2,i的取值范围为1,ΛM3。在优化过程中若满足式(17),切换到正交小波常数模盲均衡方法,若不满足此条件,则在进行M3次混合优化后,切换到正交小波常数模盲均衡方法。
实施实例
为了验证本发明CWTCMA的有效性,用水声信道进行仿真研究,并与CMA、WT-CMA进行比较。
仿真实验中,采用水声信道[0.3132,-0.104,0.8908,0.3134],信噪比为25dB,均衡器的权长为16。
【实施例1】发射信号为16QAM,CMA、WT-CMA、CWTCMA中步长因子μ分别为0.00001、0.0002、0.0001,M1、M2、M3的值分别为500、800、20,N为20;都采用第4个抽头系数为1,其余的全为0;ci的值都为0、di的值都为1;混沌初始化时采用均衡器输入数据的前500点对权向量进行初始化,初始化切换条件ζ为10-5;蒙特卡罗仿真次数为5000次的仿真结果,如图3所示。
由图3(a)可知,本发明CWTCMA收敛后的均方误差比CMA约小2dB,比WT-CMA约小0.5dB;本发明CWTCMA的收敛速度比CMA快约5000步,比WT-CMA快约1000步;由图3(c)~(e)可知,本发明CWTCMA均衡后的星座图明显比CMA、WT-CMA清晰。
【实施例2】发射信号为16PSK,CMA、WT-CMA、CWTCMA中步长因子μ分别为0.001、0.002、0.001,M1、M2、M3的值分别为300、800、20,N为20;都采用第4个抽头系数为1,其余的全为0;ci的值都为0、di的值都为1;混沌初始化时采用均衡器输入数据的前300点对权向量进行初始化,初始化切换条件ζ为10-5;蒙特卡罗仿真次数为5000次的仿真结果,如图4所示。图4(a)表明,在稳态误差上,本发明CWTCMA比CMA减小约5dB,与WT-CMA基本相同;在收敛速度上,本发明CWTCMA,比CMA快了近4200步比WT-CMA快了约1500步;由图4(c)~(e)可知,本发明CWTCMA均衡后的星座图明显比CMA、WT-CMA清晰。