Using spread pulses in wireless communication

Abstract

Methods, apparatus, systems for digital communication are described. One example method includes generating, by a transmitter apparatus, a transmission waveform by applying a discrete filter to an input signal in a delay-Doppler domain, wherein the discrete filter is periodic in delay domain and in Doppler domain; and controlling a transceiver to transmit the transmission waveform over a channel.

Description

PCT Patent Application Attorney Docket No.: 119314.8123.WO00 USING SPREAD PULSES IN WIRELESS COMMUNICATION CROSS-REFERENCE TO RELATED APPLICATION

[0001] This application is claims priority to U.S. Provisional Application No.63 / 567,776, filed on March 20, 2024, the disclosure of which is hereby incorporated by reference herein in its entirety. TECHNICAL FIELD

[0002] The present document relates to digital communication. BACKGROUND

[0003] Due to an explosive growth in the number of wireless user devices and the amount of wireless data that these devices can generate or consume, current wireless communication networks are fast running out of bandwidth to accommodate such a high growth in data traffic and provide high quality of service to users.

[0004] Various efforts are underway in the telecommunication industry to come up with next generation of wireless technologies that can keep up with the demand on performance of wireless devices and networks. Many of those activities involve situations in which a large number of user devices may be served by a network. SUMMARY

[0005] This document discloses techniques that may be used by wireless networks to achieve several operational improvements. In particular, methods, apparatus, and systems that use spread pulse based digital communication are disclosed.

[0006] In one example aspect, a method of digital communication is disclosed. The method includes generating, by a transmitter apparatus, a transmission waveform by applying a discrete filter to an input signal in a delay-Doppler domain, wherein the discrete filter is periodic in delay domain and in Doppler domain; and controlling a transceiver to transmit the transmission waveform over a channel.

[0007] In another aspect, another method of digital communications is disclosed. The method includes generating a transmission waveform comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain, and wherein L is a positive integer; and controlling a transceiver to transmit the transmission waveform over a channel. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0008] In yet another aspect, a method of wireless communication is disclosed. The method includes generating K transmission waveforms, each of the K transmission waveforms comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain, and wherein K and L are positive integers and K > L; and transmitting, the each of the K transmission waveforms using a corresponding transmission antenna.

[0009] In yet another aspect, a method of digital communication is disclosed. The method includes receiving a received signal comprising L layers, each of the L layers comprising corresponding data and a corresponding spread pilot over a channel, wherein L>=1, and wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain; and extracting data from at least some of the L layers of the received waveform by subtracting, for each extracted layer, the corresponding spread pilot from the received signal and by performing a channel equalization based on a channel sensing using the corresponding spread pilot.

[0010] In yet another aspect, an apparatus is disclosed. The apparatus is configured to implement an above-described method.

[0011] In yet another example aspect, a wireless system in which one or more of the above-described methods are implemented is disclosed.

[0012] These, and other, features are described in this document. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] FIG.1 shows an example communication network.

[0014] FIG.2 shows a simplified example of a wireless communication system in which uplink and downlink transmissions are performed.

[0015] FIG.3 shows examples of future applications of wireless technologies.

[0016] FIGS.4A and 4B show examples of pulses that are used as basis signals for wireless communication.

[0017] FIG.4C is a visual depiction example of rotations undergone by a delay-Doppler domain pulse.

[0018] FIG.5 is an example of a quasi-periodic delay Doppler domain pulse.

[0019] FIG.6 depicts an example of channel distortions undergone by a pilot signal.

[0020] FIG.7 is a block diagram of an example of a Zak transform based wireless system including a transmitter and a receiver.

[0021] FIG.8 is a pictorial depiction of effects of a doubly spread channel.

[0022] FIG.9 is a graphical example of predictability achieved in a crystalline region of wireless communications. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0023] FIG.10 is a graph showing example results for model-free and model-dependent implementations.

[0024] FIG.11A-B are flowcharts of different example implementations of an integrated sensing and communication with point pulse-tone (SCPP) and spread pulse-tone.

[0025] FIG.12 shows an example heat map of results obtained when using a pilot signal without a guard band.

[0026] FIG.13 shows an example heat map of results obtained when using a pilot signal with a guard band.

[0027] FIG.14 is a graph of example results obtained for an SCPP implementation.

[0028] FIG.15 graphically depicts examples of peak-to-average power ratio (PAPR).

[0029] FIG.16 shows an example of Zak-OTFS transceiver processing.

[0030] FIG.17 shows an example of a Zak-OTFS carrier waveform.

[0031] FIG.18 graphically depicts examples where chirp filters rotate the period lattice.

[0032] FIG.19A graphically depicts an example time domain (TD) realization of a point pulse-tone.

[0033] FIG.19B graphically depicts a magnitude of the cross-ambiguity between an example transmitted point pilot signal and an example received point pilot signal.

[0034] FIG.19C graphically depicts a magnitude of the cross-ambiguity between an example transmitted point pilot signal and an example received signal.

[0035] FIG.19D graphically depicts a magnitude of the cross-ambiguity between an example transmitted pilot signal and a data component of an example received signal.

[0036] FIG.19E graphically depicts energy profiles of example delay-Doppler domain pilot signals.

[0037] FIG.19F graphically depicts a TD realization of an example spread pilot signal.

[0038] FIG.19G is a plot showing the complementary CDF of the instantaneous-to-average-power ratio (IAPR) for example spread and point pilots.

[0039] FIG.19H-19J graphically depicts the support of a self-ambiguity function of example spread pulse-tones.

[0040] FIG.19K graphically depicts example heat maps for an effective discrete delay-Doppler domain channel filter.

[0041] FIG.19L graphically depicts example heat maps for a cross-ambiguity.

[0042] FIG.19M graphically depicts example heat maps for a cross-ambiguity.

[0043] FIG.19N graphically depicts example heat map

[0044] FIG.20A is an example plot of bit error rate (BER) of uncoded 4-QAM as a function of increasing Doppler spread.

[0045] FIG.20B is an example plot of BER of uncoded 4-QAM as a function of increasing pilot to data power ratio. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0046] FIG.20C is an example plot of normalized mean squared error (NMSE) as a function of pilot to data power ratio.

[0047] FIG.20D is an example plot of NMSE as a function of increasing pilot to data power ratio.

[0048] FIG.20E is an example plot of BER of uncoded 4-QAM as a function of increasing pilot to data power ratio.

[0049] FIG.20F is an example plot of signal to interference ratio (SIR) as a function of increasing pilot to data power ratio.

[0050] FIG.20G shows a graph of an example BER performance.

[0051] FIG.20H shows a graph of an example BER performance.

[0052] FIG.20I is a graph showing example effective throughput results.

[0053] FIG.20J are graphs showing examples of spreading operations in time domain (TD) and delay- Doppler (DD) domain.

[0054] FIG.20K is a graph showing PAPR result examples.

[0055] FIG.21 is a graphical example of a sensing waveform.

[0056] FIG.22 graphically depicts examples of effect of spreading.

[0057] FIG.23 graphically depicts examples of effect of spreading.

[0058] FIG.24 shows an example relationship among filters in discrete DD domain.

[0059] FIG.25 shows an example of multiplexing data with a spread pilot signal.

[0060] FIG.26 is a graphical depiction of an example of sensing accuracy as a function of pulse dynamic range (PDR).

[0061] FIG.27 is a graph showing example interference from a residual spread pilot.

[0062] FIG.28 shows a graph of an example bit error rate (BER) performance.

[0063] FIG.29 shows a graph of an example BER performance as a function of maximum Doppler spread.

[0064] FIG.30 is a graph showing examples of geometric criteria for filter design.

[0065] FIG.31 is a graph showing example performance of using spread pilots.

[0066] FIG.32 is a graph showing example results obtained in a turbo-aided communication and sensing implementation.

[0067] FIG.33 shows an example of a hardware platform.

[0068] FIGS.34A-34D are flowcharts for example methods of facilitating digital communication. DETAILED DESCRIPTION

[0069] To make the purposes, technical solutions and advantages of this disclosure more apparent, various embodiments are described in detail below with reference to the drawings. Unless otherwise 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 noted, embodiments and features in embodiments of the present document may be combined with each other.

[0070] Section headings are used in the present document to improve readability of the description and do not in any way limit the discussion or the embodiments to the respective sections only. Furthermore, certain standard-specific terms are used for illustrative purpose only, and the disclosed techniques are applicable to any wireless and / or digital communication systems.

[0071] 1. Introduction – wireless communication environment examples

[0072] The wireless or time-variant nature of the communication channel poses several challenges in design a transmission protocol suitable for wireless communication scenarios. These days, users expect their wireless devices to work everywhere and in a variety of mobile or stationary situations. Presently, the communication industry is working on a next generation of digital communication technology, sometimes called 6G (sixth generation) wireless technology.

[0073] At the same time, Artificial Intelligence AI and Communications: Machine Learning has revolutionized image and natural language processing - data-driven discovery has revolutionized bioinformatics. If Machine Learning is to revolutionize wireless then it needs to work at the speed of wireless.

[0074] It is expected that the next generation wireless systems will adopt Integrated Sensing and Communication (ISAC) in which signals used to transmit data interfere with signals used to sense Input / Output relations.

[0075] 2. Example wireless systems

[0076] FIG.1 shows an example of a wireless communication system 100 in which a transmitter device 102 transmits signals to a receiver 104. The signals may undergo various wireless channels and multipaths, as depicted. Some reflectors such as buildings and trees may be static, while others such as cars, may be moving scatterers. The transmitter device 102 may be, for example, a user device, a mobile phone, a tablet, a computer, or another Internet of Things (IoT) device such as a smartwatch, a camera, and so on. The receiver device 104 may be a network device such as the base station. The signals transmitted from the base station to the transmitter 102 may experience similar channel degradations produced by static or moving scatterers. The techniques described in the present document may be implemented by the devices in the wireless communication system 100. The terms “transmitter” and “receiver” are simply used for convenience of explanation. As further described herein, depending on the direction of transmission (uplink or downlink), the network station may be transmitting or receiving, and / or the user device may be receiving or transmitting.

[0077] FIG.2 shows a simplified wireless network to highlight certain aspects of the disclosed technology. A transmitter transmits wireless signals to a receiver in the wireless network. Some transmissions in the network, variously called as downlink or downstream transmissions, a network-side node such as a base station acts as a transmitter of wireless signals and one or more user devices act as 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 the receiver of these wireless signals. For some other transmissions, as depicted in FIG.2, the direction of transmission may be reversed. Such transmissions are often called uplink or upstream transmissions. For such transmissions, one or more user devices act as transmitters of the wireless signals and a network-side node such as a base station acts as the receiver of these signals (as depicted in FIG.2). Other type of transmissions in the network may include device-to-device transmissions, sometimes called direct or sideband transmissions. While the present document primarily uses the terms “downlink” and “uplink” for the sake of convenience, similar techniques may also be used for other situations in which transmissions in two directions are performed - e.g., inbound, or incoming transmissions that are received by a wireless device and outbound or outgoing transmissions that are transmitted by a wireless device. For example, downlink transmissions may be inbound transmissions for a user device, while outbound transmissions for a network device. Similarly, uplink transmission may be inbound transmissions for a network device while outbound transmissions from a wireless device. Therefore, for some embodiments, the disclosed techniques may also be described using terms such as “inbound” and “outbound” transmission without importing any 3GPP-specific or other wireless protocol-specific meaning to the terms “uplink” and “downlink.”

[0078] In frequency division multiplexing (FDM) networks, the transmissions to a base station and the transmissions from the base station may occupy different frequency bands (each of which may occupy continuous or discontinuous spectrum). In time division multiplexing (TDM) networks, the transmissions to a base station and the transmissions from the base station occupy a same frequency band but are separated in time domain using a TDM mechanism such as time slot-based transmissions.

[0079] 3. Brief introduction to use of Spread pulses and delay Doppler representations

[0080] FIG.3 shows examples of future applications of wireless technologies such as with Artificial Intelligence (AI) and Communication along with Integrated Sensing and Communication (ISAC) for 6G and beyond. As described above, with AI & Communication, machine learning algorithms have revolutionized image and natural language processing, but if they are to revolutionize wireless then they need to learn at the speed of wireless. With ISAC, coexistence in the same subframe can increase effective throughput but requires minimizing interference between sensing and data transmission. These are some of the research challenges for 6G, which the present document aims to address.

[0081] Huge demand for high-speed data led to the transition from code division multiple access CDMA to orthogonal frequency division multiplexing OFDM, and to the opportunity to think about the benefits of measuring and adapting to instantaneous channel gains.

[0082] Recently, new communication technology, called orthogonal time frequency space (OTFS) was mainly investigated in the context of communication – e.g., using the OTFS waveform for conducting robust and error-free communication through all kinds of complicated channel conditions. The present document discloses use of the OTFS waveform in a broader context of joint communication and sensing. The present document not only analyzes the advantages of OTFS as a way for doing communication, but 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 also as a way for doing radar imaging of the environment and going even beyond that and analyzing it as a way for doing radar and communication jointly on the same band in a very efficient and precise manner.

[0083] In 6G, there is a possibility that radar theory will converge with communication. OTFS technology offers this advantage, as is further disclosed in the present document. As further disclosed in the present document, a waveform that combines desirable properties, called pulse-tone, may be used as a basic building block for digital communication.

[0084] In particular, a spread waveform called “spread pulse-tone” is disclosed in the present document. This waveform exhibits all the good properties of OTFS and can also be used like a radar signal for channel sensing. Furthermore, the present document also discloses ways to achieve a control over the peak to average power ratio (PAPR), to allow implementation of the disclosed technology in a practically efficient manner.

[0085] In one advantageous aspect, spread pulse-tone signals have weak correlation with each other. So, they interfere with one another in a code division multiple access CDMA kind of manner. Systems can actually put spread pulse-tone signals on top of one another, conduct whatever thing you want to do with each of them, and then separate them in the receiver like in the case of CDMA. One way to look at the spreading is as a kind of code which is associated with the pulse-tone.

[0086] The disclosed waveforms can offer at least the following three distinct operational advantages. First, the waveform is amenable to joint communication and sensing. Second, the waveform represents a large family of waveforms without sacrificing anything on the level of performance, while opening up new capabilities of multiplexing. Third, because this waveform behaves like separable coding (e.g., CDMA) because of their weak correlation with each other, they can be used simultaneously for data transmission and sensing, without causing any significant interference to each other.

[0087] In the next generation of wireless technology, low latency of communication may be an important feature. Low latency also enables moving functions like scheduling towards the cloud. Also, it would be desirable for a machine learning application and AI application to operate on a physical layer that is predictable and has a long-time horizon. The disclosed technology offers such a solution because a representation of the propagation environment in delay Doppler can provide such features. You cannot get a representation that changes any slower because it's change captured at the speed at which physics changes. And this ability allows to on a 1 millisecond frame to acquire the environment and acquire the data; that is an enabler of machine learning.

[0088] 4. Introduction to pulse-tone waveforms

[0089] FIGS.4A and 4B show examples of pulses that can be used as basis signals for wireless communication.

[0090] In FIG.4A, an example of separating sensing and communications is depicted with a pilot signal (e.g., point pulse-tone x=xs), where the Cross-Ambiguity can be represented as Ay,x[k,l]. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0091] In some embodiments, separate Zak-OTFS subframes may be dedicated to sensing and data transmission / communication. A transmitter may transmit a pilot, the receiver may read off the I / O relation for the entire subframe from the response to a single pilot, and then use this estimate to equalize the data in a subsequent data subframe. For example, as depicted in FIG.4A, the discrete input / output (I / O) can be read off from the response to a single pilot signal transmitted through a linear time variant (LTV) channel, and the channel estimate can be used to recover data in a subsequent Zak-OTFS subframe.

[0092] How does this work? - embodiments can use linear time invariant LTI channel models as depicted in FIG.4B. A time domain TD pulse is a geometric mode of the LTI channel – the I / O relation is made up of delays that move the pulse about in time. The I / O relation can be read off from the response to a single TD pulse. After waiting longer than the delay spread, embodiments can read off the I / O relation, and equalize the data. Because delay-Doppler profile of a channel typically is static over long time periods, this approach can provide predictability, which means that the I / O relation does not change with time. And predictability can enable a model-free mode of operation such as under the crystallization condition where delay spread < delay period.

[0093] FIG.4C is a visual depiction example of rotations undergone by a delay-Doppler domain pulse. In FIG.4C, a pulse in the delay-Doppler domain is shown, which can, for example, have a delay spread ^p= 50 ^s and a Doppler spread ^p= 20 KHz.

[0094] Defining a pulse in the delay-Doppler domain runs into the problem that the Heisenberg Uncertainty Principle implies that it is not possible to simultaneously localize a signal in delay and Doppler (DD). However, it is possible to “cheat” or work around the Heisenberg Uncertainty Principle by starting with a pulse in the DD domain and extending it quasi-periodically.

[0095] The DD domain pulse is in fact a configuration of infinitely many pulses which repeat at integer multiples of the delay period ^palong the delay axis and at integer multiples of the Doppler period ^palong the Doppler axis. Refer to the box of width ^pand height ^pas the fundamental period, and take^^p^^p= 1. [e.g., the fundamental domain can be defined by delay period ^^^and the Doppler period ^^^.]

[0096] The phase of the pulse changes when the pulse location shifts by an integer multiple of ^palong the delay axis, but there is no change in phase when the pulse location shifts by an integer multiple of ^palong the Doppler axis.

[0097] Accordingly, the DD realization of a TD signal is a quasi-periodic function.

[0098] FIG.5 is an example of a quasi-periodic delay-Doppler domain pulse.

[0099] How to get from a DD domain pulse to a TD pulse-tone waveform? – This may be achieved by applying an (inverse time) Zak transform. That is, an (inverse time) Zak transform may be applied to get a TD pulse-tone from a Quasi-Periodic DD Domain Pulse. There is also an inverse (frequency) Zak 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 transform that transforms a DD domain pulse to a frequency domain FD pulse-tone. Additional details are provided in following section 5.

[0100] Here a DD domain pulse is located at (^0, ^0) within the fundamental period. The DD domain pulse is spread along the delay axis over a length 1 / B (i.e., B-1) and is spread along the Doppler axis over a length 1 / T (i.e., T-1). The TD realization x(t) can be obtained by applying the Zak transform to the DD domain pulse ^^^^^^(^^, ^^). The TD realization x(t) is a pulse train of finite duration T, with each pulse in the train spread over time duration 1 / B. Consecutive pulses of the pulse train are separated by the delay period ^p– moving the pulse location ^0along the delay axis displaces the TD pulse-tone in time. The pulse train is modulated by a sinusoid of frequency ^0– moving the pulse location ^0along the Doppler axis displaces the modulating tone in frequency.

[0101] TDM (time division multiplexing) is a limiting case. As the delay period ^pgrows, the TD pulses grow further apart, the tone structure disappears, and only a single TD pulse at ^0remains.

[0102] FDM (frequency division multiplexing) is a limiting case. As the delay period ^pshrinks, the pulses grow closer and closer together in time, and in the limit, the pulse structure disappears and only the tone remains. Pulse-tones parametrized by ^pand ^pinterpolate between TDM and FDM; that is, they live on the hyperbola ^p^^p= 1. What value of ^poptimizes predictability?

[0103] Pulse-tones are optimal as time- and band- limited signals because there is almost no overlap between two DD pulses whose delay domain locations differ by 1 / B or whose Doppler locations differ by 1 / T.

[0104] In the example depicted in FIG.5, B= 1 / 100ns = 10MHz, T = 1 / 1kHz = 1ms, N = 50μs / 100ns = 500, and M = 20kHz / 1kHz = 20, with ^p= 20 KHz (Doppler Spread) and ^p= 50 ms (Delay Spread).

[0105] FIG.6 depicts an example of channel distortions undergone by a pilot signal.

[0106] The input-output response can be given by twisted convolution – like regular convolution, but with a twist – see section 5 below for details.

