Zak-OTFS and Turbo Signal Processing for Joint Sensing and Communication

Jinu Jayachandran, Muhammad Ubadah, Saif Khan Mohammed, Ronny Hadani, Ananthanarayanan Chockalingam and Robert Calderbank J. Jayachandran, M. Ubadah and S. K. Mohammed are with Department of Electrical Engineering, Indian Institute of Technology Delhi, India (E-mail: [email protected], [email protected], [email protected]). S. K. Mohammed is also associated with Bharti School of Telecom. Technology and Management (BSTTM), IIT Delhi. R. Hadani is with Department of Mathematics, University of Texas at Austin, TX, USA (E-mail: [email protected]). A. Chockalingam is with Department of Electrical Communication Engineering, Indian Institute of Science Bangalore, India (E-mail: [email protected]). R. Calderbank is with Department of Electrical and Computer Engineering, Duke University, USA (E-mail: [email protected]).The work of S. K. Mohammed was supported in part by the Jai Gupta Chair at I.I.T. Delhi, and also in part by a project at BSTTM funded by Bharti Airtel Limited.This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible.
Abstract

The Zak-OTFS input/output (I/O) relation is predictable and non-fading when the delay and Doppler periods are greater than the effective channel delay and Doppler spreads, a condition which we refer to as the crystallization condition. The filter taps can simply be read off from the response to a single Zak-OTFS pilot pulsone, and the I/O relation can be reconstructed for a sampled system that operates under finite duration and bandwidth constraints. In previous work we had measured BER performance of a baseline system where we used separate Zak-OTFS subframes for sensing and data transmission. In this Letter we demonstrate how to use turbo signal processing to match BER performance of this baseline system when we integrate sensing and communication within the same Zak-OTFS subframe. The turbo decoder alternates between channel sensing using a noise-like waveform (spread pulsone) and recovery of data transmitted using point pulsones.

Index Terms:
Zak-OTFS, Integrated Sensing and Communication, Turbo processing.

I Introduction

6G propagation environments are changing the balance between time-frequency methods focused on OFDM signal processing and delay-Doppler methods (OTFS) [1, 3, 4, 5, 6]. OFDM is configured to prevent inter-carrier-interference (ICI) whereas OTFS is configured to embrace ICI [7, 8, 9, 10, 11, 12]. In OFDM, once the I/O relation is known, equalization is relatively simple at least when there is no ICI. However, acquisition of the I/O relation is non-trivial and dependent on the accuracy of the assumed propagation model [13]. In contrast, equalization is more involved in OTFS due to inter-symbol-interference (ISI), however acquisition of the I/O relation is simple and model free [14, 15, 16, 17]. Acquisition becomes more critical when Doppler spreads measured in KHz make it more and more difficult to estimate channels. The most challenging situation is the combination of unresolvable paths and high channel spreads. In this Letter we present simulation results for a Veh-A channel [18] where the first three paths are not separable and cannot be estimated accurately.

In previous work [14, 15] we have described how to design a parametric family of pulsone waveforms that can be matched to the delay and Doppler spreads of different propagation environments. A pulsone is a signal in the time domain which realizes a quasi-periodic localized function on the DD domain. The prototypical structure of a pulsone is a train of pulses modulated by a tone. We have analyzed performance in the situation when the pulsone parameters matches the environment channel parameters, in the sense that, the delay period of the pulsone is greater than the delay spread of the channel, and the Doppler period of the pulsone is greater than the Doppler spread of the channel. We refer to this condition as the crystallization condition. We start from this baseline system where we dedicate separate Zak-OTFS subframes to channel estimation and data transmission, abbreviated as S|C¯conditional𝑆¯𝐶S|\overline{C}italic_S | over¯ start_ARG italic_C end_ARG (sensing/channel estimation in the absence of communication signal) and C|S¯conditional𝐶¯𝑆C|\overline{S}italic_C | over¯ start_ARG italic_S end_ARG (communication in the absence of sensing signal).

The characteristic structure of a pulsone is a train of pulses modulated by a tone, a signal with unattractive peak-to-average power ratio (PAPR). In more recent work [19] we have introduced the notion of filtering in the discrete delay-Doppler domain. We have described how to construct spread waveforms with desirable characteristics by applying a chirp filter in the discrete DD domain to a point pulsone. One desirable characteristic is low PAPR, about 6dB for the exemplar spread pulsone, compared with about 15dB for the point pulsone. A second desirable characteristic is the ability to read off the I/O relation of the sampled communication system provided a second crystallization condition is satisfied. This work demonstrates how to integrate sensing and communication within a single Zak-OTFS subframe by combining a spread pulsone for channel sensing with point pulsones for data transmission. The filter in the discrete DD domain enables coexistence by minimizing interference between sensing and data transmission. We have demonstrated that sharing DD domain resources in this way increases effective throughput compared with traditional approaches that use guard bands to divide DD domain resources between sensing and communication.

In this Letter we demonstrate that turbo signal processing is able to close the performance gap between separate sensing and communication (S|C¯conditional𝑆¯𝐶S|\overline{C}italic_S | over¯ start_ARG italic_C end_ARG and C|S¯conditional𝐶¯𝑆C|\overline{S}italic_C | over¯ start_ARG italic_S end_ARG) and joint sensing and communication in the same Zak-OTFS subframe (which we abbreviate as S|Cconditional𝑆𝐶S|Citalic_S | italic_C and C|Sconditional𝐶𝑆C|Sitalic_C | italic_S). The turbo principle of iterating between functional blocks in communication receivers has proven to be a powerful method of improving performance. For example, turbo iterations between channel equalizer and channel decoder have been shown to yield tremendous improvements in bit error performance [20], [21]. Likewise, iterations between channel estimation and turbo equalizer have been shown to improve the channel estimate over the iterations by using soft information fed back from the decoder from the previous iteration to generate extended training sequences between the actual transmitted training sequences [22].

Refer to caption
Figure 1: Signal processing for proposed Zak-OTFS based iterative joint sensing and communication.

In joint sensing and communication, we estimate the effective channel, then estimate the received spread pilot, then recover the data after subtracting our estimate for the received pilot from the received signal (see Fig. 1). The residual pilot (after cancellation) interferes with data transmission. In the turbo iteration, we take the estimated data, then estimate the received data signal, then improve our estimate for the effective channel by subtracting our estimate for the received data signal from the received signal. We show that five turbo iterations suffice to match the performance of separate sensing and communication across a wide range of Doppler shifts.

