[0072] 下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
[0073] 如图1所示,常模盲均衡方法的原理图图中,a(k)为发射信号,h(k)为传输信道,s(k)为经过信道后的输出信号,n(k)为信道噪声,x(k)为均衡器的接收信号,c(k)为均衡器的权向量,e(k)为误差信号,z(k)为经过均衡后的输出信号。
[0074] 均衡器的接收信号为
[0075] x(k)=s(k)+n(k)=a(k)h(k)+n(k)
[0076] 将输入信号经过均衡器得到输出信号为
[0077] z(k)=c(k)x(k)
[0078] 误差信号为
[0079] e(k)=RCM-|z(k)|2
[0080] 式中,RCM为统计模值,定义为,
[0081]
[0082] CMA的代价函数定义为
[0083] J(k)=E{[RCM-|z(k)|2]2}
[0084] 因此,常模盲均衡器的权向量c(k)的迭代公式为,
[0085] c(k+1)=c(k)-μe(k)x(k)
[0086] 式中,μ迭代步长。
[0087] 如图2所示,为基于DNA萤火虫优化的小波盲均衡方法的原理图,首先发射信号a(k)经过传输信道h(k)后加入信道噪声n(k),得到CMA盲均衡器输入信号,将其作为DNA萤火虫优化方法的输入,并且把CMA的代价函数经适当变换后作为DNA萤火虫优化方法的适应度函数,利用DNA萤火虫优化方法的寻优能力来寻找盲均衡器初始最优权向量,具体步骤如下:
[0088] 步骤1,将发射信号a(k)经过传输信道h(k)后得到信号s(k),再加入信道噪声n(k),得到均衡器的接收信号x(k),再把均衡器的接收信号经小波变换WT后作为DNA萤火虫优化方法的输入信号y(k),;
[0089] 其中,s(k)=h(k)a(k),x(k)=h(k)a(k)+n(k),y(k)=Qx(k),k为非负整数且表示时间序列,Q为小波变换矩阵。
[0090] 步骤2,通过DNA萤火虫优化方法得到萤火虫种群的最优位置向量,并将其作为均衡器的初始最优权向量c(0),c(0)是k=0时c(k)的值;
[0091] c(k)的更新公式为
[0092] 其中,e(k)为误差信号,是由误差函数生成,μ为小波盲均衡方法的迭代步长,上标“*”表示复共轭;
[0093] 为小波变换的平均功率估计,diag[]表示对角矩阵, 和 分别表示对尺度变换
系数mj,k(k)和小波变换系数yj,k(k)的平均功率估计,kJ=L/2J-1,J为小波分解的最大尺度,L为均衡器的权长,c(k+1),c(k)分别为k+1、k时刻的最优权向量。
[0094] 步骤3,将步骤1所述的输入信号y(k)经过均衡器得到输出信号z(k);
[0095] 其中,z(k)=y(k)c(k)。
[0096] 如图3所示,通过DNA萤火虫优化方法得到萤火虫种群的最优位置向量的步骤如下:
[0097] 步骤2-1,设置DNA萤火虫优化方法的初始种群及参数;
[0098] 在一个n维搜索空间中,创建一个规模为D的初始萤火虫群的初始位置向量W=[W1,W2,···,WD];其中,Wi表示第i个萤火虫的位置向量,0<i≤D;
[0099] 每只萤火虫都具有相同的初始荧光素值、初始动态决策范围以及一个随机分配的n维位置向量Wi=(wi1,wi2,···,win),每个初始位置向量对应着一组均衡器权向量系数;
[0100] 定义最大进化代数sg1,定义最大迭代次数sg2;
[0101] 步骤2-2,确定适应度函数;
[0102] 将常模盲均衡方法CMA的代价函数的倒数定义为DNA萤火虫优化方法的适应度函数,
[0103]
[0104] 其中,J(Wi)为CMA方法的代价函数,RCM为CMA的统计模值;
[0105] 步骤2-3,计算萤火虫群中每个萤火虫位置向量的适应度函数值并将适应度函数值从大到小排列;
[0106] 其中,前一半对应的是优质萤火虫种群,后一半则为劣质萤火虫种群;选取适应度函数值最大的位置向量为当前最优位置向量,对应的萤火虫个体为最优萤火虫个体;
[0107] 步骤2-4,编码;
[0108] 将萤火虫个体的位置向量转换为十进制位置向量,从而得到了DNA序列位置向量;
[0109] 步骤2-5,执行交叉操作和变异操作,获得新一代萤火虫群;
[0110] 产生一个随机数rand1∈(0,1),与置换交叉概率p1比较,若rand1<p1,则执行置换交叉操作;再次产生一个随机数rand2∈(0,1),与转位交叉概率p2比较,若rand2<p2,则执行转位交叉操作;再产生一组与萤火虫个体DNA序列位置向量维数相同的(0,1)上的随机数,这组随机数中的元素与DNA序列位置向量中的元素一一对应,将所有随机数分别与变异概率pm比较,若随机数小于pm,则执行变异操作;
[0111] 步骤2-6,判断是否达到最大进化代数sg1,如果否,则转至步骤2-7,否则,转至步骤2-8;
[0112] 步骤2-7,将新一代萤火虫群中所有萤火虫个体的DNA序列位置向量解码,得到所有萤火虫个体的位置向量,转至步骤2-3;
[0113] 步骤2-8,获得最优萤火虫群;
[0114] 步骤2-9,计算最优萤火虫群中所有萤火虫位置向量的适应度函数值,选取适应度函数值最大的位置向量为当前最优位置向量;
[0115] 步骤2-10,更新萤火虫荧光素值,确定邻域,更新萤火虫位置和动态决策范围,进行适应度函数值计算,选取适应度函数值最大的位置向量为当前最优位置向量;
[0116] 步骤2-11,比较步骤2-10和步骤2-9中最优位置向量的适应度函数值,取适应度函数值大的位置向量为全局最优位置向量;
[0117] 步骤2-12,判断是否达到最大迭代次数sg2,如果是,则转至步骤2-13,如果否,则转至步骤2-10;
[0118] 步骤2-13,获得全局最优位置向量,及对应的最优萤火虫个体。
[0119] 步骤2-4中的编码过程如下:
[0120] 步骤2-4-1,由第i个萤火虫的位置向量Wi=(wi1,wi2,···,win)计算得到十进制位置过渡向量Bi=(bi1,bi2,···,bin);
[0121] 其中,wig表示第i个萤火虫的位置向量Wi中第g个位置值,big表示十进制位置过渡向量中第g个位置值,1≤g≤n且g为整数, d为编码长度,Dmaxg和Dming分别为第i个萤火虫的位置向量Wi中第g个位置的最大值和最小值;
[0122] 步骤2-4-2:将十进制位置过渡向量中第g个位置值big转换成一串四进制数sig,则第i个萤火虫的DNA序列位置向量 由n串四进制数组成;
[0123] 其中,长度为d, 表示第i个萤火虫的DNA序列位置向量Si中第g个子整数串中第t′位的数字,1≤t′≤d且t′为整数,d为串的长度,串的长度与编码长度一致。
[0124] 步骤2-5中交叉操作和变异操作的具体过程为:
[0125] 置换交叉操作:从DNA序列位置向量中任意选取两个作为父体,从两个父体中随机选取相同位置的交叉点,并将交叉点间的碱基串的位置进行互换,在交叉完成后,从两个新的序列位置向量中分别随机选取一段碱基数目相等的转座子序列位置向量,并进行置换,最终产生两个新的序列位置向量代替父体;选取父体时,从优质种群选取;
[0126] 转位交叉操作:从DNA序列位置向量中任意选取一个作为父体,从父体中截取一段任意碱基数目的转座子序列位置向量插入截取后父体的任意位置,产生一个新的序列位置向量代替父体;选取父体时,从优质种群选取;
[0127] 变异操作:从DNA序列位置向量中任意选取一个序列位置向量,将该序列位置向量中任一碱基以概率pm被其它碱基代替,从而产生一个新的序列位置向量。
[0128] 步骤2-7中DNA序列位置向量解码过程如下:
[0129] 步骤2-7-1,将第i个萤火虫的DNA序列位置向量解码为十进制位置过渡向量Bi=(bi1,bi2,···,bin),
[0130] 步骤2-7-2:通过以下公式按比例将big转换成第i个萤火虫的位置向量Wi中第g个位置值wig;
[0131]
[0132] 步骤2-10中萤火虫的邻域确定的具体过程为:
[0133] 步骤2-10-1,首先通过以下公式确定第i个萤火虫在第t次迭代寻优时的荧光素值li(t);
[0134] li(t)=(1-ρ)li(t-1)+βF(Wi(t))
[0135] 其中,ρ∈(0,1)为荧光素挥发率,β表示荧光素的更新率,F(Wi(t))表示适应度函数值,Wi(t)表示第i个萤火虫在第t次迭代寻优时的位置向量;
[0136] 步骤2-10-2,其次通过以下公式确定第i个萤火虫在第t次迭代寻优时的邻域Ni(t);
[0137] Ni(t)={i′:||Wi′(t)-Wi(t)||<Ri(t);li′(t)>li(t)}
[0138] 其中,||||为范数,0<i′≤D,Ri(t)为第i个萤火虫在第t次迭代寻优时的动态决策域范围。
[0139] 步骤2-10中萤火虫的动态决策域确定的具体过程为,
[0140] 动态决策域的半径按以下公式进行调整;
[0141] Ri(t+1)=min{Rs,max{0,Ri(t)+γ(nt-|Ni(t)|)}}
[0142] 其中,Ri(t+1)为第i个萤火虫在第t+1次迭代寻优过程中的动态决策域的半径,Rs为萤火虫的可视范围,γ为邻域变化率,nt为邻居阈值。
[0143] 上述方法将均衡器的输入信号作为DNA遗传优化方法的输入,并把CMA的代价函数经适当变换后作为DNA遗传优化方法的适应度函数,利用DNA遗传优化方法的寻优能力来寻找盲均衡器的最优初始权向量;利用DNA遗传优化方法与萤火虫优化方法的优点,将DNA遗传优化方法引入到萤火虫优化方法中,得到了一种基于DNA遗传优化的萤火虫方法,利用这个新型方法来优化萤火虫种群的初始值,从而优化了萤火虫种群的最优位置向量,优化后的最优位置向量作为小波盲均衡的初始化最优权向量,从而加快了其的收敛速度,减小了均方误差;与现有技术相比,本发明具有收敛速度最快、均方误差最小和全局最优的性能,在通信技术领域有很强的实用价值。
[0144] 为了验证DNA-GSO-WTCMA的有效性及稳定性,分别对CMA、WTCMA、GSO-WTCMA、DNA-GSO-WTCMA进行比较,仿真结果如图4所示。各参数为:信道h1=[0.005 0.009 -0.024 0.854 -0.218 0.049 -0.016],发射信号为16QAM信号,信道均衡器权长为16,信号采样点均为5000点,信噪比为25dB;在CMA中,将第9个抽头系数设置为1,其余为0,其步长μCMA=
0.00004;在WTCMA中,取第7个抽头系数为1,其余为0,其步长μWCMA=0.0004;在GSO-WTCMA中,其步长为μGSO-WTCMA=0.00008;而在DNA-GSO-WTCMA中,设置的步长则为μDNA-GSO-WTCMA=
0.00001。输入信号采用DB2小波进行分解,分解层次为2层,功率初始值设置为8,遗忘因子β=0.999。经过1000次的蒙特卡诺仿真后过得到的结果如图4所示。
[0145] 由图4和6可知,均方误差方面,DNA-GSO-WTCMA的均方误差要比GSO-WTCMA的小1dB左右,比WTCMA的小近3dB,比CMA减小了4dB左右,均方误差得到了有效降低;在收敛速度方面,DNA-GSO-WTCMA和GSO-WTCMA比WTCMA快了1500步左右,比CMA快了近1000步,收敛速度明显加快。由图7和8可知,在星座图方面DNA-GSO-WTCMA比其他两个更加紧凑,恢复出的传输信号更准确。
[0146] 可见,使用本发明方法DNA-GSO-WT-CMA可以显著提高盲均衡方法的收敛速度和减少均方稳态误差,在通信技术领域有很强的实用价值。
[0147] 以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。