[0107] Embodiments will want the action of a doubly spread wireless channel on a pulse-tone to be predictable and geometric. Geometric may mean that for a channel path delay, the DD domain pulse is simply translated along the delay axis by an amount equal to the path delay; and it may be advantageous that for channel path Doppler shift, the DD domain pulse is simply translated along the Doppler axis by an amount equal to the Doppler shift.

[0108] Pulse-tone waveforms live on the hyperbola ^p^p= 1, and accordingly the action of the doubly spread channel on a pulse-tone is predictable if the fundamental period captures the delay and Doppler spreads of the channel.

[0109] The fundamental period captures the channel spreads when the delay domain period ^pis greater than the channel path delay spread and the Doppler domain period ^pis greater than the 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 channel path Doppler spread. When these conditions are met, we are in the crystalline regime and the action of a doubly spread channel on a pulse-tone is predictable. In some embodiments, such as FIG.6, this can allow the predictability of LTV channels, where wireless channel dynamics can lead to replicas.

[0110] In the example of FIG.6, the input-output response can be given by the twistedconvolution: ^^(^^,^^)=ℎ(^^, ^^) ∗^^ ^^(^^, ^^), where ℎ^^^, ^^^ ൌ ∑ସ ^ୀ^ ℎ^^^^^^ െ ^^^^^^^^^ െ ^^^^.

[0111] system including aof FIG.7 supports signal processing in Zak-OTFS.

[0112] The continuous output is the twisted convolution of the input with the effective DD channel filter. Quasi-periodic signals are shown in gray scale.

[0113] The discrete output is obtained by sampling the continuous output on the information grid.

[0114] The discrete output is the discrete twisted convolution of the discrete input with the discrete effective DD channel filter.

[0115] For the example system of FIG.7, the Zak-OTFS I / O relation can be represented as: ^^௪^^ ^^^ = ^^^ ∗ఙ ℎ^^^, ^^^ ∗ఙ ^^௧௫^^^, ^^^ ∗ఙ ^^ௗௗ^^^, ^^^ ൌ ℎௗௗ^^^, ^^^ ∗ఙ ^^ௗௗ^^^, ^^^.is a pictorial depiction of effects of a doubly spread channel (e.g. doubly spread pulse-tones)interaction of a doubly spread channel with a TD pulse-tone can be predictable and geometric.

[0118] In the crystalline regime, the delay domain period ^pis greater than the channel path delay spread, and the Doppler domain period ^pis greater than the path Doppler spread: ^p> delay spread and ^p> Doppler spread.

[0119] When communicating with pulse-tone (e.g., PulsoneTM) in the crystalline regime, the input-output relations can be governed by twisted convolution, in which the I / O relation for the entire subframe can be learned directly, without learning the channel, by reading off the response to a single transmitted pilot.

[0120] Unpredictability results from aliasing in the Delay-Doppler domain. Aliasing occurs when the channel spreads are bigger than the Delay-Doppler periods of the pulse-tone.

[0121] FIG.9 is a graphical example of predictability achieved in a crystalline region / regime of wireless communications.

[0122] In this example, a doubly spread channel is provided and M=N=4 is fixed to have 16 bins.

[0123] The average power of the received discrete DD domain signal for varying ^p, ^pis plotted on the hyperbola ^p^p= 1. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0124] As the delay domain period shrinks, aliasing in delay and frequency selectivity in FDM results.

[0125] As the Doppler domain period shrinks, aliasing in Doppler and time selectivity in TDM results.

[0126] There is a sweet spot in the middle of the hyperbola. In the sweet spot, the fundamental period can capture the channel spread. When the fundamental period captures the channel spread, there is no fading. Here, the received power profile is flat like the surface of a crystalline solid and this is called the crystalline regime.

[0127] In the crystalline regime, there is no fading and the input-output (I / O) relation is predictable.

[0128] As described above, DD domain aliasing controls predictability – as the delay domain period shrinks, aliasing in delay and frequency selectivity in FDM results; as the Doppler domain period shrinks, aliasing in Doppler and time selectivity in TDM results; and in the crystalline regime, there is no fading and the input-output (I / O) relation is predictable.

[0129] FIG.10 is a graph showing example results for model-free and model-dependent implementations (e.g., Model Free vs Model-Dependent). Here BER performance in the crystalline regime is explored where the I / O relation is predictable and communications can take place effectively even when it is not possible to estimate the channel accurately. In various embodiments, pulse-tone waveforms can support model-free operation in the crystalline regime when it is not possible to learn the channel. Although not shown in the BER performance example results, additional advantages such as improvements in filtering, root raised cosine versus sinc, and extending the region of reliable operation are also achievable in the crystalline regime.

[0130] Channels are becoming infinitely complicated, and it is becoming more and more challenging to learn them in the time available. Also, importantly, the effect of the channel on the modulation may matter more than the channel itself. This effect can be captured in the I / O relation, and it is predictable.

[0131] Operating model-free can comprise estimating the I / O relation everywhere in the Zak- OTFS subframe from the response to a single pilot waveform.

[0132] Model-free performance can be very close to performance with perfect channel state information CSI.

[0133] FIG.11A is an example flowchart of different implementations of an integrated sensing and communication with point pulse-tones (SCPP) [e.g., integrated sensing and communication with point pulse-tones].

[0134] The present document considers the Veh-A channel model comprising 6 channel paths first because it is representative of real propagation environments, second because it is very difficult to make the model-dependent mode of operation work. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0135] The model-dependent approach limits the number of paths and constrains delay and Doppler shifts.

[0136] The model-free approach was described in the previous section and it can be similar in spirit to estimating the taps of an effective LTI channel.

[0137] When the crystallization condition holds, the taps of the effective channel filter heff[k, l] can be read off from the response to a single point pilot pulse-tone.

[0138] This has been shown to work better than the model-dependent approach. However, in some embodiments, it may not be optimal and in this section the (maximum-likelihood) ML estimator is introduced.

[0139] Given the effective channel estimate, an estimate for the received pilot can be formed, the estimate can be subtracted from the received signal, and the data can be recovered.

[0140] It will be shown that when point pulse-tones are used for sensing and for data transmission then it may be necessary to divide DD domain resources between the two functions. A guard band may need to be introduced around the point pilot, and data is not transmitted within this guard band.

[0141] Also, an example flowchart of an implementation of an integrated sensing and communication with point and spread pulse-tones is provided in FIG.11B.

[0142] FIG.12 shows an example heat map of results obtained when using a pilot signal without a guard band.

[0143] A point pilot xpand a data signal xdare transmitted, and a signal y comprising a contribution ypfrom the pilot and a contribution ydfrom the data can be received.

[0144] The ML estimator can be the cross-ambiguity between yp and xp observed within the guard band which suppresses interference to yp from yd.

[0145] If a guard band is not included, then interference from data could completely obscure the effective channel taps preventing channel sensing.

[0146] In some embodiments, in the crystalline regime, the ML estimate for the tap ℎ^^^^^^, ^^^ canbe given by the Cross-ambiguity function ^^௬,௫^^^^, ^^^^^௬,௫^^^^, ^^^ ൌ ^^௬^,௫^^^^, ^^^ ^ ^^௬^,௫^^^^, ^^^

[0147] When there is no guard band, data interference could prevent channel sensing bycompletely obscuring the target term ^^௬^,௫^^^^, ^^^.

[0148] Hence, guard bands may be needed for examples such as this example for FIG.12.

[0149] In the example of FIG.12, ^^^௬,௫^^^^, ^^^^ is represented for Veh-A channel with ^^^^௫ ൌ 815Hz, Doppler period 30 KHz, ^^ ൌ 31,^^ ൌ 37, and RRC pulse shaping ^^^ఛ ൌ ^^ఔ ൌ 0.6^.

[0150] FIG.13 shows an example heat map of results obtainedusing a pilot signal with a guard band. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0151] Here, the guard band is a 7x7 rectangle centered at the pilot location.

[0152] It creates a little window around the pilot location where there is very little interference from the data signal yd.

[0153] Hence, data interference preventing channel sensing can be suppressed in the example of FIG.13. That is, guard bands can suppress data interference.

[0154] For the example of FIG.13, ^^௬^,௫^^^^, ^^^ represents cross-ambiguity between the receiveddata and the transmitted pilot.

[0155] In the example of FIG.13, ^^^௬^,௫^^^^, ^^^^ is represented for Veh-A channel with ^^^^௫ ൌ 815Hz, Doppler period 30 KHz, ^^ ൌ 31,^^ ൌ 37, and RRC pulse Shaping ^^^ఛ ൌ ^^ఔ ൌ 0.6^.

[0156] FIG.14 is a graph of example results obtained for an SCPP implementation.

[0157] The graph shows example results for Uncoded 4-QAM BER vs. ^^^^௫for different sensing and communication cases with a Point pulse-tone based on a Veh-A channel with Doppler period 30 KHz,^^ ൌ 31,^^ ൌ 37].

[0158] The graph provides example results for the following cases: ^^|^^ & ^^|^^ v^^. ^^|^^̅ & ^^|^^ ̅.

[0159] S stands for sensing, and C stands for data communications.

[0160] The notation S|C means that sensing takes place with interference from data, and the notation S|C bar means that sensing takes place without interference from data.

[0161] The BER for uncoded 4-QAM as a function of Doppler spread is plotted in FIG.14.

[0162] The case S|C bar & C|S bar is the baseline considered in the previous section, where separate Zak-OTFS subframes are dedicated to sensing and to data transmission. These are the two bottom curves with the gray circle and triangle markers. Note that, in this example, the RRC filter is better than the sinc filter at confining energy to the guard band.

[0163] The case S|C & C|S is integrated sensing and communications with point pulse-tones.

[0164] In this example, it is shown that when ^max> 1KHz, the channel response to a point pilot extends beyond the guard band regardless of the choice of filter. Here (when ^max> 1KHz), data interference can severely compromise BER performance because the channel response to a point pilot (e.g., pulse-tone) extends beyond the guard band regardless of the choice of filter.

[0165] FIG.15 graphically depicts examples of peak-to-average power ratio (PAPR).

[0166] An example TD realization of a point pulse-tone with Doppler period 30 KHz, ^^ ൌ 31,^^ ൌ37, and RRC pulse shaping ^^^ఛ ൌ ^^ఔ ൌ 0.6^ is depicted.

[0167] The examples of PAPR for different cases are shown in the provided Complementary CDF (CCDF) [complementary cumulative distribution function] of the Instantaneous to Average PowerRatio (IAPR) plot with RRC pulse shaping ^^^ఛ ൌ ^^ఔ ൌ 0.6^.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0168] Peak to Average Power Ratio (PAPR) can be an unattractive feature of the pulse-tone waveform. In some embodiments, it can be about 15dB and may require highly linear power amplifiers which are typically power inefficient. For high PAPR, highly linear power amplifiers are typically required.

[0169] 5. Zak-OTFS and Integration of Sensing and Communication (ISAC)

[0170] 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 is referred to as the crystallization condition. The filter taps can be read off from the response to a single Zak-OTFS point (impulse) pulse-tone waveform, and the I / O relation can be reconstructed for a sampled system that operates under finite duration and bandwidth constraints. Predictability opens up the possibility of a model-free mode of operation.

[0171] The time-domain realization of a Zak-OTFS point pulse-tone is a pulse train modulated by a tone, hence the name, pulse-tone, with PulsoneTMbeing an example implementation. The Peak-to- Average Power Ratio (PAPR) of a Pulse-tone is about 15 dB, and a general method for constructing a spread Pulse-tone for which the time-domain realization has a PAPR of about 5~6 dB is described. The spread Pulse-tone can be constructed by applying a type of discrete spreading filter to a Zak-OTFS point pulse-tone.

[0172] The self-ambiguity function of the point pulse-tone is supported on the period lattice Λp. By applying a discrete chirp filter, a spread pulse-tone with a self-ambiguity function that is supported on a rotated lattice can be obtained. If the channel satisfies the crystallization conditions with respect to the rotated lattice then the effective DD domain filter taps can be read off from the cross-ambiguity between the channel response to the spread pulse-tone and the transmitted spread pulse-tone.

[0173] If, in addition, the channel satisfies the crystallization conditions with respect to the period lattice Λp, then in an OTFS frame consisting of a spread pilot pulse-tone and point data pulse-tones, after cancelling the received signal corresponding to the spread pulse-tone, the channel response to any data pulse-tone can be recovered. This translates integration of communication and sensing within a single OTFS frame into geometric properties of a lattice Λp used for data transmission and a rotated lattice used for sensing. The spread pilot pulse-tone looks like noise to the point data pulse-tones, and it is this incoherence that makes it possible to integrate communications and sensing without time-sharing delay- Doppler resources. Integrated sensing and communication can increase effective throughput.

[0174] Pulses in the delay-Doppler domain can enable sensing, since they are geometric modes that are moved around by the Linear Time-Variant (LTV) channel.

[0175] Brief Introduction for Zak-OTFS and ISAC

[0176] In 6G propagation environments, as we encounter Doppler spreads measured in KHz, especially for non-terrestrial networks, it is becoming more and more difficult to estimate channels, and the standard channel-model-dependent approach to wireless communication is starting to break down. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 The present document describes a physical layer that is designed to integrate sensing and communications (ISAC) in challenging propagation environments, typical of 6G networks.

[0177] 4G and 5G wireless communication networks use OFDM. Here, a cyclic prefix is used to create shared eigenfunctions of the group of time shifts, and we have shared eigenfunctions because the time shift group is commutative. Note that pulses in the time domain serve as geometric modes that are moved around by the Linear Time-Invariant (LTI) channel. In the present document, we discuss and describe doubly-spread channels, where there are no shared eigenfunctions because delay shifts and Doppler shifts do not commute. Nevertheless, pulses in the delay-Doppler domain enable sensing, since they are geometric modes that are moved around by the Linear Time-Variant (LTV) channel. These pulses are the Orthogonal Time Frequency Space (OTFS) carrier waveforms, their time domain (TD) realizations are called pulse-tones as mentioned above.

[0178] A pulse in the delay-Doppler DD domain is a quasi-periodic localized function, defined by a delay period τpand a Doppler period νp. We have shown that the Zak-OTFS Input / Output (I / O) relation is predictable and non-fading when the delay period τpis greater than the effective channel delay spread, and the Doppler period νpis greater than the effective Doppler spread. We refer to this condition as the crystallization condition with respect to the period lattice and speak of operating in the crystalline regime. We have explained how non-predictability and fading result from aliasing in the DD domain, and why the crystallization condition prevents aliasing. In the crystalline regime we have shown that the effective DD domain channel filter taps can simply be read off from the response to a single Zak-OTFS point pulse- tone.

[0179] Why OTFS rather than OFDM? The first reason is that 6G propagation environments are changing the balance between time-frequency methods characteristic of OFDM and delay-Doppler methods. In OFDM, once the I / O relation is known, equalization is relatively simple, at least when there is no inter-carrier interference. However, acquisition of the I / O relation is non-trivial and model-dependent. In contrast, equalization is more involved in OTFS, due to intersymbol interference, but acquisition of the I / O relation can be simple and model-free (it can be read off from the response to a single point pulse- tone). Acquisition becomes more critical in 6G, as Doppler spreads measured in KHz make it more and more challenging to estimate channels.

[0180] The second reason is the integration of sensing and communication. The present document describes how to spread an OTFS point pulse-tone so that sensing is simple, and so that pilot waveform and data can coexist in the same OTFS subframe.

[0181] Sensing is introduced through a discussion of practical limitations on pilot design with a point pulse-tone. As mentioned in other sections, a pulse-tone is a pulse train modulated by a tone, and this is a waveform with high PAPR. Also, when we use the same lattice for communication and sensing we end up dedicating delay Doppler resources to either sensing or communication. For example, when we combine point data symbols with a point (impulse) pilot in a single OTFS subframe, we need to avoid 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 interference between data symbols and pilot. The standard approach is to introduce a guard band, but this is an overhead that reduces spectral efficiency, and the greater the channel spread, the greater is the overhead. This tradeoff motivates our use of different lattices for communications and sensing.

[0182] We also need to solve the issue of high PAPR which plagues multicarrier modulation. For example, OFDM signals are generated by adding many subcarrier components, and they can have high peak values in the time domain. As a consequence, OFDM systems suffer from high PAPR compared with single carrier systems. Many techniques have been proposed for reducing PAPR in OFDM but none avoid degrading BER performance. Over the past several years several variants of OTFS have been reported in the literature. A multicarrier approximation to Zak-OTFS, which we refer to as MC-OTFS, has been the focus of most research attention so far. MC-OTFS suffers from high PAPR, techniques for PAPR reduction have been presented, and again none avoid degrading BER performance.

[0183] Zak-OTFS modulation starts with a quasi-periodic Dirac-delta DD domain pulse, and we apply a DD domain pulse shaping filter. The time-domain (TD) representation of the filtered signal is a TD pulse-tone with time duration T and bandwidth B inversely proportional to the Doppler and delay spread respectively of the filter. We transmit information using non-overlapping DD domain pulses spaced 1 / T apart along the delay axis and 1 / B apart along the Doppler axis. Since each pulse repeats quasi- periodically, there are M = τp / (1 / B) = Bτppulse locations along the delay axis and N = νp / (1 / T ) = Tνppulse locations along the Doppler axis. The number of distinct non-overlapping information carriers is the time- bandwidth product BT = MN, and the pulse locations are the points in the information lattice.

[0184] A discrete quasi-periodic DD domain signal is periodic along both delay and Doppler axes with period MN (e.g., M*N). In a later section, we define a discrete DD domain filter to be a discrete periodic DD domain function with period MN along both delay and Doppler axes (M2N2degrees of freedom). A later section develops the fundamentals of filtering in the discrete DD domain. We describe how we obtain a spread pilot by applying a discrete spreading filter to a DD domain signal localized at a point in the information lattice. This is a general method. In the present document, however, we also discuss and describe using a discrete chirp filter which distributes energy equally to all lattice points. Different lattice points correspond to pulse-tones with different delays and Doppler shifts; hence we observe a TD signal with almost constant amplitude. As it will be further discussed, the PAPR of the spread pulse-tone is about 5 dB, significantly less than the PAPR of the point pulse-tone, which is about 15 dB.

[0185] A later section describes how we recover the effective channel filter by de-spreading the received signal. We can estimate the taps of the discrete effective channel filter from the discrete cross- ambiguity function between the received DD signal and the transmitted spread pilot signal. The self- ambiguity function of a point pulse-tone is a rectangular lattice Λz,p, and we show that the self-ambiguity function of the spread pulse-tone is obtained by rotating Λp to obtain e.g., a lattice Λ∗. If the channel 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 satisfies the crystallization conditions defined by Λp then we show that we are able to accurately estimate the discrete effective channel taps using a point pulse-tone. If the channel satisfies the crystallization conditions defined by e.g. Λ∗ then we show that we are able to accurately estimate the discrete effective channel taps using the spread pulse-tone. Embodiments described in the present document for spreading the pilot are connected to the dynamics of the sensing environment.

[0186] Integration of sensing and communication (ISAC) is central to 6G; hence we consider embodiments involving transmission of pilot and data within a single OTFS subframe with no division of DD domain resources between sensing and communication. Data can be carried by point pulse-tones, the data signal interferes with channel estimation using a spread pulse-tone, and this interference can adversely affects the estimation of certain discrete effective channel filter taps. We describe how the significance of this interference depends on the ratio of pilot power to data power. After equalization, the spread pulse-tone interferes with data carried by the point pulse-tones, and again the degree of interference depends on the ratio of pilot power to data power. The present document identifies embodiments on how to choose this ratio to accomplish integration of sensing and communication within a single OTFS subframe and avoid time-sharing of delay-Doppler resources.

