Definition and Motivation

What Transform Do We Need?

Chapter 1 told us the channel lives naturally in the delay-Doppler plane. It did not tell us how to map an ordinary time-domain signal x(t)x(t) into that plane. This chapter builds that map. The answer is the Zak transform.

A naive guess would be the 2D Fourier transform of xx, but x(t)x(t) is a function of one variable, not two — we cannot simply 2D-Fourier-transform it. The Zak transform is, cleverly, a 1D transform of xx whose output is a function of two variables (t,ν)(t, \nu). It does this by folding xx onto itself at period T0T_0 and simultaneously collecting a Fourier series in that folding. The result lives on a 2D torus and has exactly the quasi-periodicity structure that the delay-Doppler plane needs.

The Zak transform was introduced by Joshua Zak in 1967 in the context of solid-state physics (Bloch functions in electron band theory). It was rediscovered in signal processing in the 1980s by Janssen, as the right tool for Gabor analysis. For OTFS, the Zak transform is how a continuous-time waveform meets the 2D DD grid.

Definition:

Zak Transform

Let f:RCf: \mathbb{R} \to \mathbb{C} be a sufficiently regular function (e.g., in L2(R)L^2(\mathbb{R})). For a parameter T0>0T_0 > 0 (the Zak period), the Zak transform of ff is Zf(t,ν)  =  kZf(tkT0)ej2πkνT0,Z_f(t, \nu) \;=\; \sum_{k \in \mathbb{Z}} f(t - k\,T_0)\,e^{-j 2\pi k\,\nu\,T_0}, where tRt \in \mathbb{R} and νR\nu \in \mathbb{R} are the time and Doppler variables on the delay-Doppler plane. The dual frequency period is ν0=1/T0\nu_0 = 1/T_0.

Different authors use different sign conventions. Hadani–Rakib (2017) and Raviteja–Viterbo (2018) use ej2πkνT0e^{-j 2\pi k\,\nu\,T_0} — we adopt that convention throughout this book to match the OTFS literature. Some mathematical references (e.g., Folland 1989) use the opposite sign, which flips the role of the Doppler axis.

, ,

The Zak Transform Is a Fourier Series in Disguise

Fix tt. The sequence {f(tkT0)}kZ\{f(t - k T_0)\}_{k \in \mathbb{Z}} consists of samples of ff spaced by T0T_0. The Zak transform in ν\nu is exactly the Fourier series of this sequence. So Zf(t,ν)Z_f(t, \nu) takes a continuous function f(t)f(t) and encodes, at each time tt, the Fourier series of the infinite samples {f(tkT0)}\{f(t - k T_0)\}.

The reason this is useful: the sparse DD channel acts by delay (ttτi)(t \mapsto t - \tau_i) and Doppler (ffej2πνit)(f \mapsto f\,e^{j 2\pi \nu_i t}) — exactly the two operations that the Zak transform naturally respects, as we now prove.

Theorem: Zak Transform of a Delay-and-Doppler-Shifted Signal

Let πτ0,ν0\pi_{\tau_0, \nu_0} denote the delay-and-Doppler-shift operator: (πτ0,ν0f)(t)=f(tτ0)ej2πν0t(\pi_{\tau_0, \nu_0} f)(t) = f(t - \tau_0)\,e^{j 2\pi \nu_0 t}. Then Zπτ0,ν0f(t,ν)  =  ej2πν0tZf(tτ0,νν0).Z_{\pi_{\tau_0, \nu_0} f}(t, \nu) \;=\; e^{j 2\pi \nu_0 t}\,Z_f(t - \tau_0, \nu - \nu_0). That is: a delay and Doppler shift of ff becomes a translation in the DD plane of ZfZ_f, up to a phase factor.

This is the defining "covariance property" of the Zak transform — and the reason it is the correct signal space for the DD domain. A channel that delays and Doppler-shifts the transmit signal by (τ0,ν0)(\tau_0, \nu_0) simply translates the DD-domain picture by (τ0,ν0)(\tau_0, \nu_0). The input-output relation becomes a 2D convolution, as promised in Chapter 1.

Key Takeaway

The Zak transform is covariant under delay-Doppler translation. A time-domain signal that has been delayed by τ0\tau_0 and Doppler-shifted by ν0\nu_0 appears, in the DD domain, as the same DD picture shifted by (τ0,ν0)(\tau_0, \nu_0). This single property is what makes the Zak transform the right signal space for the delay-Doppler channel. The rest of OTFS is, loosely, bookkeeping around this fact.