II System model

In previous work [19] we have introduced a framework for joint sensing and communication where the pilot/sensing signal and the data signal are embedded in a single Zak-OTFS subframe. We recall [19] that Zak-OTFS carrier waveforms are quasi-periodic functions in the discrete DD domain with period M𝑀Mitalic_M along the discrete delay axis and period N𝑁Nitalic_N along the discrete Doppler axis. The superposition of pilot and data signals in the discrete DD domain is given by

xdd[k,l]subscript𝑥dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{dd}}}[k,l]italic_x start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== Edxd,dd[k,l]+Epxs,dd[k,l]subscript𝐸𝑑subscript𝑥d,dd𝑘𝑙subscript𝐸𝑝subscript𝑥s,dd𝑘𝑙\displaystyle\sqrt{E_{d}}\,x_{\mbox{\scriptsize{d,dd}}}[k,l]\,+\,\sqrt{E_{p}}% \,x_{\mbox{\scriptsize{s,dd}}}[k,l]square-root start_ARG italic_E start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT end_ARG italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] + square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] (1)

where Epsubscript𝐸𝑝E_{p}italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT denotes the energy of the pilot and Edsubscript𝐸𝑑E_{d}italic_E start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT the energy of the data signal. The ratio Ep/Edsubscript𝐸𝑝subscript𝐸𝑑E_{p}/E_{d}italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT / italic_E start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT is therefore the pilot to data power ratio (PDR). The data signal xd,dd[k,l]subscript𝑥d,dd𝑘𝑙x_{\mbox{\scriptsize{d,dd}}}[k,l]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ], k,l𝑘𝑙k,l\in{\mathbb{Z}}italic_k , italic_l ∈ blackboard_Z is given by

xd,dd[k,l]subscript𝑥d,dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{d,dd}}}[k,l]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== 1MNx[kmodM,lmodN]ej2πkMlN,1𝑀𝑁𝑥𝑘mod𝑀𝑙mod𝑁superscript𝑒𝑗2𝜋𝑘𝑀𝑙𝑁\displaystyle\frac{1}{\sqrt{MN}}\,x[k\,\mbox{\footnotesize{mod}}\,M\,,\,l\,% \mbox{\footnotesize{mod}}\,N]\,e^{j2\pi\frac{\left\lfloor\frac{k}{M}\right% \rfloor l}{N}},divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_M italic_N end_ARG end_ARG italic_x [ italic_k mod italic_M , italic_l mod italic_N ] italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG ⌊ divide start_ARG italic_k end_ARG start_ARG italic_M end_ARG ⌋ italic_l end_ARG start_ARG italic_N end_ARG end_POSTSUPERSCRIPT , (2)

where each of the MN𝑀𝑁MNitalic_M italic_N information symbols x[k,l]𝑥𝑘𝑙x[k,l]italic_x [ italic_k , italic_l ], k=0,,M1𝑘0𝑀1k=0,\cdots,M-1italic_k = 0 , ⋯ , italic_M - 1, l=0,1,,N1𝑙01𝑁1l=0,1,\cdots,N-1italic_l = 0 , 1 , ⋯ , italic_N - 1 have average energy 𝔼[|x[k,l]|2]=1𝔼delimited-[]superscript𝑥𝑘𝑙21{\mathbb{E}}\left[|x[k,l]|^{2}\right]=1blackboard_E [ | italic_x [ italic_k , italic_l ] | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] = 1. The exponential phase term in (2) renders xd,dd[k,l]subscript𝑥d,dd𝑘𝑙x_{\mbox{\scriptsize{d,dd}}}[k,l]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] quasi-periodic. For all n,m𝑛𝑚n,m\in{\mathbb{Z}}italic_n , italic_m ∈ blackboard_Z

xd,dd[k+nM,l+mN]subscript𝑥d,dd𝑘𝑛𝑀𝑙𝑚𝑁\displaystyle x_{\mbox{\scriptsize{d,dd}}}[k+nM,l+mN]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k + italic_n italic_M , italic_l + italic_m italic_N ] =\displaystyle== ej2πnlNxd,dd[k,l].superscript𝑒𝑗2𝜋𝑛𝑙𝑁subscript𝑥d,dd𝑘𝑙\displaystyle e^{j2\pi\frac{nl}{N}}\,x_{\mbox{\scriptsize{d,dd}}}[k,l].italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_n italic_l end_ARG start_ARG italic_N end_ARG end_POSTSUPERSCRIPT italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] . (3)

The data signal in (1) can also be expressed in terms of filtering (in the DD domain) a quasi-periodic DD domain pulse at the origin.

xd,dd[k,l]subscript𝑥d,dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{d,dd}}}[k,l]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== [(k=0M1l=0N1x[k,l]δ[kk]δ[ll])\displaystyle{\Bigg{[}}\,\left(\sum\limits_{k^{\prime}=0}^{M-1}\sum\limits_{l^% {\prime}=0}^{N-1}x[k^{\prime},l^{\prime}]\,\delta[k-k^{\prime}]\,\delta[l-l^{% \prime}]\right)\,[ ( ∑ start_POSTSUBSCRIPT italic_k start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_M - 1 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_x [ italic_k start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ] italic_δ [ italic_k - italic_k start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ] italic_δ [ italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ] ) (4)
σx0,dd[k,l]],\displaystyle\hskip 8.53581pt*_{\sigma}\,x_{\mbox{\scriptsize{0,dd}}}[k,l]\,{% \Bigg{]}},∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT 0,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] ] ,

where

x0,dd[k,l]subscript𝑥0,dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{0,dd}}}[k,l]italic_x start_POSTSUBSCRIPT 0,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =ΔsuperscriptΔ\displaystyle\stackrel{{\scriptstyle\Delta}}{{=}}start_RELOP SUPERSCRIPTOP start_ARG = end_ARG start_ARG roman_Δ end_ARG end_RELOP n,mδ[knM]δ[lmN]subscript𝑛𝑚𝛿delimited-[]𝑘𝑛𝑀𝛿delimited-[]𝑙𝑚𝑁\displaystyle\sum\limits_{n,m\in{\mathbb{Z}}}\delta[k-nM]\,\delta[l-mN]∑ start_POSTSUBSCRIPT italic_n , italic_m ∈ blackboard_Z end_POSTSUBSCRIPT italic_δ [ italic_k - italic_n italic_M ] italic_δ [ italic_l - italic_m italic_N ] (5)

