OTFS as a 2D SFT of QAM Symbols

OTFS = QAM Symbols Γ— 2D SFT

With all three pieces β€” the DD grid, the ISFFT, and the standard OFDM Heisenberg transform β€” we can now state what OTFS is, operationally:

OTFS maps QAM symbols placed on the delay-Doppler grid to a time-domain waveform via the composition: ISFFT β†’ OFDM modulator.

This two-step construction is due to Hadani–Rakib (2017) and is the universally-adopted OTFS modulator. An alternative, mathematically cleaner construction β€” Zak-OTFS β€” works directly with the discrete Zak transform and is treated in Chapter 6. Both constructions produce the same time-domain waveform in the idealized limit of bi-orthogonal pulse shapes.

This section examines the ISFFT-plus-OFDM construction from the DD-grid perspective and sets up the notation used in subsequent chapters.

Definition:

OTFS Modulator (Hadani–Rakib Construction)

The OTFS transmit signal is generated by three operations:

  1. Place QAM symbols XDD[β„“,k]X_{DD}[\ell, k] on the MΓ—NM \times N DD grid.
  2. Apply ISFFT to obtain XTF[n,m]=ISFFT(XDD)[n,m]X_{TF}[n, m] = \text{ISFFT}(X_{DD})[n, m].
  3. Apply OFDM modulation (Heisenberg transform) to XTFX_{TF}: the time-domain waveform is x(t)β€…β€Š=β€…β€Šβˆ‘n=0Nβˆ’1βˆ‘m=0Mβˆ’1XTF[n,m] gtx(tβˆ’nTs) ej2Ο€mΞ”f(tβˆ’nTs),x(t) \;=\; \sum_{n = 0}^{N-1}\sum_{m = 0}^{M-1} X_{TF}[n, m]\,g_{tx}(t - n T_s)\,e^{j 2\pi m \Delta f(t - n T_s)}, where gtx(t)g_{tx}(t) is the pulse shape (typically a rectangular window of width TsT_s plus cyclic prefix). The overall mapping DD grid β†’ time-domain waveform is a unitary linear map from CMN\mathbb{C}^{MN} to a finite-dimensional subspace of L2([0,T+TCP))L^2([0, T + T_{CP})).

Theorem: OTFS is a Gabor Expansion with DD Coefficients

Applying the definition of the ISFFT, the OTFS transmit signal can be written directly as a delay-Doppler Gabor expansion: x(t)β€…β€Š=β€…β€Š1MNβˆ‘β„“=0Mβˆ’1βˆ‘k=0Nβˆ’1XDD[β„“,k] ϕℓ,k(t),x(t) \;=\; \frac{1}{\sqrt{MN}}\sum_{\ell = 0}^{M-1}\sum_{k = 0}^{N-1} X_{DD}[\ell, k]\,\phi_{\ell, k}(t), where Ο•β„“,k(t)\phi_{\ell, k}(t) is a prototype pulse delayed by ℓΔτ\ell \Delta\tau and Doppler-shifted by kΔνk \Delta\nu: Ο•β„“,k(t)β€…β€Š=β€…β€Šβˆ‘n=0Nβˆ’1βˆ‘m=0Mβˆ’1ej2Ο€(kn/Nβˆ’mβ„“/M) gtx(tβˆ’nTs) ej2Ο€mΞ”f(tβˆ’nTs).\phi_{\ell, k}(t) \;=\; \sum_{n = 0}^{N-1}\sum_{m = 0}^{M-1} e^{j 2\pi (k n / N - m \ell / M)}\,g_{tx}(t - n T_s)\,e^{j 2\pi m \Delta f(t - n T_s)}.

The OTFS transmit signal is a sum of delay-Doppler shifted copies of a prototype pulse, with coefficients equal to the data symbols on the DD grid. This is exactly the Gabor expansion structure discussed in RThe Gabor Expansion Is the OTFS Modulator.

The functions Ο•β„“,k\phi_{\ell, k} are designed to occupy (approximately) the (β„“,k)(\ell, k)-th cell of the DD grid. Their mutual orthogonality (up to the ideal bi-orthogonal pulse limit) is what makes OTFS invertible at the receiver.

Key Takeaway

