The present disclosure relates to channel estimation in a wireless communication system.
The background description provided herein is for the purpose of generally presenting the context of the disclosure. Work of the presently named inventors, to the extent it is described in this background section, as well as aspects of the description which may not otherwise qualify as prior art at the time of filing, are neither expressly or impliedly admitted as prior art against the present disclosure.
Some multiple input, multiple output (MIMO) wireless communication systems can estimate channel conditions, or gains, in the communication path between the transmitting and receiving antennas. The channel estimation process can include transmitting known training symbols, receiving the known training symbols, and processing the received symbols to estimate the channel conditions. The estimation is based on differences between the known training symbols and the received symbols. Information regarding the channel conditions can then be used to program coefficients of an equalizer of the receiver. The equalizer then compensates for the channel conditions.
Referring now to
Transmitter module 12 periodically generates a plurality of long training fields (LTFs) 18-1, . . . , 18-j, referred to collectively as LTFs 18. Each LTF 18 includes a plurality of training symbols 20-1, . . . , 20-k, referred to collectively as training symbols 20. A multiplier module 22 multiplies each training symbol 20 by a corresponding column of a preamble steering matrix P. A number of rows n of matrix P corresponds with a number of transmit antennas 26-1, . . . , 26-n, collectively referred to as antennas 26. The number of columns j of matrix P corresponds with the number of LTFs 18. Matrix P assures the orthogonality of training symbols 20 as they are transmitted from antennas 26. Matrix P has a condition number of 1, i.e. cond(P)=1.
Receiver module 14 includes receiver antennas 30-1, . . . , 30-n, collectively referred to as antennas 30, that receive the training symbols via channel 16. After receiving all of the training symbols, receiver module 14 generates matrix H based on known training symbols 20, matrix P, and the received training symbols. Receiver module 14 can then use matrix H to adjust coefficients of an internal equalization module for signals from antennas 30. It is generally desirable for receiver module 14 to generate matrix H as quickly as possible.
A sample estimation of matrix H will now be described. Assume that n=3 and j=4. Transmitter module 12 then sends 4 LTFs 18 and matrix P is a 3×4 matrix
The effective MIMO channel estimated at receiver module 14 is given by Hest=HPH=HestP−1, where Hest represents an estimation of matrix H. For the data associated with each LTF 18 the transmitter-to-receiver communication model can be described by y=HestP−1x+n, where x represents transmitted data symbols. The ZF solution is applied to the matrix HestP−1.
If the matrix P were not used, y=Hx+n{circumflex over (x)}=R−1Q*y, where y represents received data symbols. With matrix P, y=HestP−1x+n. Receiver module 14 uses each LTF 18 to estimate a column of matrix Hest. Matrix H can therefore be estimated by waiting until all columns have been estimated and then estimating the matrix HestP−1. Receiver module 14 can then perform orthogonal-triangular decomposition (QR) on HestP−1, i.e QR(HestP−1). However, the computational density increases to the order of n3, i.e. O(n3). To meet the processing latency the hardware burden would also increase based on O(n3).
The effect of matrix P will now be described to shed light on the above equations. Let Hest=QR. Without matrix P the equalized vector is given by
With matrix P the equalized vector is given by:
The equivalent matrix RP−1 is a full matrix and it is difficult to compute its inverse for an n×n communication system.
A receiver module includes an input that receives a data message from a wireless communication channel. The data message has a plurality of training fields and data. A channel estimator module recursively estimates a matrix H that represents the channel based on the plurality of training fields. The recursive estimation is performed as the plurality of training fields are being received. An equalizer module applies coefficients to the data based on the matrix H.
Further areas of applicability of the present disclosure will become apparent from the detailed description provided hereinafter. It should be understood that the detailed description and specific examples, while indicating the preferred embodiment of the disclosure, are intended for purposes of illustration only and are not intended to limit the scope of the disclosure.
The present disclosure will become more fully understood from the detailed description and the accompanying drawings, wherein:
The following description is merely exemplary in nature and is in no way intended to limit the disclosure, its application, or uses. For purposes of clarity, the same reference numbers will be used in the drawings to identify similar elements. As used herein, the term module, circuit and/or device refers to an Application Specific Integrated Circuit (ASIC), an electronic circuit, a processor (shared, dedicated, or group) and memory that execute one or more software or firmware programs, a combinational logic circuit, and/or other suitable components that provide the described functionality. As used herein, the phrase at least one of A, B, and C should be construed to mean a logical (A or B or C), using a non-exclusive logical or. It should be understood that steps within a method may be executed in different order without altering the principles of the present disclosure.
Referring now to
Communication system 50 will now be described in pertinent part. A baseband module 58 generates data messages based on m streams of incoming data. Baseband module 58 communicates the data messages (shown in
The transmitted data streams propagate through communication channel 54. Communication channel 54 perturbs the transmitted data streams due to phenomena such as reflections, signal attenuation, and so forth. The perturbations can be represented by matrix H.
Receiver module 52 includes n RF receivers that are represented by antennas 70-1, . . . , 70-n. The RF receivers receive the transmitted data streams and communicate the perturbed time domain signals to a channel estimator module 72. Channel estimator module 72 estimates matrix H based on matrix P and the long training fields that are included in the received data streams. In some embodiments channel estimator module 72 includes a processor 73 and associated memory 75 for storing and/or executing the recursive channel estimation methods that are described below.
Channel estimator module 72 communicates the n received data streams to n respective fast-Fourier transform (FFT) modules 74 and adjusts gains of an equalizer module 76. FFT modules 74 convert the time-domain data streams to frequency-domain data streams and communicate them to equalizer module 76. Equalizer module 76 compensates the respective data streams based on the gains and communicates the compensated gains to a Viterbi decoder module 78. Viterbi decoder module 78 decodes the n data streams to generate received data streams ym.
Referring now to
Referring now to
Referring now to
Control enters at a block 122 and proceeds to decision block 124. In decision block 124 control determines whether an LTF 104 is being received. If not then control returns to block 122. If an LTF 104 is being received then control branches from decision block 124 to block 126. In block 126 control receives a training symbol 106 that is associated with the current LTF 104. Control then proceeds to block 128 and updates a matrix Hest, which is described below in more detail, based on the current training symbol 106. Control then proceeds to decision block 130 and determines whether the current training symbol 106 was the last training symbol 106 of the present LTF 104. If not then control branches to block 132 and waits for the next training symbol 106 of the current LTF 104. When the next training symbol 106 is received control returns to block 126 and repeats the aforementioned steps for the new training symbol 106. On the other hand, if the training symbol 106 in decision block 130 was the last training symbol 106 of the present LTF 104 then control branches to decision block 134.
In decision block 134 control determines whether the current LTF 104 was the last LTF 104 (i.e. LTF 104-x) of the current group of LTFs 104. If not then control branches to block 136 and waits for the next LTF 104 to begin before returning to block 126. On the other hand, if the current LTF 104 was the last LTF 104-x then control branches from decision block 134 to block 138. In block 138 control generates matrix H based on matrix Hest and matrix P. Control then proceeds to block 140 and adjusts the gains of equalizer module 76 based on matrix H. Control then returns to other processes via termination block 142.
As channel estimator module 72 executes method 120 it performs distributed QR across LTFs 104. That is, QR(Hest). Computational density therefore increases as O(n2) and the processing latency of associated hardware and/or processor 73 would increase by about O(n2). This represents an improvement, e.g. reduced need for processing power, over the prior art. Example estimations of matrix Hest will now be provided for various MIMO dimensions of communications system 50.
In 2×2 and 2×3 MIMO cases the equalized vector is given by
It can be seen that
For a general n×n MIMO communication system 50 the equalized vector is given by
where Hest,n=QnRn.
Methods are known in the art for recursively solving the Q*ny term of the above equation. Recursive computation of the PnRn−1 and wll,n=1./diag(PnRn−1Rn−*PnT) terms of the above equation will now be described.
Let wll,n=1./diag(Pn×nRn×n−1Rn×n−*Pn×nT). It can be seen that for 1≦j<n the jth element of the wll vector for n streams can be recursively computed as follows:
For j=n,
where λn=P(n,1:n−1)Rn-1−1Rn-1−*P(n,1:n−1)T.
A proof of the immediately preceding equations will now be provided.
Compute R2×2.
Compute
Compute 1/wll,1=(r222+∥r11+r12∥2)/r112r222.
Compute 1/wll,2=(r222+∥r11−r12∥2)/r112r222.
Let λ=1/wll,2.
Update the inverse of the triangular matrix based on
Update the substream SNRs based on
1/wll,1→1/wll,1+∥ρ1−ρ2+ρ3∥2
1/wll,2→1/wll,2+∥ρ1−ρ2+ρ3∥2
1/wll,3→λ+∥ρ1−ρ2+ρ3∥2.
Update the lamda factor based on
λ=1/wll,1+1/wll,1−1/wll,3+8real(ρ2)ρ3.
Update the inverse of the triangular matrix based on
where
ρ1=(r33(r12r24−r14r22)−(r12r23−r13r22)r34)/r11r22r33r44,
ρ2=(r34r23−r24r33)/r22r33r44,ρ3=−r33/r33r44,ρ4=1/r44.
Update the substream SNR based on
1/wll,2→1/wll,2+∥ρ1+ρ2−ρ3+ρ4∥2
1/wll,3→1/wll,3+∥ρ1+ρ2−ρ3+ρ4∥2
1/wll,4→λ+∥ρ1+ρ2−ρ3+ρ4∥2
1/wll,4→λ+∥ρ1+ρ2−ρ3+ρ4∥2.
Compute their inverses and store them in a memory that can be included in channel estimator module 72.
The 3×3 MIMO case employs a non-square matrix P. For 3 streams transmitter module 56 sends 4 LTFs 104 and employs matrix P of
Channel estimator module 72 estimates channel matrix Hest3×4 and the real matrix is H=Hest3×4P3×4⊥. The received vector can be based on y=Hest3×4P3×4⊥x=Hx.
A distributed solution may also be employed by working directly with Hest3×4 without forming H. In the distributed solution let
Hest3×4=QR3×4 and
Hest3×4=[Hest3×3h4].
Then Q*Hest3×4=[RQ*h4] where QR=Hest3×3.
In some embodiments the equalized vector can be based on
{circumflex over (x)}=(Hest,3×4P3×4⊥)⊥y.
It can be seen that
{circumflex over (x)}=P1−1z−μv where u=R−1Q*h4, and μ=[1 −1 1]z.
The vector v is based on
v=kP1−1u where u=R−1Q*h4.
The scalar k is based on
k=1/(1+[1 −1 1]u).
The matrix P1−1 is based on
A proof of the immediately preceding equations will now be provided.
A solution is given by {circumflex over (x)}=(Hest,3×4P3×4⊥)−1y.
and write the matrix as
where v=P1−1R−1Q*h4.
Channel estimator module 72 computes
Wll=1/((Hest,3×4P3×4⊥)*(Hest,3×4P3×4⊥))−1
based on terms of v=kP1−1R−1Q*h4 and R−1R−*.
Let this matrix be
and during LTFs 104 compute
S1=4(I11−2I21+I22)
S2=4(I11−2I31+I33)
S3=4(I22+2I32+I33)
S=(I11+I22+I33)−2(I21+I31−I32)
S4=I11−2I21+I22−I31+I32+j(Q11+Q22+Q31−Q32)
S5=I11I21+I32−2I31+I33+j(Q11+Q21+Q32−Q33)
S6=I21−I22+I31−2I32+I33+j(Q21+Q22+Q31−Q33)
It can be seen that
Wll1=1/(S1+S∥v1∥2−4real(S4v*1))
Wll2=1/(S2+S∥v2∥2−4real(S5v*2))
Wll3=1/(S3+S∥v3∥2−4real(S6v*3))
and v=kP1−1u.
Proof of the above solution for Wll will now be provided. Let
H=[QRh4]P3×4⊥Q*H=[RQ*h4]P3×4⊥=R(P1+up1),u=R−1Q*h4
Then
Channel estimator module 72 still needs to recursively update R−1R−* after determining v=kP1−1u as described above. At a time n let
Rn−1=[Rn-1−1ρ]
and apply the following identity
Rn−1Rn−*=Rn-1−1Rn-1−*+ρρ*.
For example,
Channel estimator module 72 performs first column nulling by computing
R1−1=1/r11.
Channel estimator module 72 performs second column QR processing based on
where ρ1=−r12/r11r22, and ρ2=1/r22.
Channel estimator module 72 then computes ∥ρ1∥2,∥ρ2∥2,ρ1,ρ2* and updates R−1R−* based on
Channel estimator module 72 performs third column QR processing by recursively updating
where ρ1=(r12r23−r13r22)/r11r22r33 ρ2=−r23/r22r33 ρ3=1/r33.
Channel estimator module 72 can then compute
∥ρ1∥2,∥ρ2∥2,∥ρ2∥2,ρ1ρ*2,ρ1ρ*3,ρ2ρ*3
and update R−1R−* based on
R3×3−1R3×3−*=R2×2−1R2×2−*+ρρ*
Channel estimator module 72 can then compute the sums
S1=4(I11−2I21I22)
S2=4(I11−2I31+I33)
S3=4(I22+2I32+I33)
S=(I11+I22+I33)−2(I21+I31−I32)
S4=I11−2I21+I22−I31+I32+j(Q11+Q22+Q31−Q32)
S5=I11−I21+I32−2I31+I33+j(Q11+Q21+Q32−Q33)
S6=I21−I22+I31−2I32+I33+j(Q21+Q22+Q31−Q33)
Channel estimator module 72 computes
u=R−1Q*h4,
k=1/(1+[1 −1 1]u), and
v=kP1−1u.
Channel estimator module can store v in memory 75.
Channel estimator module 72 computes substream SNR based on
Wll1=1/(S1+S∥v1∥2−4real(S4v*1))
Wll2=1/(S2+S∥v2∥2−4real(S5v*2))
Wll3=1/(S3+S∥v3∥2−4real(S6v*3))
and store the substream SNR in memory 75.
Channel estimator module 72 can compute
z=R−1Q*y and μ=[1 −1 1]z
and read v from memory 75. Channel estimator module 72 can then compute
{circumflex over (x)}=P1−1z−μv
and read the substream SNR from memory 75. Equalizer module 76 can scale the equalized vector based on the SNRs.
Referring now to
The HDTV 420 may communicate with mass data storage 427 that stores data in a nonvolatile manner such as optical and/or magnetic storage devices. Mass data storage 427 may include at least one hard disk drive (HDD) and/or at least one digital versatile disk (DVD) drive. The HDD may be a mini HDD that includes one or more platters having a diameter that is smaller than approximately 1.8″. The HDTV 420 may be connected to memory 428 such as RAM, ROM, low latency nonvolatile memory such as flash memory and/or other suitable electronic data storage. The HDTV 420 also may support connections with a WLAN via WLAN network interface 429. The HDTV 420 also includes a power supply 423.
Referring now to
The receiver module may also be implemented in other control systems 440 of the vehicle 430. The control system 440 may likewise receive signals from input sensors 442 and/or output control signals to one or more output devices 444. In some implementations, the control system 440 may be part of an anti-lock braking system (ABS), a navigation system, a telematics system, a vehicle telematics system, a lane departure system, an adaptive cruise control system, a vehicle entertainment system such as a stereo, DVD, compact disc and the like. Still other implementations are contemplated.
The powertrain control system 432 may communicate with mass data storage 446 that stores data in a nonvolatile manner. Mass data storage 446 may include at least one HDD and/or at least one DVD drive. The HDD may be a mini HDD that includes one or more platters having a diameter that is smaller than approximately 1.8″. The powertrain control system 432 may be connected to memory 447 such as RAM, ROM, low latency nonvolatile memory such as flash memory and/or other suitable electronic data storage. The powertrain control system 432 also may support connections with a WLAN via a WLAN network interface 448. The control system 440 may also include mass data storage, memory and/or a WLAN interface (all not shown). Vehicle 430 may also include a power supply 433.
Referring now to
The cellular phone 450 may communicate with mass data storage 464 that stores data in a nonvolatile manner. Mass data storage 450 may include at least one HDD and/or at least one DVD drive. The HDD may be a mini HDD that includes one or more platters having a diameter that is smaller than approximately 1.8″. The cellular phone 450 may be connected to memory 466 such as RAM, ROM, low latency nonvolatile memory such as flash memory and/or other suitable electronic data storage. The cellular phone 450 also may support connections with a WLAN via the WLAN network interface 468. The cellular phone 450 may also include a power supply 453.
Referring now to
The set top box 480 may communicate with mass data storage 490 that stores data in a nonvolatile manner. Mass data storage 490 may include at least one HDD and/or at least one DVD drive. The HDD may be a mini HDD that includes one or more platters having a diameter that is smaller than approximately 1.8″. The set top box 480 may be connected to memory 494 such as RAM, ROM, low latency nonvolatile memory such as flash memory and/or other suitable electronic data storage. The set top box 480 also may support connections with a WLAN via the WLAN network interface 496. The set top box 480 may include a power supply 483.
Referring now to
The media player 500 may communicate with mass data storage 510 that stores data such as compressed audio and/or video content in a nonvolatile manner. In some implementations, the compressed audio files include files that are compliant with MP3 format or other suitable compressed audio and/or video formats. Mass data storage 510 may include at least one HDD and/or at least one DVD drive. The HDD may be a mini HDD that includes one or more platters having a diameter that is smaller than approximately 1.8″. The media player 500 may be connected to memory 514 such as RAM, ROM, low latency nonvolatile memory such as flash memory and/or other suitable electronic data storage. The media player 500 also may support connections with a WLAN via the WLAN network interface 516. The media player 500 may also include a power supply 513. Still other implementations in addition to those described above are contemplated.
Those skilled in the art can now appreciate from the foregoing description that the broad teachings of the disclosure can be implemented in a variety of forms. Therefore, while this disclosure includes particular examples, the true scope of the disclosure should not be so limited since other modifications will become apparent to the skilled practitioner upon a study of the drawings, the specification and the following claims.
The present disclosure is a continuation of U.S. patent application Ser. No. 13/725,167 (now U.S. Pat. No. 8,699,556), filed Dec. 21, 2012, which is a continuation of U.S. patent application Ser. No. 12/902,394 (now U.S. Pat. No. 8,340,169), filed on Oct. 12, 2010, which is a continuation of U.S. patent application Ser. No. 11/521,182 (now U.S. Pat. No. 7,813,421), filed on Sep. 14, 2006, which claims the benefit of U.S. Provisional Application No. 60/759,453, filed on Jan. 17, 2006. The entire disclosures of the above referenced applications are incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
6829312 | Laurila et al. | Dec 2004 | B1 |
7050513 | Yakhnich | May 2006 | B1 |
7239672 | Su | Jul 2007 | B2 |
7346129 | Gaikwad | Mar 2008 | B2 |
7978759 | Sarrigeorgidis | Jul 2011 | B1 |
20020027985 | Rashid-Farrokhi | Mar 2002 | A1 |
20040081074 | Piechocki | Apr 2004 | A1 |
20040091057 | Yoshida | May 2004 | A1 |
20050141657 | Maltsev et al. | Jun 2005 | A1 |
20050185738 | Gaikwad | Aug 2005 | A1 |
20050201454 | Chaudhuri et al. | Sep 2005 | A1 |
20050249303 | Takano | Nov 2005 | A1 |
20060105767 | Kim | May 2006 | A1 |
20060182206 | Coon et al. | Aug 2006 | A1 |
20070127608 | Scheim et al. | Jun 2007 | A1 |
20080112504 | Jiang et al. | May 2008 | A1 |
Number | Date | Country |
---|---|---|
1 566 917 | Aug 2005 | EP |
WO-2006004706 | Jan 2006 | WO |
Entry |
---|
Summary of Notice of Reasons for Rejection dated May 15, 2012 for Application No. 2008-;550466; 1 page. |
The International Preliminary Report on Patentability and Written Opinion of the International Searching Authority for corresponding International Application No. PCT/US2007/001233. |
Lo Tiing Tie et al., “Decoupled Maximum Likelihood Multiuser Channel Estimation for Asynchronous DS-CDMA System,” IEEE International Symposium on Spread Spectrum Techniques and Applications, Aug. 30, 2004, pp. 92-96, IEEE: Piscataway, NJ, ISBN: 0-;7803-8408-3. |
The Office Action dispatched Oct. 9, 2011, in corresponding Chinese; Application No. 200780002519.1. |
802.11 n; IEEE P802.11-0410889r6; Wireless LANs, TGn Sync Proposal Technical; Specification; 131 pages. May 2005. |
Number | Date | Country | |
---|---|---|---|
20140226703 A1 | Aug 2014 | US |
Number | Date | Country | |
---|---|---|---|
60759453 | Jan 2006 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 13725167 | Dec 2012 | US |
Child | 14252860 | US | |
Parent | 12902394 | Oct 2010 | US |
Child | 13725167 | US | |
Parent | 11521182 | Sep 2006 | US |
Child | 12902394 | US |