is the data pulsone corresponding to a DD pulse at the origin. DD domain filtering is implemented through twisted convolution σsubscript𝜎*_{\sigma}∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT and the filter taps in (4) depend on the information symbols. The pilot signal in (1) is given by

xs,dd[k,l]subscript𝑥s,dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{s,dd}}}[k,l]italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== ws[k,l]σxp,dd[k,l],subscript𝜎subscript𝑤𝑠𝑘𝑙subscript𝑥p,dd𝑘𝑙\displaystyle w_{s}[k,l]\,*_{\sigma}\,x_{\mbox{\scriptsize{p,dd}}}[k,l],italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT p,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] , (6)

where the spreading filter ws[k,l]subscript𝑤𝑠𝑘𝑙w_{s}[k,l]italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT [ italic_k , italic_l ] acts on the point pilot xp,dd[k,l]subscript𝑥p,dd𝑘𝑙x_{\mbox{\scriptsize{p,dd}}}[k,l]italic_x start_POSTSUBSCRIPT p,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] by twisted convolution. We follow [19] by considering chirp filters, where for k=0,1,,M1𝑘01𝑀1k=0,1,\cdots,M-1italic_k = 0 , 1 , ⋯ , italic_M - 1, l=0,1,,N1𝑙01𝑁1l=0,1,\cdots,N-1italic_l = 0 , 1 , ⋯ , italic_N - 1

ws[k,l]subscript𝑤𝑠𝑘𝑙\displaystyle w_{s}[k,l]italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== 1MNej2πq(k2+l2)MN,1𝑀𝑁superscript𝑒𝑗2𝜋𝑞superscript𝑘2superscript𝑙2𝑀𝑁\displaystyle\frac{1}{MN}\,e^{j2\pi\frac{q(k^{2}+l^{2})}{MN}},divide start_ARG 1 end_ARG start_ARG italic_M italic_N end_ARG italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_q ( italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_M italic_N end_ARG end_POSTSUPERSCRIPT , (7)

and we refer to q𝑞q\in{\mathbb{Z}}italic_q ∈ blackboard_Z as the slope parameter. The point pilot xp,dd[k,l]subscript𝑥p,dd𝑘𝑙x_{\mbox{\scriptsize{p,dd}}}[k,l]italic_x start_POSTSUBSCRIPT p,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] appearing in (6) is a discrete DD domain quasi-periodic impulse located at (kp,lp)subscript𝑘𝑝subscript𝑙𝑝(k_{p},l_{p})( italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) and is given by

xp,dd[k,l]subscript𝑥p,dd𝑘𝑙\displaystyle x_{\mbox{\scriptsize{p,dd}}}[k,l]italic_x start_POSTSUBSCRIPT p,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== n,mej2πnlpNδ[kkpnM]δ[llpmN].subscript𝑛𝑚superscript𝑒𝑗2𝜋𝑛subscript𝑙𝑝𝑁𝛿delimited-[]𝑘subscript𝑘𝑝𝑛𝑀𝛿delimited-[]𝑙subscript𝑙𝑝𝑚𝑁\displaystyle\sum\limits_{n,m\in{\mathbb{Z}}}e^{j2\pi\frac{nl_{p}}{N}}\,\delta% [k-k_{p}-nM]\,\delta[l-l_{p}-mN].∑ start_POSTSUBSCRIPT italic_n , italic_m ∈ blackboard_Z end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_n italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG end_POSTSUPERSCRIPT italic_δ [ italic_k - italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT - italic_n italic_M ] italic_δ [ italic_l - italic_l start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT - italic_m italic_N ] . (8)

We then lift the discrete DD domain signal xdd[k,l]subscript𝑥dd𝑘𝑙x_{\mbox{\scriptsize{dd}}}[k,l]italic_x start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] given by (1) to obtain a continuous DD domain signal

xdd(τ,ν)=k,lxdd[k,l]δ(τkτp/M)δ(νlνp/N)subscript𝑥dd𝜏𝜈subscript𝑘𝑙subscript𝑥dd𝑘𝑙𝛿𝜏𝑘subscript𝜏𝑝𝑀𝛿𝜈𝑙subscript𝜈𝑝𝑁\displaystyle x_{\mbox{\scriptsize{dd}}}(\tau,\nu)=\sum\limits_{k,l\in{\mathbb% {Z}}}x_{\mbox{\scriptsize{dd}}}[k,l]\delta(\tau-k\tau_{p}/M)\delta(\nu-l\nu_{p% }/N)italic_x start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT ( italic_τ , italic_ν ) = ∑ start_POSTSUBSCRIPT italic_k , italic_l ∈ blackboard_Z end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] italic_δ ( italic_τ - italic_k italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT / italic_M ) italic_δ ( italic_ν - italic_l italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT / italic_N ) (9)

that is quasi-periodic with delay period τpsubscript𝜏𝑝\tau_{p}italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT and Doppler peirod νp=1/τpsubscript𝜈𝑝1subscript𝜏𝑝\nu_{p}=1/\tau_{p}italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 1 / italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. We then apply a pulse shaping filter wtx(τ,ν)subscript𝑤𝑡𝑥𝜏𝜈w_{tx}(\tau,\nu)italic_w start_POSTSUBSCRIPT italic_t italic_x end_POSTSUBSCRIPT ( italic_τ , italic_ν ) to limit the transmitted TD signal to the duration and bandwidth of the Zak-OTFS subframe. This TD signal is obtained by applying the inverse Zak transform to the filtered DD domain signal wtx(τ,ν)σxdd(τ,ν)subscript𝜎subscript𝑤𝑡𝑥𝜏𝜈subscript𝑥dd𝜏𝜈w_{tx}(\tau,\nu)\,*_{\sigma}\,x_{\mbox{\scriptsize{dd}}}(\tau,\nu)italic_w start_POSTSUBSCRIPT italic_t italic_x end_POSTSUBSCRIPT ( italic_τ , italic_ν ) ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT ( italic_τ , italic_ν ).