[0187] Passive radar is an important instance of integrated sensing and communication. Passive radar exploits readily available, non-cooperative sources of radio energy to measure reflections from the environment and targets of interest. When a suitable illuminator is available, covert surveillance becomes possible, without the need for deployment and operation of a dedicated transmitter. For example, terrestrial digital television transmissions (DVB–T) provide an especially attractive opportunity for radar. We describe how to design the illuminator in a passive radar to be a geometric mode of the radar scene (a spread pulse-tone). We describe how this choice can simplify acquisition of the radar scene, since the scene can be simply read off from the received signal, unlike the DVB-T system. The radar application also benefits from the noise-like characteristics of the spread pulse-tone which increase the fraction of energy on target.

[0188] Pilot Design: In the present document, we describe constructing a spread pulse-tone, where the self-ambiguity function of the spread pulse-tone is supported on e.g. a lattice Λ∗ obtained by rotating the period lattice Λp associated with the data-bearing pulse-tones. We show that if the channel satisfies the crystallization conditions with respect to Λ∗ then the effective DD domain filter taps can simply be read off from the channel response. If the channel also satisfies the crystallization conditions with respect to the Λp then given the I / O response at one point in the OTFS subframe, it is possible to predict the I / O response at all other points in the subframe. We also show that the spread pulse-tone is a noise-like waveform with excellent PAPR.

[0189] Integrating Sensing and Communication (ISAC): We have translated integration of communication and sensing into geometric properties of a lattice Λp used for data transmission and a rotated lattice e.g., Λ∗ used for sensing. The data pulse-tones look like noise to the sensing pulse-tones, 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 and we have demonstrated that this incoherence makes it possible to integrate sensing and communication without the loss in spectral efficiency associated with time-sharing delay-Doppler resources. We have also demonstrated other advantages that ISAC can offer such as increase in effective throughput by avoiding time-sharing delay-Doppler resources. For example, an important application is the design of a passive radar system.

[0190] System Model

[0191] A time-domain (TD) pulse is an ideal waveform for delay-only channels (where paths induce zero Doppler shift) since it is possible to separate signals received along different paths according to their path length / distance. Similarly, a frequency-domain (FD) pulse is an ideal waveform for Doppler- only channels since it is possible to separate signals received along different paths according to the Doppler shift induced on the transmitted signal. However, neither a TD pulse nor a FD pulse is suited for doubly-spread channels where paths induce both delay and Doppler shift.

[0192] Here, it is described how a pulse in the delay-Doppler (DD) domain is matched to doubly- spread channels. A pulse in the DD domain is a quasi-periodic localized function, defined by a delayperiod τp and Doppler period νp = 1 / τp. In the period lattice ^p ^^ ^n^ p, m ^ p ^ n , m ^^ ^ , there is onlyone pulse within the fundamental region ^0 ^^ ^^ ,v ^ 0 ^ ^are infinitelymany replicas along the delay and Doppler(1)for all n ,in the time domain, this function can be realized as a pulse train modulated by a tone (see e.g. FIG.17), hence the name Pulsone™. The DD domain pulse is the carrier waveform for Zak-OTFS modulation.

[0193] 5.1.A. Zak-OTFS Modulation

[0194] Zak-OTFS modulation with modulation parameters (τp, νp), νp= 1 / τpis considered here. The transmitted TD Zak-OTFS frame is limited to a time duration T=Nτpand bandwidth B=Mνp. FIG. 16 illustrates Zak-OTFS transceiver processing.

[0195] Let x ^ k, l ^ , k ^ 0,1,^ , M –1, l ^ 0,1, ^ , N –1 , denote the BT = MN information, each having unit average energy, i.e., ^ 2symbols ^ k, l ^ ^^ ^ ^ 1. These MN information symbols areencoded into a discrete DD domain informationgiven by (2)180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00for all k ^ 0,1,^ , M –1, l ^ 0,1, ^ , N –1, n , m ^ ^ . From this it follows that(3)for all k , l , n ,alongthe delay axis and N along the Doppler axis. The encoding of the MN information symbols to xdd[k, l] is carried out as in (2), so that xdd[k, l] is quasi-periodic, since only quasi-periodic DD functions have TD realizations.

[0196] The information signal xdd[k, l] can then be converted to a continuous quasi-periodic DDsignal by lifting it to the information lattice ^dd ^^ ^k^p M , l ^ p N ^ k , l ^^ ^ , i.e.. (4)

[0197] along the delay and Doppler axis, i.e., for all n , m^^. (5)Dirac- a a a domain.

[0199] Pulse shape filtering in the DD domain can be implemented by twisted convolution. Twisted convolution of xdd(τ, ν) with a pulse shaping filter wtx(τ, ν) gives . (6)

[0200] dd ,Also, appropriate pulse shaping guarantees that the transmit TDsignal has time duration T and bandwidth B.

[0201] Note, in the present document, twisted convolution is denoted by ∗σ. Twisted convolution between two DD functions a(τ, ν) and b(τ, ν) can be given byc ^^, ^ ^^ a ^ ^ , ^ ^ ^^ b ^ ^ , ^ ^ ^^ ^ a ^ ^^ , ^ ^ ^ b ^ ^ ^ ^ ^ , ^ ^ ^ ^ ^ e j2^^^ ^ ^^ ^ ^ ^ d ^ ^ d ^ ^ . Twisted convolutionit is not commutative, i.e., a(τ, ν) ∗σb(τ, ν) ≠ b(τ, ν) ∗σa(τ, ν). 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[00202] The inverse Zak-transform of xwtxdd ^^ , v ^ gives the transmit TD Zak-OTFS modulated signal(7)

[0203] Equivalently,(8) where ζk,l(t) is the carrierof the Zak-OTFS carrier waveform. In the DD domain it is a quasi-periodic pulse. In both the TD / FD domains, the carrier waveform is a pulse train modulated by a tone and is therefore called a TD / FD pulse- tone.

[0204] In FIG. 17, the (k, l)-th quasi-periodic DD domain pulse located at ^^0,v 0 ^ ^ k ^^ pM , lv pN ^ and its TD / FD realizations referred to as TD / FD pulse-tone. The TD pulse-tonepulse train modulated by a TD tone. The FD pulse-tone comprises of a finitebandwidth pulse train modulated by a FD tone. The location of the pulses in the TD / FD pulse train and the frequency of the modulated TD / FD tone is determined by the location of the DD domain pulse (τ0, ν0). The time duration and bandwidth of a pulse-tone are inversely proportional to the characteristic width of the DD domain pulse along the Doppler axis and the delay axis, respectively. As τp→ ∞, the TD pulse-tone approaches a single TD pulse which is suited for delay-only channels. Similarly, as νp→ ∞, the FD pulse- tone approaches a single FD pulse which is suited for carrying information in Doppler-only channels. Zak- OTFS is therefore a family of modulations parameterized by τpthat interpolates between TD pulse modulation (TDM) and FD pulse modulation (FDM).

[0205] 5.1.B. Zak-OTFS Receiver

[0206] The received TD signal can be given by (9)DD representation ydd(τ, ν) of rtd(t) is given by ydd ^^, v ^ ^ ^ t ^ r td ^ t ^ ^ . Matched filtering with the receive DDpulse wrx(τ, ν) results in the signal 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 . (10)

[0207] periodicity. Sampling information lattice Λdd gives(11)for all k , l^^ .

[0208] 5.1.C. Zak-OTFS I / O Relation

[0209] Combining all equations above, from (4) to (11), gives the I / O relation , (12) where the

[0210] Inconvolution between two discrete DD domain functions heff[k, l] and xdd[k, l] which is explicitly given by

[0211] x ^ k, l ^ , k ^ 0,1,^ , M –1, l ^ 0,1, ^ , N –1 may be detected using the symbolsydd ^k , l ^ , k ^ 0,1,^ , M –1, l ^ 0,1, ^ , N –1. This yields the matrix-vector form of the I / O relation^ 1 vector of these MN received symbols is given by the product of an MN ^ MN effective channel matrix with the MN ^ 1 vector of information symbols. It is possible to detect transmitted Zak- OTFS information symbols using techniques developed for MIMO equalization since the I / O relations share the same form.

[0212] In the present document, for DD domain pulse shaping only sinc and root raised cosine (RRC) pulses are considered which are respectively given by 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 , (15) and, (16) where for 0 ≤ β ≤ 1. (17)

[0213] the main lobe and expansion in time / bandwidth. The receiver pulse shaping filter for both sinc and RRC waveforms can be given by . (18)

[0214] 5.2

[0215] The Zak-OTFS I / O relation is predictable and non-fading when the delay period τpis greater than the effective channel delay spread (the spread of heff(τ, ν) along the delay axis) and the Doppler period is greater than the effective channel Doppler spread spread of heff(τ,ν) along the Doppler axis). Thisis the crystallization condition. In the crystalline (when the crystallization condition holds), it is possible to read off the taps of the effective DD domain channel filter heff[k,l] from the response to a single Zak-OTFS pulse-tone. This pilot waveform can be referred to as a point pilot since it's DD domain realization is a localized quasi-periodic pulse. While simple and effective, this is not the maximum likelihood (ML) estimate. In this section, the ML method for estimating the taps of heff[k, l] in terms of ambiguity functions is described.

[0216] To begin, the model-dependent and model-free approaches to acquiring the I / O relation in the crystalline regime are reviewed. In the model-dependent approach, a model is imposed on the DD spreading function hphy(τ, ν), typically by prescribing a finite number of paths and constraining their delay and Doppler shifts. Given this model, the receiver estimates hphy(τ,ν). The accuracy of this estimate is limited by the time and bandwidth constraints on the pilot signal and by any mismatch between the channel model and the physical channel.

[0217] In the model-free approach, the receiver can estimate the taps of the effective DD domain channel filter heff[k,l] from the response to a pilot waveform. As mentioned above, this approach is similar in spirit to estimating the taps of an effective discrete-time LTI channel instead of estimating the continuous- time impulse response of the underlying physical channel. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0218] Here and throughout this present document, a Veh-A channel model which consists of six channel paths is considered. The delay-Doppler spreading function can be given by . (19) where hi, τi, and νiIlists an example power- In anis^max ^ max i ^ i ^ min i ^ i ^ 2.5 ^ s . For a maximum Doppler shift of νmax = 815 Hz, the Doppler spreadis 2 ^ 0.815 = Doppler shift of the i-th path can be modeled as νi=νmaxcos(θi), where thevariables^i , i ^ are independent and distributed uniformly in the interval [0, 2π).TABLE I: Power Delay Profile of Doubly-spread Veh-A Channel. Path no. i1 2 3 4 5 6Rel Dela τ ( s)0 031 071 109 173 251

[0219] eal propagationenvironments, and second because it can be very difficult to make the model-dependent mode of operation work. As an example, for a channel bandwidth B = 0.96 MHz, the delay domain resolution is 1 / B ^ 1.04µs, and the first three paths introduce delay shifts in the interval [0, 0.71]µs which is less than the delay domain resolution. These paths are therefore not separable, and so cannot be estimated accurately. One can choose a higher bandwidth so that the resolution is smaller than the path difference between any two paths thereby making the paths separable. However, this would only work for a particular path delay profile. In real scenarios, the path delays can change, however bandwidth cannot be increased indefinitely.

[0220] By contrast, model-free operation is always feasible irrespective of whether the paths are separable or not. FIG. 10 compares BER performance of model-dependent and model-free modes of operation. We consider the Veh-A channel with νp = 15 KHz, B = 0.96 MHz, T = 1.6 ms. Separate Zak- OTFS subframes can be dedicated to channel sensing and data transmission in order to focus on the difference between model-dependent and model-free modes of operation. It is emphasized that this is NOT integrated sensing and communication. In the model-free mode, the taps of heff[k,l] are estimated directly from the received DD domain symbols received in response to a point pulse-tone situated in the middle of^ 0. In the model-dependent mode, the receiver first estimates the parameters (τi, νi, hi) of the underlyingphysical channel, which is then used to estimate the taps of heff[k,l]. MMSE equalization of the matrix-vector form of the Zak-OTFS I / O relation can be employed to detect information symbols at the receiver. It is 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 observed that BER performance of the model-dependent mode exhibits a high error floor. This is due to inaccurate estimates of the parameters of the underlying physical channel, which is a consequence of insufficient subframe bandwidth and duration. On the other hand, BER performance of the model-free mode of Zak-OTFS can be considerably better, only slightly worse than performance with perfect knowledge of the Zak-OTFS I / O relation. For example, in the crystalline regime, BER performance of the model-free mode of Zak-OTFS is only slightly worse than performance with perfect knowledge of the I / O relation.

[0221] FIG.10 provides BER performance of Zak-OTFS in model-dependent and model-free mode of operation. In the crystalline regime (νp = 15 KHz), BER performance of the model-free mode of Zak- OTFS is only slightly worse than performance with perfect knowledge of the I / O relation.

[0222] Herein, this present document focuses on the model-free mode of operation in the crystalline regime. It is now described how to integrate sensing and communication within a single Zak-OTFS subframe. The discrete DD domain point pulse-tone indexed by (kp, lp) can be given by . (20)

[0223] ofdata, the received response can be given by

[0224] The support set ^ of heff[k, l] can consist of all pairs (k, l) for which heff[k, l] ≠ 0. Thecrystallization condition is satisfied when the delay spread of ^ is less than the delay period M and theDoppler spread is less than the Doppler period N. In the crystalline regime, the taps of heff[k, l] can be read off from (n, m) = (0, 0) term in (21). A previously presented method was simply to read off the point pulse- ^ j 2^ kp ^ l ^ lp ^MN tone response within the support set of heff[k–kp, l–lp] and multiply . This method can give a much better estimate than the model-dependent approach,estimation method. The ML estimator can be given by the samples of the cross-ambiguity between the received point pulse- 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 tone yp,dd[k, l] and the point pulse-tone xp,dd[k,l], when observed only within the support set ^ of heff[k, l].Given ^ k, l ^ ^^ we havea[k, l] and b[k, l] can be defined as in section 5.6 (under The Discrete Ambiguity Function). It follows from Theorem 6 in section 5.6 (under The Discrete Ambiguity Function) that the noise-free cross-ambiguity Ayp, x p ^ k , l ^ is given byp Self-ambiguity Function of the Point Pilot Signal) it follows that the self-ambiguity function of a point pulse- tone can be given by. (25). RHS do not overlap, and the (0, 0) term is simply equal to Ayp, x p ^ k , l ^ for ^ k, l ^ ^^ .180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0227] FIG.18 shows points on the information lattice Λdd, the period lattice Λpand the dual of the information lattice^ ^dd . If the delay and Doppler coordinates of these lattices are scaled down by τp / Mand νp / N then Λddmay be identified with^ 2 and ^ ^dd with MN^^ MN ^ . Now considerany shift (nMN, mMN) of a quasi-periodic signal xdd[k, l] in DD domain. It follows from (3) that

[0228] lattice^ ^dd are shown by black, light gray and dark gray dots respectively. Adjacent points on theinformation lattice are separated by ∆τ = τp / M and ∆ν = νp / N along the delay and Doppler axis respectively. Adjacent points on the dual of the information lattice are separated by MN∆τ=Nτpand MN∆ν=Mνpalong the delay and Doppler axis respectively.

[0229] Signals that are quasi-periodic with respect to the information lattice are periodic with respect to the dual of the information lattice. In section 5.3, filters in the discrete DD domain are introduced where this interplay will be important.

[0230]

[0042] We conclude this Section by discussing two practical issues that limit the use of point sensing pulse-tones.

[0231] 5.2.A High Peak-to-Average-Power-Ratio (PAPR)

[0232] FIG 19A depicts a TD realization of a point pulse-tone with Doppler period νp= 30 KHz, M = 31, N = 37. Point pulse-tone is located at (kp, lp) = ((M + 1) / 2, (N + 1) / 2). Only part of the entire TD pulse-tone is shown here, with samples taken every14Bseconds where B = Mνp= 930 KHz. RRC pulse shaping with βτ= βν= 0.6.

[0233] FIG.19A displays the magnitude of the TD realization of a point sensing pulse-tone. We observe a train of narrow pulses, exhibiting sharp peaks at the pulse locations, leading to high PAPR. This requires the use of highly linear power amplifiers which are typically power-inefficient.

[0234] 5.2.B Interference Between Data and Sensing pulse-tones

[0235] Data and sensing pulse-tones can combine at the transmitter to produce the noise-free DD domain signal given by 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 . (28)

[0236] DD domain filter (29)the received pilot pulse-tone signal results from the action of a second discrete DD domain filter (30) on the same quasi-tones will interfere at the receiver. It is emphasized that it is only possible to express this interference in terms of operators because signal processing takes place in the DD domain where twisted convolution is associative.

[0237] An equivalent expression for the DD domain information signal in (2) is that the data symbol x^ k^, l ^ ^ is carried by a discrete quasi-periodic DD domain pulse at ^k^, l ^ ^ so that[k, l]shifted by ^k^, l ^ ^ , i.e.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0239] that

[0240] , (34)data transmission. The ML estimate for the taps of heff[k, l] is simply the cross-ambiguity function of the received pilot with the transmitted pilot (see (129)), as shown in FIG.19B. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0242] FIG.19B graphically depicts magnitude of the cross-ambiguity between the transmitted point pilot signal and the received point pilot signal, i.e., heff ^ k, l ^^^ Axp , x p ^ k , l ^ for EVA channel with νmax =815 Hz. M = 31, N = 37, νp= 30 KHz; RRC lp) = ((M + 1) / 2, (N + 1) / 2). The support set of heff[k, l] is approximately the

[0243] that we integrate channel sensing and data transmission within the same Zak-OTFS subframe the received discrete noise-free DD signal is given by (35). The ML estimates for the taps of heff[k, l] are now given by the cross-ambiguity functiondata. FIG.19C illustrates how the first term, which represents interference from the data signal, can obscure the second term. Within the first term, Ax 0,x p ^ k,l ^ is the cross-ambiguity function between the quasi-periodicpulse at the origin and the quasi-periodic pilot signal.

[0245] FIG.19C graphically depicts the magnitude of the cross-ambiguity between the transmitted point pilot signal and the received signal which includes contributions from pilot and data ( Ay, xp ^ k , l ^ ).In FIG.19C: EVA channel with νmax= 815 Hz. M = 31, N = 37, νp= 30 KHz; RRC pulse= 0.6; (kp, lp) = ((M + 1) / 2, (N + 1) / 2). In the example of FIG.19C, interference prevents channel sensing.

[0246] The cross-ambiguity functionsAx 0 , x p and Axp , x p are both supported on the period lattice.This leads to interference between data and pilot signals at the receiver which prevents channel sensing by completely obscuring the second term. FIG.19D illustrates how the second term can be made visible by introducing a guard band ^ around the pilot signal, in this case a 7^7 rectangle. Guard bands cannotbe avoided when the cross-ambiguity functions Ax 0,x p ^ k,l ^ and Axp, x p ^ k , l ^ are both supported on thesame lattice.

[0247] In the example of FIG.19D, data / information symbols are not transmitted in a guard region ^(e.g., a 7^7 rectangle around the location (kp, lp) = ((M + 1) / 2, (N + 1) / 2) of the pilot signal) which is the180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 location of the point sensing signal. FIG. 19D graphically depicts the magnitude of the cross-ambiguity between the transmitted pilot signal and the data component of the received signal, i.e., Ay, xp ^ k, l ^^ E p h eff ^ k , l ^ ^^ A x p , x p ^ k , l ^ (see (36)). In FIG.19D: Veh-A channel with νmax = 815 Hz.βτ= βν= 0.6. In the example of FIG.19D, interference

[0248] In section 5.4, it will be described how to filter the point pilot signal so that Axp, x p ^ k , l ^ canbe supported on a lattice different from the period lattice. It will be shown that it is possible to choose thefilter so that the cross-ambiguity function Ax 0,x p ^ k,l ^ appears as noise to the self-ambiguity functionAxp, x p ^ k , l ^ of the filtered pilot signal. We seek a filter that results in cross-ambiguity functions that areor incoherent. This can eliminate the need for a guard band at the cost of reducing the SNR of sensing. In addition, the energy of the filtered pilot signal will be spread uniformly in the discrete DD domain, leading to a TD realization with excellent PAPR.

