[发明专利]一种基于DNA萤火虫优化的小波盲均衡方法

权利要求书

1.一种基于DNA萤火虫优化的小波盲均衡方法，其特征在于：包括以下步骤，

步骤1，将发射信号a(k)经过传输信道h(k)后得到信号s(k)，再加入信道噪声n(k)，得到均衡器的接收信号x(k)，再把均衡器的接收信号经小波变换WT后作为DNA萤火虫优化方法的输入信号y(k)；

其中，s(k)＝h(k)a(k)，x(k)＝h(k)a(k)+n(k)，y(k)＝Qx(k)，k为非负整数且表示时间序列，Q为小波变换矩阵；

步骤2，通过DNA萤火虫优化方法得到萤火虫种群的最优位置向量，并将其作为均衡器的初始最优权向量c(0)，c(0)是k＝0时c(k)的值；

所述c(k)的更新公式为

其中，e(k)为误差信号，是由误差函数生成，μ为迭代步长，上标“*”表示复共轭；

为小波变换的平均功率估计，diag[]表示对角矩阵，和分别表示对尺度变换系数m_j,k(k)和小波变换系数y_j,k(k)的平均功率估计，k_J＝L/2^J-1，J为小波分解的最大尺度，L为均衡器的权长，c(k+1),c(k)分别为k+1、k时刻的最优权向量；

步骤3，将步骤1所述的输入信号y(k)经过均衡器得到输出信号z(k)；

其中，z(k)＝y(k)c(k)；

通过DNA萤火虫优化方法得到萤火虫种群的最优位置向量的步骤如下，

步骤2-1，设置DNA萤火虫优化方法的初始种群及参数；

在一个n维搜索空间中，创建一个规模为D的初始萤火虫群的初始位置向量W＝[W₁,W₂,…,W_D]；其中，W_i表示第i个萤火虫的位置向量，0＜i≤D；

每只萤火虫都具有相同的初始荧光素值、初始动态决策范围以及一个随机分配的n维位置向量W_i＝(w_i1,w_i2,…,w_in)，每个初始位置向量对应着一组均衡器权向量系数；

定义最大进化代数sg1，定义最大迭代次数sg2；

步骤2-2，确定适应度函数；

将常模盲均衡方法CMA的代价函数的倒数定义为DNA萤火虫优化方法的适应度函数，

其中，J(W_i)为CMA方法的代价函数，R_CM为CMA的统计模值；

步骤2-3，计算萤火虫群中每个萤火虫位置向量的适应度函数值并将适应度函数值从大到小排列；

其中，前一半对应的是优质萤火虫种群，后一半则为劣质萤火虫种群；选取适应度函数值最大的位置向量为当前最优位置向量，对应的萤火虫个体为最优萤火虫个体；

步骤2-4，编码；

将萤火虫个体的位置向量转换为十进制位置向量，从而得到了DNA序列位置向量；

步骤2-5，执行交叉操作和变异操作，获得新一代萤火虫群；

产生一个随机数rand1∈(0,1)，与置换交叉概率p₁比较，若rand1＜p₁，则执行置换交叉操作；再次产生一个随机数rand2∈(0,1)，与转位交叉概率p₂比较，若rand2＜p₂，则执行转位交叉操作；再产生一组与萤火虫个体DNA序列位置向量维数相同的(0，1)上的随机数，这组随机数中的元素与DNA序列位置向量中的元素一一对应，将所有随机数分别与变异概率p_m比较，若随机数小于p_m，则执行变异操作；

步骤2-6，判断是否达到最大进化代数sg1，如果否，则转至步骤2-7，否则，转至步骤2-8；

步骤2-7，将新一代萤火虫群中所有萤火虫个体的DNA序列位置向量解码，得到所有萤火虫个体的位置向量，转至步骤2-3；

步骤2-8，获得最优萤火虫群；

步骤2-9，计算最优萤火虫群中所有萤火虫位置向量的适应度函数值，选取适应度函数值最大的位置向量为当前最优位置向量；

步骤2-10，更新萤火虫荧光素值，确定邻域，更新萤火虫位置和动态决策范围，进行适应度函数值计算，选取适应度函数值最大的位置向量为当前最优位置向量；

步骤2-11，比较步骤2-10和步骤2-9中最优位置向量的适应度函数值，取适应度函数值大的位置向量为全局最优位置向量；

步骤2-12，判断是否达到最大迭代次数sg2，如果是，则转至步骤2-13，如果否，则转至步骤2-10；

步骤2-13，获得全局最优位置向量，及对应的最优萤火虫个体。