Physical MIMO Channel Modeling

Why Physical Models Matter

The i.i.d. Rayleigh model (H\mathbf{H} with independent CN(0,1)\mathcal{CN}(0,1) entries) is analytically convenient but unrealistic: real antenna arrays have finite spacing, real scattering environments have limited angular spread, and real channels exhibit spatial correlation. This section develops physically motivated channel models that capture these effects while remaining tractable for analysis.

Understanding correlation is crucial for system design: it determines whether a MIMO deployment achieves the promised capacity gains or falls far short.

Definition:

Spatial Correlation

The spatial correlation of a MIMO channel describes the statistical dependence between different entries of H\mathbf{H}. For a channel with zero-mean entries, the full correlation is captured by

R\ntnch=E[vec(H) vec(H)H]∈CnrntΓ—nrnt\mathbf{R}_{\ntn{ch}} = \mathbb{E}[\mathrm{vec}(\mathbf{H})\, \mathrm{vec}(\mathbf{H})^H] \in \mathbb{C}^{n_r n_t \times n_r n_t}

where vec(H)\mathrm{vec}(\mathbf{H}) stacks the columns of H\mathbf{H} into a single vector. When R\ntnch=Inrnt\mathbf{R}_{\ntn{ch}} = \mathbf{I}_{n_r n_t}, the channel entries are uncorrelated (i.i.d. model).

The full correlation matrix has nr2nt2n_r^2 n_t^2 entries, making it unwieldy for large arrays. Structured models like Kronecker and Weichselberger reduce this to a manageable number of parameters.

Definition:

Kronecker Correlation Model

The Kronecker model assumes that transmit and receive correlations are separable:

H=Rr1/2 Hw Rt1/2\mathbf{H} = \mathbf{R}_{r}^{1/2}\, \mathbf{H}_{w}\, \mathbf{R}_{t}^{1/2}

where:

  • Hw∈CnrΓ—nt\mathbf{H}_{w} \in \mathbb{C}^{n_r \times n_t} has i.i.d. CN(0,1)\mathcal{CN}(0,1) entries
  • Rr∈CnrΓ—nr\mathbf{R}_{r} \in \mathbb{C}^{n_r \times n_r} is the receive correlation matrix, [Rr]ij=E[hikhjkβˆ—][\mathbf{R}_{r}]_{ij} = \mathbb{E}[h_{ik} h_{jk}^*]
  • Rt∈CntΓ—nt\mathbf{R}_{t} \in \mathbb{C}^{n_t \times n_t} is the transmit correlation matrix, [Rt]ij=E[hkihkjβˆ—][\mathbf{R}_{t}]_{ij} = \mathbb{E}[h_{ki} h_{kj}^*]

The full correlation becomes R\ntnch=RtTβŠ—Rr\mathbf{R}_{\ntn{ch}} = \mathbf{R}_{t}^{T} \otimes \mathbf{R}_{r}, which is the Kronecker product of the two marginal correlations.

The Kronecker assumption is that the scattering environment "seen" from the transmit side is independent of that seen from the receive side. This holds well for moderate array sizes and rich scattering but can fail for keyhole channels or very large arrays.

,

Definition:

Weichselberger Model

The Weichselberger model generalises the Kronecker model by allowing coupling between transmit and receive eigenmodes:

H=Ur(Ξ©βŠ™G)UtH\mathbf{H} = \mathbf{U}_r (\boldsymbol{\Omega} \odot \mathbf{G}) \mathbf{U}_t^H

where:

  • Ur\mathbf{U}_r and Ut\mathbf{U}_t are the eigenvector matrices of Rr\mathbf{R}_{r} and Rt\mathbf{R}_{t} respectively
  • G\mathbf{G} has i.i.d. CN(0,1)\mathcal{CN}(0,1) entries
  • Ω∈Rβ‰₯0nrΓ—nt\boldsymbol{\Omega} \in \mathbb{R}_{\geq 0}^{n_r \times n_t} is the coupling matrix with [Ξ©]ij[\boldsymbol{\Omega}]_{ij} specifying the average power coupling between receive eigenmode ii and transmit eigenmode jj
  • βŠ™\odot denotes element-wise (Hadamard) product

