[go: up one dir, main page]

CN100385824C - An Adaptive Channel Estimation Method for MIMO-OFDM System - Google Patents

An Adaptive Channel Estimation Method for MIMO-OFDM System Download PDF

Info

Publication number
CN100385824C
CN100385824C CNB2004100218648A CN200410021864A CN100385824C CN 100385824 C CN100385824 C CN 100385824C CN B2004100218648 A CNB2004100218648 A CN B2004100218648A CN 200410021864 A CN200410021864 A CN 200410021864A CN 100385824 C CN100385824 C CN 100385824C
Authority
CN
China
Prior art keywords
channel
ofdm symbol
symbol
pilot symbols
ofdm
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
CNB2004100218648A
Other languages
Chinese (zh)
Other versions
CN1658528A (en
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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CNB2004100218648A priority Critical patent/CN100385824C/en
Publication of CN1658528A publication Critical patent/CN1658528A/en
Application granted granted Critical
Publication of CN100385824C publication Critical patent/CN100385824C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Radio Transmission System (AREA)

Abstract

The present invention discloses a channel estimation method for MIMO-OFDM systems, which comprises the transmission of data frame structures and channel estimation. Transmitted data is divided into training OFDM symbols and data transmission OFDM symbols, wherein the training symbols are composed of training pilot symbols and 0 symbols and are used for establishing the initial parameter of channel estimation, and the data transmission OFDM symbols are composed of reference pilot symbols and data symbols; an adaptive algorithm is adopted for the channel estimation. The method of the present invention has the advantage of high spectrum efficiency, and an algorithm of a receiver has the advantage of low calculating complexity.

Description

一种MIMO-OFDM系统的自适应信道估计方法 An Adaptive Channel Estimation Method for MIMO-OFDM System

技术领域:Technical field:

本发明属于无线通信领域,它特别涉及MIMO-OFDM系统中的信道估计。The invention belongs to the field of wireless communication, and particularly relates to channel estimation in MIMO-OFDM system.

背景技术:Background technique:

未来无线通信系统要求具有更高的频谱效率,以提供高速可靠的数据传输服务。多入多出(MIMO)无线通信系统能显著改善通信系统容量或者通信可靠性,亦可提供二者的折衷。(见文献[L.Zheng,et al,“Diversity and Multiplexing:A Fundamental Tradeoff inMultiple-Antenna Channels,IEEE Trans.On Inform.Theory,vol.49,No.5,pp.1073-1096,May2003])。同时,未来移动通信系统应为宽带系统。宽带系统的典型问题是多径效应,为了简化接收机设计,OFDM调制被认为是一种有竞争力的解决方案,而OFDM用于MIMO系统即MIMO-OFDM系统可以同时获得MIMO和OFDM技术的优势,从而被认为在未来无线通信应用中有广泛的应用前景。(见文献[D.Agrawal,et al.,“Space-time coded OFDM for highdata-rate wireless communication over wideband channels,”in IEEE Vehi.Tech.Conference,vol.3,pp.2232--2236,1998.]和[Y.(G.)Li,et.al.,“Transmit diversity for OFDM systems and itsimpact on high-rate data wireless networks,”IEEE J.Select.Areas Commun.,vol.17,pp.1233-1243,July 1999.]Future wireless communication systems require higher spectral efficiency to provide high-speed and reliable data transmission services. A Multiple Input Multiple Output (MIMO) wireless communication system can significantly improve the communication system capacity or communication reliability, and can also provide a compromise between the two. (See literature [L. Zheng, et al, "Diversity and Multiplexing: A Fundamental Tradeoff in Multiple-Antenna Channels, IEEE Trans. On Inform. Theory, vol.49, No.5, pp.1073-1096, May2003]). At the same time, the future mobile communication system should be a broadband system. The typical problem of broadband systems is multipath effect. In order to simplify the receiver design, OFDM modulation is considered as a competitive solution, and OFDM is used in MIMO systems, namely MIMO- The OFDM system can obtain the advantages of MIMO and OFDM technology at the same time, so it is considered to have a wide range of application prospects in future wireless communication applications. (See literature [D.Agrawal, et al., "Space-time coded OFDM for highdata-rate wireless communication over wideband channels," in IEEE Vehi.Tech.Conference, vol.3, pp.2232--2236, 1998.] and [Y.(G.) Li, et.al., "Transmit diversity for OFDM systems and its impact on high-rate data wireless networks,” IEEE J.Select.Areas Commun., vol.17, pp.1233-1243, July 1999.]

MIMO-OFDM系统的频域信号模型为(如图1,图2所示):The frequency domain signal model of the MIMO-OFDM system is (as shown in Figure 1 and Figure 2):

Y i [ n , k ] = Σ q = 1 M T H i , q [ n , k ] S q [ n , k ] + V i [ n , k ] i=1,2,…,MR,k=0,1,…,N-1    (1) Y i [ no , k ] = Σ q = 1 m T h i , q [ no , k ] S q [ no , k ] + V i [ no , k ] i=1, 2, ..., M R , k = 0, 1, ..., N-1 (1)

其中Sq[n,k]为第q发射天线在第n个OFDM符号内第k子载波发送的信号;N为正整数,表示OFDM子载波数;令hi,q[n,l]为第n个OFDM符号内第q发射天线至第i接收天线间第l条多径分量的路径衰落系数: H i , q [ n , k ] = Σ l = 0 L - 1 h i , q [ n , l ] e - j 2 πkl N 为信道频率响应其中L为正整数表示信道阶数;Vi[n,k]为零均值加性白高斯噪音的,满足Among them, S q [n, k] is the signal sent by the qth transmitting antenna in the kth subcarrier in the nth OFDM symbol; N is a positive integer, indicating the number of OFDM subcarriers; let h i, q [n, l] be The path fading coefficient of the lth multipath component between the qth transmit antenna and the ith receive antenna in the nth OFDM symbol: h i , q [ no , k ] = Σ l = 0 L - 1 h i , q [ no , l ] e - j 2 πkl N is the channel frequency response where L is a positive integer representing the channel order; V i [n, k] is zero-mean additive white Gaussian noise, satisfying

EE. {{ VV ii ′′ [[ nno ′′ ,, kk ′′ ]] VV ii ** [[ nno ,, kk ]] }} == σσ VV 22 δδ [[ ii ′′ -- ii ]] δδ [[ nno ′′ -- nno ]] δδ [[ mm ′′ -- mm ]] -- -- -- (( 22 ))

其中δ[m-n]在m-n为0时取值为1,其他情况为0。σV 2表示噪音方差。Among them, δ[mn] takes the value of 1 when mn is 0, and is 0 in other cases. σ V 2 represents the noise variance.

在各种MIMO-OFDM系统中,宽带MIMO时间频率双选择性衰落信道的估计和跟踪是实现接收机最优解码的前提和关键。传统的单发单收(SISO)OFDM系统的信道估计技术取得了丰硕的成果(见Y.(G.)Li,“Pilot-Symbol-Aided Channel Estimation for OFDM inWireless Systems”,IEEE Trans On Vehicular Technology vol..49,no.4,pp1207-1215,July 2000及其参考文献)。近来,MIMO-OFDM系统的信道估计也引起广泛的关注,目前提出的主要的MIMO-OFDM系统信道估计方法主要有:In various MIMO-OFDM systems, the estimation and tracking of wideband MIMO time-frequency dual-selective fading channels is the premise and key to realize the optimal decoding of the receiver. The channel estimation technology of the traditional single-send-single-receive (SISO) OFDM system has achieved fruitful results (see Y. (G.) Li, "Pilot-Symbol-Aided Channel Estimation for OFDM in Wireless Systems", IEEE Trans On Vehicular Technology vol ..49, no.4, pp1207-1215, July 2000 and its references). Recently, the channel estimation of the MIMO-OFDM system has also attracted widespread attention. The main channel estimation methods for the MIMO-OFDM system currently proposed mainly include:

[1]MMSE准则下适用于判决反馈或最优训练序列的信道估计方案(见文献Y.(G.)Li,“Simplified Channel Estimation for OFDM Systems with Multiple Transmit Antennas”,IEEETrans.On Wireless Comm.,vol.1,pp.67-75,Jan.2002),[1] Channel estimation scheme suitable for decision feedback or optimal training sequence under MMSE criterion (see literature Y. (G.) Li, "Simplified Channel Estimation for OFDM Systems with Multiple Transmit Antennas", IEEETrans.On Wireless Comm., vol.1, pp.67-75, Jan.2002),

[2]基于导频符号的时域内插方法(见文献K.Lee and D.B.Williamos,”Pilot-Symbol-Assisted Channel Estimation for Space-Time Coded OFDMSystems,”EURASIP Journal on Applied Signal Processing 2002:5,507-516),[2] Time-domain interpolation method based on pilot symbols (see K.Lee and D.B.Williamos, "Pilot-Symbol-Assisted Channel Estimation for Space-Time Coded OFDM Systems," EURASIP Journal on Applied Signal Processing 2002: 5, 507- 516),

[3]基于子空间等方法的盲信道估计算法(见文献H.

Figure C20041002186400061
,et al.,“Blind ChannelIdentification and Equalization in OFDM-Based Multiantenna Systems”,IEEE Trans On SignalProcessing,vol.50,pp96-109,Jan 2002)[3] Blind channel estimation algorithm based on subspace and other methods (see literature H.
Figure C20041002186400061
, et al., "Blind Channel Identification and Equalization in OFDM-Based Multiantenna Systems", IEEE Trans On Signal Processing, vol.50, pp96-109, Jan 2002)

上述方法分别适用于不同场合,但是也各具一些显著的缺点。方法[1]具有优化的性能但其计算量大,要求接收机处理能力强,成本高昂,难以实用至实际无线通信系统,其中基于训练序列的方法在快衰落应用场合中要求频繁发送训练序列,导致频谱效率损失很大。方法[2]的典型问题是由于Nyquist采样速率对导频符号间隔的限制,在快衰落场合下其频谱效率损失较大,且其频谱效率损失与发射天线数成正比。方法[3]的典型问题是运算复杂度非常大,且会存在相位模糊等问题,实用价值不高。The above methods are applicable to different occasions respectively, but each has some significant disadvantages. Method [1] has optimized performance, but it requires a large amount of calculation, requires strong receiver processing capability, and is expensive, so it is difficult to apply to practical wireless communication systems. The method based on training sequences requires frequent transmission of training sequences in fast fading applications. This leads to a great loss of spectral efficiency. The typical problem of the method [2] is that due to the limitation of the Nyquist sampling rate on the pilot symbol interval, the spectral efficiency loss is relatively large in the case of fast fading, and the spectral efficiency loss is proportional to the number of transmitting antennas. The typical problem of the method [3] is that the computational complexity is very large, and there will be problems such as phase ambiguity, and the practical value is not high.

发明内容Contents of the invention

本发明的目的是提供一种MIMO-OFDM系统信道估计的方法,它具有较高的频谱效率,其接收机算法具有较低的运算复杂度。The purpose of the present invention is to provide a method for channel estimation of MIMO-OFDM system, which has higher spectral efficiency and its receiver algorithm has lower computational complexity.

为说明方便,特定义以下术语:For the convenience of explanation, the following terms are defined:

训练导频符号:收发两端均已知的符号,由数据帧的第一个OFDM符号(即训练OFDM符号)的特定子载波传送,其作用为使接收机信道估计模块估计使用最小二乘(LS)方法估计信道初值。Training pilot symbol: a symbol known at both the sending and receiving ends, transmitted by the specific subcarrier of the first OFDM symbol (i.e. training OFDM symbol) of the data frame, and its function is to make the channel estimation module of the receiver estimate using the least squares ( LS) method to estimate the initial value of the channel.

参考导频符号:收发两端均已知的符号,由数据帧的第一个OFDM符号后的OFDM符号(即数据传输OFDM符号)的特定子载波传送,用于产生自适应信道追踪算法所需的输入信号向量。Reference pilot symbol: A symbol known to both the transmitting and receiving ends, transmitted by a specific subcarrier of the OFDM symbol (that is, the data transmission OFDM symbol) after the first OFDM symbol of the data frame, used to generate the adaptive channel tracking algorithm The input signal vector of .

主多径分量选择(STC):从阶数较大的信道时域脉冲响应估计(由较多的多径衰落系数的估计值组成)中,根据降低估计误差的原则,选择其中较少的多径衰落系数作为信道时域脉冲响应的估计结果,其目的在于减少多径衰落系数,降低参数个数并降低估计误差。其实现方法见Y.(G.)Li,Simplified Channel Estimation for OFDM Systems With MultipleTransmit Antennas,IEEE Trans.On Wirel.Comm.Vol.1,pp.67-75,Jan.2002及其参考文献。Main multipath component selection (STC): From the channel time-domain impulse response estimation with a large order (consisting of more estimated values of multipath fading coefficients), according to the principle of reducing the estimation error, select less of them. The path fading coefficient is used as the estimation result of the channel time domain impulse response, and its purpose is to reduce the multipath fading coefficient, reduce the number of parameters and reduce the estimation error. See Y.(G.) Li, Simplified Channel Estimation for OFDM Systems With MultipleTransmit Antennas, IEEE Trans.On Wirel.Comm.Vol.1, pp.67-75, Jan.2002 and its references for its implementation method.

信道脉冲响应截短:从阶数较大的信道时域脉冲响应估计(由较多的多径衰落系数的估计值组成)中,按照缩短信道脉冲响应长度的原则,从时延最短的多径衰落系数开始,提取时延短的特定数目的多径衰落系数作为信道脉冲响应的估计。Channel impulse response truncation: From the channel time-domain impulse response estimation with a larger order (composed of more estimated values of multipath fading coefficients), according to the principle of shortening the length of the channel impulse response, from the multipath with the shortest delay Starting from fading coefficients, a specific number of multipath fading coefficients with short delays are extracted as channel impulse response estimates.

最小均方算法(LMS算法):一种基于最小均方误差准则(MMSE)的典型的自适应算法,其实现可参照文献Simon Haykin Adaptive Filter Theory,Third edition,Prentice Hall1998。Least Mean Square Algorithm (LMS Algorithm): A typical adaptive algorithm based on the minimum mean square error criterion (MMSE). For its implementation, please refer to the literature Simon Haykin Adaptive Filter Theory, Third edition, Prentice Hall1998.

迭代最小二乘算法(RLS算法):一种基于最小二乘准则(LS)的典型的自适应算法,其实现可参照文献Simon Haykin Adaptive Filter Theory,Third edition,Prentice Hall 1998。Iterative Least Squares Algorithm (RLS Algorithm): A typical adaptive algorithm based on the Least Squares Criterion (LS), its implementation can refer to the literature Simon Haykin Adaptive Filter Theory, Third edition, Prentice Hall 1998.

QR-RLS算法:一种RLS算法的实现方法,其实现可参照文献Simon Haykin AdaptiveFilter Theory,Third edition,Prentice Hall 1998。QR-RLS algorithm: an implementation method of RLS algorithm, its implementation can refer to the literature Simon Haykin AdaptiveFilter Theory, Third edition, Prentice Hall 1998.

本发明的一种新的自适应MIMO-OFDM信道估计方法,其特征在于它包含发端步骤和收端步骤:A kind of new adaptive MIMO-OFDM channel estimation method of the present invention is characterized in that it comprises sending end step and receiving end step:

所述发端步骤按照步骤1和步骤2进行:The sending step is carried out according to step 1 and step 2:

步骤1:引入训练导频符号和参考导频符号以形成数据帧结构Step 1: Introduce training pilot symbols and reference pilot symbols to form a data frame structure

对任意发射天线,按照分帧方式传送数据,每帧包含Nt(Nt为正整数)个OFDM符号;帧内第一个OFDM符号为训练OFDM符号(n=0,n表示OFDM符号序号);帧内训练OFDM符号后Nt-1个OFDM符号为数据传输OFDM符号(n=1,2,…,Nt-1)。(如图3所示);按照下述方法插入训练导频符号和插入参考导频符号:For any transmitting antenna, transmit data in a frame-by-frame manner, each frame contains N t (N t is a positive integer) OFDM symbols; the first OFDM symbol in the frame is a training OFDM symbol (n=0, n represents the OFDM symbol sequence number) ; N t −1 OFDM symbols after the intra-frame training OFDM symbols are OFDM symbols for data transmission (n=1, 2, . . . , N t −1). (as shown in Figure 3); insert training pilot symbols and insert reference pilot symbols according to the following method:

插入训练导频符号:在训练OFDM符号中,对于第q(其中q=1,2,…,MT,MT表示发射天线数,为正整数)发射天线,其第一个传送训练导频符号的子载波号为q-1,其后每隔Df0个子载波插入一个训练导频符号;对该发射天线,在非训练导频符号子载波处,传送幅度为0的0符号;(如图4所示)。Insert training pilot symbol: in the training OFDM symbol, for the qth (where q=1, 2, ..., M T , M T represents the number of transmitting antennas, which is a positive integer) transmitting antenna, the first one transmits the training pilot The subcarrier number of the symbol is q-1, and a training pilot symbol is inserted every D f0 subcarriers thereafter; for the transmitting antenna, at the non-training pilot symbol subcarrier, the transmission amplitude is 0 symbols of 0; (such as Figure 4).

插入参考导频符号:在任意数据传输OFDM符号中,所有发射天线使用相同的子载波传送参考导频符号;对任意发射天线,传送第1个参考导频符号的子载波序号为k0(0≤k0<Df1,Df1为正整数),其后每隔Df1个子载波插入一个参考导频符号;对任意发射天线,在非参考导频符号子载波处,传送空时编码输出的数据符号。传送参考导频符号的子载波序号为kp,p=0,1,…,P-1,

Figure C20041002186400081
N为OFDM子载波数;Insert reference pilot symbols: In any data transmission OFDM symbol, all transmit antennas use the same subcarrier to transmit reference pilot symbols; for any transmit antenna, the subcarrier sequence number of the first reference pilot symbol is k 0 (0 ≤k 0 <D f1 , D f1 is a positive integer), and then a reference pilot symbol is inserted every D f1 subcarriers; for any transmitting antenna, at the non-reference pilot symbol subcarrier, transmit the space-time coding output data symbol. The sequence number of the subcarrier transmitting the reference pilot symbol is k p , p=0, 1, ..., P-1,
Figure C20041002186400081
N is the number of OFDM subcarriers;

步骤2:按照步骤1形成的帧结构发送。Step 2: Send according to the frame structure formed in step 1.

所述收端步骤按照下述阶段1和阶段2进行:The receiving step is carried out according to the following phase 1 and phase 2:

阶段1:在每帧的训练OFDM符号处估计信道,以获取本帧内用于数据传输期间自适应信道追踪的初始参数,按照A、B和C三个步骤顺次进行:Phase 1: Estimate the channel at the training OFDM symbol of each frame to obtain the initial parameters used for adaptive channel tracking during data transmission in this frame, and follow the three steps of A, B and C in sequence:

步骤A:用最小二乘(即LS)估计获取当前信道所有训练导频符号子载波的频率响应的估计值:Step A: use the least squares (ie LS) estimation to obtain the estimated value of the frequency response of all training pilot symbol subcarriers of the current channel:

步骤B:对步骤A的结果,使用频域->时域变换估计所有发射天线至所有接收天线间信道脉冲响应(记其信道脉冲响应中,第n个OFDM符号中第q发射天线至第i接收天线间第l条多径分量多径衰落系数的估计为

Figure C20041002186400082
其中0≤l≤N-1)。其中i=1,2,…,MR,MR为接收天线数。Step B: For the result of step A, use the frequency domain->time domain transformation to estimate the channel impulse response between all transmitting antennas and all receiving antennas (record the channel impulse response, from the qth transmitting antenna to the ith in the nth OFDM symbol The estimation of the multipath fading coefficient of the lth multipath component between receiving antennas is
Figure C20041002186400082
where 0≤l≤N-1). Wherein, i=1, 2, ..., M R , where M R is the number of receiving antennas.

步骤C:获取信道主多径分量延迟时间并提取主多径分量衰落系数:Step C: Obtain the delay time of the main multipath component of the channel and extract the fading coefficient of the main multipath component:

利用步骤B中计算结果,使用信道脉冲响应直接截短法或者主多径分量选择(STC)法,对任意发射天线q至任意接收天线i之间的信道,从其信道多径衰落系数

Figure C20041002186400083
中选择(
Figure C20041002186400085
为正整数)条主多径,存储其延迟时间(以采样周期为单位)为li,q,k和相应路径的衰落系数
Figure C20041002186400086
其中 k = 1,2 , &CenterDot; &CenterDot; &CenterDot; , L &LeftRightArrow; . Using the calculation results in step B, using the channel impulse response direct truncation method or the main multipath component selection (STC) method, for the channel between any transmitting antenna q to any receiving antenna i, from its channel multipath fading coefficient
Figure C20041002186400083
choose from (
Figure C20041002186400085
is a positive integer) main multipath, store its delay time (in sampling period) as l i, q, k and the fading coefficient of the corresponding path
Figure C20041002186400086
in k = 1,2 , &Center Dot; &Center Dot; &CenterDot; , L &LeftRightArrow; .

阶段2:数据传输期间自适应信道追踪:Phase 2: Adaptive channel tracking during data transmission:

对帧内的第n(n=1,2,…,Nt-1)个数据传输OFDM符号,对其中任一接收天线i,使用如下信道估计方法(对所有接收天线重复下述步骤A、B和C):For the nth (n=1, 2, ..., N t -1) data transmission OFDM symbol in the frame, for any receiving antenna i, use the following channel estimation method (repeat the following steps A, B and C):

步骤A:设置自适应信道跟踪算法初值:Step A: Set the initial value of the adaptive channel tracking algorithm:

设置自适应信道估计初值为第n-1个OFDM符号信道估计的终值:Set the initial value of the adaptive channel estimation as the final value of the channel estimation of the n-1th OFDM symbol:

n=1,2,…,Nt-1    (3) n=1, 2, ..., N t -1 (3)

其中:

Figure C20041002186400089
表示第n个OFDM符号内,经过p次迭代后,所有发射天线至接收天线i之间信道的估计结果,
Figure C200410021864000810
表示第n个OFDM符号中,经过P次迭代后,所有发射天线至接收天线i之间信道估计的最终的结果,
Figure C200410021864000811
表示第n个OFDM符号中,经过P次迭代后,第q个发射天线至接收天线i之间信道估计的最终的结果,其中in:
Figure C20041002186400089
Indicates the channel estimation results between all transmit antennas and receive antenna i after p iterations in the nth OFDM symbol,
Figure C200410021864000810
Indicates the final result of channel estimation between all transmit antennas and receive antenna i after P iterations in the nth OFDM symbol,
Figure C200410021864000811
Indicates the final result of channel estimation between the qth transmit antenna and the receive antenna i after P iterations in the nth OFDM symbol, where

Figure C200410021864000812
Figure C200410021864000812

Figure C20041002186400091
Figure C20041002186400091

步骤B:自适应信道跟踪:Step B: Adaptive channel tracking:

任意参考导频符号子载波均对应于自适应算法的一次迭代,对任意的p(0≤p≤P-1,P为参考导频符号个数):Any reference pilot symbol subcarrier corresponds to one iteration of the adaptive algorithm, for any p (0≤p≤P-1, P is the number of reference pilot symbols):

自适应算法的参考信号为Yi[n,kp](第n个OFDM符号期间第i接收天线第p个参考导频子载波处的接收信号),其中kp为传送第p个参考导频符号的子载波序号;The reference signal of the adaptive algorithm is Y i [n, k p ] (the received signal at the pth reference pilot subcarrier of the i receiving antenna during the nth OFDM symbol), where k p is the transmission of the pth reference pilot The subcarrier sequence number of the frequency symbol;

自适应算法的输入信号向量wp[n]为:The input signal vector wp [n] of the adaptive algorithm is:

ww pp [[ nno ]] == [[ ww 11 ,, pp TT [[ nno ]] ,, ww 22 ,, pp TT [[ nno ]] ,, &CenterDot;&Center Dot; &CenterDot;&Center Dot; &CenterDot;&Center Dot; ,, ww Mm TT ,, pp TT [[ nno ]] ]] TT -- -- -- (( 66 ))

其中

Figure C20041002186400093
Figure C20041002186400094
表示第n个数据传输OFDM内第q发射天线的第p个参考导频符号的值,接收端已知该值;in
Figure C20041002186400093
Figure C20041002186400094
Indicates the value of the p-th reference pilot symbol of the q-th transmit antenna in the n-th data transmission OFDM, which is known at the receiving end;

自适应信道估计:可采用典型的自适应算法(包括LMS算法,RLS算法和QR-RLS算法中的任何一种)调整

Figure C20041002186400095
在当前数据传输OFDM符号内,迭代共进行P次(P为每个OFDM符号中参考导频符号的个数),得到
Figure C20041002186400096
作为本OFDM符号内的时域信道估计的结果,即
Figure C20041002186400097
Adaptive channel estimation: It can be adjusted by typical adaptive algorithms (including any one of LMS algorithm, RLS algorithm and QR-RLS algorithm)
Figure C20041002186400095
In the current data transmission OFDM symbol, the iteration is performed for a total of P times (P is the number of reference pilot symbols in each OFDM symbol), and we get
Figure C20041002186400096
As a result of the time-domain channel estimation within this OFDM symbol, that is
Figure C20041002186400097

步骤C:利用步骤B自适应信道估计的结果估计所有信道子载波频率响应按照(4)和(5)式组合的次序,拆分

Figure C20041002186400098
得到所有发射天线至接收天线i间的信道各主要多径分量的估计。记第q发射天线至第i接收天线间主多径分量的估计值为
Figure C20041002186400099
使用傅立叶变换(或者FFT)后得到第n个OFDM符号期间第q发射天线至第i接收天线间信道第k子载波频域信道估计。Step C: Utilize the result of adaptive channel estimation in step B to estimate the frequency responses of all channel subcarriers according to the combination order of (4) and (5), split
Figure C20041002186400098
Estimates of each main multipath component of the channel between all transmitting antennas and receiving antenna i are obtained. Note that the estimated value of the main multipath component between the qth transmitting antenna and the ith receiving antenna is
Figure C20041002186400099
After Fourier transform (or FFT) is used, the frequency-domain channel estimation of the kth subcarrier of the channel between the qth transmitting antenna and the ith receiving antenna during the nth OFDM symbol period is obtained.

需要说明的是,本发明方案的实现是在导频符号的插入和接收端处理算法上使用本发明所提出的方法的基础上,采用典型基于导频符号辅助信道估计的MIMO-OFDM系统的方法实现。It should be noted that the implementation of the scheme of the present invention is based on using the method proposed by the present invention on the insertion of pilot symbols and the processing algorithm at the receiving end, and adopts the method of a typical MIMO-OFDM system based on pilot symbol-assisted channel estimation accomplish.

本发明设计基本原理:Basic principle of design of the present invention:

1)对于不同发射天线,在训练OFDM符号处采用不同子载波传送训练导频信号,可以将MIMO信道转换为多个独立而相互正交的SISO信道,因而借助于简单的LS方法估计训练符号处的所有天线对间的信道响应;1) For different transmit antennas, different subcarriers are used to transmit training pilot signals at the training OFDM symbols, and the MIMO channel can be converted into multiple independent and mutually orthogonal SISO channels, so the simple LS method is used to estimate the The channel response between all antenna pairs of ;

2)多径信道所有多径分量经历不同的衰落和噪音干扰,因此剔除信噪比低的多径分量能有效改善信道估计的准确性。2) All multipath components of a multipath channel experience different fading and noise interference, so eliminating multipath components with low SNR can effectively improve the accuracy of channel estimation.