[0249] Finally, we define the data signal power to noise power ratio, the pilot signal power to noise power ratio, and the pilot signal power to data signal power ratio.

[0250] From (31), it follows that the average energy of the data signal is Ed, i.e. . (37)

[0251] Since thelp), from (33) it follows that the energy of the point pilot signal is Ep, i.e. . (38)

[0252] AWGN ntdThe average total energy of the discrete DD domain noise signal ndd[k, l] (see (12)) can be given by , (39)180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 2

[0253] The effective channel gain is h k,land therefore the received energy ^ ^ eff^ ^of k ,l ^^ ^ 2data / information pulse-tones and that of the pilot pulse-tone are Edh k,lan ^ ^ eff^ ^d k , l^^ ^E2 p ^heff^ k,l ^respectively. Hence the ratio of the power of the receivedto the ^ k , l^^ ^ noise SNR) can be given by and the ratio of the power of theby

[0254] The ratio ofpilot ^ pulse-tones can therefore be given byd^ p which is subsequently referred to as the pilot power to data power ratio (PDR).

[0255] 5.3. Filtering in the Discrete Delay Doppler Domain

[0256] Here, we describe how to design a filter in the discrete DD domain to spread a point pilot signal over the MN pulse-tones (M*N pulse-tones) located on the information grid, so that each pulse-tone contributes a fraction 1 / (MN) of the spread pilot energy. The discrete DD domain filter ws[k, l] can act on the point pilot xp,dd[k, l] by twisted convolution to produce a spread pilot signal xs,dd[k,l] given by . (42) where xp,dd[k, l] is. quasi-periodicity.

[0257] More generally, we develop the fundamentals of filtering in the discrete DD domain. We start from the theory of linear time invariant (LTI) systems. Here it is well known that linear convolution of a discrete-time periodic signal with a discrete-time filter is equivalent to periodic linear convolution of the periodic signal with the periodic extension of the filter. Thus, two discrete-time filters with identical periodic 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 extensions are equivalent. If the period is L, then the discrete-time filter is specified by a vector in^ L . Inthe discrete DD domain, quasi-periodic signals are periodic with period MN along both delay and Doppler axes (see (27)). We now show that twisted convolution of a quasi-periodic discrete DD signal with a discrete filter ws[k, l] is equivalent to a MN-periodic twisted convolution [i.e., M*N-periodic twisted convolution] of that quasi-periodic signal with the MN-periodic extension of the discrete filter. Thus, two discrete DD domain filters with identical periodic extensions are equivalent. Since the period is MN along both delay M2 N 2and Doppler axes, a discrete DD domain filter can be specified by a vector in^ .

[0258] Given an arbitrary discrete DD domain filter as[k, l], we can define its MN-periodic extension a[k, l] by . (43)

[0259] a[k, l] acts on a quasi-periodic signal b[k, l] by MN-periodic twisted convolution, denoted by ⊛σ. More preciselyperiodic signal b[k, l] coincides with the MN-periodic twisted convolution of its MN-periodic extension a[k, l] with b[k, l], i.e. . (45) Proof: We startMN, then simplify using the definition of an MN-periodic extension (43). The detailed derivation of (45) is given by (46). 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00where(48) is the DD domainthe spread DD pilot signal expression given by (49)., (50)180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00for k , l^^ , and we refer to q^^ as the slope parameter. In the next section, we will show that theconstant factor 1 / (MN) in (50) results in a spread pilot signal xs,dd[k, l] with unit energy. We will observethat xs,dd ^ k,l ^ is almost constant (say xs,dd ^ k,l ^ ^ ^ for all k, l) when M and N are odd primes, andwhen q prime to both Since the total energy of the spread pilot is M^ 1 N1^^^ xs,dd^ k,l 2^ ^ MN^ 2 , we conclude ^ ^ 1 MN . Hence the peak amplitude of the pointk ^ 0 l ^ 0pilot signal (1 at location (k, l) = (kp, lp)) is about MN times higher than that of the spread pilot signal.FIG.19E compares the energy profiles of point and spread pilot signals.

[0263] In the example of FIG.19E, energy profiles of DD domain pilot signals are provided. M = 31, N = 37 and the total DD domain energy is normalized to 1 for both point pilot and spread pilot signals, that M^ 1 N ^ 1 M ^ 1 N ^ 1^^ x s,dd ^ k,l 2^ ^ ^^ x p,dd ^ k 2is, ,l ^ ^1. The point pilot takes the value 0 at all locations exceptmagnitude of each spread pilot symbol is about 1 / (MN), i.e., about 10log10(MN) = 30.6 dB below the peak squared-magnitude of a point pilot DD signal.

[0264] 5.4 Zak-OTFS with Spread Sensing pulse-tone

[0265] 5.4.A Reducing PAPR

[0266] In Section 5.3, we constructed a spread pulse-tone using a discrete chirp filter, and we compared the energy profiles of point and spread pulse-tones in the discrete DD domain. Recall that the TD realization of the point pulse-tone located at (kp, lp) consists of narrow TD pulses at time instancest ^ ^ n^p ^ k ^pp M ^ , n ^^ where the spread of each pulse is approximately τp / M, which is the inverse. FIG.19F shows the TD realization of the spread pulse-tone constructed in Section 5.3. This TD realization is the sum of all MN pulse-tones located on the information grid and the constant term 1 / (MN) appearing in (50) reduces the amplitude of each point pulse-tone in the sum. The N pulse-tones located at a given kpinterfere destructively, and the point pulse-tones at different locations kpresult in trains of TD pulses at different locations. This explains why the TD realization of the spread pulse- tone is less peaky than the TD realization of the point pulse-tone shown in FIG.19A.

[0267] FIG. 19F graphically depicts a TD realization std ^ t ^ k ^4 B ^ ^ , k ^^ of a spread pilotsignal with Doppler period ν = 30 KHz, M =p 31, N = 37. chirp filter with q = 3. Only part of the entire TD realization is shown here, with samples taken every14Bseconds where B = Mνp= 930 KHz.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 RRC pulse wtx(τ, ν) with roll-off factors, βτ= βν= 0.6. Again, observe that the TD realization is not as peaky as that of the point pulse-tone shown in FIG.19A.

[0268] provides a complementary CDF (CCDF) plot of IAPR, where M = 31, N = 37; RRC pulse wtx(τ, ν) with roll-off factors, βτ= βν= 0.6; Discrete chirp filter with q = 3; and PDR ρp / ρd= 10 dB.

[0269] FIG. 19G displays the CDF of the instantaneous-to-average-power ratio (IAPR) for example spread and point pilots. In the absence of data, the IAPR of the spread pulse-tone does not exceed 5 dB, whereas the IAPR of the point pilot is almost 15 dB. With data alone, the IAPR is similar to that of Gaussian noise. In Section 5.5 we choose a pilot to data power ratio (PDR) of 10 dB to integrate sensing and communication within a single Zak-OTFS subframe. For both spread pilot with data and point pilot with data, the IAPR rarely exceeds 7 dB. However, for a higher PDR of 25 dB, while the IAPR of spread pilot with data rarely exceeds 9 dB, the IAPR of point pilot with data can be as high as 12 dB. In various embodiments, spreading can reduce the PAPR of the transmitted signal.

[0270] 5.4.B Achieving Predictability

[0271] Here, we show that the self-ambiguity function of the spread pulse-tone is supported on a lattice Λqthat is obtained by applying a linear transformation to the period lattice Λpthat depends on the slope q of the discrete chirp. When the effective channel filter satisfies the crystallization condition with respect to the rotated lattice Λq, we will describe how to estimate the taps of heff[k, l] from the cross- ambiguity function of the received pilot and the transmitted pilot. Here, on sensing in theabsence of data and noise, but in section 5.5 we will describe methods for in the presence of data and noise.

[0272] The spread pilot signal xs,dd[k, l] = w[k, l] ⊛σxp,dd[k, l] is given by (47). Here w[k, l] is the spreading filter and xp,dd[k, l] is the point pilot signal. In the absence of data, the received pilot signal can be given by . (51)

[0273] ,cross-spread pilot and the transmitted spread pilot. It follows from (113) in section 5.6 (under The Discrete Ambiguity Function) that 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00(in of data and noise), (53)where Axs , x s ^ k , l ^

[0274] Henceforth, we restrict our attention to the spread pulse-tone (54)arising from the discrete chirp filter,. Since ambiguity functions are periodicalong both delay and Doppler axes with^ k ^ p 1 MN, l ^ p 2 MN ^ ^ A x s , x s ^ k , l ^ ,for all p 1, p 2 ^^ , it suffices to characterize the support setx s ^ k , l ^ modulo MN along both axes.Theorem 2 shows that the self-ambiguity function Axs , x s ^ k , l ^ is supported on a lattice.

[0275] Theorem 2: Given odd primes N, M and q relatively prime to both M and N

[0276] Note, for any integer a and positive integer M, [a]M denotes the unique smallest non- negative integer which is congruent to a modulo M.

[0277] Theorem 2 implies that the self-ambiguity function Axs, x ^ k , l ^ is non-zero if and only ifthere exist integers n, m such that180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 . (56)

[0278] We rewrite (56) as .

[0279] It follows thatby. (58)

[0280] Note, in

[0281] We show that the lattices Λp, Λqhave the same fundamental volume by showing that they have the same density. We show the densities are the same by showing that the number of Λqpoints in the rectangle Ω bounded by (0, 0), (MN, 0), (0, MN) and (MN, MN) is the same as the number of Λppoints. ^1

[0282] Since ^ ^[ ^ 2q ^ ^ 2 q ] 2MN , ^ 1 ^ 2 q ^ ^ ^ 4 q modulo MN. Since M and N are oddprimes and q is relatively prime to both M and N, 4q2is relatively prime to MN. Therefore, (1 − 2qθ) isrelatively prime to MN, i.e., ^1^2q^ ^ ^^ 0 modulo MN and hence, the linear transformation ^ ^1^^ ^is ^ non-singular modulo MN. Hence the number of Λ points in Ω equals the number of Λp

[0283] From (113) in section 5.6 (under The Discrete Ambiguity Function) we know that the total energy of the spread pulse-tone is simply the value of its cross-ambiguity at the origin, i.e., . (59)s, x s ,. (60) where (k, l) areiq. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00[k, l] itself. The other terms in the R.H.S. correspond to the other lattice points not at the origin. Therefore, accurate estimation of the taps of heff[k, l] may not be possible if for any non-zero lattice point (ki, li) ≠ (0, l0) the support set of the corresponding term ej^ ^^li ^ k ii h eff ^ k ^ k i,l ^ l i ^ ej2^MN overlaps with thesupport set of heff[k, l]. We now show that if condition with respectto Λq, then we can accurately estimate the taps of heff[k, l].

[0286] Example 1: The dark gray dots in FIG.19H mark the points of the lattice Λq, i.e., support of the self-ambiguity function of the spread pulse-tone with M = 11, N = 13, q = 5, (kp, lp) = (0, 0). The fundamental period / region of Λqis a parallelogram with area 143 = 11 ^ 13 = MN bounded by the points (0, 0), (3, 19), (11, 22) and (8, 3). The lattice Λqis generated by (3, 19) and (8, 3). The cross-ambiguity between the received spread pulse-tone and the transmitted spread pulse-tone is supported on the union of the gray rectangles (these rectangles correspond to the terms in the RHS of (61)). Since the rectangles do not overlap, there is no DD domain aliasing (i.e., the channel satisfies the crystallization condition with respect to Λq), and we can accurately estimate heff[k, l] from the response received within the gray rectangle with the black border. In other words, the support of heff[k, l] is limited to a rectangular region with delay spread (kmax – kmin) and Doppler spread (lmax – lmin) shown as the gray rectangle with the black border.

[0287] That is, in FIG.19H, the dark gray dots mark the support of the self-ambiguity function of the spread pulse-tone with M = 11, N = 13, q = 5, (kp, lp) = (0, 0). The support of heff[k, l] is limited to a rectangular region with delay spread (kmax– kmin) and Doppler spread (lmax– lmin) shown as the gray rectangle with the black border.

[0288] Example 2: The dark gray dots in FIG.19I mark the support of the self-ambiguity function of the spread pulse-tone with M = 11, N = 13, q = 4, (kp, lp) = (0, 0). The image of the fundamental period / region is a parallelogram bounded by (0, 0), (24, 5), (5, 7) and (29, 12). The support 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 of heff[k, l] is the same as that in FIG.19H and is again shown as the gray rectangle with the black border. This rectangle overlaps the gray rectangles located at (5, 7) and (−5, −7) and taps of heff[k, l] located in the overlapped regions cannot be estimated accurately. By contrast, taps of heff[k, l] located in the non-overlapped region can be estimated accurately.

[0289] That is, in FIG.19I, the dark gray dots mark the lattice Λqsupporting the self-ambiguity function of the spread pulse-tone with M = 11, N = 13, q = 4, (kp, lp) = (0, 0). The lattice Λqis generated by the vectors (24, 5) and (5, 7). The channel spread is the same as that in Figure 19H. The support of heff[k, l] is limited to the gray rectangle with the black border.

[0290] Note that the size and shape of the gray rectangle does not change with q since it depends only on the delay Doppler spreading function hphy(τ, ν) of the underlying physical channel, the pulse shaping filters wtx(τ, ν), wrx(τ, ν), and the information grid Λdd(which in turn depends only on T and B, and is independent of (τp, νp)).

[0291] Previously, we have emphasized that the (strong) crystallization condition is satisfied when the channel delay spread is less than the delay period and the channel Doppler spread is less than the Doppler period. We have also previously described a weaker mathematical condition that eliminates DD domain aliasing. When we use discrete chirp filters to construct spread pulse-tones, the self-ambiguity function is supported on a transformed lattice Λq. The significance of the weaker mathematical condition is that there are channels for which aliasing-free acquisition of heff[k, l] may not be possible with the period lattice Λp, but becomes possible with the transformed lattice Λq.

[0292] Consider a spread pulse-tone where the self-ambiguity function Axs, x s ^ k , l ^ is supportedon a lattice Λq. Define the support ^ ^ki, l i ^ of heff[k – ki, l – li] by. (62)

[0293]

[0294] ^ki, l i ^ ^0,0 ^ ^Λq.Note also that when (63) isdistinct taps in ^ ^0,0 ^ cannot differ by a lattice point in Λq.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0295] Lemma 1: Let ^k2^, l 2 ^ ^^ ^ k 1 ^ , l 1 ^ ^ , be two distinct taps in ^ ^0,0 ^ . If the weaker crystallizationcondition (63) is satisfied then . (64)

[0296] Proof: and^k2^ ^ k 1 ^, l 2 ^ ^ l 1 ^ ^ ^^ q when (63) is satisfied. This is because, if ^ 1 2 ^ 1 ^ were to be some^k ^ (max ^ k, l ^^^ ^0,0 ^ k ^ min ^ k , l ^ ^ ^ ^ 0,0 ^ k ) andDoppler spreadmin ^ k , l ^ ^ ^ ^^ l ) satisfy ∆k ∆l > MN = BT (MN = (Bτp)(Tνp) = BT). For a given Bto M and N such that ∆k< M and ∆l< N, so the channel does notsatisfy the (strong) crystallization condition with respect to the period lattice Λp. Here we show that aliasing-free acquisition of heff[k, l] may still be possible if the channel is underspread ( ^ ^0,0 ^ ^ MN ). The supportset of the effective channel filter is shown as the gray rectangle with the black border in FIG.19J. The dotted rectangle in FIG.19J is the smallest rectangle with axes parallel to the delay and Doppler axes thatcircumscribes ^ ^0,0 ^ . The crystallization condition is not satisfied with respect to the period lattice Λp sincethe area of the dotted rectangle ∆k∆lis greater than MN. However the (weaker) crystallization condition (63) is satisfied with respect to Λqbecause the gray rectangles are disjoint. Hence, we can accurately estimate the effective channel taps heff[k, l].

[0298] That is, in FIG. 19J, the dark gray dots mark the lattice Λqsupporting the self-ambiguityfunction of the spread pulse-tone with M = 11, N = 13, q = 3, (kp, lp) = (0, 0). The support ^ ^0,0 ^ of theeffective channel taps is shown as the gray rectangle with the black border. The dottedis thesmallest rectangle with axes parallel to the delay and Doppler axes that circumscribes ^ ^0,0 ^ .

[0299] Example 4: We consider the Veh-A channel with the power delay profile given in Table-I of section 5.2 and with path Doppler shifts ν1= νmax, ν2= –νmax, ν3= νmax / 2, ν4= –νmax / 2, ν5=νmax / 4, ν6= –νmax / 4. We now illustrate the importance180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00Ays, x s ^ k , l ^ (see (53)) between the received spread pulse-tone and the transmitted spread pulse-tone, forspread pulse-tones defined by chirp filters with q = 3 and q=36. Before we study these comparisons, through FIG.19K we illustrate the heat map for the effective discrete DD domain channel filter heff[k, l]. It is observed that the effective channel Doppler spread is higher for νmax= 12 KHz than for νmax= 815 Hz.FIG. 19K provides heat map of heff ^ k , l ^ for Veh-A channel, RRC pulse shaping filter (βτ = βν = 0.6),Doppler period νp= 30 KHz, M = 31, N

[0300] We first consider the spread pulse-tone defined by the chirp filter with q = 3. FIG.19L showsthe heat map of Ays, x s ^ k , l ^ for νmax = 815 Hz and νmax = 12 KHz [for chirp filter with q = 3 and Veh-Achannelfor both cases]. In both cases, the support sets ^ki , l i , (ki,li) ^Λq of heff[k–ki, l–li] do not overlap and the (weaker) crystallization condition is satisfied.

[0301] Next we consider the spread pulse-tone defined by the chirp filter with q = 36. FIG.19Mshows the heat map of Ays, x s ^ k , l ^ for νmax = 815 Hz and νmax = 12 KHz [for chirp filter with q = 36 andVeh-A channelFIG.19K for both cases]. When νmax= 815 Hz, the support sets do not overlap. When νmax= 12 KHz, the support sets do overlap and the (weaker) crystallization condition is not satisfied. This compromises the accuracy of channel estimation, which in turn degrades BER performance (see Section 5.5 for more details).

[0302] 5.5 Integrated Sensing and Communication (ISAC)

[0303] 5.5.A ISAC with Spread Sensing pulse-tone

[0304] When we integrate channel sensing and data transmission in the same Zak-OTFS subframe, the spread pulse-tone used for channel sensing interferes with the point pulse-tones used for data transmission. In this section, we describe how to design the spread pulse-tone so that interference is noise- like. This property can eliminate the need for a guard band at the cost of reducing the SNR for channel sensing. Noise-like interference translates to a spread pulse-tone with energy that is almost uniformly distributed over the discrete DD domain. This explains why the PAPR of the spread pulse-tone is significantly lower than that of the point pulse-tone.

[0305] The discrete DD domain transmit signal xdd[k, l] is given by (65)is the sum of a, , , givenby 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 . (66) MN information symbols,1 N ^ 1^ x 2d,d^ k,l ^ ^is^^ ^^d^ ^

[0307] It follows from. (68)

[0308] We definethat the received discrete DD domain signal can then be given bywhere heff[k, l] is the effective discrete DD domain channel filter.

[0309] Channel sensing: Here we suppose that the (weaker) crystallization condition (63) holds, and we describe how to estimate the effective channel filter heff[k, l] from the received signal ydd[k, l]. Recallfrom (52), that in the absence of data, the ML estimate hˆeff ^ k , l ^ is given by the cross-ambiguity functionbetween the received signal and the transmitted spread pulse-tone.

[0310] Theorem 3: In the presence of data, the cross-ambiguity function Ay, xs ^ k , l ^ betweenydd[k, l] and xs,dd[k, l] is given by 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 ,

