CN101217288B - 虚拟导频辅助信道估计方法 - Google Patents

虚拟导频辅助信道估计方法 Download PDF

Info

Publication number
CN101217288B
CN101217288B CN2007101732933A CN200710173293A CN101217288B CN 101217288 B CN101217288 B CN 101217288B CN 2007101732933 A CN2007101732933 A CN 2007101732933A CN 200710173293 A CN200710173293 A CN 200710173293A CN 101217288 B CN101217288 B CN 101217288B
Authority
CN
China
Prior art keywords
channel
virtual pilot
circletimes
pilot frequency
overbar
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
Application number
CN2007101732933A
Other languages
English (en)
Other versions
CN101217288A (zh
Inventor
任俊彦
李晶峰
叶凡
王雪静
刘亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Original Assignee
Fudan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fudan University filed Critical Fudan University
Priority to CN2007101732933A priority Critical patent/CN101217288B/zh
Publication of CN101217288A publication Critical patent/CN101217288A/zh
Application granted granted Critical
Publication of CN101217288B publication Critical patent/CN101217288B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

本发明属于超宽带技术领域,具体是一种虚拟导频辅助信道估计方法。本发明设计了虚拟导频辅助的信道估计器,对接收到的序列信号经过一定处理后,运用一个32点FFT和LS或LMMSE运算得到某些子载波位置的频域响应值,然后通过线性插值,得到所有子载波位置的频域响应值。本发明可适用于多带正交频分复用超宽带系统。由于采用了新的信道估计算法,信道估计器的复杂度大大降低,能有效地抵抗多径信道带来的符号间干扰,并且能够适用于不同的UWB信道。

Description

虚拟导频辅助信道估计方法
技术领域
本发明属于超宽带技术领域,具体涉及一种适用于MB-OFDM UWB系统的新的信道估计方法,该方法能有效地降低信道估计的复杂度。
背景技术
超宽带(UWB)技术作为一种极具潜力的高速、短距离的无线传输技术,近些年在学术界和工业界都引起了极大的关注。结合多带正交频分复用(MB-OFDM)技术,MB-OFDMUWB系统能有效地抵抗多径衰落和各种窄带干扰(Narrow-Band Interference),在无线手持设备,PC及外围设备以及家庭消费电子类产品等领域有较广的应用前景。
信道估计对于MB-OFDM UWB系统来说是很关键的,尤其是在系统采用了相关解调来提高系统性能的情况下。在MB-OFDM联盟(MBOA)的物理层提议中,在数据符号前插入了前导(preamble)符号来进行频域信道估计,前导符号的具体形式见图1。以前提出的大多数算法都是基于信道脉冲响应(CIR)长度小于OFDM符号中采用的循环前缀(CP)长度的假设。事实上,这种假设仅适用于UWB信道模型1(CM1)。在IEEE P802.15工作组提出的其余3种UWB信道模型(CM2-CM4)中,由于CIR通常比CP长,符号间干扰(ISI)将不能通过简单地去除CP来消除,上述的信道估计方法的估计性能迅速恶化并会出现误差平底(Error Floor)。
基于离散傅立叶变换(DFT)的信道估计方法是一种变换域的方法,该方法利用信道脉冲响应(CIR)在时域的有限长度特性和能量集中的特点,有效地抑制了噪声的影响,提高了估计的均方误差性能。基于DFT的信道估计的基本框图如图2所示。由于额外的快速傅立叶变换/快速傅立叶逆变换(FFT/IFFT)对的引入,基于DFT的信道估计的性能有明显提升,但代价是较高的计算复杂度,额外的面积和功耗。在MB-OFDM UWB系统,子载波数目为128(N=128)。考虑到128点FFT/IFFT的实现复杂度,这种基于DFT的信道估计并不适合于实际的硬件实现。
发明内容
本发明的目的在于提出一种新的信道估计方法,以改善传统的基于DFT算法的信道估计器在硬件开销和抵抗符号间干扰等方面的不足。
本发明提出的信道估计方法,只需要使用一个32点的FFT模块,能有效地抵抗多径信道和噪声带来的符号间干扰。本方法采用了虚拟导频辅助技术,故称为虚拟导频辅助信道估计方法。
本发明提出的信道估计方法,将信道脉冲响应h划分为[h1,h2]两部分,每部分的长度均为32个采样值(即CP的长度)。很明显,这样做h的长度将被限制在两倍CP长度以内,这对大部分UWB信道(如CM1,CM2和CM3)来说是合理的。对于h超过两倍CP长度的情况(如CM4),本发明采用的近似方法会带来一定的性能损失。但由于h超过两倍CP长度的抽头的能量较小,本发明做法带来性能损失是有限的,可以忽略的。
如图1所示,本发明将发送的两个相同前导符号分别划分为c1-c4四部分,每部分长度均为32个值,其相对应位置的接收信号分别为rcp0,r01-r04,rcp1和r11-r14。假设在前导符号之前收到的信号为零,则相应的接收信号(如图1所示)可表示为:
r cp 0 = h 1 ⊗ ‾ c 4 - - - ( 1 )
r 01 = h 1 ⊗ ‾ c 1 + h 1 ⊗ ‾ c 4 + h 2 ⊗ ‾ c 4 - - - ( 2 )
r 02 = h 1 ⊗ ‾ c 2 + h 1 ⊗ ‾ c 1 + h 2 ⊗ ‾ c 1 + h 2 ⊗ ‾ c 4 - - - ( 3 )
r 03 = h 1 ⊗ ‾ c 3 + h 1 ⊗ ‾ c 2 + h 2 ⊗ ‾ c 2 + h 2 ⊗ ‾ c 1 - - - ( 4 )
r 04 = h 1 ⊗ ‾ c 4 + h 1 ⊗ ‾ c 3 + h 2 ⊗ ‾ c 3 + h 2 ⊗ ‾ c 2 - - - ( 5 )
r cp 1 = h 1 ⊗ c 4 + h 2 ⊗ ‾ c 4 + h 2 ⊗ ‾ c 3 - - - ( 6 )
r 11 = h 1 ⊗ ‾ c 1 + h 1 ⊗ ‾ c 4 + h 2 ⊗ c 4 - - - ( 7 )
r 12 = h 1 ⊗ ‾ c 2 + h 1 ⊗ ‾ c 1 + h 2 ⊗ ‾ c 1 + h 2 ⊗ ‾ c 4 - - - ( 8 )
c 13 = h 1 ⊗ ‾ c 3 + h 1 ⊗ ‾ c 2 + h 2 ⊗ ‾ c 2 + h 2 ⊗ ‾ c 1 - - - ( 9 )
r 14 = h 1 ⊗ ‾ c 4 + h 1 ⊗ ‾ c 3 + h 2 ⊗ ‾ c 3 + h 2 ⊗ ‾ c 2 - - - ( 10 )
其中
Figure S2007101732933D000211
表示正向卷积。如果定义 c = a ⊗ ‾ b , 且a,b,c均为M点序列,那么有:
c ( m ) = Σ n = 0 m a ( n ) · ( m - n ) , m = 0 , . . . , M - 1 - - - ( 11 )
Figure 2007101732933_3
表示反向卷积。如果定义c=a
Figure 2007101732933_4
b,且a,b,c均为M点序列,那么有:
c ( m ) = Σ n = 0 N - 1 - m a ( N - 1 - n ) · ( m + n ) , m = 0 , . . . , M - 1 - - - ( 12 )
表示循环卷积。很明显,a与b的循环卷积可表示为: ⊗ b = a ⊗ ‾ b + a ⊗ ‾ b 对上面的一系列公式(1)-(10)进行合并,可推导出:
r11+r02+r03+e04+rcp1=(h1+h2)
Figure 2007101732933_6
(c1+c2+c3+2c4)    (13)
r11+r12+r13+r14+rcp1=(h1+h2)
Figure 2007101732933_7
(c1+c2+c3+2c4)    (14)
由于时域的循环卷积转换到频域中即为乘法,将接收到的序列通过32点FFT转化到频域后可以很简单地根据式(13)、(14)得到(h1+h2)的频域估计值。若F(·)表示FFT变换,上述两式可合并转变成:
Rp=Hp·Cp    (15)
其中Rp=F[((r11+r02+r03+r04+rcp1)+(r11+r12+r13+r14+rcp1))/2]    (16)
Hp=F(h1+h2)    (17)
Cp=F(c1+c2+c3+2c4)    (18)
可以看到,即使前导符号之前的信号不为零,或者说存在着干扰,干扰也是有限的。从时域上看,干扰的存在最多会破坏所接收到的rcp0和r01的数据,而实际上在进行信道估计时本发明提出的算法并没有使用这两部分数据,因此可以说本发明提出的基于32点FFT的信道估计方法是具有鲁棒性的,能有效地消除CIR过长带来的干扰。
下面将介绍本发明提出的基于32点FFT来估计h的频域响应H(k)的方法。该方法并不需要从(h1+h2)中分离出h1和h2,进一步降低了计算复杂度。根据最小二乘法(LS)或最小均方误差(LMMSE)信道估计方法,我们可以得到(h1+h2)的频域估计值Hp(r)。基于本文前述的h的长度小于两倍CP长度的假设,可以对做如下hFFT变换:
H ( k ) = Σ n = 0 N - 1 h ( n ) W N nk = Σ n = 0 N 2 - 1 h ( n ) W N nk + Σ n = N 2 N - 1 h ( n ) W N nk
= Σ n = 0 N 2 - 1 h ( n ) W N nk = Σ n = 0 N 4 - 1 h ( n ) W N nk + Σ n = N 4 N 2 - 1 h ( n ) W N nk
= Σ n = 0 N 4 - 1 h ( n ) W N nk + Σ n = 0 N 4 - 1 h ( n + N 4 ) W N ( n + N 4 ) k
= Σ n = 0 N 4 - 1 [ h ( n ) + h ( n + N 4 ) W N Nk / 4 ] W N nk - - - ( 19 )
对于k=4r(r=0,1,...,N/4,N=128)的情况,式(19)可进一步转化成:
Figure S2007101732933D00035
= Σ n = 0 N 4 - 1 [ h ( n ) + h ( n + N 4 ) ] W N / 4 nr
= Σ n = 0 N 4 - 1 [ h 1 ( n ) + h 2 ( n ) ] W N / 4 nr = H p ( r ) - - - ( 20 )
由式(20),可以知道h经过128点FFT得到的频域响应H(k)与(h1+h2)经过32点FFT的频域响应Hp(r)的对应关系。这样,从式(20)可以看出问题已转化到已知子载波k=4r(r=0,1,...,N/4)处的CIR来估计出其余子载波位置的CIR,这可以通过插值的方法来实现。可以认为,原来的基于前导符号的信道估计方法就转变成了利用导频(pilot)的信道估计方法,本发明称之为虚拟导频辅助的信道估计方法。
出于降低系统复杂度的考虑,本发明采用了线性内插的方法。如式(21)所示,当rD<k<(r+1)D时,非导频子载波处的信道响应值H(k)可由相邻两个导频位置所估计的信道响应来近似。
H ( k ) = H ( rD + l )
= ( H p ( r + 1 ) - H p ( r ) ) l D + H p ( r ) - - - ( 21 )
其中D表示相邻虚拟导频的位置间隔数,l表示该子载波与前一个虚拟导频位置之间的间隔,Hp(r)表示第r个虚拟导频子载波位置的信道响应值。
现将本发明提出的虚拟导频辅助的信道估计方法归纳如下:
(1)将信道脉冲响应h划分为h1和h2两个部分,每一部分的长度均为32个采样值,即为CP长度;
(2)根据(13),(14)等式左边的表达式对接收信号进行处理,得到r11+r02+r03+r04+rcp1和r11+r12+r13+r14+rcp1序列;
(3)根据式(16),对上述两序列的平均值进行32点FFT,得到其频域响应值Rp
(4)根据已知的前导序列C1-C4,由算式(18)计算得到Cp。在已知Rp与Cp的基础上,再根据式(15),选择LS或LMMSE来得到(h1+h2)的频域响应估计值Hp
(5)根据式(20)所描述的Hp(r)和H(k)的相应关系,通过式(21)的线性插值方法得到所有子载波位置的频域响应估计值H(k)。
附图说明
图1前导符号(preamble)结构图
图2基于DFT的信道估计器结构图
图3虚拟导频辅助的信道估计器结构图
图中标号:1为多路选择器(MUX),2为128点傅立叶变换模块,3为均衡模块(Equalization),4为虚拟导频辅助的信道估计器,5为32点傅立叶变换模块,6为LS或LMMSE信道估计模块,7为插值模块(Interpolation)。
具体实施方式
下面结合附图3进一步描述本文发明的算法。
如图3所示,在信道估计阶段,接收数据通过多路复选器(MUX)进入虚拟导频辅助的信道估计器。具体来说,在MUX中,接收数据按照r11+r02+r03+r04+rcp1和r11+r12+r13+r14+rcp1的方式收集并取平均,得到的32点序列即为虚拟导频辅助信道估计器的输入序列。在虚拟导频信道估计器中,输入序列按照前面介绍的式(16)进行32点FFT,得到接收信号的频域响应值Rp。由于前导序列对于接收端是已知的,相应的(c1+c2+c3+2c4)序列的离散傅立叶变化Cp可以预先计算。这样,虚拟导频位置处的信道估计值Hp可选用前面提到的最小二乘法(LS)或最小均方误差(LMMSE)估计方法来得到。在得到虚拟导频位置的信道响应值之后,出于降低系统复杂度的考虑,本发明通过线性插值(Interpolation)的方式得到其余非虚拟导频位置的信道响应值H(k)。如式(21)所示,非虚拟导频子载波处的信道响应值H(k)可由相邻两个虚拟导频位置所估计的信道响应来近似。在我们的设计中,可取D=4。这样式(21)中的乘法可以简化为移位加法的操作,大大简化了硬件的实现复杂度。

Claims (1)

1.一种虚拟导频辅助信道估计方法,其特征在于具体步骤如下:
(1)将信道脉冲响应h划分为h1和h2两个部分,每一部分的长度均为32个采样值;
(2)将发送的两个相同前导序列分别划分为c1-c4四部分,每部分长度均为32个值,其相对应位置的接收信号分别为rcp0,r01-r04,rcp1和r11-r14,按如下公式计算:
Figure FSB00000758070000011
Figure FSB00000758070000012
得到两个序列r11+r02+r03+r04+rcp1和r11+r12+r13+r14+rcp1,这里 表示循环卷积,在频域中即为乘法;
(3)将上述两个序列的平均值进行32点傅立叶变换,得到其频域响应值Rp
Rp=F[((r11+r02+r03+r04+rcp1)+(r11+r12+r13+r14+rcp1))/2]                (16)
(4)根据已知的前导序列c1-c4,按下式进行FFT,得到Cp
Cp=F(c1+c2+c3+2c4)                                    (18)
然后,对下式:
Rp=Hp·Cp                                             (15)
采用最小二乘法或最小均方误差法,计算h1+h2的频域响应估计值Hp;这里F(·)表示傅立叶变换式;
(5)最后根据下式:
Figure FSB00000758070000014
r=0,1,...,N/4,N=128,
当rD<k<(r+1)D时,采用线性插值计算得到所有子载波位置的频域响应估计值H(k):
Figure FSB00000758070000015
其中D表示相邻虚拟导频的位置间隔数,l表示该子载波与前一个虚拟导频位置之间的间隔,Hp(r)表示第r个虚拟导频子载波位置的信道响应值。 
CN2007101732933A 2007-12-27 2007-12-27 虚拟导频辅助信道估计方法 Expired - Fee Related CN101217288B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007101732933A CN101217288B (zh) 2007-12-27 2007-12-27 虚拟导频辅助信道估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007101732933A CN101217288B (zh) 2007-12-27 2007-12-27 虚拟导频辅助信道估计方法

Publications (2)

Publication Number Publication Date
CN101217288A CN101217288A (zh) 2008-07-09
CN101217288B true CN101217288B (zh) 2012-08-22

Family

ID=39623660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007101732933A Expired - Fee Related CN101217288B (zh) 2007-12-27 2007-12-27 虚拟导频辅助信道估计方法

Country Status (1)

Country Link
CN (1) CN101217288B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101364966B (zh) * 2008-09-12 2011-03-16 河南科技大学 多入多出正交频分多址下行系统的时域信道估计方法
CN102065035B (zh) * 2008-12-31 2014-03-12 宁波大学 多带正交频分复用超宽带系统的信道估计方法
CN101447969B (zh) * 2008-12-31 2011-04-20 宁波大学 一种多带正交频分复用超宽带系统的信道估计方法
CN101808053A (zh) * 2010-02-26 2010-08-18 大唐联诚信息系统技术有限公司 基于ofdm的信道估计方法及装置
US8644330B2 (en) * 2011-03-29 2014-02-04 Intel Corporation Architecture and method of channel estimation for wireless communication system
CN104301263B (zh) * 2014-01-23 2019-04-02 河南科技大学 一种多带uwb系统低复杂度信道估计方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1710894A (zh) * 2005-07-07 2005-12-21 上海交通大学 高速移动正交频分复用系统的信道估计方法
CN1747462A (zh) * 2005-10-20 2006-03-15 清华大学 一种正交频分复用系统信号发送及接收方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1710894A (zh) * 2005-07-07 2005-12-21 上海交通大学 高速移动正交频分复用系统的信道估计方法
CN1747462A (zh) * 2005-10-20 2006-03-15 清华大学 一种正交频分复用系统信号发送及接收方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Dong X,Lu W.Song A cK..Linear interpolation in pilot symbol assisted channel estimation for OFDM.《IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS》.2007,第6卷(第5期),1910-1920. *

Also Published As

Publication number Publication date
CN101217288A (zh) 2008-07-09

Similar Documents

Publication Publication Date Title
CN107438038B (zh) 一种fbmc/oqam的导频设计和同步信道估计方法
Raghavendra et al. Improving channel estimation in OFDM systems for sparse multipath channels
CN101217288B (zh) 虚拟导频辅助信道估计方法
CN101350800B (zh) 载波间干扰处理装置和方法以及无线通信系统的接收机
CN101136884B (zh) 用于tds-ofdm系统的信道估计方法
CN102404268B (zh) 高速移动环境下莱斯信道中多普勒频偏估计与补偿方法
CN101951353B (zh) 一种干扰环境下的正交频分复用系统信道估计方法
CN100385824C (zh) 一种mimo-ofdm系统的自适应信道估计方法
CN103051578B (zh) 带有ici消除的迭代差错扩散判决ofdm信道估计方法
CN103905351A (zh) 一种dvb-t2信道估计方法及系统
CN111049766A (zh) 一种用于5g系统pdsch的估计方法
CN103873406A (zh) 水声正交频分复用通信系统帧间干扰消除方法
CN102130860B (zh) 一种带相位补偿的二维离散傅里叶变换信道估计方法
CN110266627B (zh) 基于伪导频和判决反馈的cir和cfo联合估计方法
CN103166878B (zh) 一种信道估计的方法和装置
Jie et al. An improved DFT-based channel estimation algorithm for MIMO-OFDM systems
CN104683268A (zh) 基于qr分解的ofdm系统信道估计方法
CN102647372A (zh) 信道估计方法
KR100992369B1 (ko) Ofdm 시스템의 채널 추정 장치
Shi et al. Block-type pilot channel estimation for OFDM systems under frequency selective fading channels
Chen et al. Multipath effects on time reversal OFDM communications between wireless sensors
Hou et al. Time-domain MMSE channel estimation based on subspace tracking for OFDM systems
Qin et al. Channel estimation with timing offset based on PSD & LS estimation for wireless OFDM systems
Roman Advanced receiver structures for mobile MIMO multicarrier communication systems
CN100502267C (zh) Ofdm系统中信道估计的方法

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: 20120822

Termination date: 20141227

EXPY Termination of patent right or utility model