3)传统的基于训练的信道估计器只适合于信道时变速度低的场合,而借助于分布插入地训练导频符号和参考导频符号的信道估计方法可以有效跟踪时变信道。3) The traditional channel estimator based on training is only suitable for the occasions where the time-varying speed of the channel is low, and the channel estimation method with the help of distributed training pilot symbols and reference pilot symbols can effectively track the time-varying channel.

4)最优信道估计器往往需要有信道的统计特性作为先验信息,这在实际应用中往往需要训练的时间开销,而自适应信道估计可以避免这一要求。另外自适应信道估计也避免了最优信道估计在信道统计特性变化时的性能恶化问题。4) The optimal channel estimator often needs the statistical characteristics of the channel as prior information, which often requires training time overhead in practical applications, and adaptive channel estimation can avoid this requirement. In addition, the adaptive channel estimation also avoids the performance deterioration problem of the optimal channel estimation when the channel statistical characteristics change.

5)自适应算法需要有一定数量的参考信号才能收敛,同时信道响应在连续地OFDM符号之间是连续变化的,因而将上一OFDM符号的信道估计结果作为当前OFDM符号自适应估计的初值能降低所需的参考数据数,从而可以避免频谱效率的恶化。5) The adaptive algorithm needs a certain number of reference signals to converge, and the channel response changes continuously between consecutive OFDM symbols, so the channel estimation result of the previous OFDM symbol is used as the initial value of the adaptive estimation of the current OFDM symbol The required number of reference data can be reduced, so that deterioration of spectral efficiency can be avoided.