Matched filtering at the receiver using wrx(τ,ν)=wtx(τ,ν)ej2πντsubscript𝑤𝑟𝑥𝜏𝜈superscriptsubscript𝑤𝑡𝑥𝜏𝜈superscript𝑒𝑗2𝜋𝜈𝜏w_{rx}(\tau,\nu)=w_{tx}^{*}(-\tau,-\nu)\,e^{j2\pi\nu\tau}italic_w start_POSTSUBSCRIPT italic_r italic_x end_POSTSUBSCRIPT ( italic_τ , italic_ν ) = italic_w start_POSTSUBSCRIPT italic_t italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( - italic_τ , - italic_ν ) italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π italic_ν italic_τ end_POSTSUPERSCRIPT followed by sampling on the information grid results in the discrete DD domain signal given by

ydd[k,l]subscript𝑦dd𝑘𝑙\displaystyle y_{\mbox{\scriptsize{dd}}}[k,l]italic_y start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== Edheff[k,l]σxd,dd[k,l]Rx. data signalsubscriptsubscript𝜎subscript𝐸𝑑subscripteff𝑘𝑙subscript𝑥d,dd𝑘𝑙Rx. data signal\displaystyle\underbrace{\sqrt{E_{d}}\,h_{\mbox{\scriptsize{eff}}}[k,l]*_{% \sigma}x_{\mbox{\scriptsize{d,dd}}}[k,l]}_{\mbox{\small{Rx. data signal}}}under⏟ start_ARG square-root start_ARG italic_E start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT end_ARG italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] end_ARG start_POSTSUBSCRIPT Rx. data signal end_POSTSUBSCRIPT (10)
+Epheff[k,l]σxs,dd[k,l]Rx. sensing/pilot signal+ndd[k,l],subscriptsubscript𝜎subscript𝐸𝑝subscripteff𝑘𝑙subscript𝑥s,dd𝑘𝑙Rx. sensing/pilot signalsubscript𝑛dd𝑘𝑙\displaystyle+\underbrace{\sqrt{E_{p}}h_{\mbox{\scriptsize{eff}}}[k,l]*_{% \sigma}x_{\mbox{\scriptsize{s,dd}}}[k,l]}_{\mbox{\small{Rx. sensing/pilot % signal}}}\,+\,n_{\mbox{\scriptsize{dd}}}[k,l],+ under⏟ start_ARG square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] end_ARG start_POSTSUBSCRIPT Rx. sensing/pilot signal end_POSTSUBSCRIPT + italic_n start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] ,

where ndd[k,l]subscript𝑛dd𝑘𝑙n_{\mbox{\scriptsize{dd}}}[k,l]italic_n start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] represents noise in the discrete DD domain. Note that it is associativity of twisted convolution that allows us to represent a cascade of filters/channels as a single effective channel filter heff[k,l]subscripteff𝑘𝑙h_{\mbox{\scriptsize{eff}}}[k,l]italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ]. In our previous work [19] we have described how the receiver senses the channel (estimates heff[k,l]subscripteff𝑘𝑙h_{\mbox{\scriptsize{eff}}}[k,l]italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ]) from the cross-ambiguity Ay,xs[k,l]subscript𝐴𝑦subscript𝑥𝑠𝑘𝑙A_{y,x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_y , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] between ydd[k,l]subscript𝑦dd𝑘𝑙y_{\mbox{\scriptsize{dd}}}[k,l]italic_y start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] and the transmitted spread pilot xs,dd[k,l]subscript𝑥s,dd𝑘𝑙x_{\mbox{\scriptsize{s,dd}}}[k,l]italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ]. We recall (Theorem 3333 from [19]) that

Ay,xs[k,l]subscript𝐴𝑦subscript𝑥𝑠𝑘𝑙\displaystyle{A}_{y,x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_y , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== Epheff[k,l]σAxs,xs[k,l]subscript𝜎subscript𝐸𝑝subscripteff𝑘𝑙subscript𝐴subscript𝑥𝑠subscript𝑥𝑠𝑘𝑙\displaystyle\sqrt{E_{p}}\,h_{\mbox{\scriptsize{eff}}}[k,l]\,*_{\sigma}\,{A}_{% x_{s},x_{s}}[k,l]square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ]
+Edheff[k,l]σAxd,xs[k,l]+An,xs[k,l],subscript𝜎subscript𝐸𝑑subscripteff𝑘𝑙subscript𝐴subscript𝑥𝑑subscript𝑥𝑠𝑘𝑙subscript𝐴𝑛subscript𝑥𝑠𝑘𝑙\displaystyle+\sqrt{E_{d}}\,h_{\mbox{\scriptsize{eff}}}[k,l]\,*_{\sigma}\,{A}_% {x_{d},x_{s}}[k,l]\,+\,{A}_{n,x_{s}}[k,l],+ square-root start_ARG italic_E start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT end_ARG italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] + italic_A start_POSTSUBSCRIPT italic_n , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] ,

where Axs,xs[k,l]subscript𝐴subscript𝑥𝑠subscript𝑥𝑠𝑘𝑙{A}_{x_{s},x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] is the self-ambiguity function of the spread pilot signal, Axd,xs[k,l]subscript𝐴subscript𝑥𝑑subscript𝑥𝑠𝑘𝑙{A}_{x_{d},x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] and An,xs[k,l]subscript𝐴𝑛subscript𝑥𝑠𝑘𝑙{A}_{n,x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_n , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] are the cross-ambiguity functions of the transmitted data signal xd,dd[k,l]subscript𝑥d,dd𝑘𝑙x_{\mbox{\scriptsize{d,dd}}}[k,l]italic_x start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] and noise signal ndd[k,l]subscript𝑛dd𝑘𝑙n_{\mbox{\scriptsize{dd}}}[k,l]italic_n start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] respectively with the spread pilot signal.

The self-ambiguity function of the point pilot xp,dd[k,l]subscript𝑥p,dd𝑘𝑙x_{\mbox{\scriptsize{p,dd}}}[k,l]italic_x start_POSTSUBSCRIPT p,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] is supported on the period lattice ΛpsubscriptΛ𝑝\Lambda_{p}roman_Λ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT comprising integer linear combinations of (τp,0)subscript𝜏𝑝0(\tau_{p},0)( italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , 0 ) and (0,νp)0subscript𝜈𝑝(0,\nu_{p})( 0 , italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ).

