Multiple-input-multiple-output (MIMO) channel estimation

Abstract

Method for propagation channel estimation in an NM Multiple-Input-Multiple-Output (MIMO) communication system comprising a transmitter (202) comprising N transmit antennas and a receiver (204) comprising M receive antennas and a MIMO equalizer (206) comprising multiple taps, where N>1 and M>1. The method includes producing (s402) a single tap equalizer (Q) based on a multi-tap equalizer (Q). The method also includes producing (s404) an inverse effective-channel estimate (Q.sub.e) based on Q, The method also includes inverting (s406) Q.sub.e to produce an effective-channel estimate (H.sub.e), The method also includes producing (s408) H.sub.a based on H.sub.ewherein H.sub.a can be used to determine one or more performance metrics.

Claims

1. A method for propagation channel estimation in an NM Multiple-Input-Multiple-Output (MIMO) communication system comprising a transmitter comprising N transmit antennas and a receiver comprising M receive antennas and a MIMO equalizer comprising multiple taps, where N>1 and M>1, the method comprising: producing a single tap equalizer (Q) based on a multi-tap equalizer (Q); producing an inverse effective-channel estimate (Q.sub.e) based on Q; inverting Q.sub.e to produce an effective-channel estimate (H.sub.e); and producing H.sub.a based on H.sub.e, wherein H.sub.a can be used to determine one or more performance metrics.

2. The method of claim 1, wherein producing Q comprises: computing Q.sub.F=(Q), where is an element-wise Fourier transform.

3. The method of claim 2, wherein producing Q further comprises: setting $\stackrel{\_}{Q}=\left(\frac{1}{N_u-N_l+1}{.Math.}_{k=N_l}^{N_u}.Math.Q_F\left[k\right].Math.\right)e^{j\frac{1}{N_u-N_l+1}{.Math.}_{k=N_l}^{N_u}\left(Q_F\left[k\right]\right)},$ where () is the angle operation, O.sub.F is an oversampling factor, N.sub.taps is the number of tapes of the multi-tap equalizer (Q), and N.sub.l and N.sub.u are predetermined values.

4. The method of claim 1, wherein producing Q.sub.e based on Q comprises producing Q.sub.e by combining Q with phase shifters.

5. The method of claim 1, wherein producing Q.sub.e based on Q comprises calculating:
Q.sub.e=Q, where =a phase noise compensating matrix, and is the element-wise Hadamard product.