OTFS is the Gabor expansion on the critical DD lattice. Each data symbol XDD[β„“,k]X_{DD}[\ell, k] contributes a delay-Doppler shifted copy of a prototype pulse to the transmit waveform. The symbolic expansion x(t)=βˆ‘β„“,kXDD[β„“,k] ϕℓ,k(t)x(t) = \sum_{\ell, k} X_{DD}[\ell, k]\,\phi_{\ell, k}(t) is the direct implementation of Chapter 2's Gabor-lattice theory. The ISFFT-plus-OFDM formulation is the FFT-friendly way to compute this expansion, but the underlying structure is the Gabor lattice on [0,T]Γ—[0,W][0, T] \times [0, W] at critical density.

Prototype Pulses Ο•β„“,k(t)\phi_{\ell, k}(t) in OTFS

For various (β„“0,k0)(\ell_0, k_0) DD grid positions, plot the prototype pulse Ο•β„“0,k0(t)\phi_{\ell_0, k_0}(t). Observe that it is localized in time near t=β„“0Δτ(modTs)t = \ell_0 \Delta\tau \pmod{T_s} and oscillates at the Doppler rate k0Δνk_0 \Delta\nu. The collection {Ο•β„“,k}\{\phi_{\ell, k}\} tiles the time-frequency plane at critical density.

Parameters
4
2
16
16

OTFS Transmitter Chain: QAM Grid to Waveform

A 16Γ—4 DD grid populated with random QPSK symbols. We animate: (1) the DD grid display, (2) the ISFFT producing a TF-grid image, (3) the OFDM Heisenberg transform producing the time-domain waveform. Each step is labelled with the corresponding math expression. The final waveform is the OTFS transmit signal, ready for the channel.

Example: A 2Γ—2 OTFS End-to-End Example

Build the OTFS transmit signal for M=N=2M = N = 2 and BPSK data XDD[0,0]=+1,XDD[0,1]=βˆ’1,XDD[1,0]=+1,XDD[1,1]=+1X_{DD}[0, 0] = +1, X_{DD}[0, 1] = -1, X_{DD}[1, 0] = +1, X_{DD}[1, 1] = +1. Assume a rectangular pulse gtx(t)=(1/Ts) 1[0,Ts)(t)g_{tx}(t) = (1/\sqrt{T_s})\,\mathbf{1}_{[0, T_s)}(t). Compute: (a) XTFX_{TF}, (b) x(t)x(t).

πŸŽ“CommIT Contribution(2024)

A Three-Part Tutorial on OTFS

W. Yuan, R. Schober, G. Caire β€” IEEE Communications Magazine

The Yuan–Schober–Caire tutorial (three parts, 2023–2024) is the single best entry point to the mathematical structure of OTFS for readers who already know OFDM. Part I (OFDM ↔ OTFS relationship), Part II (waveform design and receiver), and Part III (ISAC applications) together give a complete picture.

The CommIT group's involvement in this tutorial and in the broader standardization discussions around OTFS is what anchors the book's cell-free (Ch. 17) and LEO-satellite (Ch. 18) CommIT contributions in the next parts of the library. Part III in particular is the conceptual bridge between the DD-domain waveform machinery developed here and the ISAC-specific analysis in Chapter 12.

otfs-tutorialwaveform-designcaire

Common Mistake: The Ideal Bi-Orthogonal Pulse Assumption

Mistake:

Many derivations assume the prototype pulse gtxg_{tx} is bi-orthogonal β€” i.e., ⟨gtx(β‹…βˆ’nTs) ej2Ο€mΞ”fβ‹…,grx(β‹…)⟩=Ξ΄n,0Ξ΄m,0\langle g_{tx}(\cdot - n T_s)\,e^{j 2\pi m \Delta f\cdot}, g_{rx}(\cdot)\rangle = \delta_{n, 0}\delta_{m, 0} for a matched receiver pulse grxg_{rx}. This is not achievable exactly for any finite-energy pulse (Balian-Low theorem).

Correction:

With the rectangular pulse gtx=grx=(1/Ts)1[0,Ts)g_{tx} = g_{rx} = (1/\sqrt{T_s})\mathbf{1}_{[0, T_s)} and a cyclic prefix of length β‰₯Ο„max⁑\geq \tau_{\max}, the bi-orthogonality condition is satisfied exactly for the discrete-time OFDM grid β€” but only approximately for the Doppler direction. The approximate nature shows up as off-diagonal entries in the DD channel matrix (Chapter 4), typically small for fDTsβ‰ͺ1f_D T_s \ll 1. For fDTsf_D T_s not small (LEO, HST at mmWave), pulse shaping matters more β€” see Chapter 20.