We have shown [19] that the self-ambiguity function of the spread pilot is supported on a lattice ΛΛ\Lambdaroman_Λ obtained by rotating ΛpsubscriptΛ𝑝\Lambda_{p}roman_Λ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. Let 𝒮𝒮{\mathcal{S}}caligraphic_S denote the support region of the effective channel in the discrete DD domain. The crystallization conditions with respect to the lattice ΛΛ\Lambdaroman_Λ are satisfied if the translates of 𝒮𝒮{\mathcal{S}}caligraphic_S by lattice points in ΛΛ\Lambdaroman_Λ are disjoint. In this case, we can obtain the effective channel tap heff[k,l]subscripteff𝑘𝑙h_{\mbox{\scriptsize{eff}}}[k,l]italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] by evaluating the first term of (II) at points (k,l)𝑘𝑙(k,l)( italic_k , italic_l ) in a fundamental domain of ΛΛ\Lambdaroman_Λ. The second term in (II) represents interference to sensing from data. We recall (Theorem 4444 from [19]) that the magnitudes |Axd,xs[k,l]|subscript𝐴subscript𝑥𝑑subscript𝑥𝑠𝑘𝑙|{A}_{x_{d},x_{s}}[k,l]|| italic_A start_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] | are essentially independent of k,l𝑘𝑙k,litalic_k , italic_l so that interference from data to sensing is noise-like (see also Fig.18181818 from [19]). The third term in (II) represents interference to sensing from noise. We recall (Appendix J from [19]) that An,xs[k,l]subscript𝐴𝑛subscript𝑥𝑠𝑘𝑙{A}_{n,x_{s}}[k,l]italic_A start_POSTSUBSCRIPT italic_n , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] is zero-mean Gaussian distributed with variance essentially independent of k,l𝑘𝑙k,litalic_k , italic_l. Our estimate of heff[k,l]subscripteff𝑘𝑙h_{\mbox{\scriptsize{eff}}}[k,l]italic_h start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] is then

h^eff[k,l]subscript^eff𝑘𝑙\displaystyle{\widehat{h}}_{\mbox{\scriptsize{eff}}}[k,l]over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] =\displaystyle== Ay,xs[k,l]Epfork,l𝒮.subscript𝐴𝑦subscript𝑥𝑠𝑘𝑙subscript𝐸𝑝for𝑘𝑙𝒮\displaystyle\frac{{A}_{y,x_{s}}[k,l]}{\sqrt{E_{p}}}\,\,\,\mbox{\small{for}}\,% \,k,l\in{\mathcal{S}}.divide start_ARG italic_A start_POSTSUBSCRIPT italic_y , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] end_ARG start_ARG square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG end_ARG for italic_k , italic_l ∈ caligraphic_S . (12)

We suppose that the crystallization conditions are also satisfied with respect to the lattice ΛpsubscriptΛ𝑝\Lambda_{p}roman_Λ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. We use the estimate (12) to cancel the contribution made by the pilot to the received DD signal. After cancellation, the signal

ydd[k,l]Eph^eff[k,l]σxs,dd[k,l].subscript𝑦dd𝑘𝑙subscript𝜎subscript𝐸𝑝subscript^eff𝑘𝑙subscript𝑥s,dd𝑘𝑙\displaystyle y_{\mbox{\scriptsize{dd}}}[k,l]-\sqrt{E_{p}}\,{\widehat{h}}_{% \mbox{\scriptsize{eff}}}[k,l]*_{\sigma}x_{\mbox{\scriptsize{s,dd}}}[k,l].italic_y start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] - square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] . (13)

is used to recover the data using the matrix-vector formulation of the Zak-OTFS I/O relation (see [15] for details).

By spreading the pilot signal, we integrate sensing and communication within the same Zak-OTFS subframe. Sharing DD domain resources in this way increases effective throughput compared with traditional approaches that use guard bands to divide DD domain resources between sensing and communication (see Fig.28282828 from [19]). Spreading also reduces the PAPR of the pilot signal to about 5555 dB compared with 15151515 dB for the original point pilot (see Fig.10101010 and Fig.11111111 from [19]).

We recall that integrated sensing and communication with spread pilots results in an uncoded 4444-QAM BER of about 102superscript10210^{-2}10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT, compared with 105superscript10510^{-5}10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT for sensing and communication in separate Zak-OTFS subframes (see Fig.26262626 from [19]). The difference is three orders of magnitude. Section III describes a turbo signal processing method that is able to close this gap.

III Iterative Sensing and Communication

In this Section we describe the iterative signal processing algorithm illustrated in Fig. 1. The first iteration is described in Section II and the tthsuperscript𝑡𝑡t^{th}italic_t start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT iteration consists of four steps.

STEP 1111: Inputs are the detected information symbols x^(t1)[k,l]superscript^𝑥𝑡1𝑘𝑙{\widehat{x}}^{(t-1)}[k,l]over^ start_ARG italic_x end_ARG start_POSTSUPERSCRIPT ( italic_t - 1 ) end_POSTSUPERSCRIPT [ italic_k , italic_l ], k=0,1,,M1𝑘01𝑀1k=0,1,\cdots,M-1italic_k = 0 , 1 , ⋯ , italic_M - 1, l=0,1,,N1𝑙01𝑁1l=0,1,\cdots,N-1italic_l = 0 , 1 , ⋯ , italic_N - 1 and estimated channel filter taps h^eff(t1)[k,l],(k,l)𝒮superscriptsubscript^eff𝑡1𝑘𝑙𝑘𝑙𝒮{\widehat{h}}_{\mbox{\scriptsize{eff}}}^{(t-1)}[k,l],(k,l)\in{\mathcal{S}}over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t - 1 ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] , ( italic_k , italic_l ) ∈ caligraphic_S from iteration t1𝑡1t-1italic_t - 1. We form y^d,dd(t)[k,l]superscriptsubscript^𝑦d,dd𝑡𝑘𝑙{\widehat{y}}_{\mbox{\scriptsize{d,dd}}}^{(t)}[k,l]over^ start_ARG italic_y end_ARG start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] (see (LABEL:eqn196402) at top of page) and subtract this estimate for the received data signal from the received signal to obtain

