MIMO Radar and ISAC Beamforming

MIMO Radar: From Phased Arrays to Waveform Diversity

Classical phased-array radar transmits a coherent waveform from all antennas, forming a narrow beam in a single direction. The beam must be mechanically or electronically steered to scan the angular domain, limiting the achievable search rate. MIMO radar fundamentally changes this paradigm by transmitting orthogonal waveforms from each antenna, illuminating the entire angular sector simultaneously. The orthogonality enables the receiver to separate the contribution of each transmit antenna and form a virtual aperture that is much larger than the physical array.

Consider a MIMO radar with NTN_T transmit antennas and NRN_R receive antennas. Each transmit antenna ii transmits a unique waveform si(t)s_i(t), designed to be orthogonal to the waveforms of all other antennas:

∫si(t)sjβˆ—(t) dt={Esi=j0iβ‰ j\int s_i(t) s_j^*(t) \, dt = \begin{cases} E_s & i = j \\ 0 & i \neq j \end{cases}

Orthogonality can be achieved through time-division (sequential transmission), frequency-division (non-overlapping subbands), code-division (orthogonal codes), or Doppler-division (slow-time coding with interleaved phase shifts).

ISAC System Architecture

ISAC System Architecture
Architecture of a dual-function MIMO ISAC system. The transmitter uses NTN_T antennas to simultaneously serve communication users in the forward (downlink) direction and illuminate radar targets. The transmitted signal x(t)=Wcd(t)+Wss(t)\mathbf{x}(t) = \mathbf{W}_{c} \mathbf{d}(t) + \mathbf{W}_{s} \mathbf{s}(t) combines a communication precoder Wc\mathbf{W}_{c} applied to data streams d(t)\mathbf{d}(t) with a sensing precoder Ws\mathbf{W}_{s} applied to radar probing waveforms s(t)\mathbf{s}(t). The radar receiver processes the backscattered echoes using NRN_R receive antennas, while communication users decode their intended data from the forward-link signal. The trade-off parameter ρ\rho controls the power allocation between sensing and communication.

Definition:

MIMO Radar

A MIMO radar system employs NTN_T transmit antennas, each transmitting an independent, orthogonal waveform, and NRN_R receive antennas. The signal model for a single point target at angle ΞΈ\theta with complex reflectivity Ξ±\alpha is:

y(t)=α aR(ΞΈ) aTT(ΞΈ) s(tβˆ’Ο„) ej2Ο€fdt+n(t)\mathbf{y}(t) = \alpha \, \mathbf{a}_R(\theta) \, \mathbf{a}_T^T(\theta) \, \mathbf{s}(t - \tau) \, e^{j2\pi f_d t} + \mathbf{n}(t)

where:

  • aT(ΞΈ)∈CNT\mathbf{a}_T(\theta) \in \mathbb{C}^{N_T}: transmit array steering vector
  • aR(ΞΈ)∈CNR\mathbf{a}_R(\theta) \in \mathbb{C}^{N_R}: receive array steering vector
  • s(t)=[s1(t),…,sNT(t)]T\mathbf{s}(t) = [s_1(t), \ldots, s_{N_T}(t)]^T: vector of orthogonal transmit waveforms
  • n(t)∼CN(0,Οƒ2I)\mathbf{n}(t) \sim \mathcal{CN}(\mathbf{0}, \sigma^2 \mathbf{I}): noise vector

After matched filtering to separate the waveforms, the NRΓ—NTN_R \times N_T virtual data matrix for the target is:

Z=α aR(ΞΈ) aTT(ΞΈ)+N\mathbf{Z} = \alpha \, \mathbf{a}_R(\theta) \, \mathbf{a}_T^T(\theta) + \mathbf{N}

Vectorising: vec(Z)=α (aT(ΞΈ)βŠ—aR(ΞΈ))+vec(N)\mathrm{vec}(\mathbf{Z}) = \alpha \, (\mathbf{a}_T(\theta) \otimes \mathbf{a}_R(\theta)) + \mathrm{vec}(\mathbf{N}), where aT(ΞΈ)βŠ—aR(ΞΈ)\mathbf{a}_T(\theta) \otimes \mathbf{a}_R(\theta) is the virtual steering vector of dimension NTNRN_T N_R.

The key distinction from phased-array radar is that MIMO radar does not form a transmit beam; instead, it illuminates the entire scene and achieves spatial resolution through the virtual aperture at the receiver.

Two MIMO radar paradigms exist: (i) widely separated antennas (statistical MIMO), where spatial diversity improves target detection against fading, and (ii) colocated antennas (waveform-diversity MIMO), where the virtual aperture improves angular resolution. This chapter focuses on the colocated case, which is most relevant to ISAC systems with compact antenna arrays.

Theorem: Virtual Aperture of Colocated MIMO Radar

A colocated MIMO radar with NTN_T transmit antennas at positions {dt(i)}i=1NT\{d_t^{(i)}\}_{i=1}^{N_T} and NRN_R receive antennas at positions {dr(j)}j=1NR\{d_r^{(j)}\}_{j=1}^{N_R} creates a virtual array with element positions:

Dv={dt(i)+dr(j):i=1,…,NT,β€…β€Šj=1,…,NR}\mathcal{D}_v = \{d_t^{(i)} + d_r^{(j)} : i = 1, \ldots, N_T, \; j = 1, \ldots, N_R\}

The virtual array has up to NTβ‹…NRN_T \cdot N_R distinct elements. For a transmit ULA with spacing dTd_T and a receive ULA with spacing dR=dT/NTd_R = d_T / N_T (or equivalently, dT=NRdRd_T = N_R d_R):

  • The virtual array is a filled ULA with NTNRN_T N_R elements and spacing dRd_R.
  • The virtual aperture is Lv=(NTNRβˆ’1)dRL_v = (N_T N_R - 1) d_R, which is NTN_T times the physical receive aperture.
  • The angular resolution is: Ξ”ΞΈβ‰ˆΞ»Lv=Ξ»(NTNRβˆ’1)dR\Delta\theta \approx \frac{\lambda}{L_v} = \frac{\lambda}{(N_T N_R - 1) d_R}

This virtual aperture provides the angular resolution of an array with NTNRN_T N_R physical elements, using only NT+NRN_T + N_R physical antennas.

Joint Communication-Sensing Beamforming

In a MIMO ISAC system, the transmit signal must simultaneously form communication beams towards users and a radar beam towards the target. Consider a system with NTN_T transmit antennas serving LL single-antenna communication users and sensing one target at angle ΞΈs\theta_s. The transmitted signal is:

x(t)=Wcd(t)⏟communication+wss(t)⏟sensing\mathbf{x}(t) = \underbrace{\mathbf{W}_{c} \mathbf{d}(t)}_{\text{communication}} + \underbrace{\mathbf{w}_s s(t)}_{\text{sensing}}

where Wc=[w1,…,wL]∈CNTΓ—L\mathbf{W}_{c} = [\mathbf{w}_1, \ldots, \mathbf{w}_L] \in \mathbb{C}^{N_T \times L} is the communication precoding matrix, d(t)=[d1(t),…,dL(t)]T\mathbf{d}(t) = [d_1(t), \ldots, d_L(t)]^T is the data vector, ws∈CNT\mathbf{w}_s \in \mathbb{C}^{N_T} is the sensing beamforming vector, and s(t)s(t) is the radar probing waveform.

The transmit covariance matrix is:

Rx=E[xxH]=WcWcH+wswsH\mathbf{R}_{x} = \mathbb{E}[\mathbf{x}\mathbf{x}^H] = \mathbf{W}_{c} \mathbf{W}_{c}^{H} + \mathbf{w}_s \mathbf{w}_s^H

subject to the total power constraint tr(Rx)≀P\mathrm{tr}(\mathbf{R}_{x}) \leq P.

Communication metric: The SINR for user ll is:

SINRl=∣hlHwl∣2βˆ‘lβ€²β‰ l∣hlHwlβ€²βˆ£2+∣hlHws∣2+Οƒc2\mathrm{SINR}_l = \frac{|\mathbf{h}_l^H \mathbf{w}_l|^2}{\sum_{l' \neq l} |\mathbf{h}_l^H \mathbf{w}_{l'}|^2 + |\mathbf{h}_l^H \mathbf{w}_s|^2 + \sigma_c^2}

Sensing metric: The radar beampattern gain in direction ΞΈ\theta is:

P(ΞΈ)=aH(ΞΈ)Rxa(ΞΈ)P(\theta) = \mathbf{a}^H(\theta) \mathbf{R}_{x} \mathbf{a}(\theta)

The CRB for angle estimation of a target at ΞΈs\theta_s is inversely proportional to P(ΞΈs)P(\theta_s) and to the second derivative of the beampattern at ΞΈs\theta_s.

The joint optimisation problem is:

max⁑Wc,wsβ€…β€Šβˆ‘l=1Llog⁑2(1+SINRl)\max_{\mathbf{W}_{c}, \mathbf{w}_s} \; \sum_{l=1}^L \log_2(1 + \mathrm{SINR}_l) s.t.aH(ΞΈs)Rxa(ΞΈs)β‰₯Ξ“s\text{s.t.} \quad \mathbf{a}^H(\theta_s) \mathbf{R}_{x} \mathbf{a}(\theta_s) \geq \Gamma_s tr(Rx)≀P\mathrm{tr}(\mathbf{R}_{x}) \leq P

where Ξ“s\Gamma_s is the minimum required sensing beampattern gain. The Pareto frontier of this problem β€” the set of all achievable (rate, CRB) pairs β€” characterises the fundamental communication-sensing trade-off.

A common parametric approach uses a trade-off weight ρ∈[0,1]\rho \in [0, 1]:

max⁑Rxβͺ°0β€…β€Š(1βˆ’Ο)β‹…fcomm(Rx)+ρ⋅fsense(Rx)\max_{\mathbf{R}_{x} \succeq \mathbf{0}} \; (1-\rho) \cdot f_{\text{comm}}(\mathbf{R}_{x}) + \rho \cdot f_{\text{sense}}(\mathbf{R}_{x}) s.t.tr(Rx)≀P\text{s.t.} \quad \mathrm{tr}(\mathbf{R}_{x}) \leq P

When ρ=0\rho = 0, all resources are devoted to communication (maximum rate, no sensing constraint); when ρ=1\rho = 1, all resources are devoted to sensing (optimal radar beampattern, no communication). Intermediate values trace out the Pareto frontier.

Structure of the Optimal ISAC Beamformer

The optimal solution to the ISAC beamforming problem has a revealing structure. When the communication channel vectors {hl}\{\mathbf{h}_l\} and the sensing steering vector a(ΞΈs)\mathbf{a}(\theta_s) are not aligned (the typical case), the optimal transmit covariance can be decomposed as:

Rxβˆ—=βˆ‘l=1LplwlwlH⏟comm.Β beams+pswswsH⏟sensingΒ beam\mathbf{R}_{x}^* = \underbrace{\sum_{l=1}^L p_l \mathbf{w}_l \mathbf{w}_l^H}_{\text{comm. beams}} + \underbrace{p_s \mathbf{w}_s \mathbf{w}_s^H}_{\text{sensing beam}}

where wl\mathbf{w}_l is a regularised zero-forcing beamformer for user ll (balancing multi-user interference suppression with sensing beam compatibility) and ws\mathbf{w}_s concentrates energy towards the target direction in the null space of the communication channels.

When the target happens to lie in the same direction as a communication user (ΞΈsβ‰ˆΞΈl\theta_s \approx \theta_l), the communication beam to that user simultaneously serves as the sensing beam, and the trade-off vanishes β€” both functions benefit from the same spatial resources. This synergy between communication and sensing directions is a key benefit of ISAC in scenarios like V2X, where the vehicle communicates with another vehicle that is also a radar target.

πŸŽ“CommIT Contribution(2025)

Deterministic-Random Tradeoff for ISAC

F. Liu, G. Caire β€” IEEE Trans. Inform. Theory / IEEE Trans. Commun.

Liu and Caire established the fundamental information-theoretic tradeoff between sensing and communication in ISAC systems. Their key insight is that the communication function requires the transmitted signal to carry random information (entropy), while the sensing function benefits from deterministic (known) waveforms that maximise the estimation accuracy (minimise the CRB).

The deterministic-random tradeoff characterises the Pareto frontier of achievable (rate, CRB) pairs:

  • When the entire signal is deterministic (R=0R = 0), the CRB is minimised (optimal sensing).
  • When the entire signal carries random data (R=Rmax⁑R = R_{\max}), the CRB increases because the receiver must first estimate the data before extracting sensing information.
  • The intermediate region is governed by the allocation of signal dimensions between deterministic sensing pilots and random data symbols.

This work provides the missing converse for the ISAC problem: it proves that no system can simultaneously achieve the communication capacity and the optimal sensing CRB when the two functions compete for signal dimensions. The result was awarded the 2025 IEEE ComSoc/IT Joint Paper Award, recognising its fundamental contribution to the field.

isacinformation-theorycrbcapacitytradeoff
πŸŽ“CommIT Contribution(2023)

OTFS as an ISAC Waveform

W. Yuan, R. Schober, G. Caire β€” IEEE Trans. Wireless Commun.

Yuan, Schober, and Caire proposed OTFS (Orthogonal Time Frequency Space) as a waveform for ISAC systems operating in high-mobility scenarios. While OFDM is the natural ISAC waveform in low-to-moderate Doppler (Section 29.3), its performance degrades in high-Doppler environments (e.g., V2X at highway speeds, LEO satellite communication) where inter-carrier interference (ICI) corrupts both communication and sensing.

OTFS modulates data in the delay-Doppler domain using the Zak transform, where the channel is sparse and quasi-static. For sensing, OTFS provides a direct mapping between the delay-Doppler grid and the range-velocity plane, making target parameter extraction particularly natural.

Key advantages of OTFS for ISAC:

  • Doppler resilience: OTFS spreads each data symbol across the entire time-frequency plane, achieving full diversity and avoiding the ICI problem of OFDM.
  • Direct range-Doppler estimation: The delay-Doppler channel representation is sparse and directly reveals target parameters.
  • Efficient joint processing: A single Zak-domain matched filter serves both communication (channel estimation and equalisation) and sensing (target detection and parameter estimation).

Related work by Gaudio, Kobayashi, and Caire developed the DD-domain waveform design framework, optimising the input symbols in the delay-Doppler grid to balance communication and sensing metrics.

otfsisacdelay-dopplerzak-transformhigh-mobility

ISAC Beampattern with Trade-off Parameter

Visualise the ISAC transmit beampattern as a function of the sensing-communication trade-off parameter ρ\rho. A ULA with NN antennas forms a beam that balances between a communication direction and a sensing direction. When ρ=0\rho = 0, the beam points entirely towards the communication user; when ρ=1\rho = 1, it points towards the sensing target. Intermediate values produce a beampattern that serves both directions with reduced gain in each. Adjust the antenna count to observe how more antennas enable better spatial separation between the two beams, reducing the trade-off.

Parameters
16
0.5
30
-20

MIMO Radar Virtual Aperture Construction

Watch how a MIMO radar with NT=4N_T = 4 transmit and NR=8N_R = 8 receive antennas constructs a virtual array of NTΓ—NR=32N_T \times N_R = 32 elements. The virtual array provides the angular resolution of a 32-element phased array using only 12 physical antennas.
The virtual array (yellow) has NTΓ—NRN_T \times N_R elements, providing an NTN_T-fold angular resolution improvement over the physical receive array.

ISAC Beampattern Trade-off Animation

Watch the transmit beampattern evolve as the trade-off parameter ρ\rho sweeps from 0 (communication only) to 1 (sensing only). At ρ=0\rho = 0, the beam points towards the communication user (30 deg); at ρ=1\rho = 1, it shifts to the sensing target (-20 deg). Intermediate values produce a compromise beampattern.
The green line marks the communication direction, the red line marks the sensing direction. The beam transitions between the two as ρ\rho varies.

MIMO Radar Virtual Aperture

Visualise the virtual aperture created by a MIMO radar with NTXN_{TX} transmit and NRXN_{RX} receive antennas. The physical transmit and receive ULAs are shown along with the resulting virtual array. The angular beampattern of the virtual array is displayed, showing how the effective aperture of NTXΓ—NRXN_{TX} \times N_{RX} elements produces a narrow beam. Adjust the number of transmit and receive antennas to see how the virtual aperture grows, and change the target angle to verify that the virtual array beampattern peaks at the correct direction.

Parameters
4
8
15

Common Mistake: Grating Lobes in MIMO Virtual Arrays

Mistake:

Assuming that the MIMO virtual array always produces a clean beampattern without grating lobes, regardless of the physical array spacing.

Correction:

The virtual array element spacing dRd_R must satisfy dR≀λ/2d_R \leq \lambda/2 to avoid grating lobes, just as with a physical array. If the transmit spacing is dT=NRdRd_T = N_R d_R and dR>Ξ»/2d_R > \lambda/2, the virtual array has grating lobes that create ambiguous angle estimates.

For ISAC systems at high frequencies (mmWave, sub-THz), the half-wavelength spacing is very small (Ξ»/2=1.9\lambda/2 = 1.9 mm at 77 GHz), which limits the physical antenna size and the achievable element gain. Sparse array designs (nested, coprime) can mitigate this at the cost of increased sidelobe levels.

MIMO Radar

A radar system that transmits orthogonal waveforms from multiple antennas, creating a virtual aperture of NTΓ—NRN_T \times N_R elements from NT+NRN_T + N_R physical antennas. Provides improved angular resolution compared to phased-array radar of the same physical size.

Related: ISAC, DFRC

Quick Check

A colocated MIMO radar has NT=8N_T = 8 transmit antennas and NR=12N_R = 12 receive antennas, with the transmit ULA spacing chosen as dT=NRβ‹…dRd_T = N_R \cdot d_R to produce a filled virtual ULA. How many virtual array elements are created, and by what factor is the angular resolution improved compared to a phased-array radar using only the 12 receive antennas?

20 virtual elements, angular resolution improved by a factor of 20/12

96 virtual elements, angular resolution improved by a factor of 8 (equal to NTN_T)

96 virtual elements, angular resolution improved by a factor of 96

96 virtual elements, angular resolution improved by a factor of 12