[0311] When the crystallization condition (63) holds, we can obtain heff[k, l] by evaluating the firstterm on the RHS of (70) inside ^ ^0,0 ^ (see (62)). For ^k, l ^ ^^ ^0,0 ^ , the cross-ambiguity functionAy, xs ^ k , l ^ now reduces toterm in (72) which measures interference to sensing from data. This error term is given bydata to sensing contributes 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00to the mean squared

[0315] FIG. 19N shows the heat map for the quantities Axd , x s ^ k , l ^ which measure interferenceto sensing from data. We observe that the values are roughly^ 1 MN ^ ^ ^ 15 dB, in other words Axd , x s ^ k , l ^ ^ 1 MN .

[0316] The factor of 1 / (MN) is due to the fact that Axd , x s ^ k , l ^ is roughly of order 1 MN . InFIG.19N, we have plottedmap for Axd , x s ^ k , lN = 37, spread pulse-tone definedby chirp filter with q = 3 and for a random realization of the information symbols [including PDR = Ep / Ed= 0dB]. It is observed that the values are roughly of the order 10log10 ^ 1 MN ^ ^ ^ 15 dB .

[0317] The last term in (72) represents the contribution of noise to channel estimation error. Sincethe received DD domain noise samples are i.i.d. ^^ ^0, N0 ^ distributed, the samples of An, xs ^ k , l ^ areM ^ 1 N ^ 1 2zero mean. Since ^^ xs,dd^ k,l ^ = 1 , it follows from (70) that the variance is given byk ^ 0 l ^ 0. (76)

[0318] It follows from (73) that the normalized mean squared estimation error (NMSE) can be given bywhere we havexs , , ofthe data interference and noise terms in (72) from (75) and (76) respectively. In (77), ρdrefers to data SNR (see (40)) and ρprefers to the pilot SNR (see (41)). 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0319] Note, in practical implementations, while reading the estimate of heff[k, l] from Ay, xs ^ k, l ^ , ^ k , l ^ ^^ ^ 0,0 ^ , all taps of Ay, xs ^ k , l ^ in ^ ^0,0 ^ may not be genuine taps of heff[k, l] primarilydue to the randomness of noise and data interference. Therefore, in some embodiments, only those tapsexceeds a pre-determined threshold.times the standard deviation of the estimation error (heff ^ k, l ^^ hˆ eff ^ k , l ^(77)).

[0320] When the crystallization condition is satisfied, the factor 1 / (MN) in (77) significantly reduces the impact of data and noise on the accuracy of channel estimation. The factor 1 / (MN) is present because interference between the spread pulse-tone used forand the point pulse-tones used for data transmission is noise-like.

[0321] Cancellation of spread pilot followed by data detection: FIG.11B provides another example flow chart of a method proposed for integrated sensing and communication (ISAC). That is, FIG. 11B provides an example architecture of integrated sensing and communication (ISAC) with a spread pilot pulse-tone. We first estimate the received spread pilot pulse-tone as . (78) using the estimate hˆeff ^ k , l ^. the(weaker) crystallization condition is satisfied with respect to the lattice Λq. The estimate becomes more accurate for larger M, N, that is . (79)

[0322] . (80)Λp, we can reliably recover data from the almost sensing pulse-tone free received signal yd,dd[k, l] in (80).

[0324] This approach can increase effective throughput by eliminating the need to divide physical resources between sensing and communication.

[0325] The flow chart of an example spread pilot pulse-tone based joint sensing and communication is depicted in FIG.11B.

[0326] 5.5.B Numerical Simulations 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0327] As described in the present document, Zak-OTFS allows for ISAC with spread sensing pulse-tone where there is no sensing overhead and the transmit signal has low PAPR. In this section, we also see this is supported by numerical simulations of the BER performance of ISAC using point pulse- tones and spread sensing pulse-tones.

[0328] FIG.20A plots BER of uncoded 4-QAM as a function of increasing Doppler spread (2νmax) for the Veh-A channel described in Section 5.2 (see Table I and (19)). Point pulse-tones can be used for both channel sensing and data transmission.

[0329] In FIG.20A, a plot for Uncoded 4-QAM BER vs νmaxfor ISAC with point sensing pulse-tone ^for Veh-A channel is provided. Data SNR and PDR are fixed at ρd = 25 dB and p^d^ 10 dB respectively,with νp = 30 KHz, M = 31, N = 37; fixed guard region ^ (7 ^ 7 rectangle(kp, lp) = ((M + 1) / 2, (N +1) / 2)); RRC pulse shaping filter (βτ= βν= 0.6); and data detection with MMSE equalization.

[0330] Here, and throughout this Section, we consider channel sensing in the presence ofinterference from communication data ( S C ) and in the absence of interference from data (S C ).Similarly, we consider data equalization in the presence of interference from sensing (C S ) and in theabsence of interference from sensing (C S ).

[0331] S C and C S : This is the baseline (two bottom curves with the gray circle and trianglemarkers) where we dedicate separate Zak-OTFS subframes to channel sensing and to data transmission. Doppler domain aliasing increases with νmax, and BER increases as the accuracy of channel sensing degrades. We can reduce Doppler domain aliasing by substituting an RRC filter (gray triangles) for a sinc filter (gray circles).

[0332] S C and C S : When we integrate channel sensing and data transmission in the sameZak-OTFS subframe (solid black curves), the point pulse-tone used for channel sensing interferes with the point pulse-tones used for data transmission. This interference severely compromises BER performance when the Doppler spread exceeds 2 KHz (νmax> 1 KHz) because the channel response to a point pulse- tone extends beyond the guard band, regardless of the choice of filter. When the Doppler spread is less than 1 kHz, the choice of filter makes a significant difference. When we introduce a guard band, we divide DD domain resources between sensing and communication. By increasing the size of the guard band (sacrificing transmission rate) we can extend the range of reliable operation to higher Doppler spreads. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0333] S C and C S : We integrate sensing and data transmission in the same Zak-OTFSsubframe so the point pulse-tones used to transmit data interfere with the point pulse-tone used to sensethe channel (S C ). We ask what BER performance would be, if when recovering the data, there was nointerference from the sensing pulse-tone ( C S ). We estimate the channel within the guard band of theISAC subframe (S C ), then generate a data-only subframe carrying the same data as the ISAC subframe,then equalize the data using the channel estimate ( C S ). We observed little difference in BERperformance (compared to the solid black curves) and have not shown these curves in FIG.20A.

[0334] S C and C S : We ask what BER performance would be if channel estimation were notsubject to interference from data. FIG.20A illustrates BER performance when we use a separate Zak- OTFS subframe for channel sensing (tighter dashed black curves with black triangles and black circles). The improvement in BER performance is significant when compared to the solid black curves. For the tighter dashed black curves with black triangles and black circles, we observed only a small improvement in BER performance when we replaced the estimated channel by perfect channel state information (CSI), and we have not shown these curves in FIG.20A.

[0335] FIG.20B plots BER of uncoded 4-QAM as a function of increasing pilot to data power ratio (PDR) for the Veh-A channel described in Section 5.2. We fix the data SNR ρd= 25 dB, and we fix νmax= 815 Hz. By limiting the Doppler spread to 1.63 KHz, we limit the fraction of received pulse-tone energy that falls outside the guard band. Point pulse-tones are used for both channel sensing and data transmission, and we only consider RRC pulse-shaping filters.

[0336] In FIG.20B, a plot for Uncoded 4-QAM bit error rate (BER) as a function of the ratio of the ^ point sensing pulse-tone power to the total power of data pulse-tones p ^ d (PDR) for ISAC with pointsensing pulse-tones is provided: for Veh-A channel; RRC pulse (βτ= βν= 0.6); Data SNR fixed at ρd= 25 dB; Doppler period νp= 30 KHz; M = 31, N = 37; and the guard region is a 7 ^ 7 rectangle. Note the characteristic "U" shaped curves for integrated sensing and communication.

[0337] S C and C S : This is the baseline (gray curve with hollow triangles) where we dedicateseparate Zak-OTFS subframes to channel estimation and data transmission. The accuracy of channel sensing increases with PDR, but DD domain aliasing in the data subframe limits BER, resulting in an error floor.

[0338] S C and C S : The point pulse-tone used for channel sensing interferes with the pointpulse-tones used for data transmission. When PDR is small, the interference is small, and BER improves 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 as the channel estimate becomes more accurate. When the power of the pilot pulse-tone exceeds the total power of all data pulse-tones (PDR > 0 dB), interference becomes more significant than noise, and BER degrades with increasing PDR. This explains the characteristic "U" shape of the black curve with solid black triangles.

[0339] S C and C S : There is no interference from the pilot pulse-tone when recovering thedata, and BER is almost independent of PDR. This explains the flatness of the gray curve with solid gray triangles.

[0340] S C and C S : When there is a dedicated sensing subframe, the accuracy of the channelestimate improves with increasing PDR. When PDR > 0, interference from the pilot pulse-tone becomes more significant than noise, and BER degrades with increasing PDR. This explains why the black curve with hollow triangles has the same characteristic "U" shape as the black curve with solid black triangles. We observed only a small improvement in BER when we replaced the estimated channel by perfect CSI and we have not shown this curve in FIG.20B.

[0341] The normalized mean squared error (NMSE) of the channel estimate can be given by

[0342] FIG.20C plotsWe consider sinc and RRC pulse shaping filters to understand the significance of the fraction of received sensing pulse-tone energy that falls outside the guard band.

[0343] In FIG.20C, NMSE vs. PDR for ISAC with point sensing pulse-tone is provided, with NMSE of estimation of the taps of the effective channel filter heff[k, l] as a function of increasing PDR for a fixed data SNR of ρd= 25 dB, RRC filter (βτ= βν= 0.6), νp= 30 KHz.

[0344] S C : When the received pilot pulse-tone is confined within the guard band, NMSEdecreases linearly with increasing PDR (gray triangles). Otherwise NMSE decreases to a floor that depends on the fraction of the energy of the received pulse-tone that lies outside the guard band (gray circles).

[0345] S C : The limiting behavior is the same, but the floors are different (black triangles andblack circles). 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0346] FIG.20D plots NMSE as a function of increasing PDR for a spread sensing pulse-tone (q = 3) on the Veh-A channel considered in FIG.20A. The Doppler spread 2νmax= 1.63 KHz is significantly less than the Doppler period νp= 30 KHz, but even with RRC pulse shaping filters, there is residual DD domain aliasing.

[0347] In FIG.20D, NMSE vs. PDR for a spread sensing pulse-tone with q = 3 is provided, with Veh-A channel, RRC pulse shaping filter (βτ= βν= 0.6), data SNR ρd= 25 dB, νmax= 815 Hz, νp= 30 KHz, M = 31, N = 37.

[0348] S C : This is the baseline (gray curve) where we dedicate separate Zak-OTFS subframesto channel estimation and data transmission. At high PDR, the pilot power to noise power ratio (PNR) is high. NMSE saturates at high PNR due to DD domain aliasing.

[0349] S C : As PDR increases, the spread pilot becomes stronger than the combination of datainterference and noise. DD domain aliasing again causes NMSE to saturate at high PDR.

[0350] FIG.20E plots BER of uncoded 4-QAM as a function of increasing PDR. Note that with perfect CSI, the BER performance is independent of PDR (red curve).

[0351] In FIG.20E, BER vs. PDR for a spread sensing pulse-tone with q = 3 is provided, with Veh- A channel, RRC pulse shaping filter (βτ= βν= 0.6), data SNR ρd= 25 dB, νmax= 815 Hz, νp= 30 KHz, M = 31, N = 37.

[0352] S C and C S : This is the baseline (gray curve with solid gray triangles) where wededicate separate Zak-OTFS subframes to channel estimation and data transmission. BER performance tracks NMSE performance in FIG.20D, improving with increasing PDR and saturating at high PDR. Thegray curve with hollow triangles illustrates the case S C & C S where we estimate the channel in thepresence of interference from data. The gray curve with hollow triangles and gray curve with solid gray triangles differ by a small PDR offset.

[0353] S C and C S : The spread pulse-tone used for channel sensing interferes with the pointpulse-tones used for data transmission. When PDR < 10 dB, BER improves with increasing PDR. Beyond 10 dB, interference from the residual pilot pulse-tone (after cancellation) degrades BER with increasing PDR. This explains the characteristic "U" shape of the black curve with solid black triangles. When wededicate a separate Zak-OTFS subframe to sensing (S C & C S ) we observe only a small improvement180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 in BER, and we have not shown this curve in FIG. 20E. In this case, there is an optimal PDR which minimizes BER.

[0354] We use signal to interference ratio (SIR) to represent the ratio of the average power of the data pulse-tones to the average power of the residual sensing pulse-tones after cancellation. Figure 20F plots SIR as a function of increasing PDR for the Veh-A channel considered in FIGS.20D and 20E.

[0355] In FIG. 20F, Signal to Interference ratio (SIR) of data power to the power of the residual spread sensing pulse-tone (after cancellation) as a function of increasing PDR (for fixed data SNR ρd= 25 dB) is provided, with RRC pulse shaping filter (βτ= βν= 0.6). νp= 30 KHz, Chirp filter with q = 3.

[0356] S C and C S : Sensing is performed on a separate Zak-OTFS subframe (black curve withhollow triangles). Estimation accuracy improves with increasing PDR, but the energy of the residual pilot pulse-tone increases proportionally. This is because, there are always some taps of the effective channel filter that are not estimated. We sense taps within a support set, and outside this set, taps are noisy because they are subject to DD domain aliasing. This explains why SIR decreases with increasing PDR.

[0357] S C and C S : Here the channel estimate is slightly inferior to the estimate obtained froma dedicated sensing subframe. Hence the SIR ratio is smaller. The effect diminishes as PDR increases because the effect of data interference on sensing diminishes.

[0358] FIG 20G plots BER as a function of increasing νmax.

[0359] FIG 20G: BER vs νmaxfor spread sensing pulse-tone (q = 3). Veh-A channel model, RRC ^pulse shaping filter (β pτ = βν = 0.6), Doppler period νp = 30 KHz, PDR ^d^ 10 dB , data SNR ρd = 25 dB.

[0360] S C and C S : This is the baseline (gray curve with hollow gray triangles) where wededicate separate subframes to sensing and data transmission. There is no residual pilot pulse-tone to interfere with the data pulse-tones. Channel estimation is accurate even when νmaxis large. This is because the lattice Λqsupporting the self-ambiguity function of the spread sensing pulse-tone satisfies the crystallization condition (63), even for large νmax. This explains why BER is excellent and almostindependent of νmax. The gray curve with solid gray triangles (S C & C S ) differs from the gray curvewith hollow gray triangles by a small SNR offset. The black curve with solid gray triangles (C S , perfectCSI) approaches the gray curve with hollow gray triangles in the limit as νmaxapproaches 14 KHz. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0361] S C and C S : Here the spread pulse-tone used for channel sensing interferes with thepoint pulse-tones used for data transmission. Nevertheless, the variation in BER is small over a wide range of Doppler spreads. This is very different from BER performance for point sensing pulse-tones illustrated in FIG.20A (BER = 4 ^ 10−3at νmax= 300 Hz and BER = 2 ^ 10−2at νmax= 14 KHz in FIG.20G versus BER = 10−3at νmax= 300 Hz and BER = 0.4 at νmax= 14 KHz in FIG.20A). When we dedicate a separate Zak-OTFS subframe to sensing (S C and C S ), we observe a modest improvement in BER (black curve withhollow black triangles versus black curve with solid black triangles). We conclude that spread sensing pulse-tones improve upon point sensing pulse-tones by extending the range of reliable operation to a wider range of Doppler spreads.

[0362] FIG.20H illustrates the importance of filter design by comparing BER performance of chirp filters with q = 3 and q = 36. Accurate sensing requires crystallization with respect to the lattice Λq, and reliable data detection requires crystallization with respect to the period lattice Λp. The crystallization condition is satisfied for the period lattice Λpsince νmax< 15 KHz and the delay spread for the Veh-A channel is 2.5µs, which is less than τp= 1 / νp= 33.33µs. When q = 3, crystallization conditions w.r.t. both Λpand Λqare satisfied and BER performance is almost constant for a wide range of Doppler shifts (black curve with solid black triangles). When q = 36, the crystallization condition with respect to Λqis not satisfied for νmax> 2 KHz and BER degrades sharply (black curve with hollow triangles).

[0363] In FIG.20H, BER vs νmaxfor spread sensing pulse-tones defined by chirp filters with q = 3 ^and q = 36 is provided, with Veh-A channel described in FIG.20G, ν pp = 30 KHz, fixed PDR ^d^ 10 dB ,data SNR ρd= 25 dB.

[0364] Effective Throughput: This is defined to be the ratio of the number of bits reliably communicated in each subframe to the available degrees of freedom. In various embodiments, onlychannel sensing and data transmission in the same Zak-OTFS subframe (S C & C S ) is considered.When we use a point pulse-tone for channel sensing we need to introduce a guard band for various embodiments, and this division of DD domain resources between sensing and communication decreases effective throughput. When we use a spread pulse-tone for channel sensing we can avoid this overhead.

[0365] Assuming RRC pulse shaping filter with roll-off factors βτand βν, the number of degrees of freedom is BT (1 + βτ)(1 + βν). Setting BER = R, the number of reliably communicated bits in each subframe is 2MN(1 – H(R)) when using a spread sensing pulse-tone, compared with 2(MN – 49)(1 – H(R)) when using a point sensing pulse-tone (assuming a 7 ^ 7 guard band) and 2(MN – 7N)(1 – H(R)) (assuming the 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 guard region to be a strip along the Doppler domain having width 7 along the discrete delay axis). Here H denotes the binary entropy function and we note that for small R we have H(R) ^ R / ln 2.

[0366] FIG.20I plots effective throughput as a function of increasing νmax. Recall from FIG.20H that for a spread sensing pulse-tone defined by a chirp filter with q = 3, BER performance is almost constant for a wide range of Doppler spreads. This is reflected in FIG.20I. Recall from FIG.20H that when q = 36, the crystallization condition with respect to Λqis not satisfied for νmax> 2 KHz, and the BER degrades sharply. This explains the difference in effective throughput for spread sensing pulse-tones defined by chirp filters with q = 3 and q = 36.

[0367] In FIG.20I, effective throughput (bits / sec / Hz) as a function of increasing νmax, is provided forintegrated sensing and communication (S C & C S ) with νp = 30 KHz, M = 31, N = 37, RRC pulseshaping filter (βτ= βν= 0.6), Veh-A channel as described in FIG.20G.

[0368] Recall from FIG.20A that for νmax> 1 KHz, the channel response to a point sensing pulse- tone extends beyond the guard band regardless of the choice of pulse shaping filter. Interference between the point pulse-tone used for sensing and the point pulse-tones used for data transmission can severely compromise BER performance. This explains the characteristics of the effective throughput curve for a point sensing pulse-tone (dashed curve with solid triangles).

[0369] FIG. 20I illustrates the gain in effective throughput that results from sharing DD domain resources between sensing and communication.

[0370] Summary for Section 5. Zak-OTFS and ISAC

[0371] We started from a parametric family of pulse-tone waveforms that can be matched to the delay and Doppler spreads of different propagation environments. We reviewed how the (point) pulse-tone signal in the time domain realizes a quasi-periodic localized function on the delay-Doppler (DD) domain. The characteristic structure of a pulse-tone is a train of pulses modulated by a tone, a signal with unattractive peak-to-average power ratio (PAPR). We reviewed system performance in the crystalline regime where the delay period of the pulse-tone is greater than the delay spread of the channel, and the Doppler period of the pulse-tone is greater than the Doppler spread of the channel. When channel sensing and data transmission take place in separate subframes, the point pulse-tone used to sense the channel does not interfere with the point pulse-tones used to transmit data. We have reviewed how the I / O relation of the sampled communication system can be read off from the response to the point pulse-tone used for channel sensing. In the present document, we have described filtering in the discrete delay-Doppler domain. We have described and demonstrated that it is possible to construct spread waveforms with desirable characteristics by applying a chirp filter in the discrete DD domain to a point pulse-tone. One desirable characteristic is low PAPR, about 5 dB for the exemplar spread pulse-tone, compared with about 15 dB for 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 the point pulse-tone. 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. We have described and demonstrated how to integrate sensing and communication within a single Zak-OTFS subframe through the combination of a spread pulse-tone used for channel sensing and point pulse-tones for data transmission. The filter in the discrete DD domain enables coexistence by minimizing interference between sensing and data transmission. We have described and 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.