本发明的创新之处在于:The innovation of the present invention is:

1)发送数据由2部分组成,信道估计亦分为2个阶段:发送数据分为训练OFDM符号和数据传输OFDM符号(如图3所示)。其中训练OFDM符号由训练导频符号和0符号组成,用于建立信道估计的初始参数;数据传输OFDM符号由参考导频符号和数据符号组成。本发明中参考导频符号与传统的导频符号的区别在于其作用是用于产生自适应信道估计算法的输入信号,以自适应追踪时变的宽带信道(注意:而非用于和插值方法相结合以估计信道频率响应)。关于参考导频符号的另外一个显著的优点是,不同发射天线的参考导频符号占据相同的子载波,因而使参考导频符号引起的频率开销与发射天线数无关,且不需要特殊设计。其训练导频符号和参考导频符号的插入图案如图4所示。1) The transmitted data consists of two parts, and the channel estimation is also divided into two stages: the transmitted data is divided into training OFDM symbols and data transmission OFDM symbols (as shown in Figure 3). The training OFDM symbols are composed of training pilot symbols and 0 symbols, which are used to establish initial parameters for channel estimation; the data transmission OFDM symbols are composed of reference pilot symbols and data symbols. The difference between the reference pilot symbol and the traditional pilot symbol in the present invention is that its function is used to generate the input signal of the adaptive channel estimation algorithm to adaptively track the time-varying broadband channel (note: not for the sum interpolation method combined to estimate the channel frequency response). Another significant advantage about reference pilot symbols is that reference pilot symbols of different transmit antennas occupy the same subcarrier, so the frequency overhead caused by reference pilot symbols has nothing to do with the number of transmit antennas, and no special design is required. The insertion patterns of training pilot symbols and reference pilot symbols are shown in FIG. 4 .