6. The method of claim 1, wherein inverting Q.sub.e to produce the effective-channel estimate (H.sub.e) comprises setting H.sub.e equal to: (Q.sub.e*(Q.sub.eQ.sub.e*).sup.1, where ()* is the conjugate-transpose operation.

7. The method of claim 1, wherein producing H.sub.a based on H.sub.e comprises compensating for an automatic gain control (AGC) imbalance using H.sub.e, thereby producing H.sub.a.

8. The method of claim 1, wherein the transmitter comprises N automatic-gain-controllers (AGCs), the receiver comprises M AGCs, producing H.sub.a based on H.sub.e comprises setting H.sub.a equal to: P.sub.RX.sup.1H.sub.e P.sub.TX.sup.1, where P.sub.RX is a first matrix representing gains applied by the M AGCs of the receiver, and P.sub.TX is a second matrix representing gains applied by the N AGCs of the transmitter.

9. The method of claim 1, further comprising calculating one or more performance metrics using H.sub.a.

10. The method of claim 9, wherein calculating one or more performance metrics using H.sub.a comprises calculating a mean-square-error (MSE) using H.sub.a as input to the calculation.

11. The method of claim 9, wherein calculating one or more performance metrics using H.sub.a comprises determining: the singular values of H.sub.a, a condition number associated with H.sub.a, cross-polarization isolation, XPI, a power gain, and/or carrier-to-interference ratio between a MIMO main stream and one or more interference streams.

12. The method of claim 1, wherein the receiver is within the line-of-sight (LOS) of the transmitter.

13. The method of claim 1, wherein the method is performed by a baseband processor of the receiver.

14. The method of claim 1, wherein the method is performed by a channel estimator of a network node.

15. The method of claim 14, wherein the network node is located remotely from the receiver.

16. The method of claim 3, wherein $N_l=.Math.\frac{N_{\mathrm{taps}}\left(O_F-1\right)}{{\mathrm{cO}}_F}.Math.,N_u=.Math.\frac{N_{\mathrm{taps}}\left(O_F+1\right)}{{\mathrm{cO}}_F}.Math.,$ and c is a predetermined value.

17. A non-transitory computer readable storage medium storing a computer program comprising instructions which when executed by processing circuitry of a channel estimator causes the channel estimator to perform the method of claim 1.

18. A channel estimator, the channel estimator comprising: processing circuitry; and a memory, the memory containing instructions executable by the processing circuitry, whereby the channel estimator is operative to perform a method comprising: producing a single tap equalizer (Q) based on a multi-tap equalizer (Q); producing an inverse effective-channel estimate (Q.sub.e) based on Q; inverting Q.sub.e to produce an effective-channel estimate (H.sub.e); and producing H.sub.a based on H.sub.e, wherein H.sub.a can be used to determine one or more performance metrics.

19. A network node comprising the channel estimator of claim 18.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The accompanying drawings, which are incorporated herein and form part of the specification, illustrate various embodiments.

(2) FIG. 1 illustrates a 44 LOS-MIMO communication system according to an embodiment.

(3) FIG. 2 illustrates a 22 MIMO communication system according to an embodiment.

(4) FIG. 3 illustrates a MIMO communication system including a channel estimator according to an embodiment.

(5) FIG. 4 is a flowchart illustrating a process according to an embodiment.

(6) FIG. 5 illustrates a channel estimator according to an embodiment.

DETAILED DESCRIPTION

(7) FIG. 1 illustrates a LOS-MIMO communication system 100 according to some embodiments. As shown in FIG. 1, LOS-MIMO communication system 100 comprises is a 44 LO-MIMO system (i.e., the MIMO communication system has four transmitters and 4 receivers). As noted above, the performance of the MIMO link highly depends on the propagation channel, and in order to assess the performance of the MIMO link, the propagation channel (or channel for short) is estimated and represented using a matrix. This matrix can in turn be used to estimate the performance of the MIMO link, assess propagation conditions, set digital-signal-processing (DSP) parameters, etc. For example, the matrix can be sent to a central node as input for network-level traffic planning and/or routing planning.

(8) The channel can be estimated using a MIMO equalizer, which is used to de-multiplex the MIMO streams (i.e., invert the propagation channel and remove the interference from the other streams). In an NM MIMO system, the MIMO equalizer is composed of NM finite impulse response filters (FIR) and NM phase shifters. The FIR filters compensate for the MIMO channel and for various hardware imperfections such as frequency domain ripple. The phase shifters compensate for the combined phase noise and frequency offset coming from the transmit (TX) and receive (RX) local oscillators.

(9) The FIR filters and the phase shifters are dynamically adjusted during transmission such that they adapt to the propagation channel. They are typically adjusted based on an error function calculated at a slicer in the DSP chain based on various algorithms, such as, least mean square gradient descent, Kalman filters, etc. In most of the cases, the LOS-MIMO propagation channel does not experience any multipath, therefore the FIR filters of the MIMO equalizer can be single-tap filters, but they are chosen to be multitap in order to compensate for hardware imperfections.

(10) FIG. 2 shows the schematic for an example 22 MIMO communication system 200 having a transmitter 202, a receiver 204, a MIMO equalizer 206 coupled to the output of receiver 204, and two slicers. The transmitter 202 comprises a first AGC (AGC TX1), a second AGC (AGC TX2), a first radio RX unit coupled to AGC TX, and a second radio RX unit coupled to AGC TX2. Similarly, the receiver 204 comprises a first AGC (AGC RX1), a second AGC (AGX RX2), a first radio RX unit coupled to AGC RX1, and a second radio RX unit coupled to ACX RX2; 3) The purpose of the AGCs is to ensure that the transmitted/received signal power level is kept constant at a certain level. In this example, receiver 204 is within the LOS of transmitter 202.