Example: Zak Transform of a Rectangular Pulse

Let g(t)=1[0,T0)(t)g(t) = \mathbf{1}_{[0, T_0)}(t) be the rectangular pulse of unit height supported on [0,T0)[0, T_0), i.e., the indicator of the fundamental interval. Compute Zg(t,ν)Z_g(t, \nu) and interpret the result.

Zak Transform of a Rectangular Pulse of Variable Width

Plot Zg(t,ν)|Z_g(t, \nu)| for a rectangular pulse g(t)=1[0,W0](t)g(t) = \mathbf{1}_{[0, W_0]}(t) of variable width W0W_0 on a Zak grid with period T0T_0. When W0=T0W_0 = T_0 the Zak transform is uniform and flat; when W0T0W_0 \neq T_0 it exhibits fringes that signal aliasing in the DD grid. Slide the pulse width through one period to see the clean case and the aliased cases.

Parameters
1
4

Example: Zak Transform of a Dirac Impulse

Compute Zδ(t,ν)Z_\delta(t, \nu) where δ\delta is the Dirac impulse at the origin.

Building Zf(t,ν)Z_f(t, \nu) by Folding

The Zak transform is built by taking a time-domain signal f(t)f(t), replicating it at integer multiples of T0T_0, and summing each copy with a phase rotation along the Doppler axis. The animation builds the first few terms of the sum and shows the resulting function on the 2D DD plane.

The Zak Transform is a Unitary Map

One can show that the Zak transform extends to a unitary isomorphism between L2(R)L^2(\mathbb{R}) and L2(T2)L^2(\mathbb{T}^2), where T2=[0,T0)×[0,ν0)\mathbb{T}^2 = [0, T_0) \times [0, \nu_0) is the fundamental rectangle (the "quasi-periodic cell"). Concretely, Rf(t)2dt  =  T00T0 ⁣ ⁣0ν0Zf(t,ν)2dνdt.\int_{\mathbb{R}} |f(t)|^2\,dt \;=\; T_0 \int_0^{T_0}\!\!\int_0^{\nu_0} |Z_f(t, \nu)|^2\,d\nu\,dt. Unitarity is important because it means the Zak transform preserves signal energy (no free information loss). We will need unitarity in Chapter 6 when designing the OTFS transmitter.

⚠️Engineering Note

Choosing the Zak Period

In OTFS, the Zak period T0T_0 is tied to the OTFS frame parameters: T0=T/NT_0 = T/N where TT is the frame duration and NN is the number of Doppler bins. Equivalently, T0T_0 is the OFDM symbol duration after Heisenberg embedding. Choosing T0T_0 determines the Doppler resolution Δν=1/T\Delta\nu = 1/T (since the DD plane has NN Doppler bins over [0,ν0)=[0,1/T0)[0, \nu_0) = [0, 1/T_0)).

Key rule: T0×fD,maxT_0 \times f_{D,\max} must satisfy T0fD,max<1/2T_0 f_{D,\max} < 1/2 to avoid Doppler aliasing (the Zak-plane analog of Nyquist). For vehicular channels (fD500f_D \approx 500 Hz) this requires T0<1T_0 < 1 ms; LEO satellite channels (fD30f_D \approx 30 kHz) require T0<17μsT_0 < 17\,\mu\text{s}.

Practical Constraints
  • T0fD,max<1/2T_0\, f_{D,\max} < 1/2 to avoid Doppler aliasing

  • Frame duration T=NT0T = N T_0 determines Doppler resolution Δν=1/T\Delta\nu = 1/T

  • Practical trade-off: smaller T0T_0 gives more Doppler slack but shorter frames

📋 Ref: Hadani et al. 2017, §II-B

Common Mistake: When the Zak Sum Converges

Mistake:

Treating the Zak transform as "always well-defined" — the infinite sum kf(tkT0)ej2πkνT0\sum_k f(t - k T_0)\,e^{-j 2\pi k \nu T_0} can fail to converge for signals that do not decay, such as constant or slowly-varying signals.

Correction:

Strictly speaking, the Zak transform converges absolutely for fL2(R)f \in L^2(\mathbb{R}) in the L2L^2-sense over the fundamental domain T2\mathbb{T}^2. Pointwise convergence for L1L^1 signals is guaranteed almost everywhere. For rigor in OTFS, treat the Zak transform as a map L2(R)L2(T2)L^2(\mathbb{R}) \to L^2(\mathbb{T}^2) — the discrete version (Section 4) sidesteps these subtleties entirely.