[0372] 5.6 Zak-OTFS and ISAC Additional Information

[0373] Zak Transform and Quasi-periodicity

[0374] The definition of a DD domain pulse depends on a delay period τp, and a Doppler period νp,where the two periods are reciprocal, that is νp τp = 1. The Zak transform, denoted ^ t , provides a unitaryequivalence between TD signals and a subclass of quasi-periodic DD domain signals. The Zak transform of a TD signal x(t) is given by

[0375] The. (83) We observe thatalong the delay axis with period τp.

[0376] The inverse Zak transform, denoted by^ ^1t, of a quasi-periodic DD domain signal xdd(τ, ν)is given by

[0377]

[0378] The following theorem describes the Zak-OTFS carrier waveforms (pulse-tones). 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0379] Theorem 5: The transmit TD signal std(t) is given by where std,k,l(t) is theby whereis the Fourier. where180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 ,

[0381] Remarks: From (89) andthat the DD representation of the Zak-OTFS carrier waveform is given by^ 0 , is located at (kτp / M, lνp / N) (see FIG.17 that represents a Zak-OTFS carrier waveform).

[0382] We now make fundamental observations about the TD and FD realization of the Zak-OTFS carrier waveform.

[0383] Carrier waveforms having unlimited time and bandwidth: For wtx(τ, ν) = δ(τ) δ(ν), it follows from (86) that . (92)

[0384] For (k, l) = (0, 0),(93) which is simply an infinite(k, l)-th carrier waveform std,k,l(t) = xk,l(t) is nothing but std,0,0(t) with a Doppler shift of lνp / N and a delay shift of kτp / M, since it can be checked that .180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0385] Note that the (k, l)-th carrier waveform is the TD realization of a quasi-periodic DD domain k^^ pM , lv pimpulse at N ^ (see RHS of (90)). The (k, l)-th carrier waveform, which we refer to as a TD pulse-tone, is simply a train of impulses where the n-th impulse is located at t=nτp+ kτp / M, modulated by a ej 2 ^lvpN ( t^ k ^pM )complex tone / sinusoid (see (86)).

[0386] FD representation of the (k, l)-th TD pulse-tone is an FD impulse ppe^j 2^fk^p Mtrain (with m-th impulse at f = mν + lν / N) modulated by a FD tone . We refer to this waveform as a FD pulse-tone. When wtx(τ, ν) = δ(τ) δ(ν), it follows from (86) and (88) that the corresponding carrier waveforms have infinite time duration and bandwidth.

[0387] Carrier waveforms with limited time and bandwidth: The carrier waveforms can be limited in time duration and bandwidth by choosing an appropriate wtx(τ, ν). If the delay spread of wtx(τ, ν) is approximately 1 / B, and the Doppler spread is approximately 1 / T, then the TD pulse-tones have time duration T and bandwidth B. It follows from (86) that std,k,l(t) is simply xk,l(t) spread by 1 / B in the time domain. Given the structure of xk,l(t), it is clear that the TD pulse-tone std,k,l(t) consists of a train of narrow pulses modulated by a tone, with each narrow pulse having spread 1 / B and adjacent pulses separated by τpseconds. Similarly, from the integral expression in (87) it follows that the FD representation of the (k, l)- th carrier waveform, i.e., sfd,k,l(f), is simply Xk,l(f) spread by 1 / T in the frequency domain, since the Doppler domain spread of wtx(τ, ν) is 1 / T. Given the structure of Xk,l(f) in (88), it follows from the integral expression in (87) that the FD representation of the carrier waveform, i.e., sfd,k,l(f) consists of a train of narrow FD pulses modulated by a FD tone, with each narrow pulse having spread 1 / T and adjacent pulses separated by νpHz. Therefore, std,k,l(t) also has a pulse-tone structure in the FD, and is referred to as a FD pulse- tone. FIG.17the TD and FD representations of the Zak-OTFS carrier waveform.

[0388] Time and Band-limited pulse-tones

[0389] For factorizable pulse shaping waveform wtx(τ, ν) = w1(τ) w2(ν), it follows from (86) that the (k, l)-th carrier waveform is given by180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 where (108) is the inverse Fourier. (109)

[0390] Then, x2,k,lwe obtain where denoteswhere W1(f) and X2,k,l

[0391] It is clearofw1(∙). For example, choosing w1 ^^ ^^ Bsinc ^ B ^ ^ limits the bandwidth of the carrier waveforms toexactly B Hz. In general, for a given bandwidth constraint B, w1(τ) must have a spread of approximately 1 / B along the delay domain. Similarly, from (109) it follows that the duration of x2,k,l(t) can be limited to approximately T seconds by choosing the factor pulse w2(ν) to have a spread of approximately 1 / T alongthe Doppler domain (for example, with w2 ^^ ^^ T sinc ^ ^ T ^ , x2,k,l (t) is limited exactly to the TDinterval [−T / 2, T / 2]). Therefore, x2,k,lonly T / τp= N number of Dirac-delta 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 impulses of the infinite impulse train xk,l(t). Also, due to the TD convolution in (110), in the carrier waveform std,k,l (t) these N impulses are spread over a duration of approximately 1 / B seconds since the spread of w1(τ) is 1 / B. Therefore, under a finite duration and bandwidth constraint, the (k, l) carrier waveform consists of narrow pulses at t = nτp+ kτp / M where the width of each pulse is 1 / B = τp / M (since M = Bτp). Note that τp / M is also the time between the location of the n-th pulses of the (k, l)-th and the (k+ 1, l)-th carrier waveforms.

[0392] The Discrete Ambiguity Function

[0393] The discrete cross-ambiguity function Aa,b[k, l] between two discrete quasi-periodic DD domain functions a[k, l] and b[k, l] can be given byperiod N.

[0395] Theorem 6: Let a[k, l] and b[k, l] be quasi-periodic DD functions related by . (114)

[0396] Then, thel] is given by the twisted convolution between g[k, l] and the cross ambiguity between b[k, l] and c[k, l], i.e.[k, l] and b[k, l] is invariant w.r.t. the period over which the sum in the R.H.S. of (113) is computed, i.e., for any k0, l 0 ^^180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00(118) and . (119)

[0400] Substituting (118) gives180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00the summation variables k′ and l′ span a full delay and Doppler period respectively, it follows that as k′ varies from k0to k0+ M − 1, k′ mod M takes all possible values from 0 to M − 1. Similarly, as l′ varies from l0to l0+ N − 1, l′ mod N takes all possible values from 0 to N – 1. Therefore, replacing summation variable k′and l′ with k^^ k ^ mod M and l^ ^ l ^ mod N respectively, we get

[0402] Lemma 3: For any two quasi-periodic discrete DD domain functions a[k, l] and b[k,l] . (122)

[0403] b,180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00where in the second step the summation variables k′,l′, are changed to k^ ^ ^ k^ ^ k ^ and l^ ^ ^ l^ ^ l ^respectively.

[0404] Self-ambiguity Function of the Point Pilot Signal

[0405] It follows from the definition of the ambiguity function (113) that the self-ambiguity of the point pilot signal is given by (124),(a), from the location of the Dirac-delta impulses in the R.H.S., it follows that k′ ≡ kpmod M and (k′ − k) ≡ kpmod M. This is only possible if k is an integer multiple of M. Similarly, l′ ≡ lpmod N and (l′ − l) ≡ lpmod N. This is only possible if l is an integer multiple of N, in other words the self-ambiguity function Axp, x p ^ k , l ^ is non-zero for only the DD points ^nM, mN ^ , n , m^^ , which is simply the period lattice.Again, step (a) it is also clear that, for (k, l) = (nM, mN), the only non-zero contribution in the sum in the R.H.S. is from the term corresponding to n1= (kpmod M – kp) / M, n2= n1– n, m1= (lpmod N – lp) / N and j2 ^ n1 l^^ j 2 ^ n 2 ^ l ^^ l ^ ^ j l ^ k ^ ^ k ^m = m − m. T p p e N e N e 2 ^ MN2 1 he value of x , x ^ nM , mN ^ isevaluatedthe values of n1, n2, l′,k′ mentioned above. This leads p, x p , l in step (b).

[0407] Symmetry of Cross-ambiguity

[0408] First

[0409] 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO002 toobtain

[0413] The expression of the received signal ydd[k, l] in (69) consists of the sum of the received data signal, pilot signal and noise. The cross-ambiguity of the received signal with the spread pilot signal 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 xs,dd[k,l] is therefore the sum of, (i) the cross-ambiguity of the received spread pilot signal^^^^ℎ^^^^^^, ^^^ ∗ఙ ^^^,ௗௗ^^^, ^^^ with xs,dd[k,l]. Applying Theorem 6 (with g[k,l] = heff[k,l], a[k,l] =^^^^ℎ^^^^^^, ^^^ ∗ఙ ^^^,ௗௗ^^^, ^^^, b[k,l] = c[k,l] = xs,dd[k,l]) gives the first term in the R.H.S. of (70), (ii) the cross-of the received ^^^ ∗ఙ ^^ௗ,ௗௗ^^^, ^^^ with xs,dd[k,l]. Applying,^^^ ∗^^ ^^^^,^^^^^^^, ^^^, b[k,l] = xd,dd[k,l], c[k,l] = xs,dd[k,l])gives the second term in the R.H.S. of (70), and (iii) the cross-ambiguity between the noise signal ndd[k,l] and xs,dd[k,l], gives the third term in the R.H.S. of (70).

[0414] Derivation of the Maximum Likelihood (ML) Estimate of hEFF[k,l]

[0415] The maximum likelihood (ML) estimate of heff[k,l] can be given byit minimizes the energy of the estimation error signal edd[k, l] = (ys,dd[k, l] – h[k, l] *σxs,dd[k, l]). Since both ys,dd[k, l] and xs,dd[k, l] are quasi-periodic and twisted convolution conserves quasi-periodicity, it follows that the error signal is also quasi-periodic. The expression in (129) then follows from the fact that the energy of M^ 1 N ^ 12any discrete quasi-periodic DD domain signal edd[k, l] is given by ^^ edd^ k,l ^ .k ^ 0 l ^ 0

[0416] The expression for h[k, l] *σxs,dd[k, l] follows from the definition of twisted convolution and is given by (130)180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 Further analysis of (129) givesto the minimization. Next, inin the crystalline regime, the term M^ 1 N ^ 1^^ h ^ k,l ^ ^^ xs,dd^ k,l 2^ is equal to the total energy of the taps of h[k, l].180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 Using this in (131) gives (133). simplicity, we consider (kp, lp) = (0, 0), but the result is the same for any choice of point pilot. From (49), the spread pilot corresponding to (kp, lp) = (0, 0) is given by^ ^ ^ ,MN – 1 and l’ to l’ = 0,1, ^ ^ ^, MN-1, i.e.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0422] Substituting for xs,dd[k’, l’] using (134) gives

[0424] In (138), the inner summation over k’ vanishes unless (2qk − l) ≡ 0 (mod M) so that (139)^ , N − 1, there is a unique n2∈ {0, 1, ^ ^ ^ ,N − 1} such that (140)180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 that is given by The inner(142) and (n1, n2) satisfy (141). When these conditions are satisfied, the inner summation equals 1.

[0426] We need only evaluate the inner summation in (138) over l′, for (k, l) satisfying (142) and (n1, n2) satisfying (141). This summation vanishes unless (143) It follows from (141) that Setting(145) so that

[0427] Note that (142) and (146) are the two arithmetic conditions appearing inas (55). Thus the inner summation over l′ vanishes unless (k, l) satisfy (55). In addition, it follows from (144) that , (147) and it follows from (145) that(148)

[0428] appear . assume (k, l) satisfy (55), and rewrite (138) as the product of two terms in (149). We simplify the sum over n1by substituting for n2M using (144) to obtain (150). We simplify the sum over m1by substituting for 2qm2N using (145) to obtain (151). Combining (150) and (151) in (149) we obtain Theorem 2. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0429] Average Energy of DD Noise Samples

[0430] We follow the noise through the Zak-OTFS signal processing architecture described in Section 5. At the receiver, the Zak transform ^^௧converts AWGN ntd(t) with PSD N0to its DD representation ndd(τ, ν) = ^^௧(ntd(t)). Pulse shaping with wrx(τ,ν) results in ^^ௗ^ௗ^௫(τ,ν) = wrx(τ, ν)∗σ180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 ndd(τ, ν) which is then sampled on the information grid Λdd to give the discrete DD domain quasi-periodic noise signal ndd[k, l].

[0431] Here, we show that the expected energy of each noise sample ndd[k, l] satisfies^^^|^^ௗௗ^^^, ^^^|ଶ^ ൌ ^^^. We consider receive pulse shaping filters wrx(τ, ν) that factor as wrx(τ, ν) = w1(τ )w2(ν) where w1(τ), w2(ν) are unit energy pulses along the delay and Doppler axes respectively. We assume w1(τ), w2(ν) satisfy the Nyquist no-ISI criterion for information spacing 1 / B and 1 / T along the delay and Doppler domain respectively (^|^^^^^^^|ଶ^^^^ =^|^^ଶ^^^^|ଶ^^^^ = 1).

[0432] The DD representation of AWGN ntd(t) is given by

[0434] 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 where

[0435] is the tosubstitute for ndd(τ,ν). Step (b) follows from the identity T = Nτpafter changing the integration variable from τ′ to τ .

[00436] It follows from (154) that ndd [k, l] is zero mean and that ^^^|^^ௗௗ^^^, ^^^|ଶ^ is given by(156).^^^^^^τ^. Next, we observe that for large M, the spread of w1(τ) along the delay axis is significantly less than τp(since M = Bτp, the spread 1 / B = τp / M ≪ τp). Hence for n1≠ n2, the delay domain supports of w1(τ1) and ^^^∗^^^^ ^ ^^^ଶ െ ^^^^^^^^ do not overlap. Step (b) in (156) now follows from the fact that we canto the summation for which n1≠ n2.

[0438] The spread of the pulse w2(ν) along the Doppler domain is roughly 1 / T . Therefore the spread of its inverse Fourier transform W2(t) is roughly T = Nτp. Since w2(ν) satisfies the Nyquist no-ISI criterion along the Doppler domain with symbol spacing 1 / T, it follows that its inverse Fourier transform (i.e., W2(t)) must satisfy180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0439] for all t. The sinc and RRC pulses used in this paper are examples of Nyquist pulses. Using (157), we show that for all tindex n by (mN + q) where q ≡ n (mod N). Combining (156) and (158) we obtainT = Nτp.

[0442] Proof of Theorem 4

[0443] In this Appendix we calculate the expected interference from data transmission tosensing. We start by combining (66) and (71) to give (160). The information symbols ^^^^^, ^^^ have unitenergy and are statistically independent. Hence the mean squared interference energy is given by (161). Through (162), we simplify the inner summation in (161) using Lemma 2 (the self-ambiguity function^^௫ೞ,௫ೞ^^^, ^^^ does not depend on the period over which the summation is carried out). We now obtain(163) by substituting (162) in (161). Step (a) in (163) follows from Lemma 1 when the (weak) crystallization condition is satisfied. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 6. Examples of filters in the discrete DD domain and coexistence of sensing and communication in the same Zak-OTFS subframe

[0444] In this section, examples of symplectic transform operations are provided. In various embodiments, these example transform operations can be used for filters in the discrete delay-Doppler domain.

[0445] Radar engineers construct waveform libraries for tracking applications by manipulating self-ambiguity functions.

[0446] Taking a waveform ^ , displacing it in delay and Doppler, and taking the inner product with the original waveform gives the self-ambiguity function in various embodiments. When the absolute value of the self-ambiguity function is high at some location, then targets at this location can be obscured. Radar engineers have learned that chirping the waveform ^ produces a new waveform ^qand a new self-ambiguity function can be obtained by rotating the original self-ambiguity function. For example, chirping transforms the self-ambiguity surface, making it possible to track a greater diversity of radar targets.

[0447] The volume under the self-ambiguity function is fixed by Moyal’s Identity (as the square of the signal energy) but it is possible to make more targets visible by rotating the self-ambiguity function.

[0448] Chirping or Linear Frequency Modulation (LFM) is a particular example of a symplectic transformation. Example symplectic transforms in radar from Howard et al. are provided below. Displacement operator: ^^^^^, ^^^^^^^^^ൌ ^^ିఔ ఛ^^^ఔ௧^^^^^ ^ ^^^Self-ambiguity^^^ ^^^,^^^^^, ^^^,^^^ Symplectic transform: Linear Frequency Modulation (LFM)th slope ^^: ^^ ^^ : ^^ ^^ → ^^^^^ ൌ ^^^ ^మ Chirp wi^ ^ ^ ^ ^ ^ ௧^^^^^^Symplectic transforms rotate self-ambiguity functions ^^థ^^^^, ^^^ ൌ ^^^^,^^^^^, ^^^,^^^^ ൌ ^^^ , ^^^^^^ା ^^^^^, ^^^ ^^^^^^, ^^^

[0449] are described and provided below.

[0450] Communication in the discrete DD domain is revisited – the counterpart of LFM modulation is an MN-periodic discrete chirp filter [i.e., M*N-periodic discrete chirp filter].

[0451] The dimensions of the delay-Doppler grid can be chosen according to the channel characteristics and the available bandwidth and frame duration. The parameters M, N are directly related to the delay-Doppler grid periods (tau_p, nu_p) [e.g., (^p, ^p)] and can be selected such that the OTFS is in its crystalline regime. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0452] This filter can act on a point pulse-tone by twisted convolution to produce a spread pulse- tone xs.

[0453] The spread pulse-tone xsis quasi-periodic, hence periodic with period MN along the delay and Doppler axes.

[0454] In the LTI world, when a linear filter is applied to a signal with period L, the result can be given by periodic linear convolution of the periodic signal with the periodic extension of the filter.According to the notation used in the present document, w_s[k,l] (i.e., ^^^^^^, ^^^) is a discrete filter definedover the delay-Doppler grid and w[k,l] is its MN-periodic [i.e., M*N-periodic] extension of it. Note, that there are two periodicities here. The discrete filter w_s[k,l] is quasi-periodic in the delay-Doppler with respect to the dimensions M, N. On top of that, w[k,l] is an MN-periodic extension of w_s[k,l] that satisfies: w[k,l] = sum(w_s[k + n*MN, l + m*MN) for all integers n, m.

[0455] In the LTV world, there is a twist! The subscript s indicates the spreading operation.

[0456] In the LTI world there are L complex degrees of freedom and in the LTV world there are M2N2complex degrees of freedom.

[0457] In various embodiments, a ^^^^-Periodic Discrete Chirp Filter with slope ^^ (coprime to ^^^^) can be represented as [as also discussed in the previous section] ^^ ^^, ^^ ൌ1 ^൫^మା^మ൯ ^^^^^ଶగ^^^^ெே , ^^^^^^ ^^, ^^ ∈ ℤ

[0458] The spread pilot as: ^^^,ௗௗ^^^, ^^^ ൌ^^^^^^, ^^^ ∗ఙ ^^^,ௗௗ^^^, ^^^, where ^^^,ௗௗ period ^^^^ along delay andDoppler axis.

[0459] ^^^^-Periodic twisted convolution of ^^^,ௗௗ^^^, ^^^ with ^^^^^, ^^^(the ^^^^-periodic extension ofthe discrete filter) can be represented as:

[0460] A theorem

[0461] Forby a vector in ℂெమேమ.

[0462] FIG.18 graphically depicts examples where chirp filters rotate the period lattice.

[0463] Three lattices play an important role – the first is the information lattice (e.g., Λௗௗ) where the data lives, the second is the period lattice ^p, and the third is the dual of the information lattice (e.g., Λ^ௗௗ). The information lattice can be thought of as ZxZ, the period lattice as MZxNZ, and the dual of the information lattice as MNZxMNZ.

[0464] The self-ambiguity function of a point pulse-tone (e.g., point pulse-tone xp) is supportedon the period lattice ^p. For example, ^^௫^,௫^^^^, ^^^ is supported on the period lattice Λ^.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0465] The self-ambiguity function of a spread pulse-tone (e.g., spread pulse-tone xs) ^^௫ೞ,௫ೞ^^^, ^^^is supported on a rotated lattice ^q.