(11) In this example, MIMO equalizer 206 is made of 4 FIR filters and 4 phase shifters, which are adjusted based on a control signal coming from the respective slicer. This implementation adjusts _mn and q_mn based on the feedback signal only from the slicer n. This implementation allows the different receiver chains to be implemented and processed in parallel, which reduces the requirements on the FPGA/ASIC implementations.

(12) The signal before the MIMO equalizer 206 (denoted Y) can be modelled as:

(13) $Y=\left[\begin{array}{c}{y}_1\\ .Math.\\ y_M\end{array}\right]=G_{RX}{P}_{RX}H{P}_{TX}{G}_{TX}X+B,$
where

(14) $X=\left[\begin{array}{c}{x}_1\\ .Math.\\ x_N\end{array}\right]$
is the transmitted signal,

(15) $G_{\mathrm{TX}}=\left[\begin{array}{ccc}{G}_{{\mathrm{TX}}_1}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& G_{{\mathrm{TX}}_N}\end{array}\right]$
is the AGC gain at the transmitter,

(16) $P_{\mathrm{TX}}=\left[\begin{array}{ccc}{e}^{j{}_1}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& e^{j{}_N}\end{array}\right]$
is the phase noise at the transmitter,

(17) $H=\left[\begin{array}{ccc}{h}_{11}& .Math.& h_{1N}\\ .Math.& & .Math.\\ h_{M1}& .Math.& h_{\mathrm{MN}}\end{array}\right]$
is the MIMO propagation channel,

(18) $P_{\mathrm{RX}}=\left[\begin{array}{ccc}{e}^{j{}_1}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& e^{j{}_M}\end{array}\right]$
is the phase noise at the receiver,

(19) $G_{\mathrm{RX}}=\left[\begin{array}{ccc}{G}_{{\mathrm{RX}}_1}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& G_{{\mathrm{RX}}_M}\end{array}\right]$
is the AGC gain at the receiver, and

(20) $B=\left[\begin{array}{c}{b}_1\\ .Math.\\ b_M\end{array}\right]$
is the additive noise added at the receiver. The output of the equalizer is

(21) $Z=\left[\begin{array}{c}{z}_1\\ .Math.\\ z_N\end{array}\right]=\left(Q\right)*Y,$
where

(22) 0$Q=\left[\begin{array}{ccc}{q}_{11}\left[k\right]& .Math.& q_{1M}\left[k\right]\\ .Math.& & .Math.\\ q_{N1}\left[k\right]& .Math.& q_{\mathrm{NM}}\left[k\right]\end{array}\right]$
is the filter bank with the FIR filters of size N.sub.taps and k=1 . . . N.sub.taps with a sampling rate F.sub.s=F.sub.BO.sub.F, where F.sub.B is the baud rate of the signal and O.sub.F is the oversampling factor,

(23) $=\left[\begin{array}{ccc}{e}^{j{}_{11}}& .Math.& e^{j{}_{1M}}\\ .Math.& & .Math.\\ e^{j{}_{N1}}& .Math.& e^{j{}_{\mathrm{NM}}}\end{array}\right]$
is a matrix compensating for the phase noise, and is the element-wise Hadamard product.

(24) Note that both Q and converge to inverse of true channel and phase noise.

(25) Typically, the effective channel is estimated from the taps of the MIMO equalizer 206. Considering single-tap FIR filters, the channel matrix is estimated as =Q.sup.1. However, there are several issues with this approach: i) it is not straight forward how to apply the inverse operation when Q is a multi-tap filter bank; and ii) the estimate is accurate only when there is no phase noise and the gain of the AGCs is the same on all branches.

(26) This disclosure provides a channel estimator 302 (see FIG. 3) to overcome the limitations discussed above, which allows one to quantify the quality of the channel, such as: singular values, condition number, cross-polarization interference (XPI), power gain/attenuation during propagation. As shown in FIG. 3, channel estimator 302, in this embodiments, receives information indicating the gain of each of AGCs. Channel estimator 302 may be a component of a baseband processor for the receiver. In another embodiment, channel estimator 302 is a component of a network node 301, which network node 301 may be located remotely from receiver 204 and/or equalizer 206.

