Basis Expansion Models and Windowing

Two Families of Mitigation

Fractional Doppler is a fact of life. Section 2 embraced it in the input-output relation; this section develops two classical mitigations that reduce the IDI severity:

  1. Basis expansion models (BEM): represent the time-varying channel as a weighted sum of a small number of known basis functions over the OTFS frame. This is the BEM approach from FSI Ch. 11, specialized to OTFS's doubly-selective setting.
  2. Windowing: apply a tapered window (Hamming, Blackman) along the Doppler axis before the SFFT to suppress Dirichlet sidelobes, at the cost of a slightly widened main lobe.

Neither technique eliminates IDI, but both reduce it enough that the integer-Doppler assumption becomes approximately valid β€” or at least, the IDI kernel's support shrinks, lowering detection complexity.

Definition:

Basis Expansion Model (BEM)

A basis expansion model represents the time-varying path gain hi(t)h_i(t) over a block of duration TT as hi(t)β€…β€Šβ‰ˆβ€…β€Šβˆ‘q=0QBEMβˆ’1hi(q)β€‰Οˆq(t),h_i(t) \;\approx\; \sum_{q=0}^{Q_{\text{BEM}} - 1} h_i^{(q)}\,\psi_q(t), where {ψq(t)}\{\psi_q(t)\} are QBEMQ_{\text{BEM}} known basis functions (typically complex exponentials, polynomials, or DPSS) and {hi(q)}\{h_i^{(q)}\} are the BEM coefficients to be estimated.

The most common choice β€” the complex-exponential BEM (CE-BEM) β€” uses ψq(t)=ej2Ο€qt/T\psi_q(t) = e^{j 2\pi q t/T} for q=βˆ’Q/2,…,Q/2βˆ’1q = -Q/2, \ldots, Q/2 - 1 with Q=QBEMQ = Q_{\text{BEM}} an odd integer capturing the maximum Doppler rate.

,

Theorem: BEM Converts Fractional to Integer Doppler

Under the CE-BEM approximation with QBEMQ_{\text{BEM}} basis functions, the fractional-Doppler channel reduces to an effective Pβ‹…QBEMP \cdot Q_{\text{BEM}}-path channel with integer Doppler offsets. The IDI kernel collapses to Ξ΄\delta-functions at the BEM-basis Doppler shifts.

Formally: hi(t)=hi ej2πνitβ‰ˆhiβˆ‘q=0QBEMβˆ’1cq(Ξ½iT) ej2Ο€qt/T,h_i(t) = h_i\,e^{j 2\pi \nu_i t} \approx h_i \sum_{q=0}^{Q_{\text{BEM}}-1} c_q(\nu_i T)\,e^{j 2\pi q t/T}, where cq(β‹…)c_q(\cdot) are Fourier coefficients. Each BEM term contributes a single integer-Doppler path at bin qq, with amplitude hicq(Ξ½iT)h_i c_q(\nu_i T). The total effective path count is Pβ‹…QBEMP \cdot Q_{\text{BEM}} at integer Doppler positions.

BEM trades one fractional-Doppler path for QBEMQ_{\text{BEM}} integer-Doppler paths. At QBEM=3Q_{\text{BEM}} = 3, a single path becomes three effective paths β€” with total energy preserved but distributed across three Doppler bins. After BEM, the integer-Doppler detector of Chapter 8 applies directly, with the caveat that the "number of paths" is now Pβ‹…QBEMP \cdot Q_{\text{BEM}}.

Key Takeaway

BEM reduces fractional Doppler to a solved problem. Choose QBEM=3Q_{\text{BEM}} = 3 for moderate Doppler or 55 for large fractional offsets; apply Chapter 8's integer-Doppler detectors on the expanded channel. Residual error: ≲1\lesssim 1 dB for QBEM=3Q_{\text{BEM}} = 3, <0.3< 0.3 dB for QBEM=5Q_{\text{BEM}} = 5. The BEM approach is the workhorse of fractional-Doppler OTFS: it preserves the sparse-channel detector framework while accommodating real-world channel geometry.

BEM Coefficients for a Fractional-Doppler Path

For a single path with fractional Doppler Ο΅\epsilon, plot the BEM coefficients ∣cq(Ο΅)∣2|c_q(\epsilon)|^2 as a function of basis index qq, for increasing QBEMQ_{\text{BEM}}. At Ο΅=0\epsilon = 0: single spike at q=0q = 0. At Ο΅=0.4\epsilon = 0.4: energy spread across q=0,1,βˆ’1,2q = 0, 1, -1, 2. The BEM captures most of the energy in the first QBEMQ_{\text{BEM}} terms; the residual is the approximation error.

Parameters
0.3
3

Definition:

Doppler-Axis Windowing

Doppler windowing applies a tapered window w[k]w[k] to the NN OFDM symbols of the OTFS frame before the SFFT: YTF[n,m]β†’YTF[n,m]β‹…w[n],n=0,…,Nβˆ’1.Y_{TF}[n, m] \to Y_{TF}[n, m] \cdot w[n], \qquad n = 0, \ldots, N - 1. Typical windows: Hamming (w[n]=0.54βˆ’0.46cos⁑(2Ο€n/N)w[n] = 0.54 - 0.46\cos(2\pi n/N)), Blackman (w[n]=0.42βˆ’0.5cos⁑(2Ο€n/N)+0.08cos⁑(4Ο€n/N)w[n] = 0.42 - 0.5\cos(2\pi n/N) + 0.08\cos(4\pi n/N)), and Tukey (tapered cosine).