ys(t)[k,l]superscriptsubscript𝑦𝑠𝑡𝑘𝑙\displaystyle y_{s}^{(t)}[k,l]italic_y start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] =\displaystyle== ydd[k,l]y^d,dd(t)[k,l].subscript𝑦dd𝑘𝑙superscriptsubscript^𝑦d,dd𝑡𝑘𝑙\displaystyle y_{\mbox{\scriptsize{dd}}}[k,l]\,-\,{\widehat{y}}_{\mbox{% \scriptsize{d,dd}}}^{(t)}[k,l].italic_y start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] - over^ start_ARG italic_y end_ARG start_POSTSUBSCRIPT d,dd end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] . (15)

The output ys(t)[k,l]superscriptsubscript𝑦𝑠𝑡𝑘𝑙y_{s}^{(t)}[k,l]italic_y start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] is a data cancelled signal that approximates the received spread pilot.

STEP 2222: We use (12) to form the tthsuperscript𝑡𝑡t^{th}italic_t start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT estimate h^eff(t)[k,l]superscriptsubscript^eff𝑡𝑘𝑙{\widehat{h}}_{\mbox{\scriptsize{eff}}}^{(t)}[k,l]over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] for the effective channel filter taps

h^eff(t)[k,l]superscriptsubscript^eff𝑡𝑘𝑙\displaystyle{\widehat{h}}_{\mbox{\scriptsize{eff}}}^{(t)}[k,l]over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] =\displaystyle== Ays(t),xs[k,l]Epfor(k,l)𝒮.subscript𝐴superscriptsubscript𝑦𝑠𝑡subscript𝑥𝑠𝑘𝑙subscript𝐸𝑝for𝑘𝑙𝒮\displaystyle\frac{A_{y_{s}^{(t)},x_{s}}[k,l]}{\sqrt{E_{p}}}\,\,\mbox{\small{% for}}\,\,(k,l)\in{\mathcal{S}}.divide start_ARG italic_A start_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT , italic_x start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_k , italic_l ] end_ARG start_ARG square-root start_ARG italic_E start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG end_ARG for ( italic_k , italic_l ) ∈ caligraphic_S . (16)

STEP 3333: We use h^eff(t)[k,l]superscriptsubscript^eff𝑡𝑘𝑙{\widehat{h}}_{\mbox{\scriptsize{eff}}}^{(t)}[k,l]over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ], (k,l)𝕊𝑘𝑙𝕊(k,l)\in{\mathbb{S}}( italic_k , italic_l ) ∈ blackboard_S to form the tthsuperscript𝑡𝑡t^{th}italic_t start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT estimate y^s,dd(t)[k,l]superscriptsubscript^𝑦s,dd𝑡𝑘𝑙{\widehat{y}}_{\mbox{\scriptsize{s,dd}}}^{(t)}[k,l]over^ start_ARG italic_y end_ARG start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] of the received pilot signal using

y^s,dd(t)[k,l]superscriptsubscript^𝑦s,dd𝑡𝑘𝑙\displaystyle{\widehat{y}}_{\mbox{\scriptsize{s,dd}}}^{(t)}[k,l]over^ start_ARG italic_y end_ARG start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] =\displaystyle== h^eff(t)[k,l]σxs,dd[k,l].subscript𝜎superscriptsubscript^eff𝑡𝑘𝑙subscript𝑥s,dd𝑘𝑙\displaystyle{\widehat{h}}_{\mbox{\scriptsize{eff}}}^{(t)}[k,l]\,*_{\sigma}\,{% x}_{\mbox{\scriptsize{s,dd}}}[k,l].over^ start_ARG italic_h end_ARG start_POSTSUBSCRIPT eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] ∗ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT [ italic_k , italic_l ] . (17)

We subtract this estimate from the received signal to obtain

yd(t)[k,l]superscriptsubscript𝑦𝑑𝑡𝑘𝑙\displaystyle y_{d}^{(t)}[k,l]italic_y start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] =\displaystyle== ydd[k,l]y^s,dd(t)[k,l].subscript𝑦dd𝑘𝑙superscriptsubscript^𝑦s,dd𝑡𝑘𝑙\displaystyle y_{\mbox{\scriptsize{dd}}}[k,l]\,-\,{\widehat{y}}_{\mbox{% \scriptsize{s,dd}}}^{(t)}[k,l].italic_y start_POSTSUBSCRIPT dd end_POSTSUBSCRIPT [ italic_k , italic_l ] - over^ start_ARG italic_y end_ARG start_POSTSUBSCRIPT s,dd end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] . (18)

This output is a pilot cancelled signal that approximates the received data signal.

STEP 4444: We detect data/information symbols x^(t)[k,l]superscript^𝑥𝑡𝑘𝑙\widehat{x}^{(t)}[k,l]over^ start_ARG italic_x end_ARG start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] from yd(t)[k,l]superscriptsubscript𝑦𝑑𝑡𝑘𝑙y_{d}^{(t)}[k,l]italic_y start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT [ italic_k , italic_l ] using the method described in Section II. We then move to STEP 1111 of iteration t+1𝑡1t+1italic_t + 1, and the algorithm halts after a fixed number of iterations.

Section IV presents numerical simulations illustrating that multiple iterations improve the fidelity of the data cancelled and pilot cancelled signals.

TABLE I: Power-delay profile of Veh-A channel model
Path number (i𝑖iitalic_i) 1 2 3 4 5 6
τisubscript𝜏𝑖\tau_{i}italic_τ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (μ𝜇\muitalic_μs) 0 0.31 0.71 1.09 1.73 2.51
Relative power (pisubscript𝑝𝑖p_{i}italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT) dB 0 -1 -9 -10 -15 -20

IV Numerical simulations

This Section reports simulation results for the Veh-A channel model [18] which consists of six channel paths. The channel gains hisubscript𝑖h_{i}italic_h start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT are modeled as independent zero-mean complex circularly symmetric Gaussian random variables, normalized so that i=16𝔼[|hi|2]=1superscriptsubscript𝑖16𝔼delimited-[]superscriptsubscript𝑖21\sum\limits_{i=1}^{6}{\mathbb{E}}\left[|h_{i}|^{2}\right]=1∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT blackboard_E [ | italic_h start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] = 1. Table I lists the power-delay profile for the six channel paths. The Doppler shift of the i𝑖iitalic_i-th path is modeled as νi=νmaxcos(θi)subscript𝜈𝑖subscript𝜈𝑚𝑎𝑥subscript𝜃𝑖\nu_{i}=\nu_{max}\cos(\theta_{i})italic_ν start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_ν start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT roman_cos ( italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ), where νmaxsubscript𝜈𝑚𝑎𝑥\nu_{max}italic_ν start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT is the maximum Doppler shift of any path, and the variables θisubscript𝜃𝑖\theta_{i}italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, i=1,2,,6𝑖126i=1,2,\cdots,6italic_i = 1 , 2 , ⋯ , 6, are independent and distributed uniformly in [π,π)𝜋𝜋[-\pi\,,\,\pi)[ - italic_π , italic_π ).