(27) FIG. 4 is a flow chart illustrating a process 400 performed by a channel estimator (e.g., channel estimator 302) for propagation channel estimation in an NM (e.g., 22 or 44) MIMO communication system comprising a transmitter comprising N transmit antennas and a receiver comprising M receive antennas and MIMO equalizer 206 comprising multiple taps, where N>1 and M>1.

(28) Process 400 may begin in steps s402. Step s402 comprises producing a single tap equalizer (Q) based on a multi-tap equalizer (Q). Step s404 comprises producing an inverse effective-channel estimate (Q.sub.e) based on Q. Step s406 comprises inverting Q.sub.e to produce an effective-channel estimate (H.sub.e). And step s408 comprises producing H.sub.a based on H.sub.e, wherein H.sub.a can be used to determine one or more performance metrics (i.e., H.sub.a represents the propagation channel). That is, the estimate H.sub.a is the estimate of the combination between the propagation channel and phase noise from the TX and RX local oscillators (LOs); therefore, it cannot be used to examine the propagation channel (H) directly, but it shares some of the properties of H such as: singular values, condition number, XPI, power gain/attenuation during propagation, which can be used to, e.g., estimate the performance of the link, assess propagation conditions, set DSP parameters, etc.

(29) In some embodiments, prior to performing process 400, channel estimator 302 first determines whether the MIMO equalizer has reached convergence based on, for example, the bit-error-rate (BER), mean-square-error (MSE), or packer-error-rate (PER). In such an embodiment, channel estimator 302 performs process 400 in response to determining that the MIMO equalizer has reached convergence.

(30) In some embodiments, producing Q comprises: computing Q.sub.F=(Q), where is an element-wise Fourier transform, and then setting:

(31) $\stackrel{\_}{Q}=\left(\frac{1}{N_u-N_l+1}{.Math.}_{k=N_l}^{N_u}.Math.Q_F\left[k\right].Math.\right)e^{j\frac{1}{N_u-N_l+1}{.Math.}_{k=N_l}^{N_u}\left(Q_F\left[k\right]\right)},$
where () is the angle operation, O.sub.F is an oversampling factor, N.sub.taps is the number of tapes of the multi-tap equalizer (Q), and N.sub.l and N.sub.u are choses based on, for example, the sampling rate. In one embodiment,

(32) $N_l=.Math.\frac{N_{\mathrm{taps}}\left(O_F-1\right)}{{\mathrm{cO}}_F}.Math.,$

(33) $N_u=.Math.\frac{N_{\mathrm{taps}}\left(O_F+1\right)}{{\mathrm{cO}}_F}.Math.,$
where c is a predetermined constant (eg., c=2, c=3, etc.)

(34) In some embodiments, producing Q.sub.e based on Q comprises producing Q.sub.e by combining Q with phase shifters. For example, producing Q.sub.e based on Q may comprise calculating:
Q.sub.e=Q, where =a phase noise compensating matrix, and is the element-wise Hadamard product.