2)信道估计算法采用自适应算法:自适应算法具有能随着信道变化而动态跟踪信道参数的功能,具有无需信道统计特性、计算量小等特点。2) The channel estimation algorithm adopts an adaptive algorithm: the adaptive algorithm has the function of dynamically tracking channel parameters as the channel changes, and has the characteristics of no need for channel statistical characteristics and a small amount of calculation.

本发明具有以下优点:The present invention has the following advantages:

1)信道估计使用训练OFDM符号作为信道估计的初值,加快自适应信道估计算法的收敛速度;1) Channel estimation uses training OFDM symbols as the initial value of channel estimation to speed up the convergence speed of the adaptive channel estimation algorithm;

2)数据传输OFDM符号期间所有发射天线使用相同子载波传送参考导频符号;且其作用是用于产生自适应算法的输入信号和参考信号而非用作插值,具有高频谱利用率;2) During data transmission OFDM symbols, all transmitting antennas use the same subcarrier to transmit reference pilot symbols; and its function is to generate the input signal and reference signal of the adaptive algorithm instead of interpolation, with high spectrum utilization;

3)采用自适应算法跟踪信道,无需信道统计特性;3) Adaptive algorithm is used to track the channel without channel statistical characteristics;

4)信道追踪在时域进行;4) Channel tracking is performed in the time domain;

5)本发明方法可以采用主多径分量选择技术以降低信道估计初值的误差,并能减少需要追踪的信道参数的个数,从而降低运算量并提高估计准确度。5) The method of the present invention can adopt the main multipath component selection technology to reduce the error of the initial value of channel estimation, and can reduce the number of channel parameters to be tracked, thereby reducing the amount of computation and improving the estimation accuracy.

附图说明Description of drawings

图1是典型的MIMO-OFDM系统发射框图Figure 1 is a typical MIMO-OFDM system transmission block diagram

其中Sq[n,k]为第nOFDM符号期间第q发射天线第k子载波传送的信号;

Figure C20041002186400111
为第n个OFDM符号期间第q发射天线第
Figure C20041002186400112
采样时刻传送的时域信号;n表示一帧内OFDM符号的序号(n为整数,0 ≤n<NT),k表示OFDM子载波序号(k为整数,0 ≤k<N,N表示OFDM子载波总数),
Figure C20041002186400113
表示采样时刻(
Figure C20041002186400114
为整数, - N g &le; m ~ < N , 其中Ng表示OFDM循环前缀采样点数)。其中的帧形成模块完成插入训练符号和导频符号的功能,本发明中,帧形成模块需插入特殊的训练导频符号和参考导频符号。Wherein S q [n, k] is the signal transmitted by the kth subcarrier of the qth transmit antenna during the nth OFDM symbol period;
Figure C20041002186400111
is the qth transmit antenna during the nth OFDM symbol
Figure C20041002186400112
The time-domain signal transmitted at the sampling moment; n represents the sequence number of OFDM symbols in a frame (n is an integer, 0 ≤ n< NT ), k represents the sequence number of OFDM subcarriers (k is an integer, 0 ≤ k<N, N represents OFDM total number of subcarriers),
Figure C20041002186400113
Indicates the sampling time (
Figure C20041002186400114
is an integer, - N g &le; m ~ < N , Where N g represents the number of OFDM cyclic prefix sampling points). The frame forming module therein completes the function of inserting training symbols and pilot symbols. In the present invention, the frame forming module needs to insert special training pilot symbols and reference pilot symbols.

图2是典型的MIMO-OFDM接收机框图Figure 2 is a block diagram of a typical MIMO-OFDM receiver

其中

Figure C20041002186400116
表示第n个OFDM符号期间第i根接收天线第
Figure C20041002186400117
时刻接收到的采样信号,Yi[n,k]为接收机第n个OFDM符号期间第i接收天线第k子载波分离出来用于空时解码的解码器输入信号;其中i为正整数,1≤i≤MR,MR为接收天线数;其中的信道估计模块使用从接收信号中提取的导频符号估计MIMO信道频率响应。in
Figure C20041002186400116
Indicates that during the nth OFDM symbol period, the i-th receiving antenna
Figure C20041002186400117
The sampled signal received at any time, Y i [n, k] is the decoder input signal separated from the kth subcarrier of the i-th receiving antenna during the n-th OFDM symbol of the receiver for space-time decoding; where i is a positive integer, 1≤i≤M R , MR is the number of receiving antennas; the channel estimation module uses the pilot symbols extracted from the received signal to estimate the MIMO channel frequency response.

图3是本发明的数据传输帧结构Fig. 3 is the data transmission frame structure of the present invention

每一帧内共有Nt个OFDM符号,其中第一个(n=0)OFDM符号为训练OFDM符号,其后共Nt-1个OFDM符号为数据传输OFDM符号(n=1,2,…,Nt-1)。There are N t OFDM symbols in each frame, wherein the first (n=0) OFDM symbol is a training OFDM symbol, and the following N t -1 OFDM symbols are data transmission OFDM symbols (n=1, 2, ... , N t -1).

