[go: up one dir, main page]

CN1921463B - Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System - Google Patents

Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System Download PDF

Info

Publication number
CN1921463B
CN1921463B CN200510093257A CN200510093257A CN1921463B CN 1921463 B CN1921463 B CN 1921463B CN 200510093257 A CN200510093257 A CN 200510093257A CN 200510093257 A CN200510093257 A CN 200510093257A CN 1921463 B CN1921463 B CN 1921463B
Authority
CN
China
Prior art keywords
mrow
channel
impulse response
channel estimation
division multiplexing
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
CN200510093257A
Other languages
Chinese (zh)
Other versions
CN1921463A (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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CN200510093257A priority Critical patent/CN1921463B/en
Publication of CN1921463A publication Critical patent/CN1921463A/en
Application granted granted Critical
Publication of CN1921463B publication Critical patent/CN1921463B/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

本发明涉及正交频分复用移动通信系统的信道估计方法和实现装置,包括:根据前缀、或者判决反馈信息、或上一个OFDM符号的冲激响应信息获得相应符号信道的冲激响应粗估计;产生构造矩阵;根据本地的参考信号和接收到的导频信号计算互相关矢量;通过递归等式获得信道冲激响应的细估计结果;判定根据递归等式获得的细估计结果是否收敛;获得所有载波数目的信道估计结果;从结果中选择有用的数据载波信道估计结果并输出。该方法适应快衰落信道和慢衰落信道,包括导频的密度小于相关带宽的情况,不需要判断信道的类型以及信道的主要传输路径选择,具有计算量小,要求存储空间小,不用进行门限选择,实现容易等特点。

Figure 200510093257

The present invention relates to a channel estimation method and implementation device of an orthogonal frequency division multiplexing mobile communication system, including: obtaining a rough estimation of the impulse response of the corresponding symbol channel according to the prefix, or decision feedback information, or the impulse response information of the last OFDM symbol ; Generate a construction matrix; Calculate the cross-correlation vector according to the local reference signal and the received pilot signal; Obtain the fine estimation result of the channel impulse response through the recursive equation; Determine whether the fine estimation result obtained according to the recursive equation converges; Obtain Channel estimation results of all carrier numbers; select and output useful data carrier channel estimation results from the results. This method is suitable for fast-fading channels and slow-fading channels, including the case where the density of the pilot frequency is smaller than the relevant bandwidth. It does not need to judge the type of the channel and the main transmission path selection of the channel. It has a small amount of calculation, requires a small storage space, and does not need to select the threshold. , easy to realize and so on.

Figure 200510093257

Description

Channel estimation method and implementing device for orthogonal frequency division multiplexing mobile communication system
Technical Field
The present invention relates to a channel estimation method in a wireless mobile communication system, and more particularly, to a channel estimation method and an implementation apparatus in an orthogonal frequency division multiplexing mobile communication system, and more particularly, to a channel estimation method and an implementation apparatus in combination with a prefix and a pilot.
Background
Orthogonal Frequency Division Multiplexing (OFDM), a multi-carrier transmission mode, greatly reduces the sensitivity of the system to the frequency selectivity of a multipath fading channel by converting a high-speed transmission data stream into a set of low-speed parallel transmission data streams, and the introduction of a cyclic prefix further enhances the capability of the system to resist inter-symbol interference (ISI). In addition, the characteristics of high bandwidth utilization rate, simple implementation and the like enable OFDM to be applied more and more widely in the field of wireless communication, and channel estimation, which is a key technology of an OFDM system, also becomes a hot point of research in the field of communication and receives more and more attention.
There are two types of training sequence insertion forms for OFDM systems, block pilot and comb pilot. Most of the existing channel estimation methods are based on some pilot arrangement (block pilot or comb pilot).
These two types of channel estimation methods are each characterized: the channel estimation method based on the block pilot frequency is not very sensitive to the frequency selectivity of the channel, but is only suitable for the wireless channel with slow fading; the channel estimation method based on the comb-shaped pilot frequency can be suitable for a wireless channel with fast fading, but is sensitive to the frequency selectivity of the channel. However, combining the advantages of the two channel estimations, the method for estimating the channel is not common.
For an 802.16 system based on the OFDMA standard, it has both block pilots (hereinafter referred to as prefixes) and comb pilots (hereinafter referred to directly as pilots). If only the prefix is used for channel estimation, because the prefix symbol is only sent once in one frame, when the channel is a fast fading channel, the change of the channel can not be accurately tracked only by using the prefix, thereby affecting the demodulation performance of the system. If the channel estimation is performed only by using the pilot frequency, when the frequency selectivity of the channel is significant, the accuracy of the pilot frequency-based channel estimation method is also reduced, and even a "floor effect" is generated, that is, the bit error rate of the system is not increased along with the increase of the signal-to-noise ratio.
Taking an 802.16 system based on the OFDMA standard as an example, in the FUSC mode, the pilot carrier spacing is 66kHz, and for a channel whose channel model is SUI-5 (modified university of stanford temporary channel model-5), the coherence bandwidth of the channel is 70 kHz; for a channel with a channel model SUI-6 (modified Stanford university temporary channel model-6), the coherence bandwidth of the channel is 38 kHz. The pilot carrier spacing is close to the coherence bandwidth of the SUI-5 channel, and is greater than the coherence bandwidth of the SUI-6 channel, which is very prone to floor effects. Therefore, it is necessary to combine the prefix and the pilot signal to estimate the wireless channel.
The publication No. CN1505293A, publication date 6/16 in 2004, "channel estimation method and estimator in orthogonal frequency division multiplexing system", proposes a method for channel estimation by comprehensively using prefixes and pilots, which can be applied to fast fading channels, and also can be applied to channels with significant channel frequency selectivity, including the case where the pilot carrier spacing is smaller than the correlation bandwidth.
If the channel type selection method is not proper, the channel type cannot be accurately judged, the performance of the method is drastically deteriorated, and the patent does not give a judgment method of the channel type; for the main propagation path selection, the patent uses prefix impulse response to select the main path, wherein the threshold is selected, and the channel estimation performance is sensitive to the selection of the threshold, and simulation shows that if the threshold is not selected properly, the channel estimation performance is reduced sharply.
Disclosure of Invention
One of the objectives of the present invention is to provide a channel estimation method for an ofdm mobile communication system by combining a prefix and a pilot, which has excellent performance and is easy to implement, and can be adapted to various channel types, such as a frequency selective channel, a fast fading channel, and the like.
The invention also provides a device for realizing the method.
In order to achieve the purpose of the invention, the invention provides a method and a device for channel estimation by combining a prefix and a pilot frequency. Generating a pilot signal used for channel estimation by a system transmitter, wherein the pilot signal comprises a pilot insertion mode and two parts of content selected by a pilot value; the generated pilot signal is transformed into a time domain pilot signal through inverse Fourier transform; transmitting the time domain pilot signal on a wireless multi-path fading channel; and the system receiver performs Fourier transform on the time domain pilot signal to obtain a frequency domain pilot signal.
The invention is realized in such a way that:
a channel estimation method of an orthogonal frequency division multiplexing mobile communication system, comprising the processes of:
the first step, the impulse response rough estimation of the corresponding symbol channel is obtained according to the prefix, or the decision feedback information, or the impulse response information of the last OFDM symbol;
a second step of generating a construction matrix according to a local reference pilot signal, wherein the construction matrix is a conjugate transpose square matrix, the dimension of the construction matrix is the length of the rough estimation of the channel impulse response, and the element t (p) of the construction matrix is calculated according to the formula (1);
<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><msup><mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>
wherein k represents an offset of the carrier; p represents a set of pilot carrier offsets; a is a convergence factor, and satisfies 0 < a < 1; x (k) denotes a local reference pilot signal;
Figure G2005100932577D00032
is an inverse fourier transform factor;
a third step of calculating a cross-correlation vector from the local reference pilot signal and the received pilot signal, the length of the cross-correlation vector being equal to the dimension of the construction matrix, the elements r (p) of the cross-correlation vector being calculated according to equation (2);
<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
wherein y (k) represents a received pilot signal; a is a convergence factor, and satisfies 0 < a < 1;
a fourth step of obtaining a fine estimation result of the channel impulse response by a recursive equation according to the construction matrix and the cross-correlation vector, wherein the recursive equation is as follows: hm,n+1=Hm,n+R-THm,nIn the formula, Hm,nRepresents the fine estimation result of the nth channel impulse response of the mth orthogonal frequency division multiplexing symbol in a certain frame, wherein n is 0, 1, 2, …;
a fifth step of judging whether the fine estimation result obtained according to the recursion equation is converged, if so, outputting the fine estimation result of the impulse response, and if not, continuing to calculate the recursion equation;
sixthly, performing fast Fourier transform of all the carrier numbers of the system on the fine estimation result of the output impulse response to obtain a channel estimation result of all the carrier numbers;
and a seventh step of selecting and outputting useful data carrier channel estimation results from the channel estimation results of all the carriers in the system.
And the fifth step of outputting an impulse response fine estimation result by judging whether the mean square error of the results of two continuous recursive operations is less than a fixed constant or not, and if not, taking the result of the recursive operation as the input of the next recursive operation and continuing the recursive operation.
And the fifth step fixes recursion times, outputs a fine estimation result of impulse response after carrying out recursion calculation of corresponding times, and takes the fine estimation result as the impulse response information of the next orthogonal frequency division multiplexing symbol.
The length of the fine estimation of the channel impulse response is equal to the dimension of the construction matrix.
An apparatus for implementing a channel estimation method of an orthogonal frequency division multiplexing mobile communication system, comprising:
the channel impulse response rough estimator obtains the impulse response rough estimation of the channel according to the prefix, the decision feedback or the impulse response information of the last orthogonal frequency division multiplexing symbol;
a construction matrix generator, which generates a construction matrix according to the local reference pilot signal, wherein the construction matrix is a conjugate transpose square matrix, the dimension of the construction matrix is the length of the rough estimation of the channel impulse response, and the element t (p) of the construction matrix is calculated according to the formula (1);
<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><msup><mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>
wherein k represents an offset of the carrier; p represents a set of pilot carrier offsets; a is a convergence factor, and satisfies 0 < a < 1; x (k) denotes a local reference pilot signal;is an inverse fourier transform factor;
a cross-correlation vector generator for generating a cross-correlation vector by calculation based on the local reference pilot signal and the received pilot signal, wherein the length of the cross-correlation vector is equal to the dimension of the construction matrix, and the element r (p) of the cross-correlation vector is calculated according to the formula (2);
<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
wherein y (k) represents a received pilot signal; a is a convergence factor, and satisfies 0 < a < 1;
a data selector that selects either a coarse estimate of the channel impulse response or the result of a recursive equation as one input to the multiplier, the recursive equation being: hm,n+1=Hm,n+R-THm,nIn the formula, Hm,nRepresents the fine estimation result of the channel impulse response of the nth time of the mth orthogonal frequency division multiplexing symbol in a certain frame, wherein n is 0, 1, 2,…;
A multiplier for calculating and outputting the product of the construction matrix and the output of the data selector;
an adder for calculating and outputting the sum of outputs of the cross-correlation vector generator, the multiplier, and the data selector;
a termination condition decider for judging whether the recursion equation is converged, if so, terminating the recursion equation calculation and outputting channel impulse response, otherwise, continuing the recursion equation calculation;
a fast Fourier transformer for performing fast Fourier transform of the number of system carriers on the output of the termination condition decider;
and the channel estimation result selector selects useful data carrier channel estimation from the channel estimation results of all the carriers and outputs the useful data carrier channel estimation.
And the termination condition decision device judges whether the mean square error of the results of two continuous recursive calculations is smaller than a fixed constant or not, if so, outputs an impulse response fine estimation result, otherwise, takes the result of the recursive calculation as the input of the next recursive calculation and continues the recursive calculation.
The termination condition decider fixes recursion times, outputs a fine estimation result of impulse response after carrying out recursion calculation of corresponding times, and takes the fine estimation result as the impulse response information of the next orthogonal frequency division multiplexing symbol.
Compared with the existing channel estimation method, the method is suitable for fast fading channels and slow fading channels, including the condition that the density of pilot frequency is less than the related bandwidth, does not need to judge the type of the channel and the main transmission path selection of the channel, and has the characteristics of small calculated amount, small required storage space, no need of threshold selection, easy realization and the like.
Drawings
FIG. 1 is a schematic diagram of a receiver architecture incorporating the present invention;
fig. 2 is a schematic diagram of an internal implementation of the channel estimator proposed by the present invention;
fig. 3 is a flow chart corresponding to the internal implementation schematic diagram of the channel estimator proposed by the present invention.
Detailed Description
The invention is described in detail below with reference to the attached drawing figures:
in order to achieve the first object, the method provided by the invention mainly comprises the following steps:
1. obtaining the impulse response rough estimation H of the channel of the mth symbol according to the prefix or the decision feedback information or the impulse response information of the last OFDM symbolm,0
2. A construction matrix T is generated according to equation (1), the dimension of T being equal to the length of the coarse estimate of the impulse response. The construction matrix T is a matrix satisfying T ═ THSquare matrix of relationships, here THRepresents the conjugate transpose of the construction matrix T, which, when the dimension of the construction matrix T is L, is of the form:
T = t ( 0 ) , t * ( 1 ) , . . . t * ( L - 1 ) t ( 1 ) , t ( 0 ) , . . . t * ( L - 2 ) . . . . . . . . . t ( L - 1 ) , t ( L - 2 ) , . . . t ( 0 ) - - - ( 1 )
wherein,
<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><msup><mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
k represents an offset of the carrier;
p represents a set of pilot carrier offsets;
a is a convergence factor, and satisfies 0 < a < 1;
x (k) denotes a local reference pilot signal;
Figure G2005100932577D00061
is an inverse fourier transform factor.
3. The cross-correlation vector R is calculated, the length of R being equal to the dimension L of the matrix T, the cross-correlation vector R being of the form:
R=[r(0),......r(L-1)]T; (3)
wherein,
<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>
y (k) denotes a received pilot signal;
a has the same meaning as above.
4. And obtaining a fine estimation result of the channel impulse response according to the formula (5).
Hm,n+1=Hm,n+R-THm,n (5)
Wherein,
Hm,ndenotes the fine estimation result of the nth channel impulse response of the mth OFDM symbol in a certain frame, where n is 0, 1, 2, …. When n is 0, Hm,0Is a coarse estimate of the channel impulse response, H1,0=H0,0,H0,0Is the channel impulse response estimate of the prefix. Hm,nIs a vector of length equal toThe dimensions of the matrix T are constructed.
5. To Hm,n+1And performing N-point FFT to obtain the channel estimation results of all N carriers, wherein N refers to the number of all the carriers of the system.
6. And selecting useful data carrier channel estimation results from the channel estimation results of the N carriers and outputting the useful data carrier channel estimation results.
The channel estimation device corresponding to the above steps comprises the following modules:
1. and a cross-correlation vector generator for generating a cross-correlation vector R represented by equation (3).
2. And the construction matrix generator is used for generating a construction matrix T according to the channel impulse response coarse estimator, wherein the dimension of the matrix T is L.
3. And the channel impulse response rough estimator calculates the impulse response rough estimation of the channel according to the prefix, the decision feedback or the impulse response information of the previous OFDM symbol.
4. A data selector operative to select either the coarse estimate of the channel impulse response or the result of recursive equation (5) as one input to the multiplier.
5. And the multiplier is used for calculating and outputting the product of the construction matrix T and the output of the data selector.
6. And an adder for calculating and outputting the sum of the outputs of the cross-correlation vector generator, the multiplier, and the data selector.
7. A termination condition determiner for determining whether the recursive equation (5) has converged, terminating the recursive calculation if it has converged, and outputting a channel impulse response Hm,n+1Otherwise, the recursive computation is continued.
8. FFT converter for outputting H to termination condition deciderm,n+1An N-point FFT is performed, N referring to the number of all carriers of the system.
9. And a channel estimation result selector for selecting useful data carrier channel estimation from the channel estimation results of the N carriers and outputting the useful data carrier channel estimation.
In the receiver in fig. 1, a front-end processing module 101 receives a radio frequency signal through a receiving antenna, then performs down-conversion, analog-to-digital conversion, cyclic prefix removal and the like on the received signal, and finally, data is output to an FFT module 102.
The FFT block 102 performs a fast fourier transform on the input data so that the data is transformed from the time domain to the frequency domain. The FFT module 102 outputs data from the data carrier to the channel compensation module 104 and data from the pilot carrier to the channel estimation module 103.
The channel compensation module 104 receives the user data output by the FFT module 102 and the data carrier channel estimation output by the channel estimation module 103, and performs amplitude and phase compensation on the user data using the channel estimation of the data carrier.
Finally, the data is output to the demodulation and decoding module 105 for subsequent processing.
Fig. 2 is a schematic diagram of an internal implementation of the channel estimator proposed by the present invention.
The channel impulse response rough estimator 201 obtains the impulse response rough estimation of the channel according to the prefix, the decision feedback or the impulse response information of the previous OFDM symbol.
The construction matrix generator 202 generates a construction matrix T of a corresponding dimension according to equation (1) based on the dimension of the matrix determined by the channel impulse response coarse estimator.
The cross-correlation vector generator 203 generates a cross-correlation vector R shown in equation (3) based on the local reference pilot signal and the pilot signal at the receiving end.
The data selector 204 is used to select the initial value of each symbol for each recursion, and if it is the first symbol, the first initial value H1,0Equal to prefix impulse response H0,0Else, each time the initial value Hm,0Are all the lastImpulse response H of OFDM symbolm-1,n
The multiplier 205 calculates and outputs the product of the outputs of the construction matrix generator and the data selector.
The adder 206 calculates the sum of the outputs of the cross-correlation vector generator, the multiplier, and the data selector, and outputs the sum.
The termination condition decider 207 is for deciding the end of the recursive computation.
The fast fourier transformer 208 is a channel estimate H for the output of the termination condition deciderm,n+1And taking N-point Fourier transform, wherein N refers to the number of all carriers of the system.
The channel estimation result selector 209 selects and outputs channel estimation of a useful data carrier from the channel estimation results of the N carriers.
Fig. 3 corresponds to a flow chart of a channel estimation implementation.
Step 301 calculates the rough estimation of the impulse response of the channel according to the information such as the prefix, the decision feedback or the impulse response of the previous OFDM symbol, and the rough estimation is equivalent to the initial value H calculated in each recursionm,0If it is the first OFDM symbol of a frame, the initial value H of the recursive computation1,0Selected prefix impulse response estimate H0,0Otherwise, selecting the impulse response estimation H of the last symbolm-1,n. And determining the dimension of the matrix T according to the length of the channel impulse response rough estimation.
Step 302 determines the dimension of matrix T according to the pilot frequency, prefix or channel impulse response information of the previous OFDM symbol, and generates a constructed matrix T according to equation (1). Taking the OFDMA system based on the 802.16 protocol as an example, the length of the rough estimation of the channel impulse response can be obtained from the prefix information, and the dimension of the construction matrix T is equal to this length. The present invention does not limit the method of determining the dimension of the matrix T using the pilot, the prefix, or the channel impulse response information of the last OFDM symbol. Simulation shows that the method is insensitive to the accuracy of the channel impulse response length.
Step 303 generates a cross-correlation vector R according to the cross-correlation R (p) between the local reference pilot signal and the pilot signal at the receiving end, as shown in equation (3), where the length of the vector R is equal to the dimension of the matrix T.
Step 304 computes recursive equation (5) from the construction matrix T and the cross-correlation vector R.
Step 305 determines whether the recursive computation is finished according to the result of each recursive computation, and the termination condition may be determined in two ways:
one is to determine whether the average difference of the results of two consecutive recursion operations is less than a fixed constant, such as 0.01. Is, an impulse response H is outputm,n+1Step 306 is to get the last recursion result Hm,n+1As input H for the next recursionm,nThe recursive computation is continued.
The other is to fix the recursion number n, recursively calculate equation (5) n times, and then output the impulse response Hm,n+1
Step 307 is Hm,n+1As the initial value H for the recursive calculation of the next OFDM symbolm+1,0
Step 308 is the impulse response H to the channelm,n+1An N-point FFT is performed, resulting in channel estimation for all N carriers, where N refers to the number of all carriers in the system.
Step 309 is to select useful data carrier channel estimates from the channel estimates of N carriers and output them, where N refers to the number of all carriers in the system.
The foregoing is a description of one embodiment of the present invention and those skilled in the art will appreciate that various modifications and changes may be made to the embodiment of the present invention without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (7)

1. A channel estimation method of an orthogonal frequency division multiplexing mobile communication system, comprising the processes of:
the first step, the impulse response rough estimation of the corresponding symbol channel is obtained according to the prefix, or the decision feedback information, or the impulse response information of the last OFDM symbol;
a second step of generating a construction matrix according to a local reference pilot signal, wherein the construction matrix is a conjugate transpose square matrix, the dimension of the construction matrix is the length of the rough estimation of the channel impulse response, and the element t (p) of the construction matrix is calculated according to the formula (1);
<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><msup><mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>
wherein k represents an offset of the carrier; p represents a set of pilot carrier offsets; a is a convergence factor, and satisfies 0 < a < 1; x (k) denotes a local reference pilot signal;is an inverse fourier transform factor;
a third step of calculating a cross-correlation vector from the local reference pilot signal and the received pilot signal, the length of the cross-correlation vector being equal to the dimension of the construction matrix, the elements r (p) of the cross-correlation vector being calculated according to equation (2);
<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
wherein y (k) represents a received pilot signal; a is a convergence factor, and satisfies 0 < a < 1;
a fourth step of obtaining a fine estimation result of the channel impulse response through a recursion equation according to the construction matrix and the cross-correlation vector, wherein the recursion equation is; hm,n+1=Hm,n+R-THm,nIn the formula, Hm,nRepresents the fine estimation result of the nth channel impulse response of the mth orthogonal frequency division multiplexing symbol in a certain frame, wherein n is 0, 1, 2, …;
a fifth step of judging whether the fine estimation result obtained according to the recursion equation is converged, if so, outputting the fine estimation result of the impulse response, and if not, continuing to calculate the recursion equation;
sixthly, performing fast Fourier transform of all the carrier numbers of the system on the fine estimation result of the output impulse response to obtain a channel estimation result of all the carrier numbers;
and a seventh step of selecting and outputting useful data carrier channel estimation results from the channel estimation results of all the carriers in the system.
2. The channel estimation method of the orthogonal frequency division multiplexing mobile communication system as claimed in claim 1, wherein:
and in the fifth step, whether the mean square error of the result calculated by the two recursion equations is smaller than a fixed constant or not is judged, if so, a fine estimation result of impulse response is output, otherwise, the result of the recursion calculation is used as the input of the next recursion calculation, and the recursion calculation is continued.
3. The channel estimation method of the orthogonal frequency division multiplexing mobile communication system as claimed in claim 1, wherein:
and the fifth step fixes recursion times, outputs a fine estimation result of impulse response after carrying out recursion calculation of corresponding times, and takes the fine estimation result as the impulse response information of the next orthogonal frequency division multiplexing symbol.
4. The channel estimation method of the orthogonal frequency division multiplexing mobile communication system as claimed in claim 1, 2 or 3, wherein:
the length of the fine estimation of the channel impulse response is equal to the dimension of the construction matrix.
5. An apparatus for implementing a channel estimation method of an orthogonal frequency division multiplexing mobile communication system, comprising:
the channel impulse response rough estimator obtains the impulse response rough estimation of the channel according to the prefix, the decision feedback or the impulse response information of the last orthogonal frequency division multiplexing symbol;
a construction matrix generator, which generates a construction matrix according to the local reference pilot signal, wherein the construction matrix is a conjugate transpose square matrix, the dimension of the construction matrix is the length of the rough estimation of the channel impulse response, and the element t (p) of the construction matrix is calculated according to the formula (1);
<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><msup><mrow><mo>|</mo><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>
wherein k represents an offset of the carrier; p represents a set of pilot carrier offsets; a is a convergence factor, and satisfies 0 < a < 1; x (k) denotes a local reference pilot signal;is an inverse fourier transform factor;
a cross-correlation vector generator for generating a cross-correlation vector by calculation according to a local reference pilot signal and a received pilot signal, wherein the length of the cross-correlation vector is equal to the dimension of the construction matrix, and the element r (p) of the cross-correlation vector is calculated according to the formula (2);
<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Subset;</mo><mi>P</mi></mrow></munder><mi>Y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mrow><mo>-</mo><mi>kp</mi></mrow></msubsup><mo>,</mo><mi>P</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
wherein y (k) represents a received pilot signal; a is a convergence factor, and satisfies 0 < a < 1;
a data selector that selects either a coarse estimate of the channel impulse response or the result of a recursive equation as one input to the multiplier, the recursive equation being: hm,n+1=Hm,n+R-THm,nIn the formula, Hm,nRepresents the fine estimation result of the nth channel impulse response of the mth orthogonal frequency division multiplexing symbol in a certain frame, wherein n is 0, 1, 2, …;
a multiplier for calculating and outputting the product of the construction matrix and the output of the data selector;
an adder for calculating and outputting the sum of outputs of the cross-correlation vector generator, the multiplier and the data selector;
a termination condition decider for judging whether the recursion equation is converged, if so, terminating the recursion equation calculation and outputting channel impulse response, otherwise, continuing the recursion equation calculation;
a fast Fourier transformer for performing fast Fourier transform of the number of system carriers on the output of the termination condition decider;
and the channel estimation result selector selects useful data carrier channel estimation from the channel estimation results of all the carriers and outputs the useful data carrier channel estimation.
6. The apparatus for implementing a channel estimation method of an orthogonal frequency division multiplexing mobile communication system as claimed in claim 5, wherein:
and the termination condition decision device judges whether the mean square error of the results of two continuous recursive calculations is smaller than a fixed constant or not, if so, outputs an impulse response fine estimation result, otherwise, takes the result of the recursive calculation as the input of the next recursive calculation and continues the recursive calculation.
7. The apparatus for implementing a channel estimation method of an orthogonal frequency division multiplexing mobile communication system as claimed in claim 5, wherein:
the termination condition decider fixes recursion times, outputs a fine estimation result of impulse response after carrying out recursion calculation of corresponding times, and takes the fine estimation result as the impulse response information of the next orthogonal frequency division multiplexing symbol.
CN200510093257A 2005-08-23 2005-08-23 Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System Expired - Fee Related CN1921463B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200510093257A CN1921463B (en) 2005-08-23 2005-08-23 Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200510093257A CN1921463B (en) 2005-08-23 2005-08-23 Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System

Publications (2)

Publication Number Publication Date
CN1921463A CN1921463A (en) 2007-02-28
CN1921463B true CN1921463B (en) 2010-05-05

Family

ID=37779038

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200510093257A Expired - Fee Related CN1921463B (en) 2005-08-23 2005-08-23 Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System

Country Status (1)

Country Link
CN (1) CN1921463B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101257471B (en) * 2008-03-31 2012-03-28 浙江大学 Mobile multimedia digital broadcast OFDM system channel estimation method
US8649456B2 (en) 2009-03-12 2014-02-11 Futurewei Technologies, Inc. System and method for channel information feedback in a wireless communications system
US8675627B2 (en) 2009-03-23 2014-03-18 Futurewei Technologies, Inc. Adaptive precoding codebooks for wireless communications
EP2451202B1 (en) 2009-06-29 2017-08-23 Huawei Technologies Co., Ltd. Method, device and system for reference signal processing
JP5032538B2 (en) * 2009-08-28 2012-09-26 株式会社東芝 Transmission line response estimator
CN102036264B (en) 2009-09-30 2013-12-04 华为技术有限公司 Method and device for measuring channel of cell
CN102056220B (en) 2009-10-28 2014-02-19 华为技术有限公司 Method and device for realizing channel measurement
CN101815042B (en) * 2010-04-13 2012-12-05 新邮通信设备有限公司 Orthogonal frequency division multiplexing (OFDM) system channel estimation method and device
CN101895311B (en) * 2010-06-18 2013-08-07 华为技术有限公司 Signal-receiving method and receiver
CN101951353B (en) * 2010-09-30 2013-02-13 电子科技大学 Channel estimation method for orthogonal frequency division multiplexing (OFDM) system under interference environment
WO2015143598A1 (en) * 2014-03-24 2015-10-01 华为技术有限公司 Method and apparatus for channel estimation processing

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1437338A (en) * 2003-03-08 2003-08-20 华中科技大学 Channel estimation method for orthogonal frequency-division multiplexing communication system
WO2004040813A1 (en) * 2002-10-29 2004-05-13 Qualcomm Incorporated Channel estimation for ofdm communication systems
CN1505293A (en) * 2002-12-02 2004-06-16 北京三星通信技术研究有限公司 Wireless channel estimation method and evaluation device for orthogonal frequency division multiplexing communication system
CN1514557A (en) * 2002-12-31 2004-07-21 上海贝尔阿尔卡特股份有限公司 Method and device of equalized fast fading channel in orthogonal frequency dirision multiplex system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004040813A1 (en) * 2002-10-29 2004-05-13 Qualcomm Incorporated Channel estimation for ofdm communication systems
CN1505293A (en) * 2002-12-02 2004-06-16 北京三星通信技术研究有限公司 Wireless channel estimation method and evaluation device for orthogonal frequency division multiplexing communication system
CN1514557A (en) * 2002-12-31 2004-07-21 上海贝尔阿尔卡特股份有限公司 Method and device of equalized fast fading channel in orthogonal frequency dirision multiplex system
CN1437338A (en) * 2003-03-08 2003-08-20 华中科技大学 Channel estimation method for orthogonal frequency-division multiplexing communication system

Also Published As

Publication number Publication date
CN1921463A (en) 2007-02-28

Similar Documents

Publication Publication Date Title
US7693036B2 (en) Decision feedback channel estimation and pilot tracking for OFDM systems
US7327812B2 (en) Apparatus and method for estimating a plurality of channels
US8064328B2 (en) Channel estimation device
CN101981879B (en) Broadband pilot channel estimation using a reduced order FFT and a hardware interpolator
CN103107969B (en) Incremental iterative time-varying channel evaluation and inter carrier interference (ICI) elimination method of fast orthogonal frequency division multiplexing (OFDM) system
CN103051578B (en) With the iteration error propagation judgement OFDM channel estimation method that ICI eliminates
JP2009532957A (en) Channel estimation for rapidly dispersive fading channels
US7218693B2 (en) Method and apparatus for channel estimation in a packet based transmission system having reduced complexity
CN101577692A (en) Channel estimating method of orthogonal frequency division multiplexing system and device thereof
CN103873422A (en) Method for eliminating multipath interference in system symbol in underwater sound orthogonal frequency-division multiplexing system
CN1921463B (en) Channel Estimation Method and Implementation Device for Orthogonal Frequency Division Multiplexing Mobile Communication System
CN1738300B (en) Method for estimating maximum likelihood frequency offset in mobile communication system
CN1937598A (en) Channel estimation method in orthogonal frequency-division multiplexing system and channel estimation device
CN101309243A (en) A New Parametric Channel Estimator for OFDM
CN101064571B (en) Apparatus and method for enhanced channel estimation in OFDM receiver
CN105656816A (en) Channel estimation method for comb-shaped pilot frequency OFDM system
US7619964B2 (en) High doppler channel estimation for OFD multiple antenna systems
CN101075829B (en) A Channel Estimation Method Applicable to Orthogonal Frequency Division Multiplexing System
CN104717168B (en) Orthogonal frequency division multiplexing (OFDM) ultra wide band system anti-multipath regular synchronization scheme
JP2007037151A (en) Apparatus, method, and computer program for estimating channel of OFDM transmission system
CN1988520B (en) A Channel Estimation Method for Orthogonal Frequency Division Multiplexing System
CN1909528B (en) Channel estimation method and device for crossing frequency division multiplexing system
CN113055318B (en) Channel estimation method
US20080219360A1 (en) System and method for performing channel estimation using interpolation
RU2298286C1 (en) Method for evaluating channel in multi-frequency radio communication systems with several transmitting and receiving antennas

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

Termination date: 20200823