(35) In some embodiments, inverting Q.sub.e to produce the effective-channel estimate (H.sub.e) comprises setting H.sub.e equal to: (Q.sub.e*(Q.sub.eQ.sub.e*).sup.1, where ()* is the conjugate-transpose operation.

(36) In some embodiments, producing H.sub.a based on H.sub.e comprises compensating for an automatic gain control (AGC) imbalance using H.sub.e, thereby producing H.sub.a.

(37) In some embodiments, producing H.sub.a based on H.sub.e comprises setting H.sub.a equal to: P.sub.RX.sup.1H.sub.eP.sub.TX.sup.1, where P.sub.RX is a first matrix representing gains applied by the M AGCs of the receiver, and P.sub.TX is a second matrix representing gains applied by the N AGCs of the transmitter. In one embodiment,

(38) $P_{\mathrm{TX}}^{-1}=\left[\begin{array}{ccc}\frac{1}{G_{{\mathrm{TX}}_1}}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& \frac{1}{G_{{\mathrm{TX}}_N}}\end{array}\right],P_{\mathrm{RX}}^{-1}=\left[\begin{array}{ccc}\frac{1}{G_{{\mathrm{RX}}_1}}& .Math.& 0\\ .Math.& & .Math.\\ 0& .Math.& \frac{1}{G_{{\mathrm{RX}}_M}}\end{array}\right],$ Gtx.sub.i is the gain of the i-th AGC of the transmitter, and Grx.sub.i is the gain of the i-th AGC of the receiver.

(39) In some embodiments, process 400 further includes calculating one or more performance metric using H.sub.a. In some embodiments, calculating one or more performance metric using H.sub.a comprises calculating a mean-square-error (MSE) for a slice using H.sub.a as input to the calculation (i.e., the mean squared error between H.sub.a.sup.1*yx, where y are the received samples and x are the closest constellation points to H.sub.a.sup.1*y. In some embodiments, calculating one or more performance metric using H.sub.a comprises determining one or more of: the singular values of H.sub.a, a condition number associated with H.sub.a, cross-polarization isolation, XPI, a power gain, a carrier-to-interference ratio. In one embodiment, calculating XPI comprises calculating: the ratio between abs(H.sub.a(n_x, m_x)).sup.2 and abs(H.sub.a(n_y, m_y)).sup.2, where H.sub.a(n_x, m_x) is the estimated channel gain from TX antenna n and RX antenna m on polarization x, whereas H.sub.a(n_y, m_y) is the estimated channel gain from TX antenna n and RX antenna m on polarization y.

(40) FIG. 5 is a block diagram of channel estimator 302, according to some embodiments. As shown in FIG. 5, channel estimator 302 may comprise: processing circuitry (PC) 502, which may include one or more processors (P) 555 (e.g., one or more general purpose microprocessors and/or one or more other processors, such as an application specific integrated circuit (ASIC), field-programmable gate arrays (FPGAs), and the like), which processors may be co-located in a single housing or in a single data center or may be geographically distributed (i.e., channel estimator 302 may be a distributed computing apparatus); and a local storage unit (a.k.a., data storage system) 508, which may include one or more non-volatile storage devices and/or one or more volatile storage devices. Channel estimator may also include at least one network interface 548 comprising a transmitter (Tx) 545 and a receiver (Rx) 547 for enabling channel estimator 302 to transmit data to and receive data from other nodes connected to a network 110 (e.g., an Internet Protocol (IP) network) to which network interface 548 is connected.

(41) In embodiments where PC 502 includes a programmable processor, a computer program product (CPP) 541 may be provided. CPP 541 includes a computer readable medium (CRM) 542 storing a computer program (CP) 543 comprising computer readable instructions (CRI) 544. CRM 542 may be a non-transitory computer readable medium, such as, magnetic media (e.g., a hard disk), optical media, memory devices (e.g., random access memory, flash memory), and the like. In some embodiments, the CRI 544 of computer program 543 is configured such that when executed by PC 502, the CRI causes channel estimator 302 to perform steps described herein (e.g., steps described herein with reference to the flow charts). In other embodiments, channel estimator 302 may be configured to perform steps described herein without the need for code. That is, for example, PC 502 may consist merely of one or more ASICs. Hence, the features of the embodiments described herein may be implemented in hardware and/or software.

(42) While various embodiments are described herein, it should be understood that they have been presented by way of example only, and not limitation. Thus, the breadth and scope of this disclosure should not be limited by any of the above-described exemplary embodiments. Moreover, any combination of the above-described elements in all possible variations thereof is encompassed by the disclosure unless otherwise indicated herein or otherwise clearly contradicted by context.

(43) Additionally, while the processes described above and illustrated in the drawings are shown as a sequence of steps, this was done solely for the sake of illustration. Accordingly, it is contemplated that some steps may be added, some steps may be omitted, the order of the steps may be re-arranged, and some steps may be performed in parallel.

REFERENCES

(44) [1] Larsson, P., Lattice array receiver and sender for spatially orthonormal MIMO communication, in 2005 IEEE 61st Vehicular Technology Conference, 2005, pp. 192-196 Vol. 1. [2] Sjodin, M., et. al., A 40.2 bps/Hz Single Polarization 44 Line-of-Sight MIMO Link With Unsynchronized Oscillators, in 2019 IEEE Radio and Wireless Symposium (RWS), 2019, pp. 1-3.