We consider Zak-OTFS modulation with Doppler period νp=30subscript𝜈𝑝30\nu_{p}=30italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 30 KHz, delay period τp=1/νp=33.3subscript𝜏𝑝1subscript𝜈𝑝33.3\tau_{p}=1/\nu_{p}=33.3italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 1 / italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 33.3 ms, M=31𝑀31M=31italic_M = 31 and N=37𝑁37N=37italic_N = 37. The channel bandwidth B=Mνp=0.93𝐵𝑀subscript𝜈𝑝0.93B=M\nu_{p}=0.93italic_B = italic_M italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 0.93 MHz and the subframe duration T=Nτp=1.2𝑇𝑁subscript𝜏𝑝1.2T=N\tau_{p}=1.2italic_T = italic_N italic_τ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 1.2 ms. Note that the first three channel paths introduce delay shifts in the interval [0,0.71]00.71[0,0.71][ 0 , 0.71 ] ms and each is less than the delay resolution. We consider matched filtering using root raised cosine (RRC) pulse shaping filters with roll-off factors βν=βτ=0.6subscript𝛽𝜈subscript𝛽𝜏0.6\beta_{\nu}=\beta_{\tau}=0.6italic_β start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT = italic_β start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT = 0.6. This increases the subframe duration to T=(1+βν)T=1.92superscript𝑇1subscript𝛽𝜈𝑇1.92T^{\prime}=(1+\beta_{\nu})T=1.92italic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = ( 1 + italic_β start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) italic_T = 1.92 ms and the bandwidth to B=(1+βτ)B=1.40superscript𝐵1subscript𝛽𝜏𝐵1.40B^{\prime}=(1+\beta_{\tau})B=1.40italic_B start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = ( 1 + italic_β start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT ) italic_B = 1.40 MHz. We consider a spread pilot pulsone constructed using a chirp filter in the discrete delay-Doppler domain ((7) with q=3𝑞3q=3italic_q = 3) as described in [19]. We implement the turbo signal processing pipeline illustrated in Fig. 1 using MMSE equalization to recover information symbols using the effective channel matrix.

For a fixed data SNR of 25252525 dB we set the PDR to 10101010 dB so that the pilot power is sufficient to start the turbo process. Fig. 2 illustrates that in the crystalline regime, five turbo iterations (dashed curve with red triangles) are sufficient to match the BER performance of separate sensing and communications (dashed curve with black squares). When νmax>12subscript𝜈𝑚𝑎𝑥12\nu_{max}>12italic_ν start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT > 12 KHz the estimate of the effective channel becomes less accurate because of Doppler domain aliasing, and BER performance degrades because interference from the residual pilot becomes more significant than noise. We focus on the role of channel estimation by designing a reference system (solid blue curve with diamonds) where sensing takes place in a separate subframe (S|C¯conditional𝑆¯𝐶S|\overline{C}italic_S | over¯ start_ARG italic_C end_ARG) but data transmission is still subject to interference from the residual pilot (C|Sconditional𝐶𝑆C|Sitalic_C | italic_S). Fig. 2 illustrates that BER performance with five turbo iterations differs from that of the reference system by a small SNR offset.

Next we consider how BER performance depends on PDR, and we set νmax=6subscript𝜈𝑚𝑎𝑥6\nu_{max}=6italic_ν start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 6 KHz so that we are operating deep in the interior of the crystalline regime. Fig. 3 illustrates that when PDR <10absent10<-10< - 10 dB the pilot power is not sufficient to start the turbo process, that the initial estimate of the effective channel is not sufficiently accurate. When PDR >25absent25>25> 25 dB interference from the residual pilot (after cancellation) is more significant than noise, and BER degrades as PDR increases. For intermediate values of the PDR, BER performance improves with increasing PDR as estimates of the effective channel become more accurate. This explains the characteristic “U” shape of the intermediate curves in Fig. 3. Again, we focus on the role of channel estimation by considering the reference system described above (S|C¯conditional𝑆¯𝐶S|\overline{C}italic_S | over¯ start_ARG italic_C end_ARG and C|Sconditional𝐶𝑆C|Sitalic_C | italic_S), where the residual pilot has very little energy and interference offered to the transmitted data is inconsequential (channel estimation is very accurate since it is based on a separate sensing-only subframe). We also consider a second reference system (green curve with diamonds), where sensing is subject to interference from data (S|Cconditional𝑆𝐶S|{C}italic_S | italic_C), but data recovery is not subject to interference from a residual pilot (C|S¯conditional𝐶¯𝑆C|\overline{S}italic_C | over¯ start_ARG italic_S end_ARG). Data is transmitted in a separate subframe and there are no turbo iterations for this second reference system. Fig. 3 confirms that interference from the residual pilot is responsible for the degradation in BER performance with increasing PDR.

Refer to caption
Figure 2: Uncoded 4444-QAM BER as a function of increasing νmaxsubscript𝜈𝑚𝑎𝑥\nu_{max}italic_ν start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT. Veh-A channel, data SNR =25absent25=25= 25 dB, PDR =10absent10=10= 10 dB, Doppler period νp=30subscript𝜈𝑝30\nu_{p}=30italic_ν start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 30 KHz, M=31,N=37formulae-sequence𝑀31𝑁37M=31,N=37italic_M = 31 , italic_N = 37, RRC pulse shaping filter (βν=βτ=0.6subscript𝛽𝜈subscript𝛽𝜏0.6\beta_{\nu}=\beta_{\tau}=0.6italic_β start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT = italic_β start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT = 0.6).
Refer to caption
Figure 3: Uncoded 4444-QAM BER as a function of increasing PDR. Veh-A channel considered in Fig. 2.