The Kronecker model is the special case where Ξ©=Ξ»rΞ»tT\boldsymbol{\Omega} = \boldsymbol{\lambda}_r \boldsymbol{\lambda}_t^T (rank-1 coupling matrix).

Definition:

Keyhole Channel

A keyhole (pinhole) channel arises when all propagation passes through a single spatial bottleneck (e.g., a tunnel, a narrow gap between buildings). The channel has the structure

H=aratH\mathbf{H} = \mathbf{a}_r \mathbf{a}_t^H

where ar∈Cnr\mathbf{a}_r \in \mathbb{C}^{n_r} and at∈Cnt\mathbf{a}_t \in \mathbb{C}^{n_t} are random vectors with possibly rich distributions.

Despite potentially rich scattering at both the transmitter and receiver, the channel has rank 1 because all energy funnels through a single spatial mode. This means rank(H)=1\mathrm{rank}(\mathbf{H}) = 1 regardless of the number of antennas, and spatial multiplexing is impossible.

Keyhole channels are pathological for spatial multiplexing but still benefit from array gain (nt+nrn_t + n_r in dB) and diversity gain. They demonstrate that spatial multiplexing requires both multiple antennas AND a rich scattering environment.

Theorem: Capacity Loss from Spatial Correlation (Kronecker Model)

Under the Kronecker model, the ergodic capacity satisfies

Ccorr≀CiidC_{\mathrm{corr}} \leq C_{\mathrm{iid}}

with equality if and only if Rt=Int\mathbf{R}_{t} = \mathbf{I}_{n_t} and Rr=Inr\mathbf{R}_{r} = \mathbf{I}_{n_r}.

More precisely, for equal power allocation at high SNR:

Ccorrβ‰ˆCiidβˆ’log⁑2det⁑(Rt)βˆ’log⁑2det⁑(Rr)+constC_{\mathrm{corr}} \approx C_{\mathrm{iid}} - \log_2 \det(\mathbf{R}_{t}) - \log_2 \det(\mathbf{R}_{r}) + \text{const}

The capacity loss is governed by the "effective rank" of the correlation matrices: the more singular the correlation matrices, the greater the loss.

Spatial correlation reduces the effective number of independent channel coefficients. In the extreme case of full correlation (Rt\mathbf{R}_{t} or Rr\mathbf{R}_{r} rank-1), the MIMO channel degenerates to a single spatial stream regardless of the number of antennas.

Kronecker Correlated MIMO Channel

Visualise how transmit and receive spatial correlation affect the MIMO channel eigenvalue distribution and capacity. Compare with the i.i.d. baseline at various correlation levels.

Parameters
4

Number of transmit antennas

4

Number of receive antennas

0.5

Transmit correlation coefficient (exponential model)

0.5

Receive correlation coefficient (exponential model)

15

Keyhole Channel vs. Full-Rank Channel

Compare the eigenvalue spread and capacity of a keyhole (rank-1) channel with a full-rank i.i.d. Rayleigh channel. Observe how the keyhole channel collapses all spatial degrees of freedom into a single stream.

Parameters
4
4
20

Example: Exponential Correlation Model

A 4Γ—44 \times 4 MIMO system uses a uniform linear array (ULA) at both transmitter and receiver. The spatial correlation follows the exponential model:

[Rt]ij=ρt∣iβˆ’j∣,[Rr]ij=ρr∣iβˆ’j∣[\mathbf{R}_{t}]_{ij} = \rho_t^{|i-j|}, \qquad [\mathbf{R}_{r}]_{ij} = \rho_r^{|i-j|}

with ρt=0.7\rho_t = 0.7 and ρr=0.3\rho_r = 0.3.

(a) Compute the eigenvalues of Rt\mathbf{R}_{t} and Rr\mathbf{R}_{r}.

(b) How many "effective" spatial modes does each side contribute?