图4是本方案设计的训练导频符号、参考导频符号图案Figure 4 is the pattern of training pilot symbols and reference pilot symbols designed in this scheme

图中给出了发射天线q(其中q=1,2,…,MT,MT为发射天线数)的训练导频符号和参考导频符号的构成图案;图中第一列表示训练OFDM符号,其余共Nt-1列表示数据传输OFDM符号;n表示OFDM符号的序号,0≤n≤Nt-1),k表示子载波序号,0≤k≤N-1,Nt为每帧OFDM符号个数;Df0为训练OFDM符号中,任意发射天线二相邻训练导频符号子载波间隔;Df1为任意发射天线数据传输OFDM符号中二相邻参考导频符号的子载波间隔;

Figure C20041002186400118
表示训练导频符号,
Figure C20041002186400119
表示参考导频符号,
Figure C200410021864001110
表示数据符号,
Figure C200410021864001111
表示0符号。The figure shows the composition patterns of the training pilot symbols and reference pilot symbols of the transmitting antenna q (where q=1, 2, ..., M T , M T is the number of transmitting antennas); the first column in the figure represents the training OFDM symbol, and the remaining N t -1 columns represent OFDM symbols for data transmission; n represents the sequence number of OFDM symbols, 0≤n≤N t -1), k represents the sequence number of subcarriers, 0≤k≤N-1, N t is each The number of frame OFDM symbols; D f0 is the subcarrier spacing of two adjacent training pilot symbols in any transmitting antenna in the training OFDM symbol; D f1 is the subcarrier spacing of two adjacent reference pilot symbols in the data transmission OFDM symbol of any transmitting antenna ;
Figure C20041002186400118
Denotes training pilot symbols,
Figure C20041002186400119
Indicates the reference pilot symbol,
Figure C200410021864001110
represents the data symbol,
Figure C200410021864001111
Represents the 0 symbol.

图5是本发明信道自适应追踪方案示意图Fig. 5 is a schematic diagram of the channel adaptive tracking scheme of the present invention

根据当前OFDM符号的类型,当为帧内第一个OFDM符号即训练OFDM符号时,使用最小二乘(LS)算法估计当前信道响应并选择主多径分量,这位接收端处理的第一个阶段;当当前OFDM符号为数据传输OFDM符号时,进入接收端信道估计第二阶段,即使用自适应信道跟踪算法跟踪信道。According to the type of the current OFDM symbol, when the first OFDM symbol in the frame is the training OFDM symbol, use the least squares (LS) algorithm to estimate the current channel response and select the main multipath component, which is the first one processed by the receiver stage; when the current OFDM symbol is an OFDM symbol for data transmission, enter the second stage of channel estimation at the receiving end, that is, use an adaptive channel tracking algorithm to track the channel.

图6是使用本发明方案信道估计的系统性能仿真举例Fig. 6 is a system performance simulation example using channel estimation of the scheme of the present invention

图中给出了使用本发明方案的一个特例的误比特率(BER)计算机方针性能曲线,其中SNR表示信噪比,esti表示使用本发明方法估计的结果,ideal表示理想信道估计的结果。自适应算法采用的是LMS算法。图中给出了下述具体实施方式中参数选择下的系统误比特率性能。其中的虚线为理想信道估计性能,实线为本文方案的性能。The figure shows a bit error rate (BER) computer target performance curve of a special case of the scheme of the present invention, wherein SNR represents the signal-to-noise ratio, esti represents the result estimated using the method of the present invention, and ideal represents the result of ideal channel estimation. The adaptive algorithm adopts the LMS algorithm. The figure shows the bit error rate performance of the system under parameter selection in the following specific embodiments. The dotted line is the ideal channel estimation performance, and the solid line is the performance of the scheme in this paper.

具体实施方式Detailed ways

下面以给出一个具体的MIMO-OFDM配置下,本专利的实现步骤。本例中的参数并不影响本发明的一般性。The implementation steps of this patent under a specific MIMO-OFDM configuration are given below. The parameters in this example do not affect the generality of the invention.

发端处理步骤:Originating processing steps:

步骤A:引入训练导频符号和参考导频符号以形成数据帧结构Step A: Introduce training pilot symbols and reference pilot symbols to form a data frame structure

对每一根发射天线,数据传输按照分帧方式传送;每帧第一个OFDM符号为训练OFDM符号;其后为数据传输OFDM符号;每一帧包含的OFDM符号个数为Nt=40。For each transmitting antenna, the data transmission is transmitted in a frame-by-frame manner; the first OFDM symbol in each frame is a training OFDM symbol; the subsequent OFDM symbols are data transmission OFDM symbols; the number of OFDM symbols contained in each frame is N t =40.

训练导频符号位于训练OFDM符号中,不同发射天线采用不同子载波传送训练导频符号。(如图(4)所示)。其第q发射天线第p训练导频符号分配的子载波为

Figure C20041002186400121
The training pilot symbols are located in the training OFDM symbols, and different transmitting antennas use different subcarriers to transmit the training pilot symbols. (As shown in Figure 4). The subcarrier allocated to the pth training pilot symbol of the qth transmit antenna is
Figure C20041002186400121

k &OverBar; q , p = q - 1 + p D f 0 , p = 0,1 , &CenterDot; &CenterDot; &CenterDot; , P &OverBar; - 1 , q=1,2,…,MT k &OverBar; q , p = q - 1 + p D. f 0 , p = 0,1 , &Center Dot; &CenterDot; &Center Dot; , P &OverBar; - 1 , q=1, 2, ..., M T

其中

Figure C20041002186400124
为训练OFDM符号处单个发射天线发射的训练导频符号个数;in
Figure C20041002186400124
The number of training pilot symbols transmitted by a single transmit antenna at the training OFDM symbol;

其中q为发射天线序号,1≤q≤MT,本例中选择发射天线数MT=2,OFDM调制载波数为N=1024;训练OFDM符号中非训练导频符号处传送幅度为0的0符号,Where q is the serial number of the transmitting antenna, 1≤q≤M T , in this example, the number of transmitting antennas M T =2 is selected, and the number of OFDM modulation carriers is N=1024; 0 symbol,

在任意数据传送OFDM符号中,所有发射天线均采用相同的子载波传送参考导频符号。用于传送参考导频符号的子载波的序号为:In any data transmission OFDM symbol, all transmit antennas use the same subcarrier to transmit reference pilot symbols. The sequence numbers of the subcarriers used to transmit the reference pilot symbols are:

kp=pDf1+k0,p=0,1,…,P-1k p = pD f1 + k 0 , p = 0, 1, . . . , P-1

其中共有参考子载波个数为

Figure C20041002186400125
其中本例中选择Df1=16,k0=0 。在其它子载波处传送空时网格编码的编码结果。The total number of reference subcarriers is
Figure C20041002186400125
In this example, D f1 =16 and k 0 =0 are selected. The coding results of space-time trellis coding are transmitted at other subcarriers.

步骤B:按照步骤A的帧结构发送;Step B: send according to the frame structure of step A;

信道估计在接收端进行;按照阶段1和阶段2两个阶段进行:Channel estimation is performed at the receiving end; it is performed in two stages: phase 1 and phase 2:

阶段1:在每帧的训练OFDM符号处估计信道,以获取本帧内用于数据传输期间自适应信道追踪的初始参数,按照A、B和C三个步骤顺次进行:Phase 1: Estimate the channel at the training OFDM symbol of each frame to obtain the initial parameters used for adaptive channel tracking during data transmission in this frame, and follow the three steps of A, B and C in sequence:

步骤A:用最小二乘(即LS)估计获取当前信道所有训练导频符号子载波的频率响应的估计值:Step A: use the least squares (ie LS) estimation to obtain the estimated value of the frequency response of all training pilot symbol subcarriers of the current channel:

Figure C20041002186400131
i=1,2,…,MR,q=1,2,…,MT p = 0,1 , &CenterDot; &CenterDot; &CenterDot; , P &OverBar; - 1 , n=0
Figure C20041002186400131
i=1, 2, ..., M R , q = 1, 2, ..., M T , p = 0,1 , &Center Dot; &Center Dot; &CenterDot; , P &OverBar; - 1 , n=0

其中

Figure C20041002186400133
为发射天线q至接收天线i处第子载波处的频率响应的估计,
Figure C20041002186400135
为第n个OFDM符号内第i接收天线第
Figure C20041002186400136
子载波接收信号;
Figure C20041002186400137
为第q发射天线在第
Figure C20041002186400138
个子载波传送的已知的训练信号;其中i为接收天线序号,1≤i≤MR,本例中选择MR=2。in
Figure C20041002186400133
is the position from the transmitting antenna q to the receiving antenna i An estimate of the frequency response at the subcarriers,
Figure C20041002186400135
is the i-th receive antenna in the n-th OFDM symbol
Figure C20041002186400136
The subcarrier receives the signal;
Figure C20041002186400137
is the qth transmitting antenna at the
Figure C20041002186400138
The known training signal transmitted by subcarriers; where i is the serial number of the receiving antenna, 1≤i≤MR , and M R =2 is selected in this example.

