This application claims the benefit of Korean Patent Application No. 10-2007-0086278, filed on Aug. 27, 2007, in the Korean Intellectual Property Office, the disclosure of which is incorporated herein in its entirety by reference.
1. Technical Field
The present disclosure relates to a method of estimating inter-carrier interference (ICI) and an ICI mitigating equalizer and, more particularly, to a method of estimating an ICI component by approximating a channel to an (M1)-th order model and an ICI mitigating equalizer.
2. Discussion of Related Art
Orthogonal frequency division multiplexing (OFDM) systems are applied to various broadcasting/communication systems. An OFDM system serves to transmit data by using a plurality of sub-carriers that are orthogonal to one another.
Referring to
In the OFDM transmitter 110, a cyclic prefix (CP) is inserted into a transmitted signal that is transformed through IFFT so as to prevent inter-symbol interference (ISI) and to estimate a channel. The transmitted signal Si(n) into which the CP is inserted is output to a wired/wireless channel via the RF transmitter 115 and the transmitting antenna 116. In the OFDM receiver 120, a received signal ri(n) from which the CP is removed is transformed into a received signal Ri(k) in a frequency domain by using the FFT block 124. The equalizer 125 outputs a transmitting signal Ŝi(k) estimated by equalizing the received signal Ri(k) in the frequency domain.
In order to accurately estimate the transmitted signal by using the OFDM receiver 120, a response characteristic of a channel has to be accurately recognized. The response characteristic of the channel may be changed according to time, however due to time-selectivity. The response characteristic of the channel may be changed according to frequency due to frequency-selectivity. Accordingly, it is difficult to accurately recognize the response characteristic of the channel. On the other hand, the time-selectivity and the frequency-selectivity of a channel damages orthogonality among sub-carriers and causes inter-carrier interference (ICI). Although the response characteristic of the channel may be accurately recognized by using pilots included in the CP, it may be inaccurate to estimate the transmitted signal due to an influence of the ICI.
The ICI increases the complexity of the accurate estimation of the transmitted signal by using the OFDM receiver 120. In order to allow the OFDM receiver 120 to accurately estimate the transmitted signal, the ICI component included in the received signal is initially estimated, and the ICI component included in the received signal has to be removed or mitigated based on the estimation result. Because a calculation amount and computational complexity increase in order to remove or mitigate the ICI component included in the received signal, the structure of the OFDM receiver 120 becomes complex. The OFDM receiver 120 is designed in consideration of a trade-off relation between accuracy of estimation and computational complexity.
Exemplary embodiments of the present invention provide a method of estimating an inter-carrier interference (ICI) component by approximating a channel to an (M1)-th order model. Exemplary embodiments of the present invention also provide an ICI mitigating equalizer for removing or mitigating an ICI component that is estimated by using the method of estimating the ICI component.
According to an exemplary embodiment of the present invention, there is provided an inter-carrier interference (ICI) mitigating equalizer comprising a channel estimator, a channel calculator, an ICI estimator, a subtracter, and an equalizer. The channel estimator estimates a channel response from a received signal. The channel calculator calculates a basic component and fluctuating components of a channel response from the estimated channel response. The ICI estimator multiplies the fluctuating components of the channel response by a received signal in a frequency domain, filters the multiplication results according to filter coefficients, and estimates an ICI component included in the received signal based on the filtering results. The subtracter subtracts the ICI component from the received signal in the frequency domain. The equalizer equalizes an output signal of the subtracter based on the basic component of the channel response.
In an exemplary embodiment of the present invention, the basic component of the channel response may be channel response component that is time-invariant for one symbol interval, and the fluctuating components of the channel response may be a channel response components that are time-variant for one symbol interval. In addition, in a case where a channel is approximated to an (M1)-th order model, the channel calculator may output first order to (M1)-th order fluctuating components of the channel response to the ICI estimator.
In addition, the ICI estimator may include: a multiplication unit multiplying the fluctuating components of the channel response by the received signal in the frequency domain and outputting the multiplication results; a filter bank filtering the multiplication result according to the filter coefficients and outputting the filtering results; and an adder adding up the filtering result and outputting the ICI component.
In addition, in a case where a channel is approximated to an (M1)-th order model, the multiplication unit may include first to (M1)-th multipliers that respectively multiply first order to (M1)-th order fluctuating components of the channel response that are respectively output from the channel calculator by the received signal in the frequency domain to output the first to (M1)-th multiplication results.
In addition, in a case where a channel is approximated to an (M1)-th order model, the filter bank may include first to (M1)-th filters that respectively filter the first to (M1)-th multiplication results that are respectively output from the first to (M1)-th multipliers of the multiplication unit according to first to (M1)-th filter coefficients to output the first to (M1)-th filtering results. In addition, the first to (M1)-th filters may be linear, time-invariant filters. In addition, the first to (M1)-th filters may be finite impulse response (FIR) filters.
In addition, the equalizer may be embodied as a one-tap equalizer.
Furthermore, the aforementioned ICI mitigating equalizer may be applied to a receiver of an orthogonal frequency division multiplexing (OFDM) system.
According to an exemplary embodiment of the present invention, there is provided a method of estimating the ICI, the method comprising: calculating first order to (M1)-th order fluctuating components of a channel response from a channel response that is estimated from a received signal; multiplying each of the first order to (M1)-th order fluctuating components of the channel response by a received signal in a frequency domain and outputting first to (M1)-th multiplication results; filtering the first to (M1)-th multiplication results respectively according to first to (M1)-th filter coefficients and outputting first to (M1)-th filtering results; and estimating an ICI component included in the received signal by adding up the first to (M1)-th filtering results.
In an exemplary embodiment of the present invention, in a case where a channel is approximated to an (M1)-th order model, the channel response that is estimated from the received signal may be divided into a basic component that is time-invariant for one symbol interval and the first order to (M1)-th order fluctuating components that are time-variant for one symbol interval.
According to an exemplary embodiment of the present invention, there is provided an ICI mitigating equalizer comprising: a channel estimator estimating a channel response from a received signal; a channel calculator calculating a basic component and a first order fluctuating component of a channel response from the estimated channel response; a first multiplier multiplying the first order fluctuating component of the channel response by a received signal in a frequency domain and outputting a first multiplication result; a first filter FIR-filtering the first multiplication result according to first filter coefficients and outputting an ICI component included in the received signal; a subtracter subtracting the ICI component from the received signal in the frequency domain; and an equalizer equalizing an output signal of the subtracter according to the basic component of the channel response.
In an exemplary embodiment of the present invention, the channel calculator may approximate a channel to a first order linear model. In addition, the channel calculator may output the estimated value of the channel response in the present symbol interval as the basic component of the channel response. In addition, the channel calculator may subtract the estimated value of the channel response in the previous symbol interval from the estimated value of the channel response in the next symbol interval, divide the subtraction result by the estimated value of the channel response in the present symbol interval, and output the division result as the first order fluctuating component of the channel response.
In addition, the channel calculator may include: a first delayer delaying the estimated value of the channel response that is output from the channel estimator by one symbol interval and outputting the delayed estimated value; a second delayer delaying an output signal of the first delayer by one symbol interval again and outputting the delayed output signal; a subtraction unit subtracting the estimated value of the channel response in the previous symbol interval that is output from the second delayer from the estimated value of the channel response in the next symbol interval that is output from the channel estimator; and a divider dividing the subtraction result that is output from the subtraction unit by the estimated value of the channel response in the present symbol interval that is output from the first delayer and outputting the division result to the first multiplier as the first order fluctuating component of the channel response.
In addition, the equalizer may divide the output signal of the subtraction unit by the basic component of the channel response that is output from the first delayer to output the division result. In addition, the equalizer may be embodied as a one-tap equalizer.
The ICI mitigating equalizer according to an exemplary embodiment of the present invention requires a lesser amount of calculation and lower computational complexity as compared with a conventional technique that requires a large amount of calculation and high computational complexity in order to remove the ICI component included in the received signal.
Exemplary embodiments of the present invention will be understood in more detail from the following description taken in conjunction with the attached drawings, in which:
Hereinafter, the present invention will be described in detail by explaining exemplary embodiments with reference to the attached drawings.
Referring to
The channel estimator 240 estimates a channel response {tilde over (H)} from a received signal.
The channel calculator 250 calculates a basic component 0 and fluctuating components
The ICI estimator 260 multiplies the fluctuating components
The subtracter 270 subtracts the ICI component Îi(k) from the received signal Ri(k) in the frequency domain. The equalizer 280 equalizes an output signal of the subtracter 270 according to the basic component 0 of the channel response. The equalizer 280 outputs the equalization result as an estimated transmitted signal Ŝi(k).
The ICI mitigating equalizer 230 illustrated in
Referring to
A transmitted signal Si(n) that is output from a transmitter of an OFDM system can be represented by Equation 1 as follows:
where Ng is a size of a cyclic prefix (CP), N is a FFT size, that is, the total number of sub-carriers, and {Si(k)}k=0N−1 is an i-th quadrature amplitude modulation (QAM)-mapping symbol.
When a discrete response of a linear time variant multi-path channel is hltv(n,l), a channel response hi(n,l) for the path with a relative delay of l numbers of sampling intervals can be defined by Equation 2 as follows:
where P is a maximum delay among delays due to a multi-path.
In addition, when the transmitted signal Si(n) is represented by Equation 1 and when the channel response hi(n,l) is defined by Equation 2, the received signal ri(n) obtained by removing the CP in the receiver of the OFDM system is represented by Equation 3 as follows:
where vi(n) is additional white Gaussian noise (AWGN).
The received signal ri(n) in the time domain is transformed into the received signal Ri(k) in the frequency domain through FFT represented by Equation 4 as follows:
and wi(k) is a signal obtained by transforming v1(n) through FFT.
Referring to Equation 4, the received signal in the frequency domain can be represented in a matrix form by Equation 5 as follows:
Ri=Gi·Si+wi [Equation 5]
Ri=[Ri(0) . . . Ri(N−1)]T
Si=[Si(0) . . . Si(N−1)]T
wi=[wi(0) . . . wi(N−1)]T
Gi=(Gk,m)N×N
where the matrix Gi is a channel gain matrix (or an equalization matrix). When orthogonality among sub-carriers is maintained, all the elements except the main diagonal elements of the matrix Gi are zero. When the orthogonality among the sub-carriers is damaged and when the ICI component is included in the received signal, however, elements except the main diagonal elements of the matrix Gi are non-zero. When the elements except the main diagonal elements of the matrix Gi are non-zero, it becomes more difficult for the OFDM receiver to estimate the transmitted signal.
If the AWGN component wi(k) s excluded in Equation 5, a symbol can be estimated by Equation 6.
Ŝi=Gi−1Ri [Equation 6]
In order to simplify the equalization matrix Gi, in an exemplary embodiment of the present invention, a channel is approximated to an (M1)-th order model by Equation 8.
First, a channel vector {tilde over (h)}(l) s defined for an i-th symbol interval by Equation 7 as follows:
{tilde over (h)}(l)≡[{tilde over (h)}1(l){tilde over (h)}2(l) . . . {tilde over (h)}M(l)]T [Equation 7]
where components of the channel vector {tilde over (h)}(l) are M numbers of sample values that are estimated from the channel response hltv(n,l). The components of the channel vector {tilde over (h)}(l) correspond to the channel response {tilde over (H)} of
In Equation 8, 0≦p≦M1, and the (M1)-th order model of the channel is determined by αp(n) and a channel coefficient βpq. In an exemplary embodiment of the present invention, the channel coefficient βpq can be set to a constant for an i-th symbol interval. In addition, in an embodiment of the present invention, αp(n) can be set so as to satisfy Equation 9.
As shown in Equation 8, αp(n) is related to the fluctuating components of the channel response. In the fluctuating components of the channel response in Equation 8, an argument n is separated from an argument l. On the other hand, the basic component of the channel response includes most of the decisive energy of the channel response hi(n,l).
Equation 10 is obtained by substituting Equation 8 for hi(n,l) in the equation of Gk,m in Equation 4.
In Equation 10, δ(m) is a Kronecker delta, that is, a unit delta function,
Equation 10 can be represented as a matrix form by Equation 11 as follows:
That is, {hacek over (H)}p is an N×N diagonal matrix.
In addition, a filter coefficient matrix Fp in Equation 11 is as follows:
where the filter coefficient matrix Fp is independent of a symbol index i and has a Toplitzian property.
Equation 12 is obtained by substituting Equation 11 for Gi of Equation 6.
In Equation 12, IN×N is an N×N identity matrix, and
Equation 13 is obtained by truncating high-order components so as to simplify Equation 12.
Because the channel response is changed very slowly for the symbol interval, αp(n) is also changed very slowly. Accordingly, Fp(m) obtained by transforming αp(n) through FFT has a low-pass property. In an exemplary embodiment of the present invention, only components due to 2Qp numbers of sub-carriers among N numbers of sub-carriers are considered according to the aforementioned property. That is, only sub-carriers of which carrier indices ranges from −Qp to +Qp are considered. Equation 14 is induced from Equation 13 by reflecting the aforementioned fact.
Hereinafter, an exemplary embodiment of the present invention will be described with reference to Equation 14 and
The channel calculator 350 calculates a basic component 0 and fluctuating components
The channel calculator 350 calculates {hacek over (H)}p (0≦p≦M1) from the estimated channel response H by using Equation 15. Here, {hacek over (H)}p (p=0) is output to the equalizer 380 as the basic component 0 of the channel response. In addition, the channel calculator 350 calculates
In
In
Hereinafter, the FIR filters will be described with reference to
Referring to
Although
In
In
In
As described above, the ICI mitigating equalizer for estimating and removing the ICI component included in the received signal according to an exemplary embodiment of the present invention performs (2Σp=1M
On the other hand, in
A method of estimating ICI according to an exemplary embodiment of the present invention will be described below based on the aforementioned description.
First, first order to (M1)-th order fluctuating components
As described above, in a case where a channel is approximated to the (M1)-th order model, the channel response {tilde over (H)} or {tilde over (h)} that is estimated from the received signal can be divided into a basic component that is time-invariant for a symbol interval and first order to (M1)-th order fluctuating components that are time-variant for a symbol interval.
Referring to
The channel estimator 540 estimates a channel response {tilde over (H)} or {tilde over (h)} that is estimated from a received signal. Hereinafter, an estimation of a channel response by a channel estimator 540 in a case where a channel is approximated to a first order linear model will be described with reference to
Referring to
The i-th symbol is assumed to be a present symbol. A case where sample values of channel responses in the (i−1)-th to (i+1)-th symbol intervals are used so as to approximate a channel to the first order linear model for the i-th symbol interval is considered. In this case, the channel estimator 540 outputs an estimated value Ĥi−1 of a channel response in the previous symbol interval according to a sample value of a channel response at the center of the (i−1)-th symbol interval, outputs an estimated value Ĥi of a channel response in the present symbol interval according to a sample value of a channel response at the center of the i-th symbol interval, and outputs an estimated value {tilde over (H)}i+1 of a channel response in the next symbol interval according to a sample value of a channel response at the center of (i+1)-th symbol interval. It will be understood by those of ordinary skill in the art that Ĥi−1 corresponds to {tilde over (H)}i−1, Ĥi corresponds to {tilde over (H)}i, and {tilde over (H)}i+1 corresponds to {tilde over (H)}i+1.
In a case where a channel is approximated to a first order linear model for the i-th symbol interval by using a sample value {tilde over (h)}i−1 of the channel response at the center of the (i−1)-th symbol interval and a sample value {tilde over (h)}i+1 of the channel response at the center of the (i+1)-th symbol interval, Equation 7 that defines the channel vector for the i-th symbol interval can be represented by Equation 16. Here, M in Equation 7 is 3 in Equation 16.
{tilde over (h)}(l)(n,l)≡[{tilde over (h)}i−1(N/2−1,l){tilde over (h)}1(N/2−1,l){tilde over (h)}i+1(N/2−1,l)]T [Equation 16]
In order to approximate the channel to the first order linear model, Equation 17 is obtained by substituting 1 for M1 of Equation 8.
On the other hand, Equation 19 is obtained by substituting Equation 18 for the channel coefficient βpq included in an equation of {hacek over (h)}p(l) in Equation 8.
A minimum-mean-squared-error (MMSE) of a signal obtained by transforming {hacek over (h)}i(l) through FFT is estimated by Equation 20 with reference to an equation of {hacek over (h)}1(l) in Equation 19.
dH1(k)=Ĥi+1(k)−Ĥi−1(k) [Equation 20]
The first order fluctuating component
The basic component of the channel response can be represented by Equation 22 with reference to an equation of {hacek over (h)}0(l) in Equation 19. In this exemplary embodiment, the basic component of the channel response is a mean value of the channel responses in the (i−1)-th to (i+1)-th symbol intervals.
0=Ĥi(k) [Equation 22]
Equation 23 is obtained by rewriting Equation 14 according to Equations 16 to 22 so as to be suitable for a case where a channel is approximated to the first order linear model.
Here, F1(q) in Equation 23 is obtained by transforming α1(n) in Equation 17 through FFT. As shown in
Hereinafter, the ICI mitigating equalizer shown in
In
As shown in
In order to calculate Equations 22 and 21, the first delayer 551 included in the channel calculator 550 delays the estimated value of the channel response that is output from the channel estimator 540 by a symbol interval and outputs the delayed estimated value. The second delayer 553 included in the channel calculator 550 delays the output signal of the first delayer 551 by a symbol interval again and outputs the delayed output signal. The subtracter 555 included in the channel calculator 550 subtracts the estimated value Ĥi−1 of the channel response in the previous symbol interval that is output from the second delayer 553 from the estimated value Ĥi+1 of the channel response in the next symbol interval that is output from the channel estimator 540. The divider 557 included in the channel calculator 550 divides the subtraction result that is output from the subtracter 555 by the estimated value Ĥi of the channel response in the current symbol interval that is output from the first delayer 551 and outputs the division result to the first multiplier 561 as the first order fluctuating component
In
In
In
In
The simulation conditions of
In
In
While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it will be understood by those of ordinary skill in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. The exemplary embodiments should be considered in descriptive sense only and not for purposes of limitation. Therefore, the scope of the invention is defined not by the detailed description of the invention but by the appended claims, and all differences within the scope will be construed as being included in the present invention.
Number | Date | Country | Kind |
---|---|---|---|
10-2007-0086278 | Aug 2007 | KR | national |
Number | Name | Date | Kind |
---|---|---|---|
7120209 | Gorokhov et al. | Oct 2006 | B2 |
20040005010 | He et al. | Jan 2004 | A1 |
20060239367 | Wilhelmsson et al. | Oct 2006 | A1 |
20060251198 | Ma et al. | Nov 2006 | A1 |
Number | Date | Country |
---|---|---|
1514555 | Jul 2004 | CN |
1635724 | Jul 2005 | CN |
2005-150839 | Jun 2005 | JP |
Number | Date | Country | |
---|---|---|---|
20090060016 A1 | Mar 2009 | US |