CN101478349B - 基于正交小波包变换的双模式变步长盲均衡方法 - Google Patents
基于正交小波包变换的双模式变步长盲均衡方法 Download PDFInfo
- Publication number
- CN101478349B CN101478349B CN 200910028459 CN200910028459A CN101478349B CN 101478349 B CN101478349 B CN 101478349B CN 200910028459 CN200910028459 CN 200910028459 CN 200910028459 A CN200910028459 A CN 200910028459A CN 101478349 B CN101478349 B CN 101478349B
- Authority
- CN
- China
- Prior art keywords
- wavelet packet
- wpt
- orthogonal wavelet
- packet transform
- centerdot
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明公布了一种基于正交小波包变换的双模式变步长盲均衡方法(WPT-VDMA),本发明利用符号判决实现基于正交小波包变换的常数模盲均衡算法(WPT-CMA)与基于正交小波包变换的判决导引盲均衡算法(WPT-DD)间的切换,以减小均方误差,并以输出信号功率为自变量的变步长来加快WPT-CMA的收敛速度,进一步提高WPT-VDMA算法的收敛性能。该发明方法收敛速度快、均方误差小。
Description
技术领域
本发明涉及一种双模式变步长盲均衡算法,尤其涉及一种基于正交小波包变换的双模式变步长盲均衡方法。
背景技术
在水声数字通信中,多径效应和有限带宽会产生严重的码间干扰(inter-symbol interference,ISI),降低了水下数据传输速率和可靠性。因此,必须采用均衡技术来提高通信质量。目前,不需要发射周期性的训练序列的盲均衡技术是水声通信领域的研究热点。而收敛速度、均方误差是影响均衡器性能优劣的主要因素。因此,目前国内外学者在提高收敛速度和减小均方误差方面提出了许多方法。文献[1](CHEN S.Low complexity concurrent constant modulusalgorithm and soft decision-directed scheme for blind equalization[J].IEE Proc.Vis.Image Signal Processing,2003,137(5):312-320.)与[2](LITWIN L R,ZOLTOWSKI M D,ENDERS T J.Blended CMA:smooth,adaptive transfer fromCMA to DD-LMS[A]IEEE Wireless Communications and NetworkingConference[C].1999.797-800.)将判决导引算法(DD,Decision Directed)引入到盲均衡算法中,有效地减小了均方误差,降低了通信的误码率,但该算法没有改变均衡器输入信号的自相关性,这是影响算法收敛性的主要因素之一。文献[2](Cooklev T An Efficient Architecture for Orthogonal Wavelet Transforms[J].IEEESignal Processing Letters(S1070-9980),2006,13(2):77-79.)、[3](HUANG K,LV R.Adaptive equalization algorithm based on wavelet packet transform[J].ActaElectronic Sinica,2003,31(8):1205-1208.)、[4](REN K CH,TU Y Q.A method ofimproving frequency resolution of modified covariance spectrum estimation and itsabilities to adapt noise using multistage wavelet de-composition[J].Journal ofElectronic Measurement and Instrument,2006,20(4):15-21.)、[5](Attallah S.Thewavelet transform-domains LMS adaptive filter with partial subband coefficientupdating[J].IEEE Trans Cricuits and Systems,2006,53(1):8-12.)与[6](S Long P.Dimensional Finite wavelet Filters[J].Journal of Computation Mathernatics,2003,5:595-602.)利用正交小波变换对均衡器输入信号进行变换,再对信号进行能量归一化处理,降低了信号的自相关性,不仅提高了收敛速度,而且算法的复杂度也不高。然而,由于小波变换只对尺度空间作了分解,当信号的高频部分信息较丰富时,由于这种划分把高频都分到了一个频带,因此细节难以分辨。与小波变换相比,小波包变换对信号的高、低频分量均进行分解,能有效分辨信号细节,去相关能力更强,因而效果更好。
发明内容
本发明要解决的技术问题是针对现有技术存在的缺陷提出一种基于正交小波包变换的双模式变步长盲均衡方法。
本发明基于正交小波包变换的双模式变步长盲均衡方法,其特征在于包括如下步骤:
a.)将发射信号a(n)经过脉冲响应信道c(n)得到信道输出向量x(n),其中n为正整数表示时间序列,下同;
b.)将信道噪声w(n)与步骤a所述的信道输出向量x(n)求和得到均衡器的输入向量:y(n)=x(n)+w(n);
c.)将步骤b所述的均衡器的输入向量y(n)先经过正交小波包变换得到正交小波包变换器WPT的输出向量:R(n)=y(n)Q,其中Q为正交小波包变换矩阵;
d.)将步骤c所述的正交小波包变换器WPT的输出向量R(n)经过功率归一化后与常数模误差或判决导引误差更新均衡器权向量f(n)后得到均衡器输出信号z(n)=fH(n)R(n),其中上标H表示对均衡器权向量f(n)取共轭转置;
其中步骤d所述的常数模误差和判决导引误差的选择通过以下方法判断:
当 则将步骤e所述的判决输出信号与步骤d所述的均衡器输出信号z(n)所构成的误差信号更新基于正交小波包变换的判决导引盲均衡算法WPT-DD权向量:其中是功率归一化平均能量构成的对角矩阵,n+1为当前时刻n的后一时刻,下同,R*(n)是R(n)的共轭矩阵,下同,μ1为基于正交小波包变换的判决导引盲均衡算法WPT-DD的迭代步长;
当 则将步骤d所述的均衡器输出信号z(n)和发射信号a(n)的模R2所构成误差信号|z(n)|2-R2更新基于正交小波包变换的常数模盲均衡算法WPT-CMA权向量:
本发明利用符号判决实现基于正交小波包变换的常数模盲均衡算法(WPT-CMA)与基于正交小波包变换的判决导引盲均衡算法(WPT-DD)间的切换,以减小均方误差;并以输出信号功率为自变量的变步长来进一步提高算法的收敛性能,具有收敛速度快、均方误差小的特点。
附图说明
图1:基于正交小波变包变换的盲均衡器结构图;
图2:本发明结构图;
图3:本发明实施例1仿真图:
(a)均方误差曲线图,(b)CMA输出结果图,(c)WPT-CMA输出结果图,
(d)WPT-DMA输出结果图,(e)WPT-VDMA输出结果图;
图4:本发明实施例2仿真图:
(a)均方误差曲线图,(b)CMA输出结果图,(c)WPT-CMA输出结果图,
(d)WPT-DMA输出结果图,(e)WPT-VDMA输出结果图。
具体实施方式
如图1所示。基于正交小波变包变换的盲均衡器结构,图中,n∈Z+(表示时间序列),a(n)是零均值独立同分布发射信号;c(n)是信道脉冲响应;w(n)是加性高斯白噪声;y(n)为均衡器的输入向量;R(n)是y(n)经过正交小波包变换后的信号;f(n)是均衡器权向量且长度为L;ψ(·)是无记忆非线性函数,用来表示无记忆非线性估计器;z(n)是均衡器输出信号;是判决装置对z(n)的判决输出信号。
根据小波包理论,有限冲击响应(FIR)均衡器权向量f(n)可用一簇正交小波包基函数来表示,即
式中,rj,k,m(n)=<f(n),ψj,k,m(n)>(表示均衡器f(n)与小波包基ψj,k,m(n)取内积),ψj,k,m(n)=2-j/2ψm(2-jn-k)(表示小波包基函数),n=0,1,…,2N-1(表示时间序列),J为小波包分解的最大尺度,j=1,2,…J(表示分解层数),KJ=2N/2J表示在尺度J下小波函数的最大平移,k表示小波函数的平移,m表示小波函数的尺度,rj,k,m(n)为均衡器的权系数。此时,均衡器的输出为
式中, l表示长为2N的均衡器权向量的第l个抽头,0≤l≤2N-1。该式表明,输入向量y(n)需与每一个尺度上的小波包基函数ψj,k,m(n)作卷积,即相当于对输入y(n)作离散正交小波包变换,Wj,k,m(n)为相应的变换系数。
设N=2J,均衡器的长度为2N,经过推导,可以得到i(1≤i≤J-log2 2N+1)级小波包分解对应的2J-i+1×2J-i+1矩阵
所以,若令Q为J级小波包分解的2J×2J矩阵,则Q可表示为
式中,Hj和Gj分别为由小波滤波器系数h(n)和尺度滤波器系数g(n)所构成的矩阵。
设 其中表示第j层分解中的第kj组的第m个信号,表示与信号相对应的权向量抽头系数。j∈(0,J)表示正交小波包变换的分解层数,J表示最大分解层数;kj∈(0,kJ)表示j层分解中第kj组,最大组数是kJ;m∈(1,2N)表示每一组中的第m个信号,2N表示均衡器的长度。
根据最小均方误差准则,可以得到基于正交小波包变换的常数模盲均衡算法(orthogonal Wavelet Packet Transform based CMA,WPT-CMA)为
R(n)=y(n)Q,
z(n)=fH(n)R(n),
e(n)=|z(n)|2-R2,
f(n+1)=f(n)+μiR-1(n)e(n)R*(n)z(n),
式中,μ是迭代步长,R2=E{|a(k)|4}/E{|a(k)|2}是一个依赖于信源序列高阶统计量的实常数, 且 式中,diag[]表示对角矩阵,0<β<1为遗忘因子,σj,k,m 2(n)表示的是对信号rj,k,m(n)的功率估计值。信号经过正交小波包变换后,又对其进行了能量归一化处理,使得收敛速度得到进一步的提高。
如图2所示。本发明基于正交小波包变换的双模式变步长盲均衡方法,其特征在于包括如下步骤:
a.)将发射信号a(n)经过脉冲响应信道c(n)得到信道输出向量x(n),其中n为正整数表示时间序列,下同;
b.)将信道噪声w(n)与步骤a所述的信道输出向量x(n)求和得到均衡器的输入向量:y(n)=x(n)+w(n);
c.)将步骤b所述的均衡器的输入向量y(n)先经过正交小波包变换得到正交小波包变换器WPT的输出向量:R(n)=y(n)Q,其中Q为正交小波包变换矩阵;
d.)将步骤c所述的正交小波包变换器WPT的输出向量R(n)经过功率归一化后与常数模误差或判决导引误差更新均衡器权向量f(n)后得到均衡器输出信号z(n)=fH(n)R(n),其中上标H表示对均衡器权向量f(n)取共轭转置;
e.)将步骤d所述的均衡器输出信号z(n)经过判决装置得到判决输出信号
其中步骤d所述的常数模误差和判决导引误差的选择通过以下方法判断:
当 则将步骤e所述的判决输出信号与步骤d所述的均衡器输出信号z(n)所构成的误差信号更新基于正交小波包变换的判决导引盲均衡算法WPT-DD权向量:是功率归一化平均能量构成的对角矩阵,n+1为当前时刻n的后一时刻,下同,R*(n)是R(n)的共轭矩阵,下同,μ1为基于正交小波包变换的判决导引盲均衡算法WPT-DD的迭代步长;
当 则将步骤d所述的均衡器输出信号z(n)和发射信号a(n)的模R2所构成误差信号|z(n)|2-R2更新基于正交小波包变换的常数模盲均衡算法WPT-CMA权向量:
其中,是发射信号a(n)的模,E[]为数学期望,sign[]为符号函数,μ2为基于正交小波包变换的常数模盲均衡算法WPT-CMA的迭代步长。
本发明在迭代初始阶段,由于算法未收敛,sign[|z(n)|2-R2]与不等的概率较大,算法多工作在WPT-CMA模式,同时利用WPT-CMA算法稳定性,保证了算法收敛;而当算法收敛后,sign[|z(n)|2-R2]与相等的概率增大,则算法多以WPT-DD模式迭代,从而使得算法在收敛后具有较小的均方误差。
对于盲均衡算法中的步长,如果采用大步长,会产生较大的均方误差,但算法收敛速度和跟踪速度快;反之,采用小步长,每次调整权系数的幅度就小,算法收敛速度和跟踪速度慢,但当均衡器权向量接近最优值时,权向量将在最优值附近一个较小的范围内来回抖动,因而均方误差较小。而盲均衡算法收敛的过程就是均衡器输出信号功率逐渐增加,向发送信号功率逐渐逼近的过程。基于此,本发明提出一种以均衡器输出信号功率为自变量的变步长表达式,即
μ(n)=η[1-exp(-α|r(n)|)],
r(n)=E{|z(n)|2},
式中,r(n)表示均衡器输出信号的功率估值,α、η是参数。对μ(n)关于r(n)求导即
dμ(n)/dr(n)=α·ηexp(-α|r(n)|),
因为α、η是大于零的,所以dμ(n)/dr(n)>0,即步长μ(n)随r(n)的增大是单调递增的。
当算法工作在WPT-CMA模式时,采用以均衡器输出信号功率为自变量的变步长,此时算法未收敛,采用逐渐增大步长的办法能够有效的加快收敛速度;当算法工作在WPT-DD模式时,此时算法进入初始收敛阶段,采用较小的固定步长,保证收敛后的均方误差较小,从而使得该方法收敛后具有较快的收敛速度和较小的均方误差。
为了验证基于正交小波包变换的变步长双模式盲均衡方法(WPT-VDMA)的有效性,用水声信道进行仿真实验,并与CMA、WPT-CMA、WPT-DMA(基于正交小波包变换的双模式盲均衡算法)进行比较。
实施例1:复水声信道
如图3所示,信道为c=[e-0.7j 0 0 0.3e-1.8j];发射信号为4QAM;均衡器权长为16,信噪比为20dB;其它参数设置如表1,1500次蒙特卡罗仿真结果,如图3所示。
表1 仿真参数值
实施例2:最小相位水声信道
如图4所示,信道为c=[0.9656 0.0906 0.0578 0.2368];发射信号为32QAM;均衡器权长为16,信噪比为20dB;其它参数设置如表2,3000次蒙特卡罗仿真结果,如图4所示。
表2 仿真参数值
Claims (1)
1.一种基于正交小波包变换的双模式变步长盲均衡方法,其特征在于包括如下步骤:
a.)将发射信号a(n)经过脉冲响应信道c(n)得到信道输出向量x(n),其中n为正整数表示时间序列,下同;
b.)将信道噪声w(n)与步骤a所述的信道输出向量x(n)求和得到均衡器的输入向量:y(n)=x(n)+w(n);
c.)将步骤b所述的均衡器的输入向量y(n)先经过正交小波包变换得到正交小波包变换器WPT的输出向量:R(n)=y(n)Q,其中Q为正交小波包变换矩阵;
d.)将步骤c所述的正交小波包变换器WPT的输出向量R(n)经过功率归一化后与常数模误差或判决导引误差更新均衡器权向量f(n)后得到均衡器输出信号z(n)=fH(n)R(n),其中上标H表示对均衡器权向量f(n)取共轭转置;
其中步骤d所述的常数模误差和判决导引误差的选择通过以下方法判断:
当则将步骤e所述的判决输出信号 与步骤d所述的均衡器输出信号z(n)所构成的误差信号 更新基于正交小波包变换的判决导引盲均衡算法WPT-DD权向量: 其中 是功率归一化平均能量构成的对角矩阵,下同,n+1为当前时刻n的后一时刻,下同,R*(n)是R(n)的共轭矩阵,下同,μ1为基于正交小波包变换的判决导引盲均衡算法WPT-DD的迭代步长;
当则将步骤d所述的均衡器输出信号z(n)和发射信号a(n)的模R2所构成误差信号|z(n)|2-R2更新基于正交小波包变换的常数模盲均衡算法WPT-CMA权向量:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910028459 CN101478349B (zh) | 2009-01-20 | 2009-01-20 | 基于正交小波包变换的双模式变步长盲均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910028459 CN101478349B (zh) | 2009-01-20 | 2009-01-20 | 基于正交小波包变换的双模式变步长盲均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101478349A CN101478349A (zh) | 2009-07-08 |
CN101478349B true CN101478349B (zh) | 2013-05-08 |
Family
ID=40838986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910028459 Expired - Fee Related CN101478349B (zh) | 2009-01-20 | 2009-01-20 | 基于正交小波包变换的双模式变步长盲均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101478349B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101656697B (zh) * | 2009-09-18 | 2013-04-03 | 南京信息工程大学 | 基于t/2分数间隔的频域盲均衡方法 |
CN101651643B (zh) * | 2009-09-18 | 2013-01-02 | 南京信息工程大学 | 基于空间分集的小波神经网络盲均衡方法 |
CN101710884B (zh) * | 2009-12-04 | 2013-07-10 | 深圳国微技术有限公司 | 基于信道估计和均衡的qam模式识别的方法 |
JP2011234223A (ja) * | 2010-04-28 | 2011-11-17 | Toshiba Corp | 等化装置及び放送受信装置 |
CN101924718B (zh) * | 2010-08-30 | 2013-07-03 | 南京信息工程大学 | 模糊神经网络控制的混合小波神经网络盲均衡方法 |
CN102184455B (zh) * | 2011-04-15 | 2013-12-04 | 南京信息工程大学 | 基于自适应免疫克隆的正交小波超指数迭代盲均衡方法 |
CN102361475B (zh) * | 2011-06-15 | 2014-07-16 | 南京信息工程大学 | 基于混沌支持向量机优化的小波加权多模盲均衡方法 |
CN106850035A (zh) * | 2016-12-15 | 2017-06-13 | 四川九洲电器集团有限责任公司 | 一种多通道幅相失真校正的方法 |
CN111800356B (zh) * | 2020-06-16 | 2023-01-31 | 北京银河信通科技有限公司 | 并行变步长cma均衡算法、装置、电子设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395780A (zh) * | 2000-09-11 | 2003-02-05 | 福克斯数码公司 | 用于使用能够利用自适应信道均衡器中的目标权重矢量的稀疏性的自适应算法的设备和方法 |
-
2009
- 2009-01-20 CN CN 200910028459 patent/CN101478349B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395780A (zh) * | 2000-09-11 | 2003-02-05 | 福克斯数码公司 | 用于使用能够利用自适应信道均衡器中的目标权重矢量的稀疏性的自适应算法的设备和方法 |
Non-Patent Citations (9)
Title |
---|
> * |
< * |
.1997,第18卷(第2期), * |
.2002,第23卷(第6期), * |
徐金标,葛建华等.基于CMA算法的双模式盲均衡算法.< * |
徐金标,葛建华等.基于CMA算法的双模式盲均衡算法.<<通信学报>>.1997,第18卷(第2期), |
王蜂,赵俊渭等.一种常数模与判决导引相结合的盲均衡算法研究.< * |
王蜂,赵俊渭等.一种常数模与判决导引相结合的盲均衡算法研究.<<通信学报>>.2002,第23卷(第6期), |
通信学报> * |
Also Published As
Publication number | Publication date |
---|---|
CN101478349A (zh) | 2009-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101478349B (zh) | 基于正交小波包变换的双模式变步长盲均衡方法 | |
Teo et al. | Optimal MMSE finite parameter model for doubly-selective channels | |
CN101355541B (zh) | 快变信道条件下正交频分复用系统中分块均衡方法 | |
CN108712353B (zh) | 软迭代信道估计方法 | |
Han et al. | A deep learning based channel estimation scheme for IEEE 802.11 p systems | |
CN100385824C (zh) | 一种mimo-ofdm系统的自适应信道估计方法 | |
CN110266617B (zh) | 超奈奎斯特系统的多径信道估计方法 | |
CN101478350A (zh) | 基于正交小波变换的超指数迭代联合盲均衡方法 | |
CN113381951A (zh) | 时变频选衰落信道下mftn联合信道估计与均衡方法 | |
CN107508778B (zh) | 一种循环相关信道估计方法及装置 | |
CN115296970B (zh) | 基于逐元素外部信息的迭代正交时频空波形检测方法 | |
Zhang et al. | An adaptive matching pursuit algorithm for sparse channel estimation | |
CN111355677A (zh) | 一种基于滤波器组的多载波水下高速通信系统 | |
CN101656697B (zh) | 基于t/2分数间隔的频域盲均衡方法 | |
CN101447969A (zh) | 一种多带正交频分复用超宽带系统的信道估计方法 | |
Pereira et al. | Tibwb-ofdm: A promising modulation technique for mimo 5g transmissions | |
CN104135455B (zh) | 一种通信系统迭代接收方法 | |
CN105471774A (zh) | 一种相位噪声的估计方法、估计装置、接收机及通信设备 | |
CN113381954B (zh) | 一种基于广义近似消息传递的频域均衡方法 | |
Baumgartner et al. | A soft interference cancellation inspired neural network for SC-FDE | |
Oyerinde | Channel estimation for SISO and MIMO OFDM communications systems. | |
Gümüş et al. | Channel estimation and symbol demodulation for OFDM systems over rapidly varying multipath channels with hybrid deep neural networks | |
WO2008026835A1 (en) | Adaptive frequency domain equalizer and adaptive frequency domain equalization method | |
CN116633737B (zh) | 针对超奈奎斯特系统的低复杂度svd预编码方法 | |
CN116633734B (zh) | 适用于高阶调制的超奈奎斯特系统svd预编码方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130508 Termination date: 20170120 |