步骤B:使用频域->时域变换估计所有发射天线至所有接收天线间信道脉冲响应:Step B: Estimate the channel impulse response from all transmit antennas to all receive antennas using frequency domain -> time domain transform:

Figure C20041002186400139
i=1,2,…,MR,q=1,2,…,MT,n=0
Figure C20041002186400139
i=1, 2, ..., M R , q = 1, 2, ..., M T , n = 0

其中:in:

Figure C200410021864001311
为帧内第n个OFDM符号处,第q发射天线至第i接收天线间信道脉冲响应的估计,
Figure C200410021864001312
为相应的帧内第n个OFDM符号处,第q发射天线至第i接收天线间信道第l条路径多径系数的估计;
Figure C200410021864001311
is the estimate of the channel impulse response between the qth transmit antenna and the ith receive antenna at the nth OFDM symbol in the frame,
Figure C200410021864001312
is the estimation of the multipath coefficient of the lth path of the channel between the qth transmit antenna and the ith receive antenna at the nth OFDM symbol in the corresponding frame;

Figure C200410021864001313
Figure C200410021864001313

Figure C200410021864001314
为由每根发射天线
Figure C200410021864001315
个训练导频按照步骤A估计的频域信道(在第n个OFDM符号期间的第q发射天线至第i接收天线频域信道系数)。Wq为与第q发射天线相关的
Figure C200410021864001316
维变换矩阵,其第u行第v列元素为:
Figure C200410021864001314
for each transmit antenna
Figure C200410021864001315
The frequency-domain channel estimated by the training pilots according to step A (frequency-domain channel coefficients from the qth transmit antenna to the i-th receive antenna during the nth OFDM symbol period). W q is associated with the qth transmit antenna
Figure C200410021864001316
Dimensional transformation matrix, the elements of row u and column v are:

[[ WW qq ]] uu ,, vv == ee -- jj 22 &pi;&pi; NN (( qq -- 11 ++ (( uu -- 11 )) DD. &Integral;&Integral; 00 )) (( vv -- 11 ))

其中q为发射天线序号,1≤q≤MT u = 1,2 , &CenterDot; &CenterDot; &CenterDot; , P &OverBar; , v=1,2,…,N。Where q is the serial number of the transmitting antenna, 1≤q≤M T , u = 1,2 , &CenterDot; &Center Dot; &Center Dot; , P &OverBar; , v=1, 2, . . . , N.

步骤C:获取信道主多径分量延迟时间并提取主多径分量衰落系数利用步骤B中计算结果,按照主多径分量选择(STC)法,对任意发射天线q至任意接收天线i之间的单发单收(SISO)信道,从其信道多径衰落系数

Figure C200410021864001319
(其中0≤l≤N-1)中选择选择
Figure C200410021864001320
较大的
Figure C200410021864001321
条主多径,记录其延迟时间(以采样周期为单位)为
Figure C200410021864001322
(第q发射天线至第i接收天线之间),相应路径的衰落系数
Figure C200410021864001323
本例中主多径分量选择 L &LeftRightArrow; = 6 . Step C: Obtain the delay time of the main multipath component of the channel and extract the fading coefficient of the main multipath component. Using the calculation results in step B, according to the main multipath component selection (STC) method, for any transmission antenna q to any reception antenna i Single send single receive (SISO) channel, from its channel multipath fading coefficient
Figure C200410021864001319
(where 0≤l≤N-1) choose to choose
Figure C200410021864001320
larger
Figure C200410021864001321
main multipath, record its delay time (in the unit of sampling period) as
Figure C200410021864001322
(between the qth transmit antenna and the i-th receive antenna), the fading coefficient of the corresponding path
Figure C200410021864001323
In this example, the main multipath component selection L &LeftRightArrow; = 6 .

阶段2:数据传输期间自适应信道追踪:Phase 2: Adaptive channel tracking during data transmission:

在接收数据传输OFDM符号时(n=1,2,…,NT-1),对其中任一接收天线i,信道估计方法如下(所有接收天线均需要重复下述A、B和C三个步骤):When receiving data and transmitting OFDM symbols (n=1, 2, ..., N T -1), for any receiving antenna i, the channel estimation method is as follows (all receiving antennas need to repeat the following three steps of A, B and C step):

步骤A:设置自适应信道跟踪算法初值:Step A: Set the initial value of the adaptive channel tracking algorithm:

将所有发射天线至接收天线i之间信道估计结果(实际操作时应存储于接收端存储器)记为以下向量:The channel estimation results between all transmitting antennas and receiving antenna i (should be stored in the memory of the receiving end during actual operation) are recorded as the following vectors:

Figure C20041002186400142
Figure C20041002186400142

其中

Figure C20041002186400143
为第n个OFDM符号中,第q发射天线至第i接收天线之间SISO信道主多径衰落系数的估计值组成的向量。当接收机处理,第1个数据传输OFDM符号时(n=1),
Figure C20041002186400144
由阶段1的步骤C得到;其它的数据传输OFDM符号中(n=2,3,…,Nt-1),  (7)式的值由对第n-1个OFDM符号期间信道自适应追踪的结果得到。此时,设置信道估计初值为上一个OFDM符号信道估计的终值:in
Figure C20041002186400143
is a vector composed of estimated values of the main multipath fading coefficients of the SISO channel between the qth transmit antenna and the ith receive antenna in the nth OFDM symbol. When the receiver processes the first data transmission OFDM symbol (n=1),
Figure C20041002186400144
Obtained by step C of stage 1; in other data transmission OFDM symbols (n=2, 3, ..., N t -1), the value of (7) formula is obtained by channel adaptive tracking during the n-1th OFDM symbol The result is obtained. At this point, set the initial value of the channel estimation to the final value of the channel estimation of the last OFDM symbol:

Figure C20041002186400145
Figure C20041002186400145

其中

Figure C20041002186400146
表示第n个OFDM符号内第i接收天线经过p次迭代后的估计结果,
Figure C20041002186400147
为相应的第n个OFDM符号自适应估计最终的结果,即自适应算法迭代P次后的结果。in
Figure C20041002186400146
Indicates the estimation result of the i-th receiving antenna in the n-th OFDM symbol after p iterations,
Figure C20041002186400147
The final result is adaptively estimated for the corresponding nth OFDM symbol, that is, the result after the adaptive algorithm iterates P times.

步骤B:自适应信道跟踪:Step B: Adaptive channel tracking:

这里采用LMS算法实现信道的自适应追踪。Here, the LMS algorithm is used to realize the adaptive tracking of the channel.

对任意的p(p=0,1,…,P-1)(即在任意数据传输OFDM符号内,对所有参考导频符号子载波):For any p (p = 0, 1, ..., P-1) (i.e., within any data transmission OFDM symbol, for all reference pilot symbol subcarriers):

LMS算法参考信号为Yi[n,kp](第n个OFDM符号期间第i接收天线第p个参考导频子载波处的接收信号),自适应算法的输入信号向量由已知的参考导频信号结合参考导频子载波位置、主多径分量延迟时间得到,产生方法为:The LMS algorithm reference signal is Y i [n, k p ] (the received signal at the pth reference pilot subcarrier of the ith receiving antenna during the nth OFDM symbol period), and the input signal vector of the adaptive algorithm is determined by the known reference The pilot signal is obtained by combining the position of the reference pilot subcarrier and the delay time of the main multipath component. The generation method is:

LMS算法在第p次迭代时的输入信号向量wp[n]构造为:The input signal vector wp [n] of the LMS algorithm at the pth iteration is constructed as:

ww pp [[ nno ]] == [[ ww 11 ,, pp TT [[ nno ]] ,, ww 22 ,, pp TT [[ nno ]] ,, &CenterDot;&Center Dot; &CenterDot;&Center Dot; &CenterDot;&Center Dot; ,, ww Mm TT ,, pp TT [[ nno ]] ]] TT

其中

Figure C20041002186400151
Figure C20041002186400152
表示第n个数据传输OFDM内第q发射天线的第p个参考导频符号的值,接收端已知该值;
Figure C20041002186400153
由阶段1的步骤C选择。在当前数据传输OFDM符号内,迭代共P=64次,得到
Figure C20041002186400154
作为本OFDM符号内的时域信道估计的结果。in
Figure C20041002186400151
Figure C20041002186400152
Indicates the value of the p-th reference pilot symbol of the q-th transmit antenna in the n-th data transmission OFDM, which is known at the receiving end;
Figure C20041002186400153
Selected by Step C of Phase 1. In the current data transmission OFDM symbol, a total of P = 64 iterations, get
Figure C20041002186400154
As a result of time-domain channel estimation within this OFDM symbol.

LMS算法的实现信道追踪的方法为:The method of realizing the channel tracking of the LMS algorithm is as follows:

对所有的p=0,1,…,P-1,迭代计算:For all p=0, 1, ..., P-1, iterative calculation:

ee [[ nno ,, pp ]] == YY ii [[ nno ,, kk pp ]] -- ww pp TT [[ nno ]] hh ~~ ii [[ nno ,, pp ]]

hh ~~ ii [[ nno ,, pp ++ 11 ]] == hh ~~ ii [[ nno ,, pp ]] ++ 22 &mu;&mu; ww pp ** [[ nno ]] ee [[ nno ,, pp ]]