(c) Estimate the capacity loss relative to i.i.d. at SNR = 20 dB.

Quick Check

In the Kronecker MIMO channel model H=Rr1/2 Hw Rt1/2\mathbf{H} = \mathbf{R}_{r}^{1/2}\, \mathbf{H}_{w}\, \mathbf{R}_{t}^{1/2}, what happens to the channel when Rt\mathbf{R}_{t} has rank 1?

The MIMO channel degenerates to rank 1 regardless of Rr\mathbf{R}_{r} or nrn_r

The channel is unaffected because Rr\mathbf{R}_{r} compensates

The channel becomes full-rank because the noise fills in the missing modes

The capacity doubles because all power is focused on one strong mode

Common Mistake: The Kronecker Model Can Underestimate or Overestimate Capacity

Mistake:

Blindly applying the Kronecker model to all MIMO scenarios, especially with large angular spreads or when the transmitter and receiver "see" overlapping clusters of scatterers.

Correction:

The Kronecker model assumes separability of transmit and receive correlation, which fails when the coupling between transmit and receive eigenmodes is not rank-1. For example, measured indoor channels at 5.2 GHz show 2--4 dB capacity overestimation by the Kronecker model. Use the Weichselberger model or full correlation R\ntnch\mathbf{R}_{\ntn{ch}} when accuracy matters. Always validate against measurements for your deployment scenario.

Why This Matters: Spatial Correlation in Massive MIMO

In 5G NR massive MIMO (64--256 antenna elements at the base station), spatial correlation is not a bug but a feature. The large array aperture creates highly directional beams, and the channel becomes low-rank in the angular domain. This is exploited by:

  • Codebook-based precoding (Type I/II in 5G NR) that matches the dominant eigenmodes
  • FDD channel estimation that exploits low-rank structure to reduce feedback overhead from O(nt)\mathcal{O}(n_t) to O(r)\mathcal{O}(r)
  • Pilot contamination mitigation that uses angle-of-arrival separation between users

The Weichselberger model is particularly useful for massive MIMO because the coupling matrix Ξ©\boldsymbol{\Omega} directly captures the angular power spectrum.

See full treatment in Chapter 16

Comparison of MIMO Channel Models

Propertyi.i.d. RayleighKroneckerWeichselbergerKeyhole
Correlation structureNone (uncorrelated)Separable Tx/RxFull eigenmode couplingRank-1
Parameters0nt2+nr2n_t^2 + n_r^2nt2+nr2+ntnrn_t^2 + n_r^2 + n_t n_rnt+nrn_t + n_r
Typical rankmin⁑(nt,nr)\min(n_t, n_r)≀min⁑(nt,nr)\leq \min(n_t, n_r)≀min⁑(nt,nr)\leq \min(n_t, n_r)1
Best forAnalysis baselineModerate correlationMeasured channelsTunnel/corridor
Capacity accuracyUpper boundGood (Β±2\pm 2 dB)ExcellentExact for rank-1

Kronecker correlation model

A MIMO channel model that separates transmit and receive correlation: H=Rr1/2HwRt1/2\mathbf{H} = \mathbf{R}_{r}^{1/2} \mathbf{H}_{w} \mathbf{R}_{t}^{1/2}, assuming the correlation at one end is independent of the other.

Related: Spatial correlation, Weichselberger model

Spatial correlation

The statistical dependence between channel gains at different antenna elements, arising from limited angular spread, antenna spacing, and scatterer geometry.

Related: Kronecker correlation model, Channel matrix

Weichselberger model

A generalisation of the Kronecker model that captures coupling between transmit and receive eigenmodes via a non-negative coupling matrix Ξ©\boldsymbol{\Omega}.

Related: Kronecker correlation model

Keyhole channel

A MIMO channel where all propagation passes through a single spatial bottleneck, resulting in rank-1 channel matrix H=aratH\mathbf{H} = \mathbf{a}_r \mathbf{a}_t^H regardless of the number of antennas.

Related: Channel rank, Spatial multiplexing