[0466] A theorem (e.g., Theorem 2 from section 5.4.B) provides that the ambiguity function^^௫ೞ,௫ೞ^^^, ^^^ is supported on a rotated lattice Λ^ given by^1 െ 2^^^^^Λ^ ൌ ^ ^^ 11 2^^൨ Λ^ where ^^ ൌ ^2^^^ି^ െ 2^^(mod ^^^^) .

[0467] spreading operations in time domain (TD) and delay-domain).

[0468] The effect of spreading in the DD domain is to distribute energy equally across all point pulse-tones, resulting in a unit energy signal with energy about 1 / MN at each location on the informationgrid. For example, Energy ห^^^,ௗௗ^^^, ^^^หଶ ^ ^ ெே .

[0469] The effect of spreading in the TD is to produce a noise-like waveform that is much less peaky than a point pulse-tone, where the peaks in the spread TD waveform can be an order of magnitude smaller than the peaks in the point TD waveform.

[0470] The examples of FIG.20J are based on chirp filter with slope q = 3.

[0471] FIG.20K is a graph showing peak to average power ratio (PAPR) result examples. Here PAPR is compared. Complementary CDF (CCDF) [complementary cumulative distribution function] of theInstantaneous to Average Power Ratio (IAPR) with RRC transmit filter with roll-off factors ^^^ఛ ൌ ^^ఔ ൌ 0.6^is provided for different cases. Here, RRC stands for root-raised cosine filter. The chirp can be applied as a discrete filter on the delay-Doppler symbols. After that, a non-discrete 2-D filter (RRC) can be applied in delay-Doppler, to shape the OTFS waveform in time and frequency. In the result examples of FIG.20K, it can be seen that performance of Spread pulse-tone is about 5 dB, point pulse-tone is about 15 dB, Data + spread pulse-tone is about 8 dB, and Data + point pulse-tone is about 12 dB. That is, a summary of the IAPR is as follows: Spread pulse-tone ~ 5 dB vs. Point pulse-tone ~ 15 dB, and Data + spread pulse-tone ~ 8 dB vs. Data + point pulse-tone ~ 12 dB.

[0472] FIG.21 is a graphical example of a sensing waveform. Here the geometry of sensing with spread pulse-tones is described. The q value can affect the grid rotation directly. A q value needs to be selected that will not result in a rotation that causes aliasing and takes it out of the crystalline regime. Inthis example, M=11, N=13, q=5 with pilot location: ൫^^^, ^^^൯ ൌ ^0,0^.

[0473] The dark gray dots mark the rotated lattice ^q– this is the support of the self-ambiguity function of the spread pulse-tone.

[0474] The gray rectangle with the black border, denoted by ^, marks the support of theeffective discrete DD domain channel (e.g., ℎ^^^^^^, ^^^).

[0475] Observe that the translates of ^ by lattice points in ^qare disjoint (e.g., translates of ^ by ^qdo not overlap), which means that there is no DD domain aliasing. Hence, it is possible to 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 accurately estimate the channel on the information grid (the lattice ^ddor ^p) from the response received within ^; here, the crystallization condition holds with respect to the rotated lattice ^q.

[0476] FIG.22 graphically depicts examples of effect of spreading. The depicted example heatmaps of |ℎ^^^^^^, ^^^| are represented for Veh-A channel with RRC pulse shaping filter ^^^ఛ ൌ ^^ఔ ൌ 0.6^,Doppler period ^^^ ൌ 30 KHz, ^^ ൌ 31, ^^ ൌ 37. Here, the spread pulse-tone (e.g., Pulsone) is defined bythe chirp filter with q=3. The bottom heat maps are example heat maps of thepulse-tone defined by the chirp filter with q=3.

[0477] In these examples, xsis sent and ysis received and the cross-ambiguity between ysand xs is used to visualize the effective channel.

[0478] Good geometry means that the translates of the support set do not overlap and the crystallization condition is satisfied with respect to the rotated lattice ^q. For example, good geometry is indicated in the examples of FIG.22 with a chirp filter with q=3.

[0479] FIG.23 graphically depicts examples of effect of spreading. The depicted example heatmaps of |ℎ^^^^^^, ^^^| are represented for Veh-A channel with RRC pulse shaping filter ^^^ఛ ൌ ^^ఔ ൌ 0.6^,Doppler period ^^^ ൌ 30 KHz, ^^ ൌ 31, ^^ ൌ 37.Here, the spread pulse-tone (e.g., Pulsone) is defined bythe chirp filter with q=36. The bottom heat maps are example heat maps of the received spread pulse- tone defined by the chirp filter with q=36.

[0480] In an example of FIG.23, when ^max= 815 Hz, the geometry is good - the translates of the support set do not overlap and the crystallization condition is satisfied with respect to the rotated lattice ^q.

[0481] In the example of FIG.23, when ^max= 12 KHz, the geometry is bad - the translates of the support set overlap, the crystallization condition is not satisfied with respect to the rotated lattice ^q, and DD domain aliasing compromises accuracy of channel sensing.

[0482] FIG.24 shows an example relationship among filters in discrete DD domain. The diagram captures the different components in Zak-OTFS transceiver processing. When separate Zak-OTFS subframes are dedicated to sensing and communication, transceiver signal processing can involve the components shown in the continuous DD domain and the time domain.

[0483] In this Section the discrete DD domain was focused on. It was described how filters in the discrete DD domain define spread waveforms with low PAPR – about 5 dB compared with about 15 dB for a point pulse-tone (e.g., Pulsone).

[0484] It was shown that it is possible to acquire the effective channel from the response to a single pilot when the crystallization condition is satisfied with respect to the rotated lattice ^q.That is, it was described that it is possible to construct spread waveforms with desirable characteristics by applying a chirp filter in the discrete DD domain to a point pulse-tone. It was also described that it is possible to read off the I / O relation provided a second crystallization condition is satisfied with respect to ^q. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0485] 7. Examples of ISAC in the same Zak-OTFS subframe

[0486] Simulation results that illustrate BER performance as a function of key systems parameters such as pilot power, data power and maximum Doppler shift are provided.

[0487] FIG.25 shows an example of multiplexing data with a spread pilot signal. This is an example of integrating sensing and communication with a spread pulse-tone. A spread pilot xs and a point data signal xdare transmitted. Here, the data can interfere with the spread pilot.

[0488] In the discrete DD domain, a signal ^^ௗௗ^^^, ^^^ is received, ^^ௗௗ^^^, ^^^ ൌ ^^ௗ,ௗௗ^^^, ^^^ ^ ^^^,ௗௗ^^^, ^^^ ^ ^^ௗௗ^^^, ^^^.

[0489] The signal ^^ௗௗ^^^, ^^^ comprises a contribution ^^^,ௗௗ^^^, ^^^ from the pilot and a contribution^^ௗ,ௗௗ^^^, ^^^ from the data and a noise term (^^ௗௗ^^^, ^^^) which can be neglected for simplicity.

[0490] In the absence of noise, ^^ௗௗ^^^, ^^^ can be represented as^^ௗௗ^^^, ^^^ ൌ ^^^ௗ൫ℎ^^^^^^, ^^^ ∗ఙ ^^ௗ,ௗௗ^^^, ^^^൯ ^ ^^^^൫ℎ^^^^^^, ^^^ ∗ఙ ^^^,ௗௗ^^^, ^^^൯.that can be used to make thefrom the data signal.

[0492] The cross-ambiguity ^^௬,௫ೞ^^^, ^^^ between ^^ௗௗ^^^, ^^^ and ^^^,ௗௗ^^^, ^^^ can be represented as^^௬,௫ೞ^^^, ^^^ ൌ ^^^^൫ℎ^^^^^^, ^^^ ∗ఙ ^^௫ೞ,௫ೞ^^^, ^^^൯ ^ ^^^ௗ൫ℎ^^^^^^, ^^^ ∗ఙ ^^௫^,௫ೞ^^^, ^^^൯.noise like, is described below.

[0494] The crystallization condition holds with respect to the rotated lattice ^q.

[0495] The translates of the channel support S(0, 0) by lattice points in ^qare disjoint. This means that there is no DD domain aliasing, and it is possible to accurately estimate the channel on the information grid (the lattice ^ddor ^p) from the response received within S(0, 0). ^^ ^ห^^^^^^^^^^^^^^^^, ^^^ ∗^^ ^^^^^^,^^^^^^^, ^^^ห^^ ^ =^^^^^^^^∑^^^,^^^∈^^^^^,^^^ |^^^^^^^^^^^, ^^^|^^

[0496] is independent of the location [k, l]. This contribution to the mean squared error appears as an SNR (signal to noise ratio) hit, and the factor MN in the denominator makes this SNR hit less significant.

[0497] The channel can be estimated. For example, estimate ℎ^^^^^^, ^^^ by reading off ^^௬,௫ೞ^^^, ^^^ in^^^^,^^, ^^௬,௫ೞ^^^, ^^^ ൌ ^^^^ℎ^^^^^^, ^^^ ^ ^^^ௗ൫ℎ^^^^^^, ^^^ ∗ఙ ^^௫^,௫ೞ^^^, ^^^൯

[0498] for the received pilot from^^ௗௗ^^^, ^^^ to obtain an estimate for ^^ௗ,ௗௗ^^^, ^^^. (The estimate for ^^ௗ,ௗௗ^^^, ^^^ = y^^,^^^k, l^).

[0499] Then recover the data from this estimate y^^,^^^k, l^.180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0500] FIG.26 is a graphical depiction of an example of sensing accuracy as a function of pulse dynamic range (PDR).

[0501] In the subsequent figures, a standard Veh-A channel will be referred to. This standardchannel will be described as a standard Veh-A channel with Doppler period ^^^ ൌ 30 KHz, M = 31, N = 37,and RRC pulse shaping filter (^^ఛ ൌ ^^ఔ ൌ 0.6).

[0502] Here, in this FIG.26, sensing accuracy (NMSE) is plotted as a function of the ratio of pilot energy to data energy (PDR). NMSE can be given by ∑ห ^ ^ ^ ^ ^మ ^ೖ ^^^^ ^,^ ି^^^^ ୩,୪ หNMSE:,^^∈ೄ∑^ೖ,^^∈ೄ |^^^^^^,^^|మ.

[0503] Here, ^^^^௫ ൌ 815 Hz.

[0504] The communication (S| C bar or S|C), and as PDR increases the pilot overwhelms the noise. In this example (S|C), there is no interference from data.

[0505] The black curve is sensing and communications in the same Zak-OTFS subframe (S|C), and as PDR increases the pilot (e.g., spread pilot) overwhelms the combination of data interference and noise.

[0506] NMSE saturates at high PDR because of DD domain aliasing.

[0507] FIG.27 is a graph showing example interference from a residual spread pilot. Here, theinterference is measured from the residual spread pilot. In this example, Data SNR ^^ௗ ൌ 25 dB and^^^^௫ ൌ 815 Hz for a standard Veh-A channel. And, in this example, SIR represents the ratio of data powerto residual spread pilot power (after cancellation). When the effective channel is estimated, it is estimated within a support region – taps outside this region are noisy and subject to DD domain aliasing. That is, not all effective channel taps are estimated: channel taps within a support region are read off as taps outside this region are subject to DD domain aliasing and are noisy.

[0509] Sensing accuracy can improve with increasing PDR but so does the power of the interference from the residual pilot. That is, sensing accuracy can improve with PDR, but residual spread pilot energy increases proportionately.

[0510] FIG.28 shows a graph of an example bit error rate (BER) performance. Here, BER performance of uncoded 4-QAM is plotted as a function of the ratio of pilot energy to data energy (PDR).In this example, Data SNR ^^ௗ ൌ 25 dB and ^^^^௫ ൌ 815 Hz for a standard Veh-A channel.

[0511] When CSI is perfect, so is subtraction of the received pilot, and BER performance is independent of PDR [e.g., C | ^^, perfect CSI: BER performance independent of PDR].

[0512] The dashed curve with the light gray triangle markers [middle curve] shows the baseline - separate sensing and communications (S|C bar & C|S bar) where performance tracks sensing accuracy (NMSE) [e.g., S | C & C | S : BER performance tracks sensing accuracy (NMSE)]. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0513] The dashed curve with the black triangle markers [top curve] is sensing and communications in the same Zak-OTFS subframe (S|C & C|S). Here, sensing accuracy (NMSE) improves as PDR increases, but so does the interference from the residual pilot pulse-tone after cancellation. The result is a “U” shaped curve with an optimal PDR [e.g., S | C & C | S : Characteristic ``U” shape and optimal PDR].

[0514] In example FIG.28, when PDR < 10 dB, BER improves as PDR increases. When, PDR > 10 dB, interference from residual spread pilot degrades BER as PDR increases.

[0515] FIG.29 shows a graph of an example BER performance as a function of maximum Doppler spread. Here, BER performance of uncoded 4-QAM is plotted as a function of the maximumDoppler spread. In this example, Data SNR ^^ௗ ൌ 25 dB and PDR =10 db for a standard Veh-A channel.

[0516] The gray curve [bottom curve] shows the baseline - separate sensing and communications (S|C bar & C|S bar) where there is no residual pilot pulse-tone to interfere with data transmission.

[0517] The black curve (black triangles) is sensing and communications in the same Zak-OTFS subframe (S|C & C|S) – the residual pilot interferes with data transmission but the variation in BER is small over a wide range of Doppler spreads.

[0518] This is very different from the performance of point pilots as the variation is much greater for point pilot pulse-tones.

[0519] Spread pilots can extend the range of reliable operation to a wider range of Doppler spreads.

[0520] FIG.30 is a graph showing examples of geometric criteria for filter design. In thisexample, Data SNR ^^ௗ ൌ 25 dB and PDR =10 db for a standard Veh-A channel.

[0521] The slope q of the chirp filter that defines the spread filter is very important. Here, slopes q=3 and q=36 are considered and BER performance is compared.

[0522] Reliable data detection requires crystallization with respect to the lattice ^p, and that isthe case in the example of FIG.30 – Independent of slope ^^ : ^^^^௫ ^ ^^^ ൌ 30 KHz , ^^^^௫ ൌ 2.5 ^^^^ ^ ^^^

[0523] Accurate sensing requires crystallization with respect to the lattice ^q. In the example of FIG.30, it is satisfied for q=3, but not for q=36 when ^max> 2KHz.

[0524] FIG.31 is a graph showing example performance of using spread pilots. In this example,Data SNR ^^ௗ ൌ 25 dB and PDR =10 db for a standard Veh-A channel.

[0525] Effective throughput is the ratio of the number bits reliably communicated in each subframe to the available degrees of freedom.

[0526] In various embodiments, only sensing and communications in the same Zak-OTFS subframe (S|C & C|S) is considered. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0527] DD domain resources are divided between sensing and communications when point pilots are used and DD domain resources are shared when spread pilots are used. It will be appreciated by skilled artisans that sharing increases effective throughput. For example, sharing DD domain resources can increase effective throughput compared with traditional approaches that use guard bands to divide DD domain resources between sensing and communications.

[0528] Hence, in various embodiments, the use of spread pilots can help maximize effective throughput.

[0529] FIG.32 is a graph showing example results obtained in a turbo-aided communication andsensing implementation. In this example, Data SNR ^^ௗ ൌ 25 dB and PDR =10 db for a standard Veh-Achannel.

[0530] In various embodiments, the effective channel is estimated, and then the received spread pilot is estimated. Then the data can be recovered after subtracting the estimate for the received pilot from the received signal.

[0531] Now take the estimated data, then estimate the received data signal, and then improve the estimate for the effective channel by subtracting the estimate for the received data signal from the received signal.

[0532] This can then be repeated as many times as is useful. This is turbo aided sensing and communication; and, in various embodiments, two turbo iterations suffice to match the performance of separate sensing and communication across a wide range of Doppler shifts.

[0533] In some embodiments, subtracting the received spread pilot from the received signal can include computing the twisted convolution of the received spread pilot and a result of the channel sensing.

[0534] 8. Examples and implementations of the disclosed technology

[0535] FIG.33 is a block diagram representation of a hardware platform 3300 which may be used to implement the various methods described in the present document. The hardware platform 3300 may be incorporated within a base station or a user device. The hardware platform 3300 includes at least one processor 3302, an optional memory 3304 (this may be optional because in some cases the memory may be internal to the processor) and a transceiver circuitry 3306. The processor(s) may execute instructions, e. g., by reading from the memory 3304, and control the operation of the transceiver circuitry 3306 and the hardware platform 3300 to perform the methods described herein. The transceiver circuitry 3306 may be used, for example, for transmitting or receiving signals, as disclosed herein. In some embodiments, the memory 3304 and / or the transceiver circuitry 3306 may be partially or completely contained within the processor 3302 (e.g., same semiconductor package). The hardware platform 3300 may further include one or more antennas to perform multi-layer communication as disclosed herein. The hardware platform 3300 may be used to implement the transmitter apparatus or the receiver apparatus described herein. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0536] In some embodiments, method of digital communication includes generating, by a transmitter apparatus, a transmission waveform comprising a data portion and a pilot portion and operating a transceiver of the transmitter apparatus to transmit the transmission waveform over the channel. The pilot portion comprises a single pulse-tone pilot signal that allows a receiver to estimate a channel on which the transmission waveform is transmitted. Additional details of such a method are provided with reference to the point-pulse pilot technology disclosed herein.

[0537] In some embodiments, a digital communication method includes receiving, by a receiver apparatus, over a channel, a received signal corresponding to a transmission waveform transmitted over the channel, the transmission waveform comprising a data portion and a pilot portion wherein the pilot portion comprises a single pulse-tone pilot signal; and estimating the channel based on a cross-ambiguity function between the received signal and the single pulse-tone pilot signal.

[0538] In the disclosed embodiments, the single pulse-tone pilot signal comprises a time domain signal that is represented in a delay-Doppler domain as a quasi-periodic pulse, wherein the delay-Doppler domain is defined by a delay period ^^^and the Doppler period ^^^.

[0539] The above-disclosed embodiments may use a point pulse pilot, as disclosed in the present document.

[0540] In some embodiments, a method of data communications, implemented by a transmitter apparatus, includes generating a first intermediate signal by combining information symbols in a delay- Doppler domain with a quasi-periodic impulse train in the Delay-Doppler domain; generating a second intermediate signal by performing a twisted convolution of the first intermediate signal with a transmit pulse shaping filter; generating a transmission waveform by applying inverse Zak transform to the second intermediate signal; and causing a transmitter electronics to transmit the transmission waveform over a channel. Additional features of such a method are described with respect to the Zak OTFS implementations disclosed in the present document.

[0541] In some embodiments, a method of digital communications includes, receiving, by a receiver apparatus from a transmitter apparatus a transmission waveform over a channel, generating a first intermediate signal by applying a Zak transform to the transmission waveform, generating a second intermediate signal by performing a twisted convolution between the first intermediate signal and a receiver pulse shaping filter, recovering information signal bits from the second intermediate signal.

[0542] In the above-disclosed methods, the twisted convolution may be performed in the delay- Doppler domain. Some examples are provided with reference to, e.g., FIG.7.