其中:μ为步长因子,选择依据为使算法收敛且具有较小的估计误差,本例中LMS算法步长因子取为μ=0.02。Among them: μ is the step size factor, the basis for selection is to make the algorithm converge and have a small estimation error. In this example, the step size factor of the LMS algorithm is taken as μ = 0.02.

步骤C:利用步骤B自适应信道估计的结果估计所有信道子载波频率响应Step C: Estimate the frequency response of all channel subcarriers using the result of adaptive channel estimation in step B

Figure C20041002186400157
拆分得到所有发射天线至接收天线i间的信道各主要多径分量的估计为
Figure C20041002186400158
其中
Figure C20041002186400159
为该主多径延迟时间的估计值(以采样周期为单位), l &LeftRightArrow; = 1,2 , &CenterDot; &CenterDot; &CenterDot; , L &LeftRightArrow; . 则经傅立叶变换后得到第n个OFDM符号期间第q发射天线至第i接收天线间信道第k子载波频域信道估计为:Depend on
Figure C20041002186400157
Split to obtain the estimation of each main multipath component of the channel between all transmitting antennas and receiving antenna i as
Figure C20041002186400158
in
Figure C20041002186400159
is the estimated value of the main multipath delay time (in sampling period), l &LeftRightArrow; = 1,2 , &Center Dot; &Center Dot; &CenterDot; , L &LeftRightArrow; . Then after Fourier transform, the frequency domain channel estimation of the kth subcarrier of the channel between the qth transmit antenna and the ith receive antenna during the nth OFDM symbol period is obtained for:

Figure C200410021864001512
Figure C200410021864001512

其中k表示子载波的序号(其取值范围为0≤k<N-1)。本例中频谱效率损失为 &xi; = N + ( D i - 1 ) P D i N , 即8.6%。图6给出了使用本发明方案的一个特例的误比特率(BER)计算机方针性能曲线,仿真结果显示,在较高信噪比时,本实施例中信噪比损失为1dB左右。Where k represents the serial number of the subcarrier (the value range is 0≤k<N-1). In this example, the spectral efficiency loss is &xi; = N + ( D. i - 1 ) P D. i N , That is 8.6%. Fig. 6 has provided and used the bit error rate (BER) computer target performance curve of a special example of the scheme of the present invention, and simulation result shows, when higher SNR, SNR loss is about 1dB in the present embodiment.

Claims (1)

1. the adaptive channel estimation method of a MIMO-OFDM system is characterized in that it comprises make a start step and receiving end step:
The described step of making a start is carried out according to step 1 and step 2:
Step 1: introduce training frequency guide symbol and reference pilot symbols to form data frame structure
To any transmitting antenna, transmit data according to a minute frame mode, every frame comprises N tIndividual OFDM symbol, N tBe positive integer; First OFDM symbol is a training OFDM symbol in the frame, and n=0, n represent OFDM symbol sequence number; N behind the accent white silk OFDM symbol in the frame t-1 OFDM symbol is a data transmission OFDM symbol, n=1, and 2 ..., N t-1; Insert training frequency guide symbol and insert reference pilot symbols according to following method:
Insert training frequency guide symbol: in training OFDM symbol, for q transmitting antenna, q=1,2 ..., M T, M TThe expression number of transmit antennas, M TBe positive integer, first subcarrier number that transmits training frequency guide symbol is q-1, thereafter every D F0Individual subcarrier inserts a training frequency guide symbol; To q transmitting antenna, at non-training frequency guide symbol subcarrier place, the transmission amplitude is 00 symbol;
Insert reference pilot symbols: in any data transmission OFDM symbol, all transmitting antennas use identical subcarrier to transmit reference pilot symbols; To any transmitting antenna, the subcarrier sequence number that transmits the 1st reference pilot symbols is k 0, 0≤k 0<D F1, D F1Be positive integer, thereafter every D F1Individual subcarrier inserts a reference pilot symbols; To any transmitting antenna, at non-reference pilot symbols subcarrier place, transmit the data symbol of Space Time Coding output: the subcarrier sequence number that transmits reference pilot symbols is k p, p=0,1 ..., P-1,
Figure C2004100218640002C1
N is the OFDM sub-carrier number;
Step 2: the frame structure that forms according to step 1 sends.
Described receiving end step was carried out according to following stage 1 and stage 2:
Stage 1: estimate channel at the training OFDM symbol place of every frame, be used for the initial parameter that adaptive channel is followed the trail of during the transfer of data in this frame, carry out in turn according to A, B and three steps of C to obtain:
Steps A: the estimated value of obtaining the frequency response of all training frequency guide symbol subcarriers of current channel with least-squares estimation;
Step B: to the result of steps A, use frequency domain->spatial transform estimates that all transmitting antennas are to channel impulse response between all reception antennas, if in the channel impulse response, being estimated as of l bar multipath component multipath fading coefficient between q transmitting antenna to the i reception antenna in n OFDM symbol 0≤l≤N-1, i=1,2 ..., M R, M RBe the reception antenna number;
Step C: obtain channel master multipath component time of delay and extract main multipath component fading coefficients:
Utilize result of calculation among the step B, use direct truncation method of channel impulse response or main multipath component back-and-forth method, to the channel between the extremely any reception antenna i of any transmitting antenna q, from its channel multi-path fading coefficients
Figure C2004100218640002C3
The middle selection
Figure C2004100218640002C4
Bar master multipath,
Figure C2004100218640002C5
Be positive integer, storage delay time l I, q, kFading coefficients with respective paths
Figure C2004100218640002C6
k = 1,2 , &CenterDot; &CenterDot; &CenterDot; , L &infin; ;
Stage 2: adaptive channel is followed the trail of during the transfer of data:
To n data transmission OFDM symbol in the frame, n=1,2 ..., N t-1, wherein arbitrary reception antenna i, use following channel estimation methods:
Steps A: adaptive channel track algorithm initial value is set:
Adaptive channel is set estimates that initial value is the final value that n-1 OFDM symbol estimated:
Figure C2004100218640003C1
n=1,2,…,N t-1 (1)
Wherein:
Figure C2004100218640003C2
Represent in n the OFDM symbol that through after p iteration, all transmitting antennas are to the estimated result of channel between the reception antenna i,
Figure C2004100218640003C3
Represent in n the OFDM symbol that through after P iteration, all transmitting antennas are to the final result of channel estimating between the reception antenna i,
Figure C2004100218640003C4
Represent in n the OFDM symbol that through after P iteration, q transmitting antenna is to the final result of channel estimating between the reception antenna i, wherein
Figure C2004100218640003C5
Step B: adaptive channel is followed the tracks of:
Arbitrarily the reference pilot symbols subcarrier is all corresponding to an iteration of adaptive algorithm, and to p arbitrarily, 0≤p≤P-1, P are the reference pilot symbols number, and the reference signal of adaptive algorithm is Y i[n, k p], the received signal at p reference pilot subcarrier of i reception antenna place during n OFDM symbol, k pFor transmitting the subcarrier sequence number of p reference pilot symbols;
The input signal vector w of adaptive algorithm p[n] is:
w p [ n ] = [ w 1 . p T [ n ] , w 2 . p T [ n ] , . . . , w M T . p T [ n ] ] T - - - ( 4 )
Wherein
Figure C2004100218640003C9
The value of representing p reference pilot symbols of q transmitting antenna in n the data transmission OFDM, receiving terminal is known should to be worth;
Adaptive channel is estimated: adopt the adjustment of LMS algorithm or RLS algorithm or QR-RLS algorithm
Figure C2004100218640003C10
In the current data transmission OFDM symbol, iteration is carried out P time altogether, and P is the number of reference pilot symbols in each OFDM symbol, obtains As the time domain channel results estimated in the OFDM symbol, promptly
Figure C2004100218640003C12
Step C: utilize step B adaptive channel results estimated to estimate the order of all channel sub-carrier frequencies responses, split according to the combination of (2) and (3) formula
Figure C2004100218640003C13
Obtain of the estimation of all transmitting antennas to each the main multipath component of channel between reception antenna i; The estimated value of remembering main multipath component between q transmitting antenna to the i reception antenna is
Figure C2004100218640004C1
Use to obtain during n the OFDM symbol channel k sub-carrier frequency domain channel estimating between q transmitting antenna to the i reception antenna behind the Fourier transform,
All reception antennas are repeated above-mentioned steps A, B and C.
CNB2004100218648A 2004-02-20 2004-02-20 An Adaptive Channel Estimation Method for MIMO-OFDM System Expired - Fee Related CN100385824C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2004100218648A CN100385824C (en) 2004-02-20 2004-02-20 An Adaptive Channel Estimation Method for MIMO-OFDM System

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2004100218648A CN100385824C (en) 2004-02-20 2004-02-20 An Adaptive Channel Estimation Method for MIMO-OFDM System

Publications (2)

Publication Number Publication Date
CN1658528A CN1658528A (en) 2005-08-24
CN100385824C true CN100385824C (en) 2008-04-30

Family