V Conclusions

We started by observing that the uncoded 4-QAM BER performance of sensing and communications in separate Zak-OTFS subframes is three orders of magnitude better than that of integrated sensing and communication with spread pilots. In joint sensing and communication, we estimate the effective channel, then estimate the received spread pilot, then recover the data after subtracting our estimate for the received pilot from the received signal. The residual pilot (after cancellation) interferes with data transmission, and we showed that this is the reason for the three orders of magnitude gap in performance. We described a turbo signal processing algorithm that alternates between channel sensing using a spread pilot and data recovery. We showed that five turbo iterations suffice to match the performance of separate sensing and communication across a broad range of Doppler shifts.

References

  • [1] H. Tataria, M. Shafi, A. F. Molisch, M. Dohler, H. Sjöland, and F. Tufvesson, “6G wireless systems: Vision, requirements, challenges, insights, and opportunities,” Proceedings of the IEEE, vol. 109, no. 7, pp. 1166-1199, Jul. 2021.
  • [2]
  • [3] C. -X. Wang, X. You, X. Gao, X. Zhu, Z. Li, C. Zhang, H. Wang, Y. Huang, Y. Chen, H. Haas, J. S. Thompson, E. G. Larsson, M. Di Renzo, W. Tong, P. Zhu, X. Shen, H. V. Poor, and L. Hanzo, “On the road to 6G: Visions, requirements, key technologies, and testbeds,” IEEE Commun. Surveys & Tuts., vol. 25, no. 2, pp. 905-974, 2023.
  • [4] R. Hadani et al., “Orthogonal time frequency space modulation,” Proc. IEEE WCNC’2017, pp. 1-6, Mar. 2017.
  • [5] Z. Wei, W. Yuan, S. Li, J. Yuan, G. Bharatula, R. Hadani and L. Hanzo, “Orthogonal time-frequency space modulation: a promising next-generation waveform,” IEEE Wireless Commun. Mag., vol. 28, no. 4, pp. 136-144, Aug. 2021.
  • [6] W. Yuan et al., “Best readings in orthogonal time frequency space (OTFS) and delay Doppler signal processing,” Jun. 2022. https://www.comsoc.org/publications/best-readings/orthogonal-time-frequency-space-otfs-and-delay-doppler-signal-processing.
  • [7] A. Monk, R. Hadani, M. Tsatsanis, and S. Rakib, “OTFS - orthogonal time frequency space: a novel modulation meeting 5G high mobility and massive MIMO challenges,” arXiv:1608.02993 [cs.IT] 9 Aug. 2016.
  • [8] P. Raviteja, K. T. Phan, Y. Hong, and E. Viterbo, “Embedded pilot aided channel estimation for OTFS in delay-Doppler channels,” IEEE Trans. Veh. Tech., vol. 68, no. 5, pp. 4906-4917, May 2019.
  • [9] P. Raviteja, Y. Hong, E. Viterbo, and E. Biglieri, “Practical pulse-shaping waveforms for reduced-cyclic-prefix OTFS,” IEEE Trans. Veh. Tech., vol. 68, no. 1, pp. 957-961, Jan. 2019.
  • [10] C. Shen, J. Yuan, and H. Lin, “Error performance of rectangular pulse-shaped OTFS with practical receivers,” IEEE Wireless Commun. Lett., vol. 11, no. 12, pp. 2690-2694, Dec. 2022.
  • [11] S. K. Mohammed, “Derivation of OTFS modulation from first principles,” IEEE Trans. Veh. Tech., vol. 70, no. 8, pp. 7619-7636, Aug. 2021.
  • [12] S. K. Mohammed, “Time-domain to delay-Doppler domain conversion of OTFS signals in very high mobility scenarios,” IEEE Trans. Veh. Tech., vol. 70, no. 6, pp. 6178-6183, Jun. 2021.
  • [13] E. Panayirci, H. Senol, and H. V. Poor, “Joint channel estimation, equalization, and data detection for OFDM systems in the presence of very high mobility,” IEEE Trans. Signal Process., vol. 58, no. 8, pp. 4225-4238, Aug. 2010.
  • [14] S. K. Mohammed, R. Hadani, A. Chockalingam, and R. Calderbank, “OTFS – A mathematical foundation for communication and radar sensing in the delay-Doppler domain,” IEEE BITS the Information Theory Magazine, vol. 2, no. 2, pp. 36-55, 1 Nov. 2022.
  • [15] S. K. Mohammed, R. Hadani, A. Chockalingam, and R. Calderbank, “OTFS – Predictability in the delay-Doppler domain and its value to communications and radar sensing,” IEEE BITS the Information Theory Magazine, IEEE early access, doi: 10.1109/MBITS.2023.3319595, Sep. 2023.
  • [16] S. Li, W. Yuan, Z. Wei, J. Yuan, B. Bai, and G. Caire, “On the pulse shaping for delay-Doppler communications,” IEEE GLOBECOM’2023, pp. 4909-4914, Dec. 2023.
  • [17] S. Gopalam, I. B. Collings, S. V. Hanly, H. Inaltekin, S. R. B. Pillai, and P. Whiting, “Zak-OTFS implementation via time and frequency windowing,” IEEE Trans. Commun., IEEE early access, doi: 10.1109/TCOMM.2024.3366403.
  • [18] ITU-R M.1225, “Guidelines for evaluation of radio transmission technologies for IMT-2000,” International Telecommunication Union Radio communication, 1997.
  • [19] M. Ubadah, S. K. Mohammed, R. Hadani, S. Kons, A. Chockalingam, and R. Calderbank, “Zak-OTFS for integration of sensing and communication,” available online: arXiv:2404.04182v1 [eess.SP] 5 Apr 2024.
  • [20] C. Douillard et al., “Iterative correction of intersymbol interference: turbo-equalization,” Eur. Trans. Telecommun., vol. 6, pp. 507-511, Sep.-Oct. 1995.
  • [21] M. Tüchler, R. Koetter, and A. C. Singer, “Turbo equalization: principles and new results,” IEEE Trans. Commun., vol. 50, no. 5, pp. 754-767, May 2002.
  • [22] R. Otnes and M. Tuchler, “Iterative channel estimation for turbo equalization of time varying frequency selective channels,” IEEE Trans. Wireless Commun., vol. 3, pp. 1918-1923, Nov. 2004.