CN107018103A - 一种基于自适应步长猴群优化的小波常模盲均衡方法 - Google Patents
一种基于自适应步长猴群优化的小波常模盲均衡方法 Download PDFInfo
- Publication number
- CN107018103A CN107018103A CN201710222810.5A CN201710222810A CN107018103A CN 107018103 A CN107018103 A CN 107018103A CN 201710222810 A CN201710222810 A CN 201710222810A CN 107018103 A CN107018103 A CN 107018103A
- Authority
- CN
- China
- Prior art keywords
- monkey
- group
- vector
- value
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 102
- 241000282693 Cercopithecidae Species 0.000 title claims abstract description 101
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 44
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 72
- 230000009466 transformation Effects 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims description 17
- 230000009194 climbing Effects 0.000 claims description 13
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000008054 signal transmission Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000007306 turnover Effects 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000011664 nicotinic acid Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03082—Theoretical aspects of adaptive time domain methods
- H04L25/03089—Theory of blind algorithms, recursive or not
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03159—Arrangements for removing intersymbol interference operating in the frequency domain
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03165—Arrangements for removing intersymbol interference using neural networks
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03433—Arrangements for removing intersymbol interference characterised by equaliser structure
- H04L2025/03535—Variable structures
- H04L2025/03541—Switching between domains, e.g. between time and frequency
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03592—Adaptation methods
- H04L2025/03598—Algorithms
- H04L2025/03611—Iterative algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03592—Adaptation methods
- H04L2025/03598—Algorithms
- H04L2025/03681—Control of adaptation
- H04L2025/03687—Control of adaptation of step size
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Theoretical Computer Science (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种基于自适应步长猴群优化的小波常模盲均衡方法,本发明将自适应步长引入猴群优化方法,得到了一种自适应步长猴群优化方法,这个新方法具有很好的全局寻优能力,利用这个新方法得到猴群在搜索空间的最优位置向量,将此向量作为基于小波变换的常模盲均衡方法中的初始权向量,能使盲均衡系统更接近期望的理想系统,从而收敛速度加快,均方误差降低。与现有技术相比,本发明参数少、复杂度低、收敛速度快、稳态误差小,有很好的实用价值。
Description
技术领域:
本发明涉及数据处理技术领域,具体讲是一种基于自适应步长猴群优化的小波常模盲均衡方法。
背景技术:
信号在水声数字系统中高速传输时,会因带宽有限和多径传播等因素的影响产生码间干扰(Inter-Symbol interference,ISI),导致严重失真,若在接收端引入盲均衡技术,可有效消除和减少ISI,提高通信质量。常模盲均衡方法(Constant ModulusAlgorithm,CMA)结构简单、性能稳定,目前被广泛采用,但存在收敛速度慢、稳态误差大等问题;将小波变换(WT)引入了CMA,利用WT可降低信号与噪声的自相关性,WT-CMA虽较传统CMA性能有了一定的提升,但方法中的代价函数仍是多模态的,寻优方法也仍沿用了CMA的梯度思想,所以未能从根本上解决无法获得全局最优解的问题,效果有限。自适应步长猴群优化方法(LMA)作为一种新型仿生群智能优化方法,将猴群爬山过程中的爬、望、跳、翻四个动作模拟设计成了三个搜索过程,并在爬过程中采用自适应步长,用猴群代表一组目标函数的潜在解,用猴群所能达到的最高位置作为目标函数的全局最优解,当猴群寻找到搜索空间的最高位置时,其适应度函数取得最大值。
发明内容:
本发明要解决的技术问题是,提供一种基于自适应步长猴群优化的小波常模盲均衡方法,该方法利用自适应步长猴群优化方法良好的全局搜索能力来获取猴群在搜索空间的最优位置向量,并将这个最优位置向量作为小波常模盲均衡方法中的初始权向量,结合小波变换降低信号自相关性,有效改善了常模盲均衡方法容易陷入局部最小值、收敛速度慢、稳态误差大等问题,具有能获取全局最优解,收敛速度快,稳态误差小的优点。
本发明的技术解决方案是,提供一种基于自适应步长猴群优化的小波常模盲均衡方法,包括基于小波变换的常模盲均衡方法,其中,还包括自适应步长猴群优化方法,具体步骤如下:
步骤①初始化阶段:先设置基于自适应步长猴群优化的小波常模盲均衡方法中的相关参数,接着随机产生初始猴群,确定自适应步长猴群优化方法的适应度函数,然后生成信号y(k),将其作为自适应步长猴群优化方法的输入,最后确定猴群全局最优位置向量的初始值;
步骤②迭代寻优阶段:通过自适应步长猴群优化方法获得猴群的全局最优位置向量,并将其作为基于小波变换的常模盲均衡方法(WT-CMA)的初始权向量;
步骤③信号均衡输出阶段:将上述步骤①中的输入信号y(k)通过基于小波变换的常模盲均衡方法均衡输出。
优选地,本发明所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其中,初始化阶段的具体步骤如下:
步骤a设置基于自适应步长猴群优化的小波常模盲均衡方法中所有相关参数;
步骤b生成信号y(k):信号经信道传输后将发生一定变化。
y(k)=cT(k)a(k)+b(k) (1)
式(1)中,a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列。
步骤c随机产生初始猴群:在n维空间创建一个规模为m的猴群X=[X1,X2,…,Xm],则第i只猴子的位置可以用一个n维的向量Xi=(xi1,xi2,…,xin),i=1,2,…,m来表示,每只人工猴的位置向量与基于小波变换的常模盲均衡方法中的权向量设置为相同形式。猴群的初始化采用的是随机形式,过程如下:
xij=xmin,j+(xmax,j-xmin,j)rand (2)
式(2)中,i=1,2,…,m,j=1,2,…,n,xij为第i只猴子在第j维的实际位置,xmin,j和xmax,j分别表示搜索空间第j维的下界和上界,rand可随机产生一个在区间[0,1]上的实数;
步骤d确定适应度函数:自适应步长猴群优化方法最终是取得其适应度函数的最大值,将自适应步长猴群优化方法的适应度函数f(Xi)设置为基于小波变换的常模盲均衡方法的代价函数J(Xi)的倒数。当f(Xi)取得最大值时,J(Xi)取得最小值,当代价函数取得全局最小值时,盲均衡系统成为期望中的理想系统。
步骤e确定猴群全局最优位置向量的初始值:将y(k)作为混合猴群优化方法的输入信号,根据适应度函数,分别计算每只猴子的初始位置向量对应的适应度函数值,并比较结果,将该群体中的最大适应度函数值定义为全局最优适应度函数值f(X*)的初始值,该值对应的位置向量定义为全局最优位置向量X*的初始值。
优选地,本发明所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其中,迭代寻优阶段的具体步骤如下:
步骤a爬过程,具体如下:
1)在第t次的爬行中,随机生成一个向量ΔXi=(Δxi1,Δxi2,…,Δxin),i=1,2,…,m,分量Δxij以相同的概率0.5取爬步长λ(t)(λ(t)>0)或-λ(t),爬步长将随着爬过程迭代次数的增加而减小
式(3)中,λmin、λmax分别为最小爬步长和最大爬步长,xmin和xmax为搜索空间的上界和下界,t=1,2,…tmax,tmax为爬过程的最大迭代次数。
2)计算伪梯度
式(4)中,i=1,2,…,m,j=1,2,…,n,向量f'i(xi)=(f'i1(xi),f'i2(xi),…f'in(xi))称为适应度函数在点Xi处的伪梯度;
3)设向量Y=(y1,y2,…yn),向量中各分量为
yj=xij+λ(t)·sign(f'ij(xi)) (5)
式(5)中,i=1,2,…,m,j=1,2,…,n,
4)若Y满足约束条件,则更新位置,否则,保持Xi不变;
5)返回爬过程的步骤1),继续进行迭代,直到达到设置的最大爬次数或者前后两次迭代过程中适应度函数值无明显变化,爬过程结束,转入望-跳过程;
步骤b望-跳过程具体步骤如下:
1)每只人工猴经过以上的攀爬,会停下来,在视野范围(xij-γ,xij+γ)(i=1,2,…,m,j=1,2,…,n)内多次瞭望,探索周围是否有更好的位置,每次瞭望得到的位置向量设为Y=(y1,y2,…yn),这个位置向量是在视野范围内通过随机方式产生的,若Y满足约束条件并有f(Y)>f(Xi),则更新位置,否则,保持位置不变,重复望-跳过程的步骤1)直到有满足条件的Y产生或满足一定的运行次数终止;
2)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
3)返回爬过程,进行迭代寻优,达到设定次数结束,进入翻过程;
步骤c翻过程
经过爬过程和望-跳过程的迭代,每只猴子都已到达自己探索区域内的最好位置。翻过程的目的是迫使猴群到新的区域进行探索,避免陷入局部最优。步骤如下:
1)在翻区间[c,d]内随机产生一个实数θ,令
yj=xij+θ(pj-xij) (6)
式(6)中,i=1,2,…,m,j=1,2,…,n,支点每只猴子沿当前位置指向支点的方向(当θ>0时)或反方向(当θ<0时)进行探索,得到Y=(y1,y2,…yn)。若Y满足约束条件并有f(Y)>f(Xi),则人工猴进行翻跳,更新位置,否则,重复翻过程步骤1),直到有满足条件的Y产生或者达到最大设置次数时终止;
2)计算所有猴子当前位置向量对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
检验是否满足结束条件,若满足则迭代寻优阶段结束,即为基于小波变换的常模盲均衡方法中权向量W的初始值,输出;否则,转到爬过程,继续进行迭代寻优。
优选地,本发明所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其中,对信号均衡输出阶段的具体步骤如下:
步骤a对信号进行小波变换
R(k)=y(k)L (7)
式(7)中,L为小波变换矩阵,R(k)是y(k)经过小波变换后的信号序列;
步骤b权向量的迭代
式(8)中,μ为权向量的迭代步长,R*(k)为R(k)的共轭,
分别表示对尺度变换系mI,n(k)、小波变换系数ri,n(k)的平均功率估计,其迭代公式为:
式(9)和(10)中,α(0<α<1)为平滑因子;
步骤c信号的均衡输出
z(k)=WT(k)R(k) (11)
e(k)=RCM-|z(k)|2 (13)
J(k)=E{[RCM-|z(k)|2]2} (14)
式(11)、(12)和(13)中,WT(k)为权向量W(k)的转置,a(k)为独立同分布且均值为零的发射信号,RCM为发射信号的统计模值,e(k)为常模误差函数,z(k)为经过小波变换后均衡器的输出信号。式(14)为WT-CMA的代价函数。
本发明的有益效果是:利用自适应步长猴群优化方法良好的全局搜索能力来获取猴群在搜索空间的最优位置向量,将其作为常模盲均衡方法的初始权向量并结合小波变换降低信号自相关性,有效改善了常模盲均衡方法容易陷入局部最小值、收敛速度慢、稳态误差大等问题,具有收敛速度快,均方误差小,能获得全局最优解的优点。
附图说明:
图1是本发明的原理图;
图2是本发明的流程图;
图3是使用本发明方法后的收敛曲线图;
图4是使用本发明方法后的CMA输出星座图;
图5是使用本发明方法后的WT-CMA输出星座图;
图6是使用本发明方法后的LMA-WT-CMA输出星座图。
具体实施方式:
下面结合附图和具体实施方式对本发明一种基于自适应步长猴群优化的小波常模盲均衡方法作进一步说明:
将小波变换WT引入传统的常模盲均衡方法CMA,对携带噪声的均衡器输入信号实行预处理操作,降低信号间以及信号与噪声间的相关性,能加快方法收敛速度,但WT-CMA陷入局部极值的可能性较大。为解决这一问题,本发明在WT-CMA中引入自适应步长猴群优化方法(LMA),利用其良好的全局寻优特性来寻找WT-CMA的初始权向量(该向量即为猴群的全局最优位置向量)。本发明提出的这种将自适应步长猴群优化方法(LMA)和小波变换理论(WT)与常模盲均衡方法(CMA)有机融合,称为基于自适应步长猴群优化的小波常模盲均衡方法(LMA-WT-CMA)。
如图1所示,图中a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列,R(k)是y(k)经过小波变换后的信号序列,W(k)为均衡器权向量,Ψ(·)为无记忆的非线性函数,表示无记忆非线性估计器,z(k)是均衡器的输出信号,e(k)为常模误差函数。
如图2所示,本发明一种基于自适应步长猴群优化的小波常模盲均衡方法的具体流程如下:
一、初始化阶段
具体步骤如下:
步骤a设置基于自适应步长猴群优化的小波常模盲均衡方法中所有相关参数;
步骤b生成信号y(k):信号经信道传输后将发生一定变化。
y(k)=cT(k)a(k)+b(k) (1)
式(1)中,a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列。
步骤c随机产生初始猴群:在n维空间创建一个规模为m的猴群X=[X1,X2,…,Xm],则第i只猴子的位置可以用一个n维的向量Xi=(xi1,xi2,…,xin),i=1,2,…,m来表示,每只人工猴的位置向量与基于小波变换的常模盲均衡方法中的权向量设置为相同形式。猴群的初始化采用的是随机形式,过程如下:
xij=xmin,j+(xmax,j-xmin,j)rand (2)
式(2)中,i=1,2,…,m,j=1,2,…,n,xij为第i只猴子在第j维的实际位置,xmin,j和xmax,j分别表示搜索空间第j维的下界和上界,rand可随机产生一个在区间[0,1]上的实数;
步骤d确定适应度函数:自适应步长猴群优化方法最终是取得其适应度函数的最大值,将自适应步长猴群优化方法的适应度函数f(Xi)设置为基于小波变换的常模盲均衡方法的代价函数J(Xi)的倒数。当f(Xi)取得最大值时,J(Xi)取得最小值,当代价函数取得全局最小值时,盲均衡系统成为期望中的理想系统。
步骤e确定猴群全局最优位置向量的初始值:将y(k)作为混合猴群优化方法的输入信号,根据适应度函数,分别计算每只猴子的初始位置向量对应的适应度函数值,并比较结果,将该群体中的最大适应度函数值定义为全局最优适应度函数值f(X*)的初始值,该值对应的位置向量定义为全局最优位置向量X*的初始值。
二、迭代寻优阶段
找出使适应度函数值最大的位置向量,将此向量作为基于小波变换的常模盲均衡方法的初始权向量。当代价函数取得全局最小值时,盲均衡系统成为期望中的理想系统。具体步骤如下:
步骤a爬过程。自适应步长猴群优化方法的爬过程主要是利用伪梯度思想,计算适应度函数在当前位置相邻两点的函数值,继而进行比较移动的过程。具体如下:
1)在第t次的爬行中,随机生成一个向量ΔXi=(Δxi1,Δxi2,…,Δxin),i=1,2,…,m,分量Δxij以相同的概率0.5取爬步长λ(t)(λ(t)>0)或-λ(t),爬步长将随着爬过程迭代次数的增加而减小
式(3)中,λmin、λmax分别为最小爬步长和最大爬步长,xmin和xmax为搜索空间的上界和下界,t=1,2,…tmax,tmax为爬过程的最大迭代次数。
2)计算伪梯度
式(4)中,i=1,2,…,m,j=1,2,…,n,向量f'i(xi)=(f'i1(xi),f'i2(xi),…f'in(xi))称为适应度函数在点Xi处的伪梯度。
3)设向量Y=(y1,y2,…yn),向量中各分量为
yj=xij+λ(t)·sign(f'ij(xi)) (5)
式(5)中,i=1,2,…,m,j=1,2,…,n,
4)若Y满足约束条件,则更新位置,即Xi←Y;否则,保持Xi不变。
5)返回爬过程的步骤1),继续进行迭代重新进行过程,直到达到设置的最大爬次数或者前后两次迭代过程中适应度函数值无明显变化,爬过程结束,转入望-跳过程;
步骤b望-跳过程。每只猴子经过爬过程后,在视野长度为γ的范围内,向四周进行多次眺望,寻找是否存在比当前更好的位置,如果有,则更新到更好的位置上去。具体步骤如下:
1)在视野范围(xij-γ,xij+γ)(i=1,2,…,m,j=1,2,…,n)内,在每只人工猴当前位置周围随机进行一定次数眺望,设某次得到的位置向量为Y=(y1,y2,…yn),若Y满足约束条件并有f(Y)>f(Xi),则更新位置,即Xi←Y;否则,保持位置不变,重复望-跳过程的步骤1)直到有满足条件的Y产生或满足一定的运行次数终止。
2)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*。
3)返回爬过程,重复至达到设定次数结束,进入翻过程。
步骤c翻过程
经过爬过程和望-跳过程的迭代,每只猴子都已到达自己探索区域内的最好位置。翻过程的目的是迫使猴群到新的区域进行探索,避免陷入局部最优。步骤如下:
1)在翻区间[c,d]内随机产生一个实数θ,令
yj=xij+θ(pj-xij) (6)
式(6)中,i=1,2,…,m,j=1,2,…,n,支点每只猴子沿当前位置指向支点的方向(当θ>0时)或反方向(当θ<0时)进行探索,得到Y=(y1,y2,…yn)。若Y满足约束条件并有f(Y)>f(Xi),则人工猴进行翻跳,更新位置,否则,重复翻过程步骤1),直到有满足条件的Y产生或者达到最大设置次数时终止;
2)计算所有猴子当前位置对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*。
检验是否满足结束条件,若满足则迭代寻优阶段结束,即为基于小波变换的常模盲均衡方法中权向量W的初始值,输出;否则,转到爬过程,继续进行迭代寻优。
三、信号均衡输出阶段
具体步骤如下:
步骤a对信号进行小波变换
R(k)=y(k)L (7)
式(7)中,L为小波变换矩阵,R(k)是y(k)经过小波变换后的信号序列;
步骤b权向量的迭代
式(8)中,μ为权向量的迭代步长,R*(k)为R(k)的共轭。
分别表示对尺度变换系mI,n(k)、小波变换系数ri,n(k)的平均功率估计,其迭代公式为:
式(9)和(10)中,α(0<α<1)为平滑因子。
步骤c信号的均衡输出
z(k)=WT(k)R(k) (11)
e(k)=RCM-|z(k)|2 (13)
J(k)=E{[RCM-|z(k)|2]2} (14)
式(11)、(12)和(13)中,WT(k)为权向量W(k)的转置,a(k)为独立同分布且均值为零的发射信号,RCM为发射信号的统计模值,e(k)为常模误差函数,z(k)为最终的输出信号。式(14)为WT-CMA的代价函数。
四、仿真验证
混合相位水声信道:h=[0.3132,-0.1040,0.8908,0.3134],发射信号为16QAM,均衡器抽头个数为16,信噪比为25dB。CMA中,第8个抽头系数设置为1,其余为0,步长μCMA=0.00054;WT-CMA中,第8个抽头系数设置为1,其余为0,步长μWT-CMA=0.00023;LMA-WT-CMA中,步长μLMA-WT-CMA=0.000018。输入信号采用DB4小波进行分解,分解层次2层,功率初始值设置为4,遗忘因子β=0.999;2000次蒙特卡诺仿真结果如图3、图4、图5和图6所示,从以上图中可看出:LMA-WT-CMA在收敛速度上比CMA快了大约4000步、比WT-CMA快了大约2000步。在稳态误差上比CMA减小了12dB左右,比WT-CMA减小了7dB左右,所以LMA-WT-CMA的输出星座图最清晰,紧凑,均衡效果最好。
以上所述的实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (4)
1.一种基于自适应步长猴群优化的小波常模盲均衡方法,包括基于小波变换的常模盲均衡方法(WT-CMA),其特征在于:还包括自适应步长猴群优化方法,具体步骤如下:
步骤①初始化阶段:先设置基于自适应步长猴群优化的小波常模盲均衡方法中的相关参数,接着随机产生初始猴群,确定自适应步长猴群优化方法的适应度函数,然后生成信号y(k),将其作为自适应步长猴群优化方法的输入信号,最后确定猴群全局最优位置向量的初始值;
步骤②迭代寻优阶段:通过自适应步长猴群优化方法获得猴群的全局最优位置向量,并将其作为基于小波变换的常模盲均衡方法(WT-CMA)的初始权向量;
步骤③信号均衡输出阶段:将上述步骤①中的信号y(k)通过基于小波变换的常模盲均衡方法均衡输出。
2.根据权利要求1所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其特征在于:所述初始化阶段的具体步骤如下:
步骤a设置基于自适应步长猴群优化的小波常模盲均衡方法中所有相关参数;
步骤b生成信号y(k):信号经信道传输后将发生一定变化,
y(k)=cT(k)a(k)+b(k) (1)
式(1)中,a(k)为平稳独立同分布且具有零均值的发射信号序列,c(k)为信道的脉冲响应向量,b(k)为加性高斯白噪声,y(k)为a(k)经信道传输后的信号序列;
步骤c随机产生初始猴群:在n维空间创建一个规模为m的猴群X=[X1,X2,…,Xm],则第i只猴子的位置可以用一个n维的向量Xi=(xi1,xi2,…,xin),i=1,2,…,m来表示,每只人工猴的位置向量与基于小波变换的常模盲均衡方法中的权向量设置为相同形式,猴群的初始化采用的是随机形式,过程如下:
xij=xmin,j+(xmax,j-xmin,j)rand (2)
式(2)中,i=1,2,…,m,j=1,2,…,n,xij为第i只猴子在第j维的实际位置,xmin,j和xmax,j分别表示搜索空间第j维的下界和上界,rand可随机产生一个在区间[0,1]上的实数;
步骤d确定适应度函数:自适应步长猴群优化方法最终是取得其适应度函数的最大值,将自适应步长猴群优化方法的适应度函数f(Xi)设置为基于小波变换的常模盲均衡方法的代价函数J(Xi)的倒数,当f(Xi)取得最大值时,J(Xi)取得最小值,当代价函数取得全局最小值时,盲均衡系统成为期望中的理想系统;
步骤e确定猴群全局最优位置向量的初始值:将y(k)作为混合猴群优化方法的输入信号,根据适应度函数,分别计算每只猴子的初始位置向量对应的适应度函数值,并比较结果,将该群体中的最大适应度函数值定义为全局最优适应度函数值f(X*)的初始值,该值对应的位置向量定义为全局最优位置向量X*的初始值。
3.根据权利要求1所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其特征在于:所述迭代寻优阶段的具体步骤如下:
步骤a爬过程,具体如下:
1)在第t次的爬行中,随机生成一个向量ΔXi=(Δxi1,Δxi2,…,Δxin),i=1,2,…,m,分量Δxij以相同的概率0.5取爬步长λ(t)(λ(t)>0)或-λ(t),爬步长将随着爬过程迭代次数的增加而减小,
式(3)中,λmin、λmax分别为最小爬步长和最大爬步长,xmin和xmax为搜索空间的上界和下界,t=1,2,…tmax,tmax为爬过程的最大迭代次数;
2)计算伪梯度
式(4)中,i=1,2,…,m,j=1,2,…,n,向量f'i(xi)=(f'i1(xi),f'i2(xi),…f'in(xi))称为适应度函数在点Xi处的伪梯度;
3)设向量Y=(y1,y2,…yn),向量中各分量为
yj=xij+λ(t)·sign(f'ij(xi)) (5)
式(5)中,i=1,2,…,m,j=1,2,…,n
4)若Y满足约束条件,则更新位置,否则,保持Xi不变;
5)返回爬过程的步骤1),继续进行迭代,直到达到设置的最大爬次数或者前后两次迭代过程中适应度函数值无明显变化,爬过程结束,转入望-跳过程;
步骤b望-跳过程具体步骤如下:
1)每只人工猴经过以上的攀爬,会停下来,在视野范围(xij-γ,xij+γ)(i=1,2,…,m,j=1,2,…,n)内多次瞭望,探索周围是否有更好的位置,每次瞭望得到的位置向量设为Y=(y1,y2,…yn),这个位置向量是在视野范围内通过随机方式产生的,若Y满足约束条件并有f(Y)>f(Xi),则更新位置,否则,保持位置不变,重复望-跳过程的步骤1)直到有满足条件的Y产生或满足一定的运行次数终止;
2)计算所有猴子当前位置向量对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
3)返回爬过程,重复至达到设定次数结束,进入翻过程;
步骤c翻过程
经过爬过程和望-跳过程的迭代,每只猴子都已到达自己探索区域内的最好位置。翻过程的目的是迫使猴群到新的区域进行探索,避免陷入局部最优。步骤如下:
1)在翻区间[c,d]内随机产生一个实数θ,令
yj=xij+θ(pj-xij) (6)
式(6)中,i=1,2,…,m,j=1,2,…,n,支点每只猴子沿当前位置指向支点的方向(当θ>0时)或反方向(当θ<0时)进行探索,得到Y=(y1,y2,…yn)。若Y满足约束条件并有f(Y)>f(Xi),则人工猴进行翻跳,更新位置,否则,重复翻过程步骤1),直到有满足条件的Y产生或者达到最大设置次数时终止;
2)计算所有猴子当前位置向量对应的适应度函数值,并与f(X*)比较,用最优值更新f(X*),用该值对应的位置向量更新X*;
检验是否满足结束条件,若满足则迭代寻优阶段结束,即为基于小波变换的常模盲均衡方法中权向量W的初始值,输出;否则,转到爬过程,继续进行迭代寻优。
4.根据权利要求1所述的一种基于自适应步长猴群优化的小波常模盲均衡方法,其特征在于:所述信号均衡输出阶段具体步骤如下:
步骤a对信号y(k)进行小波变换
R(k)=y(k)L (7)
式(7)中,L为小波变换矩阵,R(k)是y(k)经过小波变换后的信号序列;
步骤b权向量的迭代
式(8)中,μ为权向量的迭代步长,R*(k)为R(k)的共轭, 分别表示对尺度变换系mI,n(k)、小波变换系数ri,n(k)的平均功率估计,其迭代公式为:
式(9)和(10)中,α(0<α<1)为平滑因子;
步骤c信号的均衡输出
z(k)=WT(k)R(k) (11)
e(k)=RCM-|z(k)|2 (13)
J(k)=E{[RCM-|z(k)|2]2} (14)
式(11)、(12)和(13)中,WT(k)为权向量W(k)的转置,a(k)为独立同分布且均值为零的发射信号,RCM为发射信号的统计模值,e(k)为常模误差函数,z(k)为最终的输出信号。式(14)为WT-CMA的代价函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710222810.5A CN107018103B (zh) | 2017-04-07 | 2017-04-07 | 一种基于自适应步长猴群优化的小波常模盲均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710222810.5A CN107018103B (zh) | 2017-04-07 | 2017-04-07 | 一种基于自适应步长猴群优化的小波常模盲均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107018103A true CN107018103A (zh) | 2017-08-04 |
CN107018103B CN107018103B (zh) | 2020-02-14 |
Family
ID=59445280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710222810.5A Expired - Fee Related CN107018103B (zh) | 2017-04-07 | 2017-04-07 | 一种基于自适应步长猴群优化的小波常模盲均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107018103B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107765259A (zh) * | 2017-09-18 | 2018-03-06 | 国家电网公司 | 一种基于提升小波改进阈值的输电线路激光测距信号去噪算法 |
CN110533905A (zh) * | 2019-07-03 | 2019-12-03 | 广东工业大学 | 一种交通流预测的方法、系统及设备 |
CN111064683A (zh) * | 2019-12-23 | 2020-04-24 | 西安电子科技大学 | 一种改进的恒模均衡算法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1574812A (zh) * | 2003-06-09 | 2005-02-02 | 三星电子株式会社 | 多载波接收装置的盲均衡器及其方法 |
CN101232358A (zh) * | 2007-01-23 | 2008-07-30 | 复旦大学 | 一种在多输入多输出系统中信道盲均衡方法和系统 |
CN101304394A (zh) * | 2008-06-12 | 2008-11-12 | 上海交通大学 | 针对空时编码正交频分复用系统的免疫动态信道均衡方法 |
CN101478509A (zh) * | 2009-01-20 | 2009-07-08 | 南京信息工程大学 | 正交小波变换和时间分集技术相融合的盲均衡方法 |
CN101547173A (zh) * | 2009-05-05 | 2009-09-30 | 南京信息工程大学 | 基于平衡正交多小波变换的盲均衡方法 |
CN101902417A (zh) * | 2010-06-30 | 2010-12-01 | 南京信息工程大学 | 基于蚁群优化的正交小波变换超指数迭代盲均衡方法 |
CN102123115A (zh) * | 2011-04-15 | 2011-07-13 | 南京信息工程大学 | 基于粒子群优化的正交小波盲均衡方法 |
EP2487813B1 (en) * | 2011-02-01 | 2013-04-17 | Alcatel Lucent | Polarization de-multiplex for multilevel signals |
CN105007246A (zh) * | 2015-07-29 | 2015-10-28 | 南京信息工程大学 | 一种模因方法优化的多模盲均衡方法 |
-
2017
- 2017-04-07 CN CN201710222810.5A patent/CN107018103B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1574812A (zh) * | 2003-06-09 | 2005-02-02 | 三星电子株式会社 | 多载波接收装置的盲均衡器及其方法 |
CN101232358A (zh) * | 2007-01-23 | 2008-07-30 | 复旦大学 | 一种在多输入多输出系统中信道盲均衡方法和系统 |
CN101304394A (zh) * | 2008-06-12 | 2008-11-12 | 上海交通大学 | 针对空时编码正交频分复用系统的免疫动态信道均衡方法 |
CN101478509A (zh) * | 2009-01-20 | 2009-07-08 | 南京信息工程大学 | 正交小波变换和时间分集技术相融合的盲均衡方法 |
CN101547173A (zh) * | 2009-05-05 | 2009-09-30 | 南京信息工程大学 | 基于平衡正交多小波变换的盲均衡方法 |
CN101902417A (zh) * | 2010-06-30 | 2010-12-01 | 南京信息工程大学 | 基于蚁群优化的正交小波变换超指数迭代盲均衡方法 |
EP2487813B1 (en) * | 2011-02-01 | 2013-04-17 | Alcatel Lucent | Polarization de-multiplex for multilevel signals |
CN102123115A (zh) * | 2011-04-15 | 2011-07-13 | 南京信息工程大学 | 基于粒子群优化的正交小波盲均衡方法 |
CN105007246A (zh) * | 2015-07-29 | 2015-10-28 | 南京信息工程大学 | 一种模因方法优化的多模盲均衡方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107765259A (zh) * | 2017-09-18 | 2018-03-06 | 国家电网公司 | 一种基于提升小波改进阈值的输电线路激光测距信号去噪算法 |
CN110533905A (zh) * | 2019-07-03 | 2019-12-03 | 广东工业大学 | 一种交通流预测的方法、系统及设备 |
CN110533905B (zh) * | 2019-07-03 | 2022-02-15 | 广东工业大学 | 一种交通流预测的方法、系统及设备 |
CN111064683A (zh) * | 2019-12-23 | 2020-04-24 | 西安电子科技大学 | 一种改进的恒模均衡算法 |
Also Published As
Publication number | Publication date |
---|---|
CN107018103B (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111464465B (zh) | 一种基于集成神经网络模型的信道估计方法 | |
CN102123115B (zh) | 基于粒子群优化的正交小波盲均衡方法 | |
CN109756432B (zh) | Ofdm信道估计方法和装置 | |
CN107124379A (zh) | 一种基于改进狼群优化的正交小波常模盲均衡方法 | |
CN107018103A (zh) | 一种基于自适应步长猴群优化的小波常模盲均衡方法 | |
CN108540419A (zh) | 一种基于深度学习的抗子载波间干扰的ofdm检测方法 | |
CN105070296A (zh) | 一种活性因子集员成比例子带自适应回声消除方法 | |
CN110138459A (zh) | 基于基追踪去噪的稀疏水声正交频分复用信道估计方法及装置 | |
Tian et al. | Group sparse underwater acoustic channel estimation with impulsive noise: Simulation results based on Arctic ice cracking noise | |
CN103312422B (zh) | 一种基于人工鱼群算法的信号盲检测方法 | |
CN107070825A (zh) | 一种基于模拟退火狼群优化的小波加权多模盲均衡方法 | |
Cheng et al. | Channel estimation and equalization based on deep BLSTM for FBMC-OQAM systems | |
CN106953820A (zh) | 基于双Sigmoid复数连续神经网络的信号盲检测方法 | |
CN115499276B (zh) | 智能反射面辅助通信系统的信道估计方法、系统及产品 | |
CN115065578B (zh) | 一种基于改进自适应阈值的dft信道估计方法 | |
CN107231194B (zh) | 室内可见光通信系统中基于收敛状态的变步长均衡方案 | |
CN103841067B (zh) | 一种水声信道通信信号的均衡方法 | |
CN107018104B (zh) | 一种基于混合猴群优化的小波加权多模盲均衡方法 | |
CN102118332A (zh) | 基于免疫克隆粒子群优化的正交小波盲均衡方法 | |
CN101924718B (zh) | 模糊神经网络控制的混合小波神经网络盲均衡方法 | |
CN114611550A (zh) | 基于复数卷积模块的多特征自动调制识别方法 | |
CN103117969A (zh) | 一种分数低阶统计量模值变换小波频域多模盲均衡方法 | |
CN102185808B (zh) | 一种快速收敛的免疫克隆正交小波变换盲均衡方法 | |
CN103973617B (zh) | 多果蝇群协同优化小波常模盲均衡方法 | |
CN103346987A (zh) | 免疫人工鱼群优化频域自适应小波多模盲均衡方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210608 Address after: 230000 Laoshan South Road, Chaohu City, Hefei City, Anhui Province Patentee after: HEFEI TECHNOLOGY College Address before: 232001 No.1, Linchang Road, tianjia'an District, Huainan City, Anhui Province Patentee before: HUAINAN VOCATIONAL TECHNICAL College |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200214 |
|
CF01 | Termination of patent right due to non-payment of annual fee |