This application claims the benefit of Taiwan application Serial No. 100130179 filed Aug. 23, 2011, the disclosure of which is incorporated by reference herein in its entirety.
1. Technical Field
The disclosed embodiments relate in general to a channel estimation method with channel adaptation and a device using the same.
2. Description of the Related Art
Of the existing technologies, the orthogonal frequency division multiplexing (OFDM) technology is frequently used in various wireless communication protocols (such as worldwide interoperability for microwave access (WiMAX) or the 3rd generation partnership project-long term evolution (3GPP-LTE)) for the data transmission of physical layer. In general, the OFDM receiver executes an operation of channel estimation according to the pilot channel gains by using the interpolation technology.
However, in the occasions of high-speed movement, the Doppler effect will destroy the orthogonal properties between the sub-carriers used in the OFDM technology, and may further incur inter-carrier interference (ICI). In addition, the occasions of high-speed movement radically change the power delay profile (PDP) of the channel, so that the results obtain by the channel estimation method according to the interpolation technology are no more reliable. Consequently, in occasions of high-speed movement, the existing OFDM system may face the problem of deterioration in the efficiency of the entire system.
The disclosure is directed to a channel estimation method with channel adaptation and a device using the same. An input symbol is received, and an operation of channel estimation is executed by a Wiener filter according to the pilot channel gains corresponding to the pilot signal of the input symbol. According to the channel estimation method with channel adaptation and the device using the same related to the disclosure, an adaptive alteration for the first and second multi-path statistical characteristic parameters of a Wiener filter is further calculated according to the above channel estimation result and the pilot channel gains, and an updated Wiener filter is generated accordingly. An operation of channel estimation is executed on a next input symbol according to the updated Wiener filter, so that the Wiener filter used for channel estimation can be adjusted according to the input signal and existing channel conditions and adapted to rapid change in channel characteristics in the occasions of high-speed movement. An adaptive adjustment is executed on the operation of channel estimation in response to rapid change in channel characteristics in the occasions of high-speed movement to achieve superior efficiency of the entire system.
The disclosure provides a channel estimation method with channel adaptation for receiving an input signal having several input symbols each having several pilot signals. The channel estimation method includes the following steps. Firstly, an n-th input symbol of an input signal is received, and several pilot channel gains are obtained by calculating several pilot channels transmitting the pilot signals of the n-th input symbol, wherein n is a natural number. Next, an operation of interpolation is executed on the pilot channel gains by a Wiener filter to obtain several data channel gains, corresponding to several data channels transmitting several respective data signals of the n-th input symbol. Then, an adaptive alteration for the first and second multi-path statistical characteristic parameters is calculated according to the data channel gains and the pilot channel gains, and the first and the second multi-path statistical characteristic parameters are adjusted accordingly. Then, an updated Wiener filter is generated according to the adjusted first and second multi-path statistical characteristic parameters. Afterwards, following the step of generating the updated Wiener filter, the value of n is increased by 1 and the step of channel estimation is repeated so that an operation of channel estimation is executed on a next input symbol of the input signal.
The disclosure discloses a channel estimation device with channel adaptation for receiving an input signal having several input symbols each having several pilot signals. The channel estimation device includes a reception unit, a channel estimation unit and a channel adaptation parameter tracking unit. The reception unit receives an n-th input symbol of an input signal, wherein n is a natural number. The channel estimation unit obtains several pilot channel gains by calculating several pilot channels transmitting the pilot signals of the n-th input symbol. The channel estimation unit further includes a Wiener filter for executing an operation of interpolation on the pilot channel gains to obtain several data channel gains, corresponding to a plurality of data channels transmitting a plurality of respective data signals of the n-th input symbol. The channel adaptation parameter tracking unit calculates an adaptive alteration for the first and second multi-path statistical characteristic parameters according to the data channel gains and the pilot channel gains, and accordingly adjusts the first and the second multi-path statistical characteristic parameters. The channel estimation unit generates an updated Wiener filter according to the adjusted first and second multi-path statistical characteristic parameters to execute an operation of channel estimation on the (n+1)-th input symbol.
The above and other aspects of the disclosure will become better understood with regard to the following detailed description of the non-limiting embodiment(s). The following description is made with reference to the accompanying drawings.
In the following detailed description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the disclosed embodiments. It will be apparent, however, that one or more embodiments may be practiced without these specific details. In other instances, well-known structures and devices are schematically shown in order to simplify the drawing.
The channel estimation device with channel adaptation of the present embodiment calculates an adaptive alteration for multi-path statistical characteristic parameters of a Wiener filter according to the channel estimation result and the pilot channel gains, and accordingly generates an updated Wiener filter for executing an operation of channel estimation on a next input symbol of the input signal.
Referring to
The operation of channel estimation executed by the channel estimation device 1 according to each input symbol of the input signal Si and obtains each restored data through restoration is substantially the same for all input symbols, and such operation is exemplified below by the operation of executing channel estimation by the channel estimation device 1 according to the n-th input symbol Yn and the operation of obtaining a corresponding restored data Xn, wherein n is a natural number.
The channel estimation device 1 includes a reception unit 102, a channel estimation unit 104, a channel adaptation parameter tracking unit 106, an equalizing/quantifying unit 108 and a demapper unit 110. The reception unit 102 receives a wireless signal Sr and provides an input signal Si accordingly. For example, the reception unit 102 includes an antenna 102a, a cyclic prefix remover 102b and a fast Fourier transformer 102c. The antenna 102a receives a wireless signal Sr. The cyclic prefix remover 102b receives and outputs the wireless signal Sr. The fast Fourier transformer 102c generates and provides an input symbol Yn according to the wireless signal Sr provided by the cyclic prefix remover 102b. Thus, the channel estimation device 1 executes related operation of channel estimation in a frequency domain.
The channel estimation unit 104 obtains several pilot channel gains by calculating several pilot channels transmitting the pilot signals of the input symbol Yn. For example, the channel estimation unit 104 further includes a Wiener filter. The channel estimation unit 104 executes an operation of interpolation on the pilot channel gains according to the linear minimum mean squared error (LMMSE) algorithm to obtain several data channel gains by respectively calculating several data channels transmitting several data signals of the input symbol Yn. The channel estimation unit 104 outputs the estimated data channel gains to the channel adaptation parameter tracking unit 106. For example, the data channel estimation result outputted by the channel estimation unit 104 can be expressed by a data channel gain vector Ĥ.
Furthermore, the operation of the Wiener filter can be expressed as the following formula (1):
Ĥ=W
d
H
×{tilde over (H)}
p (1)
Wherein, WdH denotes a conjugate transpose matrix of the Wiener filter matrix Wd; {tilde over (H)}p denotes a pilot channel gain vector; Ĥ denotes an estimated data channel gain vector.
For example, the Wiener filter matrix Wd can be expressed according to the following formula (2):
W
d=(Rp+σn2I)−1rdp (2)
Wherein, Rp denotes an auto-correlation matrix E{HpHpH}; rdp denotes a cross-correlation matrix E{HpHd}; I denotes a unit matrix; σn2 denotes a noise variance.
In an example, the uniform distribution power delay profile (PDP) model is applied for estimating the auto-correlation matrix and the cross-correlation matrix related to frequency. For example, the frequency related formula following the Fourier transformation is expressed as the following formula (3):
R
f,uniform
[Δf]=e
−j2πτ
Δf sinc(TΔf) (3)
Wherein, the sinc(TΔf) satisfies formula (4):
Parameter τm and T respectively satisfy formulas (5) and (6):
Wherein parameters τm and τrms respectively denote the average value and the mean squared root of delay corresponding to respective transmission path in a multiple transmission path environment. In other words, the parameters τm and τrms can be viewed as the statistical characteristics parameters of the multi-path environment for the current transmission system.
In another example, the existing exponential distribution PDP model is applied for estimating the auto-correlation matrix and the cross-correlation matrix related to frequency. For example, the frequency related formula following Fourier transformation is expressed as the following formula (7):
Wherein τ0 satisfies formula (8):
τ0=τm−τrms (8)
In other words, the Wiener filter matrix Wd is a function using parameters τm and τrms as variables.
The data channel gain vector Ĥ provided by the channel estimation unit 104 is further provided to the equalizing/quantifying unit 108, which accordingly obtains a signal {circumflex over (X)} through restoration. For example, the equalizing/quantifying unit 108 includes an equalizer 108a and a quantifier 108b, respectively executing equalization and quantification operations.
The demapper unit 110 receives the restored signal {circumflex over (X)} obtained by the equalizing/quantifying unit 108, and accordingly generates a soft bit signal.
The channel adaptation parameter tracking unit 106 calculates an adaptive alteration for the multi-path statistical characteristic parameters τm and τrms according to the estimated channel gain vector Ĥ. For example, the channel adaptation parameter tracking unit 106 locates a channel estimation error e according to the difference between the actual data channel gain vector Hd and the estimated data channel gain vector Ĥ and accordingly generates a square absolute value function J of channel estimation error e. Since the actual data channel gain vector Hd is given, the estimated data channel gain vector Ĥ satisfies formula (1) (that is, the product of the Wiener filter matrix W multiplied by the pilot channel gain vector wherein {tilde over (H)}p), and the Wiener filter matrix W is a function of multi-path statistical characteristic parameters τm and τrms, and the square absolute value function J of channel estimation error e is also the function of the multi-path statistical characteristic parameters τm and τrms.
The channel adaptation parameter tracking unit 106 further approximates the optimum solution of the multi-path statistical characteristic parameters τm and τrms according to the square absolute value function J with respect to the maximum value of the gradient of multi-path statistical characteristic parameters τm and τrms.
Furthermore, the channel estimation error e and the square absolute value function J respectively satisfy formula (9) and formula (10):
e=H
d
−Ĥ=H
d
−W×{tilde over (H)}
P (9)
J=∥e∥
2
=∥H
d
−W×{tilde over (H)}
P∥2 (10)
For example, the channel adaptation parameter tracking unit 106 estimates the actual data channel gain vector Hd according to the following formula (11):
The gradient operator of the square absolute value function J with respect to multi-path statistical characteristic parameters τm and τrms can be expressed as formula (12):
Wherein Wm is defined as
and Wrms is defined as
After the square absolute value function J is calculated, the following formula (13) can be obtained:
Wherein Λ satisfies the following formula (14):
Wherein λm and λrms are constants. Thus, the channel adaptation parameter tracking unit 106 provides adjusted multi-path statistical characteristic parameters τm(n+1) and τrms (n+1) to the channel estimation unit 104. The channel estimation unit 104 accordingly generates an updated Wiener filter, and executes an operation of channel estimation on a next input symbol Yn+1 of the input signal Si.
Referring to
Then, the method proceeds to step (c), an adaptive alteration for multi-path statistical characteristic parameters τm(n) and rrms(n) is calculated by the channel adaptation parameter tracking unit 106 according to the data channel gain vector Ĥ and the pilot channel gain vector Ĥp, and an adjusted multi-path statistical characteristic parameters τm(n+1) and τrms(n+1) is obtained through adjustment.
Afterwards, the method proceeds to step (d), an updated Wiener filter is generated by the channel estimation unit 104 according to the adjusted multi-path statistical characteristic parameters τm(n+1) and τrms(n+1). Then, the method proceeds to step (e), the value of n is increased by 1 and steps (a) and (b) are repeated, so that an operation of channel estimation is executed on a next input symbol of the input signal Yn+1 by the channel estimation unit 104 according to the updated Wiener filter.
Referring to
Referring to
Furthermore, the ICI estimation unit 212 executes an operation of ICI estimation according to the pilot channel gain vector {tilde over (H)}p and the data channel gain vector Ĥ to locate the ICI gain I for each pilot channel and each data channel. The PIC unit 214 executes a PIC operation on each pilot channel and each data channel according to the ICI gain I corresponding to each pilot channel and each data channel to locate interference cancelled input symbol Yn(i), wherein the parameter i indicates the number of times of the iteration of ICI cancellation, and is such as a natural number ranging between 1 and a predetermined number of times i_max, which can be realized by any natural number larger than or equal to 1.
The interference cancelled input symbol Yn(i) is further provided to the channel estimation unit 204, so that the channel estimation unit 204 again executes an operation of channel estimation to locate the data channel gain vector Ĥ. Then, the ICI estimation unit 212 and the PIC unit 214 again execute a PIC operation on each pilot channel and each data channel according to the ICI gain I corresponding to each pilot channel and each data channel to locate an input symbol Yn(i+1) in which the interference is cancelled. In other words, the ICI estimation unit 212, the PIC unit 214 and the channel estimation unit 204 form an iteration path to continuously execute an iteration of ICI cancellation on the input symbol Yn. After a predetermined number of times of iteration, the above unit skips the iteration path, and outputs a data channel gain vector Ĥ in which the ICI is cancelled.
In greater details, the ICI estimation unit 212 estimates the ICI gain I[k, k+d] resulted in the k-th sub-carrier channel by the (k+d)-th sub-carrier channel according to the following formula (15):
Wherein h′[k] denotes a slope distribution of the k-th sub-carrier channel, and k and d both are a natural number.
For example, the ICI estimation unit 212 of the present embodiment of the invention simplifies formula (15) as formula (17) with reference to the formula (16) to simplify the estimation of I[k,k+d]:
For the PCI unit 214, the operation can be expressed as the following formula (18):
Wherein parameter i denotes the number of times of iteration, and the limit of the parameter i is a predetermined value, which determines the number of times of iteration, I[, k+m] denotes the estimate of ICI gain resulted in the k-th sub-carrier channel by the (k+m)-th sub-carrier, and {circumflex over (X)}[k+m] denotes a restored data.
Referring to
Furthermore, after the data channel gain vector Ĥ is located, the channel estimation method proceeds to step (g), an operation of ICI estimation is executed by the ICI estimation unit 212 according to the pilot channel gain vector {tilde over (H)}p and the data channel gain vector Ĥ to locate an ICI gain I for each pilot channel and each data channel. For example, the ICI gain I[k, k+d] resulted in the k-th sub-carrier channel by the (k+d)-th sub-carrier channel can be expressed as formula (15).
Then, the method proceeds to step (h), a PIC operation is executed on each pilot channel and each data channel by the PCI unit 214 according to the ICI gain I to locate an input symbol Yn(i) in which interference is cancelled. After that, the method proceeds to step (i), whether the iteration of interference cancellation has been executed for the predetermined number of times i_max is determined by the channel estimation unit 204 according to whether the parameter i is equal to the predetermined number of times i_max. If not, the channel estimation unit 204 again calculates the data channel gain vector Ĥ according to the input symbol Yn(i), in which interference is cancelled, and provides the data channel gains to the ICI estimation unit 212 for executing the iteration of interference cancellation again.
When the parameter i is equal to the predetermined number of times i_max, the channel estimation unit 204 provides the data channel gain vector Ĥ to the channel adaptation parameter tracking unit 206. In other words, the channel estimation method sequentially executes steps (c)˜(e) to complete the remaining operations of the channel estimation method. The details of steps (c)˜(e) are similar to corresponding steps of the channel estimation method of the first embodiment, and the similarities are not repeated here.
According to the channel estimation method with channel adaptation and the device using the same disclosed in the above embodiments, an input symbol is received, and an operation of channel estimation is executed by a Wiener filter according to the pilot channel gains corresponding to the pilot signal of the input symbol. An adaptive alteration for the first and the second multi-path statistical characteristic parameters of the Wiener filter is calculated according to the aforementioned channel estimation result and the pilot channel gains, and an updated Wiener filter is generated accordingly. An operation of channel estimation is executed on a next input symbol according to the updated Wiener filter, so that the Wiener filter used for channel estimation can be adjusted according to the input signal and existing channel conditions and adapted to rapid change in channel characteristics in the occasions of high-speed movement. Compared to the conventional channel estimation method, the channel estimation method with channel adaptation and the device using the same disclosed in above embodiments can execute adaptive adjustment on the operation of channel estimation in response to rapid change in channel characteristics in the occasions of high-speed movement to achieve superior efficiency of the entire system.
It will be apparent to those skilled in the art that various modifications and variations can be made to the disclosed embodiments. It is intended that the specification and examples be considered as exemplary only, with a true scope of the disclosure being indicated by the following claims and their equivalents.
Number | Date | Country | Kind |
---|---|---|---|
100130179 | Aug 2011 | TW | national |