The effect on the IDI kernel: sidelobes are suppressed at the cost of a wider main lobe. For Hamming, sidelobes drop from βˆ’13-13 dB (no window) to βˆ’43-43 dB.

Theorem: Windowing Sidelobe Suppression vs Main-Lobe Widening

For a window w[n]w[n] with equivalent noise bandwidth (ENBW) Ξ±ENBW\alpha_{\text{ENBW}} and peak sidelobe Ξ±PSL\alpha_{\text{PSL}} (in dB), the IDI kernel under windowing satisfies:

  • Main-lobe width (3-dB): Ξ±ENBWβ‹…(1/N)\alpha_{\text{ENBW}}\cdot(1/N) cells. For Hamming: Ξ±ENBW=1.36\alpha_{\text{ENBW}} = 1.36 (36%36\% wider than rectangular).
  • Peak sidelobe: Ξ±PSL\alpha_{\text{PSL}} dB below main lobe. Hamming: βˆ’43-43 dB (vs βˆ’13-13 dB for rectangular).

The product Ξ±ENBWβ‹…10Ξ±PSL/20\alpha_{\text{ENBW}} \cdot 10^{\alpha_{\text{PSL}}/20} is a figure of merit: smaller is better. Standard windows give values in the range [0.5,2][0.5, 2].

Windowing is a classical DSP technique (Harris 1978) for leakage suppression. In OTFS, it trades a slight increase in effective Doppler bin size for a dramatic reduction in sidelobe energy. For fractional-Doppler OTFS at high SNR, windowing is "free" because the main-lobe widening is negligible compared to the sidelobe suppression benefit.

Windowing does not eliminate IDI β€” it just concentrates the leakage in fewer neighboring bins. For kIDI=1k_{\text{IDI}} = 1 under Hamming windowing (vs kIDI=3k_{\text{IDI}} = 3 for rectangular), the extended matrix has 3Γ— fewer non-zeros, so detection is proportionally faster.

Windowing Functions for OTFS

WindowMain-lobe (ENBW)Peak sidelobe (dB)When to use
Rectangular (no window)1.0-13Theoretical baseline
Hamming1.36-43Moderate fractional (standard choice)
Blackman1.73-58Large fractional (Ο΅>0.3\epsilon > 0.3)
Blackman-Harris1.90-92Very large fractional (research)
Tukey (Ξ±=0.5\alpha = 0.5)1.22-25Light windowing, compromise

Windowed IDI Kernels: Shape vs Window

Plot the IDI kernel ∣KIDIw(k,ϡ)∣2|K_{\text{IDI}}^w(k, \epsilon)|^2 for several window choices (rectangular, Hamming, Blackman) at fixed ϡ=0.3\epsilon = 0.3. Observe: the rectangular kernel has high sidelobes that decay slowly. Hamming suppresses sidelobes by ∼30\sim 30 dB. Blackman reduces them even further but widens the main lobe.

Parameters
0.3
32
⚠️Engineering Note

Choosing Between BEM and Windowing

The two approaches are complementary:

  • BEM: explicit model of the IDI kernel. Good for high-precision detection (low BER targets). Requires estimating additional path parameters. Overhead: O(Pβ‹…QBEM)O(P \cdot Q_{\text{BEM}}) parameters vs O(P)O(P) integer case. Optimal when the fractional offsets are known or can be estimated.

  • Windowing: pre-processing step with no additional parameters. Zero-cost detection (same detector runs on the windowed grid). Suppresses sidelobes but cannot eliminate the main-lobe leakage. Best for moderate fractional offsets and when detector simplicity matters.

Combined approach: Hamming windowing + QBEM=3Q_{\text{BEM}} = 3 often gives the best trade-off: windowing suppresses far-out sidelobes; BEM captures residual main-lobe spread. Net: ∼0.5\sim 0.5 dB from ML at typical parameters.

For 5G NR-aligned OTFS: Hamming + QBEM=3Q_{\text{BEM}} = 3 is the reference implementation. For LEO satellite: Blackman + QBEM=5Q_{\text{BEM}} = 5.

Practical Constraints
  • β€’

    BEM needs path Doppler estimate (Ο΅\epsilon) to set coefficients

  • β€’

    Windowing adds no parameters but can't fully remove IDI

  • β€’

    Best-of-both: windowing for sidelobes + BEM for main lobe

Example: BER With BEM vs Integer Assumption

Compare BER at SNR = 15 dB for an OTFS system with P=4P = 4 paths, all at Ο΅=0.3\epsilon = 0.3, under (a) integer assumption (ignore IDI), (b) BEM with QBEM=3Q_{\text{BEM}} = 3.

πŸŽ“CommIT Contribution(2020)

BEM-OTFS for Fractional Doppler

G. D. Surabhi, A. Chockalingam, G. Caire β€” IEEE Trans. Vehicular Technology

Surabhi, Chockalingam, and Caire extended the diversity analysis of OTFS Chapter 9 to fractional-Doppler channels. Their key result: with BEM of order QBEMQ_{\text{BEM}}, OTFS recovers diversity order min⁑(P,min⁑(QBEM,kmax⁑))\min(P, \min(Q_{\text{BEM}}, k_{\max})) β€” the minimum of path count and the BEM's effective Doppler resolution.

The CommIT contribution refines the naive fractional-Doppler treatment by showing that BEM's effective Doppler resolution is the fundamental limit β€” beyond QBEM=2kmax⁑+1Q_{\text{BEM}} = 2 k_{\max} + 1, no gain is obtained. This informs the reference implementation's choice of QBEM=3Q_{\text{BEM}} = 3 for typical terrestrial channels.

commitbemfractional-dopplerdiversity