[0543] In some embodiments, a digital communication method includes, generating, by a transmitter apparatus, a transmission waveform comprising a data portion and / or a pilot portion in a delay-Doppler domain; wherein the pilot portion comprises a pulse-tone pilot signal; and operating a transceiver to transmit the transmission waveform over the channel. Additional features of this method are described with respect to ISAC using point pulse-tone waveforms. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0544] In some embodiments, a method of digital communication includes, receiving, by a receiver apparatus, a received signal representing a transmission waveform transmitted over a channel, wherein the transmission waveform comprises a data portion and / or a pilot portion that are combined together in a delay-Doppler domain, wherein the pilot portion comprises a pulse-tone pilot signal; determining an estimate of the channel using the pilot portion, wherein the determining is performed based on a channel sensing operation, estimating the data portion by subtracting, from the received signal, an estimate of the pilot portion determined from the channel sensing operation, and recovering information bits from the data portion. Additional features of this method are described with respect to ISAC using point pulse-tone waveforms.

[0545] In the above-disclosed methods, the channel sensing operation may be performed by calculating a cross-ambiguity function between the received signal and the pulse-tone pilot signal.

[0546] In some embodiments, the data portion and the pilot portion are separated by a guard band in the delay-Doppler domain.

[0547] In some embodiments, no guard band is used for separating the data portion and the pilot portion in the delay-Doppler domain.

[0548] In some embodiments, the transmission waveform is organized in subframes comprising one or more of: (1) a subframe for sensing in which the pilot portion is transmitted, (2) a subframe for communication in which the data portion is transmitted, or (3) a subframe for joint sensing and communication in which both the data portion and the pilot portions are transmitted.

[0549] The above-described methods may be used for joint sensing and communication in which transmissions are organized as either separate frames for sensing and communications, or joint frames for sensing and communication. Here, a point pilot may be used.

[0550] In some embodiments, a digital communication method includes generating, by a transmitter, a transmission waveform comprising a data portion and / or a pilot portion in a delay-Doppler domain and operating a transceiver to transmit the transmission waveform over the channel. In some embodiments, the pilot portion comprises a spread pulse-tone pilot signal.

[0551] In some embodiments, a digital communication method includes receiving, by a receiver apparatus, a received signal representing a transmission waveform transmitted over a channel, wherein the transmission waveform comprises a data portion and / or a pilot portion that are combined together in a delay-Doppler domain, wherein the pilot portion comprises a spread pulse-tone pilot signal, determining an estimate of the channel using the pilot portion, wherein the determining is performed based on a channel sensing operation, estimating the data portion by subtracting, from the received signal, an estimate of the pilot portion determined from the channel sensing operation, and recovering information bits from the data portion.

[0552] In the above-disclosed methods, the channel sensing operation may be performed by calculating a cross-ambiguity function between the received signal and the pulse-tone pilot signal. In 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 some embodiments, no guard band is used for separating the data portion and the pilot portion in the delay-Doppler domain. In various embodiments, the transmission waveform is organized in subframes comprising one or more of: (1) a subframe for sensing in which the pilot portion is transmitted, (2) a subframe for communication in which the data portion is transmitted, or (3) a subframe for joint sensing and communication in which both the data portion and the pilot portions are transmitted.

[0553] In some embodiments, the spread pulse-tone pilot signal is generated by applying a chirp filter to a pulse-tone signal. In some embodiments, the spread pulse-tone pilot signal is quasi periodic in the delay-Doppler domain.

[0554] In some embodiments, a digital communication method can include receiving, by a receiver apparatus, a received signal comprising multiple layers, each layer comprising corresponding data and a corresponding spread pilot over a channel, wherein, for each layer, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain; and extracting data from at least some of the layers of the received waveform by subtracting, for each extracted layer, the corresponding spread pilot from the received signal and by performing a channel equalization based on a channel sensing using the corresponding spread pilot. In some embodiments, the extracting can include performing the subtracting, the channel equalization and the channel sensing over multiple iterations. In some embodiments, subtracting the corresponding spread pilot signal from the received signal can include computing a twisted convolution of the corresponding spread pilot signal and a result of the channel sensing.

[0555] 9. Example technical solutions

[0556] The technical solutions disclosed in the present document provide a discrete filter that is periodic in M*N on both delay and Doppler dimension. One can think of the M x N basic delay-Doppler grid, as being extended N times along delay and M times along Doppler, creating an (M*N) x (N*M) grid. The filter is defined over this extended grid, which is then extended periodically to have infinite number of (M*N) x (N*M) duplicates along both axes.

[0557] For the waveforms defined in the present document, the delay-Doppler grid is infinite, but can be described using a basic grid of M x N, where signals that are quasi-periodic. This means that these signals are periodic on the Doppler dimension (every N grid elements, the signal repeats) and non- periodic on the delay dimension (repeat every M grid elements but with some phase rotation).

[0558] The present document also defines a filter in delay-Doppler. The filter is periodic with periodicity of M*N on both dimensions. This means that it repeats itself after M*N delay-Doppler grid elements.

[0559] Applying a filter to a signal (such as a delta), will result in a filtered signal that is also M*N periodic on both dimensions.

[0560] The following technical solutions may be implemented by some preferred embodiments. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0561] With respect to the use of a spread pilot signal, the following solutions may be used.

[0562] 1. A method of digital communication (e.g., method 3410 depicted in FIG.34A), comprising: generating (3402), by a transmitter apparatus, a transmission waveform by applying a discrete filter to an input signal in a delay-Doppler domain, wherein the discrete filter is periodic in delay domain and in Doppler domain; and controlling (3404) a transceiver to transmit the transmission waveform over a channel. For example, one or more processors of the transmitter apparatus may control the transceiver to cause a transmission to occur.

[0563] 2. The method of solution 1, wherein the discrete filter is M*N periodic in the delay- Doppler domain, where M and N are positive integers. In this document M*N represents a multiplication of numbers M and N.

[0564] 3. The method of solution 2, wherein M and N represent delay and Doppler dimensions of a delay-Doppler grid in the delay-Doppler domain, corresponding to the delay period and the Doppler period.

[0565] 4. The method of solutions 1-3, wherein the input signal comprises a quasi-periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N-periodic twisted convolution.

[0566] 5. The method of solution 4, wherein an intermediate signal resulting from the applying is spread over the delay-Doppler domain as an M*N-periodic signal in the delay-Doppler domain.

[0567] 6. The method of solution 5, wherein the transmission waveform is generated by applying an inverse Zak transform to the intermediate signal.

[0568] 7. The method of any of solutions 1-6, wherein the discrete filter comprises a chirp signal in the delay-Doppler domain.

[0569] 8. The method of solution 7, wherein the chirp signal is defined by a slope parameter q, and wherein the slope parameter q is configured to avoid aliasing in the delay-Doppler domain after the transmission waveform is transmitted over the channel.

[0570] Some solutions, e.g., as listed below, may be used to perform a multi-layer communication.

[0571] 9. A method of digital communication (e.g., method 3420 depicted in FIG.34B), comprising: generating (3422) a transmission waveform comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain, and wherein L is a positive integer; and controlling (3424) a transceiver to transmit the transmission waveform over a channel.

[0572] 10. The method of solution 9, wherein the data and the spread pilot are overlapping in the delay-Doppler domain. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0573] 11. The method of solution 9, wherein the data occupies entirety of the delay-Doppler domain, and the spread pilot occupies entirety of the delay-Doppler domain.

[0574] 12. The method of any of solutions 9-11, wherein the controlling comprises weighting transmission energies of the data and the spread pilot.

[0575] 13. The method of any of solutions 9-12, wherein the transmission waveform is configured for an integrated sensing and communication operation.

[0576] 14. The method of any of solutions 9-13, wherein L is greater than 1, and wherein the each of the L layers uses a corresponding different spread pilot that is generated using a different discrete filter.

[0577] 15. The method of any of solutions 9-14, wherein the filtered pulse-tone signal is generated by applying a discrete filter to an input signal in the delay-Doppler domain, wherein the input signal comprises a quasi-periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N-periodic twisted convolution, where M and N are positive integers representing delay and Doppler domain dimensions of a basic delay-Doppler grid.

[0578] The following solutions may be implemented by embodiments that use multiple antennas.

[0579] 16. A method of wireless communication (e.g., method 3430 depicted in FIG.34C), comprising: generating (3432) K transmission waveforms, each of the K transmission waveforms comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L layers, the corresponding data comprises symbols in a delay- Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain, and wherein K and L are positive integers and K > L; and transmitting (3434), the each of the K transmission waveforms using a corresponding transmission antenna.

[0580] 17. The method of solution 16, wherein spread pilots used for at least some of the K transmission waveforms are different from each other.

[0581] 18. The method of any of solutions 16-17, wherein each spread pilot is obtained by applying a corresponding discrete filter to a corresponding input signal in a delay-Doppler domain, wherein the discrete filter is M*N periodic in the delay-Doppler domain, where M and N are positive integers.

[0582] 19. The method of solution 18, wherein M and N represent delay and Doppler dimensions of a delay-Doppler grid along the delay-Doppler domain.

[0583] 20. The method of solution 19, wherein the corresponding input signal comprises a quasi- periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N- periodic twisted convolution.

[0584] 21. The method of any of solutions 16-20, wherein each of the corresponding discrete filter comprises a corresponding chirp signal in the delay-Doppler domain. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0585] 22. The method of solution 21, wherein each of the corresponding chirp signal is defined by a corresponding slope parameter q, and wherein the corresponding slope parameter q is configured to avoid aliasing after the corresponding transmission waveform is transmitted over the channel.

[0586] The following solutions may be implemented by receiver apparatus.

[0587] 23. A digital communication method (e.g., method 3440 depicted in FIG.34D), comprising: receiving (3442) a received signal comprising L layers, each of the L layers comprising corresponding data and a corresponding spread pilot over a channel, wherein L>=1, and wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain; and extracting (3444) data from at least some of the L layers of the received waveform by subtracting, for each extracted layer, the corresponding spread pilot from the received signal and by performing a channel equalization based on a channel sensing using the corresponding spread pilot.

[0588] 24. The method of solution 23, wherein the channel sensing is performed by computing a cross-ambiguity function between the extracted layer of the received signal and a priori knowledge of the corresponding spread pilot.

[0589] 25. The method of any of solutions 23-24, wherein the subtracting the corresponding spread pilot comprises computing a twisted convolution of the corresponding spread pilot and a result of the channel sensing.

[0590] 26. The method of any of solutions 23-25, wherein the extracting comprises performing the subtracting, the channel equalization and the channel sensing over multiple iterations.

[0591] 27. The method of any of solutions 23-26, further including, performing an integrated sensing and communication operation.

[0592] 28. An apparatus for digital communications, comprising at least one processor and a transceiver, wherein the transceiver is configured to transmit or receive signals under control of the at least one processor, and the at least one processor is configured to implement a method recited in any of solutions 1 to 27.

[0593] 29. A non-transitory computer readable medium having code stored thereon; the code, upon execution by one or more processors, causing the one or more processors to implement a method recited in any of solutions 1 to 27.

[0594] 10. Conclusion

[0595] It will be appreciated that the present document discloses techniques for Zak-OTFS pulse-tone shaping in delay-Doppler (Delay / Doppler) domain and Zak-OTFS pulse-tone spreading in delay-Doppler domain and heterogenous Zak-OTFS shaped pulse-tone in delay-Doppler (different pulse shaping) for Pilots & Data. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00

[0596] It will also be appreciated that the present document discloses a method of separating Heterogenous Zak-OTFS shaped pulse-tone in delay-Doppler, channel sensing based on Zak-OTFS shaped pulse-tone in delay-Doppler.

[0597] The disclosed techniques may be used to implement a passive radar based on Zak- OTFS shaped pulse-tone in delay-Doppler Pilots or a radar based on delay-Doppler Pulse shaping.

[0598] It will also be appreciated that the present document discloses Turbo receivers for Zak- OTFS.

[0599] It will be appreciated by one of skill in the art that Integrated sensing and communications reduces to geometric properties of a lattice Λp used for data transmission and a lattice Λq used for sensing.

[0600] By sharing DD domain resources between sensing andit is possible to optimize throughput without compromising BER performance achieved by separating sensing and communication.

[0601] When the channel satisfies the crystallization condition with respect to the lattice Λq the effective DD filter taps can be read off from the response to a single spread pilot.

[0602] When the channel satisfies the crystallization condition with respect to the lattice Λp, then given the I / O response at one point in a Zak-OTFS subframe it is possible to predict the response at all points in the subframe.

[0603] Filters in the discrete DD domain enable integration by minimizing interference between sensing and data transmission – the data pulse-tones look like noise to the sensing pulse-tone.

[0604] Filters in the discrete DD domain can be used to design noise-like waveforms with excellent PAPR.

[0605] The disclosed and other embodiments, modules and the functional operations described in this document can be implemented in digital electronic circuitry, or in computer software, firmware, or hardware, including the structures disclosed in this document and their structural equivalents, or in combinations of one or more of them. The disclosed and other embodiments can be implemented as one or more computer program products, i.e., one or more modules of computer program instructions encoded on a computer readable medium for execution by, or to control the operation of, data processing apparatus. The computer readable medium can be a machine-readable storage device, a machine- readable storage substrate, a memory device, a composition of matter effecting a machine-readable propagated signal, or a combination of one or more of them. The term “data processing apparatus” encompasses all apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. The apparatus can include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 operating system, or a combination of one or more of them. A propagated signal is an artificially generated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, which is generated to encode information for transmission to suitable receiver apparatus.

[0606] A computer program (also known as a program, software, software application, script, or code) can be written in any form of programming language, including compiled or interpreted languages, and it can be deployed in any form, including as a standalone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program does not necessarily correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.

[0607] The processes and logic flows described in this document can be performed by one or more programmable processors executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).

[0608] Processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer. Generally, a processor will receive instructions and data from a read -only memory or a random-access memory or both. The essential elements of a computer are a processor for performing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks. However, a computer need not have such devices. Computer readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media, and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto optical disks; and CD ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.

[0609] While this patent document contains many specifics, these should not be construed as limitations on the scope of an invention that is claimed or of what may be claimed, but rather as descriptions of features specific to particular embodiments. Certain features that are described in this document in the context of separate embodiments can also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment can also be implemented in multiple embodiments separately or in any suitable sub-combination. Moreover, 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a sub-combination or a variation of a sub- combination. Similarly, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results.

[0610] Only a few examples and implementations are disclosed. Variations, modifications, and enhancements to the described examples and implementations and other implementations can be made based on what is disclosed. 180520639.2

Claims

PCT Patent Application Attorney Docket No.: 119314.8123.WO00 CLAIMS 1. A method of digital communication, comprising: generating, by a transmitter apparatus, a transmission waveform by applying a discrete filter to an input signal in a delay-Doppler domain, wherein the discrete filter is periodic in delay domain and in Doppler domain; and controlling a transceiver to transmit the transmission waveform over a channel.

2. The method of claim 1, wherein the discrete filter is M*N periodic in the delay-Doppler domain, where M and N are positive integers.

3. The method of claim 2, wherein M and N represent delay and Doppler dimensions of a delay-Doppler grid in the delay-Doppler domain, corresponding to the delay period and the Doppler period.

4. The method of claims 1-3, wherein the input signal comprises a quasi-periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N-periodic twisted convolution.

5. The method of claim 4, wherein an intermediate signal resulting from the applying is spread over the delay-Doppler domain as an M*N-periodic signal in the delay-Doppler domain.

6. The method of claim 5, wherein the transmission waveform is generated by applying an inverse Zak transform to the intermediate signal.

7. The method of claim 1, wherein the discrete filter comprises a chirp signal in the delay-Doppler domain.

8. The method of claim 7, wherein the chirp signal is defined by a slope parameter q, and wherein the slope parameter q is configured to avoid aliasing in the delay-Doppler domain after the transmission waveform is transmitted over the channel.

9. A method of digital communication, comprising: generating a transmission waveform comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L Layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain, and wherein L is a positive integer; and controlling a transceiver to transmit the transmission waveform over a channel. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 10. The method of claim 9, wherein the data and the spread pilot are overlapping in the delay-Doppler domain.

11. The method of claim 9, wherein the data occupies entirety of the delay-Doppler domain and the spread pilot occupies entirety of the delay-Doppler domain.

12. The method of any of claims 9-11, wherein the controlling comprises weighting transmission energies of the data and the spread pilot.

13. The method of claim 12, wherein the transmission waveform is configured for an integrated sensing and communication operation.

14. The method of claim 13, wherein L is greater than 1, and wherein the each of the L layers uses a corresponding different spread pilot that is generated using a different discrete filter.

15. The method of claim 14, wherein the filtered pulse-tone signal is generated by applying a discrete filter to an input signal in the delay-Doppler domain, wherein the input signal comprises a quasi-periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N-periodic twisted convolution, where M and N are positive integers representing delay and Doppler domain dimensions of a basic delay-Doppler grid.

16. A method of wireless communication, comprising: generating K transmission waveforms, each of the K transmission waveforms comprising L layers, wherein each of the L layers comprises a corresponding data and a corresponding spread pilot, wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi- periodic delta signal over the delay-Doppler domain, and wherein K and L are positive integers and K > L; and transmitting, the each of the K transmission waveforms using a corresponding transmission antenna.

17. The method of claim 16, wherein spread pilots used for at least some of the K transmission waveforms are different from each other.

18. The method of any of claims 16-17, wherein each spread pilot is obtained by applying a corresponding discrete filter to a corresponding input signal in a delay-Doppler domain, wherein the discrete filter is M*N periodic in the delay-Doppler domain, where M and N are positive integers. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 19. The method of claim 18, wherein M and N represent delay and Doppler dimensions of a delay- Doppler grid along the delay-Doppler domain.

20. The method of claim 19, wherein the corresponding input signal comprises a quasi-periodic impulse signal in the delay-Doppler domain, and wherein the applying comprises using an M*N-periodic twisted convolution.

21. The method of claim 20, wherein each of the corresponding discrete filter comprises a corresponding chirp signal in the delay-Doppler domain.

22. The method of claim 21, wherein each of the corresponding chirp signal is defined by a corresponding slope parameter q, and wherein the corresponding slope parameter q is configured to avoid aliasing after the corresponding transmission waveform is transmitted over a channel.

23. A digital communication method, comprising: receiving a received signal comprising L layers, each of the L layers comprising corresponding data and a corresponding spread pilot over a channel, wherein L>=1, and wherein, for the each of the L layers, the corresponding data comprises symbols in a delay-Doppler domain and the corresponding spread pilot comprises a filtered pulse-tone signal that is a filtered quasi-periodic delta signal over the delay-Doppler domain; and extracting data from at least some of the L layers of the received waveform by subtracting, for each extracted layer, the corresponding spread pilot from the received signal and by performing a channel equalization based on a channel sensing using the corresponding spread pilot.

24. The method of claim 23, wherein the channel sensing is performed by computing a cross-ambiguity function between the extracted layer of the received signal and a priori knowledge of the corresponding spread pilot.

25. The method of any of claims 23-24, wherein the subtracting the corresponding spread pilot comprises computing a twisted convolution of the corresponding spread pilot and a result of the channel sensing.

26. The method of claim 23, wherein the extracting comprises performing the subtracting, the channel equalization and the channel sensing over multiple iterations.

27. The method of claim 26, further including, performing an integrated sensing and communication operation. 180520639.2PCT Patent Application Attorney Docket No.: 119314.8123.WO00 28. An apparatus for digital communications, comprising at least one processor and a transceiver, wherein the transceiver is configured to transmit or receive signals under control of the at least one processor, and the at least one processor is configured to implement a method recited in any of claims 1 to 27.

29. A non-transitory computer readable medium having code stored thereon; the code, upon execution by one or more processors, causing the one or more processors to implement a method recited in any of claims 1 to 27. 180520639.2

Images