CN101656697A - 基于t/2分数间隔的频域盲均衡方法 - Google Patents
基于t/2分数间隔的频域盲均衡方法 Download PDFInfo
- Publication number
- CN101656697A CN101656697A CN 200910183376 CN200910183376A CN101656697A CN 101656697 A CN101656697 A CN 101656697A CN 200910183376 CN200910183376 CN 200910183376 CN 200910183376 A CN200910183376 A CN 200910183376A CN 101656697 A CN101656697 A CN 101656697A
- Authority
- CN
- China
- Prior art keywords
- cma
- equalizer
- frequency
- domain
- equalization method
- 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
Images
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明公布了一种基于T/2分数间隔的频域盲均衡方法(T/2FF-CMA)。本发明方法针对波特间隔盲均衡方法稳态误差大、收敛速度慢的缺点,利用分数间隔盲均衡器来避免因欠采样引起的频谱混叠,以有效地减小盲均衡器的稳态误差;利用快速傅里叶变换技术(FFT)和重叠保留方法实现频域盲均衡方法,来改善盲均衡器的收敛性能,同时减小计算复杂度。本发明方法T/2FF-CMA的稳态误差小、收敛速度快、计算复杂度低。水声信道实施实例表明,与波特间隔时域盲均衡方法(TB-CMA)和分数间隔时域盲均衡方法(TF-CMA)相比,本发明方法T/2FF-FCMA性能更优。
Description
技术领域
发明涉及一种基于T/2分数间隔的频域盲均衡方法,属于克服水声信道多径衰落引起的码间干扰(Inter-Symbol Interference,ISI)的盲均衡方法的技术领域。
背景技术
水声通信系统中,由于多径衰落的影响以及信道失真,导致在接收端产生严重的码间干扰,降低了信息的传输质量。克服码间干扰最有效的手段是在接收端采用均衡技术。传统的均衡技术采用的是由发射机发送训练序列来辅助实现信道的估计和均衡,浪费了有限的带宽资源。盲均衡技术由于不需要发送已知训练序列,节省了带宽,提高了通信效率而成为目前研究的热点(见文献:[1]Guo Yecai,Han Yingge,Rao Wei.Blindequalization algorithms based on different error equations with exponential variable stepsize[C]//The First International Symposium on Test Automation&Instrumentation(ISTAI).Xi’an,China:World Publishing Corporation,2006:497-501;[2]Robert W.Lucky.The adaptive equalizer[J].IEEE Signal Processing.Magazine,2006:104-107;[3]J.Arenas-Garcia,A.R Figueiras-Vidal.Improved blind equalization via adaptive combination ofconstant modulus algorithms[C].2006 IEEE International Conference on Acoustics,Speechand Signal Processing,ICASSP 2006,3:756-759.)。
在众多盲均衡方法中,传统的波特间隔盲均衡方法稳态误差大、收敛速度慢。而分数间隔盲均衡器(FSE-CMA,Fractionally Spaced Equalizer based on Constant ModulusAlgorithm)具有良好的均衡效果,但FSE-CMA方法的计算复杂度较高(见文献:[4]张艳萍,赵俊谓.基于分数间隔的水声信道盲均衡算法研究[J]声学与电子工程,2005,2:21-23;[5]S.A.Elahmar,A.Djebbari M.Bouziani et al.Performance of Fractionally Spaced BlindChannel Shortening[J].International Journal of Applied Engineering Research,2007,2(1):191-199.)。利用快速傅里叶变换技术和重叠保留算法实现了频域最小均方(LMS,Least Mean Square)算法,在获得比时域更好均衡效果的同时,计算量大大减小(见文献:[6]冯存前,张永顺.变步长频域快速自适应收发隔离算法研究[J].电子对抗技术,2004,19(5):22-25)。
发明内容
本发明所要解决的技术问题是针对传统的波特间隔盲均衡方法(BSE-CMA),稳态误差大、收敛速度慢的缺点,利用分数间隔盲均衡器,来避免因欠采样引起的频谱混叠,以有效地减小盲均衡器的稳态误差;利用快速傅里叶变换技术(FFT)和重叠保留算法实现频域盲均衡方法,来改善均衡器的收敛性能,同时减小计算复杂度,从而发明了基于分数间隔的频域盲均衡方法(FF-CMA)。本发明方法稳态误差小、收敛速度快、计算复杂度低。水声信道仿真结果表明,与分数间隔时域盲均衡方法(TF-CMA)和波特间隔时域盲均衡方法(TB-CMA)相比,本发明方法T/2FF-CMA性能更优。
本发明为实现上述目的,采用如下技术方案:
本发明基于T/2分数间隔的频域盲均衡方法,其特征在于包括如下步骤:
第一步:将发射信号序列s(k)分别经过第p个子信道得到第p个子信道输出信号y(p)(k),其中k为时间序列,p=1,2,…P,P为自然数,表示子信道总数;
第二步:将第一步所述的第p个子信道输出信号y(p)(k)经过对应的L点快速傅里叶变换FFT得到第p个子均衡器输入信号Y(p)(K),其中K为正整数,表示由L点数据构成的块数;
第三步:将第二步所述的第p个子均衡器输入信号Y(p)(K)经过L点快速傅里叶反变换IFFT和重叠保留法得到基于分数间隔均衡器的时域输出信号为: 其中 为第p个盲均衡器的频域输出信号,Fk(K)为均衡器的频域权向量。
所述的基于T/2分数间隔的频域盲均衡方法,其特征在于第三步所述的均衡器权向量Fk(k)的更新方法如下:
其中,
式中,其中K为正整数,表示由L点数据构成的块数,0≤K ≤2L-1;E(·)为误差函数,L为自然数,conj(·)表示取复数共轭操作,p=1,2,…P,P为自然数,表示子信道总数。
本发明基于分数间隔的频域盲均衡方法(FF-CMA),利用分数间隔盲均衡器减小因欠采样引起的频谱混叠以减小稳态误差,利用可调滤波器的频域特性去补偿基带系统的频率特性,使包括可调滤波器在内的基带系统的中特性满足无失真传输的要求,同时利用快速傅里叶变换(FFT)技术和重叠保留法计算线性卷积,使计算量大大降低。保证了本发明方法的实时性能。通过水声信道的仿真实验,证明了本发明方法的有效性。
附图说明
图1:基于分数间隔时域盲均衡系统原理图。
图2:本发明系统模型图:基于分数间隔频域盲均衡器结构图。
图3:本发明实例1仿真图,(a)均方误差曲线;(b)均衡器的输入;(c)TB-CMA方法输出;(d)TF-CMA方法输出;(e)T/2FF-CMA方法输出。
图4:本发明实例2仿真图,(a)均方误差曲线;(b)均衡器的输入;(c)TB-CMA方法输出;(d)TF-CMA方法输出;(e)T/2FF-CMA方法输出
具体实施方式
如图1所示,基于分数间隔的时域盲均衡方法
对输入信号以T/P(T是采样间隔,P是正整数)的速率采样,盲均衡器的权向量利用经典的常数模算法进行更新,且更新过程在时域中进行,就得到基于分数间隔的时域盲均衡法(TF-CMA),其结构如图1所示。
图1中,k是时间序列;s(k)是独立同分布的发射信号序列;c(p)(k)是分数间隔均衡器的第p个子信道的脉冲响应;n(p)(k)是第p个子信道加入的高斯白噪声序列;y(p)(k)为第p个子均衡器的输入信号;f(p)(k)为子均衡器的权向量;z(p)(k)是第p个子均衡器的输出信号;z(k)是经过合并后的信号。根据分数间隔盲均衡器和时域盲均衡方法的基本理论,可得子均衡器的输入为
式中,NC为子信道长度。
子信道盲均衡器的输出为
式中,Mf为子信道长度,T表示转置操作。
分数间隔均衡器的输出为
子均衡器权向量的更新公式为
式中,R2=E[|s(k)|4]/E[|s(k)|2]是发射序列s(k)的模,该方法称为基于分数间隔的时域盲均衡算法(TF-CMA)。
如图2所示,基于分数间隔的频域盲均衡方法
频域均衡指的是均衡在频域而不是时域中完成。它的基本思想是利用可调滤波器的频率特性去补偿基带系统的频率特性,使包括可调滤波器在内的基带系统的总特性满足无失真传输的要求。
图2中,Y(p)(K)为对y(p)(k)作L点快速傅里叶变换(FFT)的输入信号,F(p)(K)为第p个子均衡器的频域权向量;Z(p)(K)为子均衡器的频域输出信号。根据傅里叶变换和重叠保留法的原理,将Y(p)(K)分割成L长的块,权系数的更新每L个样点进行一次,而每次的更新是由L个误差信号样点累加结果控制的。运用快速FFT技术和重叠保留法过程如下:
对应盲均衡器的权向量为
此块盲均衡器的频域输出为
则基于分数间隔的盲均衡器的时域输出为
权向量更新公式为
式中,
采用频域CMA方法可以减少运算量。实际上,对于实输入数据,L阶时域CMA均衡器提供L点输出需要2L2实乘。采用频域均衡器需要3个L点FFT和2L个复乘,以提供同样的输出。但对于实输入,由于对称性,有一半权可以不计算。而且L点FFT可以由L/2点FFT和L/2复乘实现。因此频域处理共需要3log2(L/2)+4L个实乘。这样,频域CMA对时域CMA的运算量之比为
当L=16时,该比值为0.41;L=32时,比值为0.25;L=256、N=256时,比值为0.049。因此,当N大时,节省的运算量是很大的。
实施实例
为验证本发明方法:基于T/2分数间隔的频域盲均衡器(T/2FF-CMA)的性能,以基于分数间隔的时域盲均衡器(TF-CMA)、基于波特间隔的时域盲均衡器(TB-CMA)为比较对象,采用水声信道c=[1,0,0.3e-0.7j,0,0,0.2e-0.8j]进行仿真分析。
【实施例1】发射信号为4QAM;方差为1,信噪比为20dB;均衡器权长均为12且均为中心抽头初始化;μTF-CMA=0.002,μFF-CMA=0.004,μT/2FF-CMA=0.004;5000次蒙特卡洛仿真结果,如图3所示。
图3(a)表明,本发明方法T/2FF-CMA的收敛速度比TF-CMA和FB-CMA分别提高了450步和300步;其稳态误差比TB-CMA和TF-CMA分别减小了5dB和3dB。图3(b)至(e)表明,本发明方法T/2FF-CMA比TF-CMA和TB-CMA输出的星座图清晰、集中。
【实施例2】发射信号为4PSK;方差为1,信噪比为20dB;均衡器权长均为12且均为中心抽头初始化;μTF-CMA=0.006,μFF-CMA=0.006,μT/2FF-CMA=0.008;5000次蒙特卡洛仿真结果,如图4所示。
图4(a)表明,本发明方法T/2FF-CMA的收敛速度比TF-CMA和TB-CMA分别提高了450步和200步;其稳态误差比TF-CMA和TB-CMA减小了3dB。图4(b)至(e)表明,本发明方法T/2FF-CMA比TF-CMA和TB-CMA输出的星座图清晰、集中。
Claims (2)
1、一种基于T/2分数间隔的频域盲均衡方法,其特征在于包括如下步骤:
第一步:将发射信号序列s(k)分别经过第p个子信道得到第p个子信道输出信号y(p)(k),其中k为时间序列,p=1,2,…P,P为自然数,表示子信道总数;
第二步:将第一步所述的第p个子信道输出信号y(p)(k)经过对应的L点快速傅里叶变换FFT得到第p个子均衡器输入信号Y(p)(K),其中K为正整数,表示由L点数据构成的块数;
第三步:将第二步所述的第p个子均衡器输入信号Y(p)(K)经过L点快速傅里叶反变换IFFT和重叠保留法得到基于分数间隔均衡器的时域输出信号为: ,其中 为第p个盲均衡器的频域输出信号,Fk(K)为均衡器的频域权向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910183376 CN101656697B (zh) | 2009-09-18 | 2009-09-18 | 基于t/2分数间隔的频域盲均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910183376 CN101656697B (zh) | 2009-09-18 | 2009-09-18 | 基于t/2分数间隔的频域盲均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101656697A true CN101656697A (zh) | 2010-02-24 |
CN101656697B CN101656697B (zh) | 2013-04-03 |
Family
ID=41710792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910183376 Expired - Fee Related CN101656697B (zh) | 2009-09-18 | 2009-09-18 | 基于t/2分数间隔的频域盲均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101656697B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441972A (zh) * | 2013-08-30 | 2013-12-11 | 四川大学 | 一种盲均衡方法及盲均衡系统 |
CN103763228A (zh) * | 2014-01-07 | 2014-04-30 | 南京信息工程大学 | 组合优化自适应频域盲均衡方法及系统 |
CN104660534A (zh) * | 2008-09-18 | 2015-05-27 | 三菱电机株式会社 | 发送装置以及接收装置 |
CN106973024A (zh) * | 2017-03-23 | 2017-07-21 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
CN111245499A (zh) * | 2020-01-08 | 2020-06-05 | 西安电子科技大学 | 基于预整形的时域并行分数间隔均衡器及均衡方法 |
CN111342905A (zh) * | 2018-12-18 | 2020-06-26 | 深圳市中兴微电子技术有限公司 | 一种信号处理方法、装置和计算机存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101083639B (zh) * | 2006-06-01 | 2010-05-12 | 上海无线通信研究中心 | 基于分数倍采样的滤波器组及其滤波方法 |
CN101478349B (zh) * | 2009-01-20 | 2013-05-08 | 南京信息工程大学 | 基于正交小波包变换的双模式变步长盲均衡方法 |
-
2009
- 2009-09-18 CN CN 200910183376 patent/CN101656697B/zh not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104660534A (zh) * | 2008-09-18 | 2015-05-27 | 三菱电机株式会社 | 发送装置以及接收装置 |
CN104660534B (zh) * | 2008-09-18 | 2018-02-09 | 三菱电机株式会社 | 发送装置以及接收装置 |
CN103441972A (zh) * | 2013-08-30 | 2013-12-11 | 四川大学 | 一种盲均衡方法及盲均衡系统 |
CN103763228A (zh) * | 2014-01-07 | 2014-04-30 | 南京信息工程大学 | 组合优化自适应频域盲均衡方法及系统 |
CN103763228B (zh) * | 2014-01-07 | 2017-01-25 | 南京信息工程大学 | 组合优化自适应频域盲均衡方法及系统 |
CN106973024A (zh) * | 2017-03-23 | 2017-07-21 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
CN106973024B (zh) * | 2017-03-23 | 2020-01-03 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
CN111342905A (zh) * | 2018-12-18 | 2020-06-26 | 深圳市中兴微电子技术有限公司 | 一种信号处理方法、装置和计算机存储介质 |
CN111342905B (zh) * | 2018-12-18 | 2022-11-04 | 深圳市中兴微电子技术有限公司 | 一种信号处理方法、装置和计算机存储介质 |
CN111245499A (zh) * | 2020-01-08 | 2020-06-05 | 西安电子科技大学 | 基于预整形的时域并行分数间隔均衡器及均衡方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101656697B (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Teo et al. | Optimal MMSE finite parameter model for doubly-selective channels | |
CN101404631B (zh) | 一种单载波系统的自适应均衡方法 | |
CN101656697B (zh) | 基于t/2分数间隔的频域盲均衡方法 | |
CN101478510B (zh) | 一种自适应均衡器及使用该均衡器的接收机系统 | |
CN107332797B (zh) | 一种电力线ofdm通信系统中的信道估计方法 | |
US20100046599A1 (en) | Apparatus and method for acquiring initial coefficient of decision feedback equalizer using fast fourier transform | |
Nissel et al. | FBMC-OQAM in doubly-selective channels: A new perspective on MMSE equalization | |
CN101309241B (zh) | 全通均衡的时间反演超宽带无线通信方法及系统 | |
CN101409694A (zh) | 一种单载波频域均衡实现方法及系统 | |
CN101656696B (zh) | 基于联合合并时空分集频域小波盲均衡方法 | |
CN101656579A (zh) | 引入小波神经网络的t/2分数间隔盲均衡方法 | |
CN103957177A (zh) | 一种ofdm信号的多级时域盲均衡方法 | |
CN109302240A (zh) | 基于双选择衰落信道的低复杂度osdm串行均衡方法 | |
Zhang et al. | A low-complexity approximate message passing equalizer for OTFS system | |
CN101330480B (zh) | 一种单载波分数阶傅立叶域均衡技术 | |
CN102164106B (zh) | 分数间隔判决反馈瑞利Renyi熵小波盲均衡方法 | |
CN103117969B (zh) | 一种分数低阶统计量模值变换小波频域多模盲均衡方法 | |
CN105119857A (zh) | 一种雷达站间低抖动、抗干扰信号通信链路技术 | |
Sirvi et al. | Wavelet based OFDM system over flat fading channel using NLMS equalization | |
CN109217954B (zh) | 基于双选择衰落信道的低复杂度osdm块均衡方法 | |
Li et al. | Sliding window-frequency domain equalization for multi-mode communication systems | |
Ysebaert et al. | Constraints in channel shortening equalizer design for DMT-based systems | |
Hossain et al. | Performance analysis of RLS and VSS-LMS channel estimation techniques for 4G MIMO OFDM systems | |
Yecai et al. | A modified T/2 fractionally spaced coordinate transformation blind equalization algorithm | |
Xie et al. | Decision feedback equalizer with frequency domain bidirectional noise prediction for MIMO-SCFDE system |
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 |
Granted publication date: 20130403 Termination date: 20150918 |
|
EXPY | Termination of patent right or utility model |