ID=35007813

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004100218648A Expired - Fee Related CN100385824C (en) 2004-02-20 2004-02-20 An Adaptive Channel Estimation Method for MIMO-OFDM System

Country Status (1)

Country Link
CN (1) CN100385824C (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101808053A (en) * 2010-02-26 2010-08-18 大唐联诚信息系统技术有限公司 Channel estimation method and device based on OFDM (Orthogonal Frequency Division Multiplexing)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007110456A (en) * 2005-10-14 2007-04-26 Hitachi Ltd Wireless communication device
CN100525275C (en) * 2005-10-27 2009-08-05 财团法人工业技术研究院 MIMO-OFDM system and precoding and feedback method thereof
DE602005010592D1 (en) * 2005-11-15 2008-12-04 Alcatel Lucent Method for transmitting channel quality information in a multi-carrier radio communication system and corresponding mobile station and base station
CN1988520B (en) * 2005-12-23 2010-09-08 中兴通讯股份有限公司 A Channel Estimation Method for Orthogonal Frequency Division Multiplexing System
CN101076001B (en) * 2006-05-15 2011-03-02 中兴通讯股份有限公司 Method for estimating channel based on orthogonal frequency division multiplexing system
CN101083646B (en) * 2006-06-01 2010-04-14 电子科技大学 A Channel Estimation Optimization Method for Clipped OFDM System
CN101155156B (en) * 2006-09-25 2011-06-08 华为技术有限公司 Channel estimation method and device and pilot frequency sequence generation method and device
CN101170531B (en) * 2006-10-24 2012-01-18 北京大学 A channel estimate method and corresponding communication method and system
CN101291312B (en) 2007-04-17 2012-01-11 电信科学技术研究院 Reference symbol processing method, system, base station and terminal
CN101383798B (en) * 2007-09-05 2012-09-26 中兴通讯股份有限公司 Channel estimation method and system based on regulating factor
CN101389050B (en) * 2007-09-14 2011-11-16 展讯通信(上海)有限公司 Information sending and receiving method using reference symbol configuration change in single frequency network
CN101222470B (en) * 2008-01-31 2010-07-14 上海交通大学 Channel Estimation Method for Dual-Antenna Discrete Fourier Spread Spectrum Generalized Multicarrier Systems
KR101689597B1 (en) 2009-02-09 2016-12-26 엘지전자 주식회사 Method for transmitting reference signals in downlink multiple input multiple output
CN101500013B (en) * 2009-03-27 2012-05-02 武汉烽火网络有限责任公司 Link Tracking Path Information Retrieval Method Based on IEEE802.1ag Standard
KR101663616B1 (en) * 2009-04-29 2016-10-07 엘지전자 주식회사 Method for generating reference signal sequence in wireless communication system and apparatus therefor
CN101605120B (en) * 2009-07-17 2011-09-14 清华大学 Method and device for decoding iterated internal code and evaluating channel for orthogonal frequency division multiplexing system
CN101621492A (en) * 2009-08-14 2010-01-06 中兴通讯股份有限公司 Resource determining method of special demodulation data reference signal
CN101848183B (en) * 2010-04-06 2013-02-06 中国人民解放军信息工程大学 Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System
CN103929380B (en) * 2013-01-10 2017-04-26 晨星软件研发(深圳)有限公司 Wireless receiving system and channel effect estimation method thereof
CN103414666B (en) * 2013-08-26 2016-06-08 电子科技大学 A kind of based on the two-dimensional adaptive OFDM channel estimation method being uniformly distributed pilot tone
WO2015154213A1 (en) * 2014-04-07 2015-10-15 Intel IP Corporation Client specific frequency offset estimation copyright notice
CN104967582B (en) * 2015-07-22 2018-02-27 浙江大学 Channel estimation methods based on odd-even alternation pilot frequency sequence in NAVDAT
CN110868367B (en) * 2018-08-27 2022-06-28 扬智科技股份有限公司 Channel estimation method and circuit
CN112260811B (en) * 2020-10-17 2022-12-06 西安交通大学深圳研究院 Pilot frequency distribution method of multi-input multi-output orthogonal frequency division multiplexing system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020181390A1 (en) * 2001-04-24 2002-12-05 Mody Apurva N. Estimating channel parameters in multi-input, multi-output (MIMO) systems
WO2003034642A2 (en) * 2001-10-17 2003-04-24 Nortel Networks Limited Synchronisation in multicarrier cdma systems
WO2003034644A1 (en) * 2001-10-17 2003-04-24 Nortel Networks Limited Scattered pilot pattern and channel estimation method for mimo-ofdm systems
WO2003049397A2 (en) * 2001-12-03 2003-06-12 Qualcomm Incorporated Iterative detection and decoding for a mimo-ofdm system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020181390A1 (en) * 2001-04-24 2002-12-05 Mody Apurva N. Estimating channel parameters in multi-input, multi-output (MIMO) systems
WO2003034642A2 (en) * 2001-10-17 2003-04-24 Nortel Networks Limited Synchronisation in multicarrier cdma systems
WO2003034644A1 (en) * 2001-10-17 2003-04-24 Nortel Networks Limited Scattered pilot pattern and channel estimation method for mimo-ofdm systems
WO2003049397A2 (en) * 2001-12-03 2003-06-12 Qualcomm Incorporated Iterative detection and decoding for a mimo-ofdm system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
BLIND CHANNEL IDENTIFICATION AND EQUALIZATIONIN OFDM-BASED MULTIANTENNA SYSTEMS. H.BOLCSKEI,ET AL.IEEE TRANS ON SIGNAL PROCESSING,Vol.50 No.1. 2002 *
PILOT-SYMBOL-ASSISTED CHANNEL ESTIMATIONFORSPACE-TIME CODED OFDM SYSTEMS. K.LEE AND D.B.WILLIAMOS.EURASIP JOURNAL ON APPLIED SIGNAL PROCESSING. 2002 *
SIMPLIFIED CHANNEL ESTIMATION FOR OFDMSYSTEMS WITH MULTIPLE TRANSMIT ANTENNAS. YE(GEOFFREY) LI.IEEE TRANSACTIONS ON WIRELESS COMM.,Vol.1 No.1. 2002 *
基于MIMO-OFDM系统的信道估计方法. 张高峰,何方白.现代通信理论与信号处理进展-2003年通信理论与信号处理年会论文集. 2003 *
多用户MIMO-OFDM系统的增强型信道估计算法. 西安电子科技大学ISN国家重点实验室,,董乐,张海林.2003年无线及移动通信委员会学术年会论文集. 2003 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101808053A (en) * 2010-02-26 2010-08-18 大唐联诚信息系统技术有限公司 Channel estimation method and device based on OFDM (Orthogonal Frequency Division Multiplexing)

Also Published As

Publication number Publication date
CN1658528A (en) 2005-08-24

Similar Documents

Publication Publication Date Title
CN100385824C (en) An Adaptive Channel Estimation Method for MIMO-OFDM System
CN101505290B (en) An Improved Frequency Offset Estimation Method for Wideband MIMO
CN101578829B (en) A channel estimation method
CN103685096B (en) A kind of MIMO-OFDM system channel estimation method based on optimal pilot
CN101557377B (en) Method, device and system for calculation of pre-filtering coefficient and interference suppression
CN101827060A (en) Adaptive modulation-demodulation method base on fractional order Fourier transform
CN115277335B (en) A Symbol Detection Method for Orthogonal Timing Multiplexing System
CN101360079A (en) Wavelet Domain Value Denoising Method Based on Maximum Likelihood Estimator of Wavelet Denoising Algorithm
CN1937598A (en) Channel estimation method in orthogonal frequency-division multiplexing system and channel estimation device
CN101355543A (en) Channel Estimation Method for MIMO-SCFDE System Based on Orthogonal Training Sequence
CN102045285A (en) Channel estimation method and device and communication system
CN109600327A (en) It is a kind of to interfere the channel estimation methods utilized based on imaginary part
CN110233688B (en) Orthogonal Space-Time Coding Transmission Method for Large-Scale Antennas in Beam Domain Based on Doppler Suppression
CN102025662B (en) Channel estimation method and device for MIMO (multiple input multiple output) OFDM (orthogonal frequency division multiplexing) system
Eldemiry et al. Overview of the orthogonal time-frequency space for high mobility communication systems
CN101267409A (en) A MIMO-OFDM Dual Selective Channel Tracking Method
CN101026606A (en) Orthogonal pilot frequency sequence design method
CN101848183A (en) Channel estimation method and device in multiple input multiple output OFDM (Orthogonal Frequency Division Multiplexing) system
CN101018219B (en) Space frequency signal processing method
CN101197796B (en) Wireless sensor network channel evaluation method based on SC-FDE and virtual multi-antenna
Ganesh et al. Channel estimation analysis in MIMO-OFDM wireless systems
CN102006250B (en) Turbo enhancement method for MIMO-SCFDE wireless communication receiver
US8369471B1 (en) Method and apparatus for improving channel estimation
CN105119856B (en) Single carrier frequency domain equalization method based on sub-block analysis
CN102487368B (en) Design method and realization device of Per-tone equalizer (PTEQ)

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080430

Termination date: 20110220