## Doppler Ambiguity Elimination Using 5G Signals in Integrated Sensing and Communication

Keivan Khosroshahi ∗† , Philippe Sehier , and Sami Mekki , † † ∗ Universit´ e Paris-Saclay, CNRS, CentraleSup´ elec, Laboratoire des Signaux et Syst` emes, Gif-sur-Yvette, France † Nokia Standards, Massy, France keivan.khosroshahi@centralesupelec.fr, { philippe.sehier, sami.mekki } @nokia.com

Abstract -The industrial point of view towards integrated sensing and communication (ISAC), the preference is to leverage existing resources and fifth-generation (5G) infrastructure to minimize deployment costs and complexity. In this context, we explore the utilization of current 5G new radio (NR) signals aligned with 3rd generation partnership project (3GPP) standards. Positioning reference signals (PRS) for sensing and physical downlink shared channel (PDSCH) for communication have been chosen to form an ISAC framework. However, PRSbased sensing suffers from Doppler ambiguity when the Doppler frequency shift is severe. To address this challenge, we introduce a novel method within the ISAC system that leverages the demodulation reference signal (DMRS) present in PDSCH to eliminate Doppler ambiguity. Furthermore, we formulate a resource allocation problem between PRS and PDSCH to achieve a Pareto optimal point between communication and sensing without Doppler ambiguity. Through simulations and analysis, we demonstrate the effectiveness of our proposed method on joint DMRS-PRS exploitation in mitigating Doppler ambiguity and the efficiency of the resource allocation scheme in achieving Pareto optimality for ISAC within a 5G NR framework.

Index Terms -PRS - PDSCH - DMRS - ISAC - Doppler ambiguity

## I. INTRODUCTION

networks, the use of the positioning reference signal (PRS) for sensing purposes holds significant promise, given its rich time-frequency resources and flexible configuration. PRS, introduced in 3rd generation partnership project (3GPP) Release 16 of the 5G specification, aims to improve the positioning accuracy of connected user equipments (UEs) due to its high resource element (RE) density and superior correlation properties compared to existing reference signals due to the diagonal or staggered PRS RE pattern. In sensing context, PRS is used in [2] to estimate the range and Doppler, and its performance is compared with other pilots available in 5G standard, namely demodulation reference signal (DMRS), channel state information reference signal (CSI-RS) and synchronization signal (SS). A two-stage joint pilot optimization, target detection, and channel estimation scheme has been introduced in [4] for ISAC. In [5], DMRS and CSI-RS are separately used for velocity and range estimation. In [6], a multi-range joint automotive radar and communication system based on pilotbased orthogonal frequency-division multiplexing (OFDM) waveform is investigated. PRS and DMRS are jointly used in [7] to remove range ambiguity in ISAC scenarios.

The adoption of integrated sensing and communication (ISAC) technology in upcoming mobile networks is increasingly gaining consensus due to the emergence of future applications in the domains of industry 4.0, automotive, internet of things (IoT) in smart cities and health care, which impose new requirements on wireless communication networks [1]. Integrating sensing services into the mobile infrastructure enables new services to be deployed at low cost, reusing the mobile coverage base already deployed and without requiring additional spectrum resources, while minimizing the impact on existing equipment.

Several reference signals have been defined in fifthgeneration (5G) new radio (NR), addressing different objectives. They have been specified to provide excellent passive detection performance, robust anti-noise capabilities, and favorable auto-correlation characteristics [2]. Moreover, using pilots in ISAC systems leads to low hardware costs, easy implementation, and low computational complexity [3]. Consequently, significant interest has emerged in designing ISAC signals based on pilots. Among the pilots available in 5G

This work was funded by the European Commission through the H2020 MSCA 5GSmartFact project under grant agreement number 956670.

However, exploiting the existing reference signals of the 5G-NR raises new issues. This is because the sensed targets might be far from the gNB, or their speed may be high. Consequently, the reflected signals may be received with a large time delay, above the cyclic prefix length, and large Doppler shifts. Periodogram-based estimation for delay and Doppler exhibit spurious peaks in these situations, creating ambiguities in distance or velocity estimates. Such fake peaks can be interpreted as real targets. While the proposed solution in [8] has focused on multiple reference signals to address the ambiguity, they are not directly applicable in ISAC scenarios due to potential interference between physical downlink shared channel (PDSCH) for data transmission and other reference signals such as PRS and CSI-RS. Additionally, the method described in [2] leads to a significant reduction in the maximum detection range, especially when employing large PRS comb sizes. Hence, to the best of our knowledge, a 3GPPcompliant solution to address range ambiguity within the ISAC framework has not been proposed.

In this work, we introduce a novel approach leveraging the current reference signals in 5G NR to enable sensing for unconnected targets with enhanced maximum detectable Doppler frequency shift without Doppler ambiguity. To provide sensing services alongside data transmission, we propose an ISAC OFDM resource grid comprising PDSCH for communication and PRS for sensing. To remove the Doppler ambiguity, we propose a method to reuse DMRS within PDSCH for sensing in addition to communication. Furthermore, we propose time and frequency resource allocation between PRS and PDSCH since sensing and communication symbols cannot share the same REs in an OFDM resource grid. Building upon this, we formulate a resource allocation problem between PDSCH and PRS to achieve Pareto optimality between communication and sensing without Doppler ambiguity. Ultimately, via extensive simulations, we prove the efficacy of our proposed method in addressing Doppler ambiguity and finding the Pareto optimal point between communication and sensing.

## II. PRS AND DMRS DESCRIPTION

The 5G-NR waveform is built on a time-frequency resource grid as described in technical specifications (TS) 38.214 [9] and the smallest time-frequency resource is a physical resource block (PRB), consisting of 12 contiguous sub-carriers and 14 consecutive symbols. In the following two sections, we explain about PRS and DMRS in 5G NR.

## A. PRS

The generation of PRS as explained in TS 38.211 [10] as follows

$$r ( m ) = \frac { 1 } { \sqrt { 2 } } ( 1 - 2 c ( 2 m ) ) + j \frac { 1 } { \sqrt { 2 } } ( 1 - 2 c ( 2 m + 1 ) ), \quad ( 1 ) \quad \text{(sq. for} } _ { \text{signals} }.$$

where c i ( ) is denoted as Gold sequence of length31 and the starting value of c i ( ) for the PRS is provided in [10]. PRS is generated from the Gold sequence and it has favorable antinoise and auto-correlation characteristics. The PRS is allocated a minimum of 24 PRBs and a maximum of 272 PRBs, depending on the carrier size, showing the flexible transmission parameters supported in 5G NR. This flexibility allows for a versatile configuration of time-frequency resources to fulfill sensing accuracy requirements across diverse applications and scenarios. As outlined in TS 38.211 [10], PRS supports four comb structures, i.e., Comb { 2 4 6 12 , , , } , in the frequency domain and five symbol number configurations, i.e., Symbol { 1 2 4 6 12 , , , , } , in the time domain. The structure of PRS with 12 symbols and two comb sizes is illustrated in Fig. 1.

## B. DMRS

Another reference signal that can potentially be used in the ISAC framework is DMRS, available in the PDSCH. DMRS replaces the cell-specific reference signal (CRS) in long-term evolution (LTE) for channel estimation, and PDSCH decoding [10]. DMRS generation is also based on (1), and the initial value of c i ( ) for DMRS can be found in [10]. DMRS generation occurs within the designated PDSCH allocation, as detailed in TS 38.211 [10]. The allocated resources for PDSCH are available in the bandwidth part (BWP) of the carrier as explained in TS 38.214 [9]. DMRS is used for channel estimation and is confined to the resource blocks (RBs) assigned for PDSCH. Its structure is adaptable to

Fig. 1: PRS allocation with starting symbol= 2 and different Comb sizes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-0bab8096-821e-4f1f-a073-3eaedfad03b4.jpg)

various deployment scenarios and use cases. The positioning of DMRS symbols depends on the mapping type, which can be slot-wise (Type A) or non-slot-wise (Type B). Additional DMRS symbol positions are determined by tables outlined in TS 38.211 [10]. Furthermore, front-load DMRS may occupy 1 to 4 OFDM symbols in the time domain. DMRS is within physical channels allocated for communication receivers and is transmitted using the same beam selected for communication receivers.

## III. ISAC SYSTEM MODEL

The existing 5G NR signals were not initially designed for radar sensing applications. This section outlines a methodology for efficiently leveraging both communication and sensing signals within the ISAC framework. Our vision involves implementing a downlink ISAC system where a gNB directs signals toward a sensing area where K point-like targets exist. The receiver captures echoes from these targets to estimate their range and velocity. We assume that the targets are within the line-of-sight (LOS) of both the transmitter and the receiver. Fig. 2 depicts an example of the considered scenario. If we consider an OFDM resource grid with N symbols in the time domain and M sub-carriers in the frequency domain, the transmitted signal in the time domain can be written in the following form: [11]

$$\underset { S \text{ with} } { \text{fucny} } \underset { n = 0 } { \text{$s(t)=\sum_{n=0}^{N-1}$ rect} } ( \frac { t - n T _ { 0 } } { T _ { 0 } } ) \sum _ { m = 0 } ^ { M - 1 } v ( m, n ) e ^ { j 2 \pi m \Delta f ( t - n T _ { 0 } ) }, \ ( 2 )$$

where rect ( t/T 0 ) denotes the rectangular pulse, ∆ f represents the subcarrier spacing. T 0 = T CP + T s is the total duration of the OFDM symbol, where T s = 1 ∆ f represents the symbol duration, and T CP denotes the duration of the cyclic prefix (CP). v m,n ( ) is the complex transmitted symbol at OFDM symbol n and sub-carrier M , where n = 0 , ..., N -1 and m = 0 , ..., M -1 within an M × N OFDM resource grid. If there are k targets present in the sensing area, the echo signals from targets received by the receiver can be written as [12]

$$\psi _ { \ } r ( t ) = \sum _ { k = 1 } ^ { K } \beta _ { k } s ( t - \tau _ { k } ) e ^ { j 2 \pi f _ { d, k } t } + u ( t ), \quad \ \ ( 3 )$$

where β k denotes the attenuation factor of the k -th target, f d,k represents the Doppler frequency shift of the k -th target, τ k

Fig. 2: ISAC system model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-a203e993-8d5a-43f7-a779-054fe802c62a.jpg)

signifies the delay of the k -th target and u t ( ) ∈ C refers to complex additive white Gaussian noise (AWGN) with zero mean and variance of 2 σ 2 .

Conventional OFDM receivers proceed by sampling the signal at each symbol time, followed by performing fast fourier transform (FFT) to extract the modulated symbols. Under normal circumstances, delay and Doppler effects are compensated, leading to mitigation of intersymbol interference (ISI). The cyclic prefix is typically leveraged to address receive timing errors and multipath spread. However, in our scenario, the path delay and Doppler may not fall within the usual range due to significant uncertainties regarding the target's position and speed. As a result, it becomes necessary to consider ISI in both the time and frequency domains. Consequently, the received modulation symbols can be expressed as:

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \text{written} \\ y ( m, n ) = & \sum _ { k = 1 } ^ { K } \beta _ { k } e ^ { j 2 \pi n T _ { 0 } f _ { d, k } } e ^ { - j 2 \pi m \Delta f \tau _ { k } } v ( m, n ) + \text{ISI} & & \text{Acco} \\ & + q ( m, n ). & \text{(4)} & \text{achievá} \\ & \text{chown} \colon$$

The term q m,n ( ) ∈ C represents the complex AWGN on the m -th subcarrier and the n -th OFDM symbol. This noise has a zero mean and a variance of 2 σ 2 , and it is obtained from the sampling and FFT process over r t ( ) .

## IV. ESTIMATION OF THE RANGE AND VELOCITY

Doppler and delay would cause a rotation on the time and frequency axis, respectively over the received samples. Determining range and Doppler necessitates a two-dimensional search. Multiple techniques are outlined in [12]; however, we focus here on a computationally efficient method operating in two stages: initially estimating the range and subsequently the Doppler, after compensating for the delay. The bistatic distance, i.e., from the transmitter to the target and then from the target to the receiver, can be derived from range assessment using the periodogram [13]. To achieve this, we eliminate the transmitted sensing symbols from the received echoes through point-wise division as described below.

̸

$$g ( m, n ) = \begin{cases} \frac { y ( m, n ) } { v ( m, n ) }, & \text{If } v ( m, n ) \neq 0, \\ 0, & \text{If } v ( m, n ) = 0. \end{cases}$$

Afterward, we compute M -point inverse fast fourier transform (IFFT) of the PRS on the n -th column of the g m,n ( )

$$r _ { n } ( l ) = & | I \text{FFT} ( g ( m, n ) ) | = | \sum _ { k = 1 } ^ { K } ( \beta _ { k } e ^ { j 2 \pi n T _ { 0 } f _ { d, k } } \sum _ { m = 0 } ^ { M - 1 } \\ & e ^ { - j 2 \pi m \Delta f \frac { R ^ { k } t \circ t } { c } } e ^ { j 2 \pi \frac { m l } { M } } + \frac { q ^ { \prime } ( m, n ) } { v ( m, n ) } e ^ { j 2 \pi \frac { m l } { M } } ) |, \quad ( 6 ) \\ - \dots \, \pi \quad \alpha \quad \pi \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \beta \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \eta \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \alpha \quad \colon$$

where l = 0 , ..., M -1 , | | . is the absolute value, q ′ ( m,n ) = q m,n ( ) + ISI, and we replaced τ k = R tot k c while R tot k is the bistatic range, c is the light speed. The absolute value is taken to reduce Doppler sensitivity. When the arguments of e -j 2 πm ∆ f R tot k c e j 2 π ml M in (6) cancel each other, the maximum value occurs. The next step is to perform IFFT over all columns of the g m,n ( ) and average over them to increase the signal-to-noise ratio (SNR):

$$\bar { \mathbf r } ( l ) = \frac { 1 } { N } \sum _ { n = 0 } ^ { N - 1 } r _ { n } ( l ).$$

Afterward, we need to find the index of the maximum value named ˆ l k in (7). Then, the bistatic distance of each target can be calculated as [2]

$$\hat { \ } \hat { R } _ { k } ^ { t o t } = \frac { \hat { l } _ { k } c } { \Delta f M }.$$

The range resolution is

$$\Delta R = \frac { c } { \Delta f M },$$

and the expression for the maximum detection range can be written as

$$R _ { m a x } = \frac { c } { \Delta f }. \text{ \quad \quad } ( 1 0 )$$

According to [2], enhancing accuracy performance is achievable by increasing the number of IFFT points m a as shown below

$$\text{GN on} \quad & \text{shown below} \\ \text{d from} \quad & \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text1$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\texts$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\$$$

In this case, (8) will change into

$$\hat { R } _ { k } ^ { t o t } = \frac { \hat { l } _ { k } c } { m _ { a } \Delta f M }.$$

As highlighted in [14], increasing m a results in range accuracy enhancement, albeit at the expense of heightened computational complexity.

After compensating the delay caused by the targets, for Doppler estimation, we perform N -points FFT on the m -th row of the g m,n ( )

$$\intertext { \text{ate the } } \text{low on one $y(m,n)$} \\ \text{through} \\ & v _ { m } ( d ) = & | \text{FFT} ( g ( m, n ) ) | = | \sum _ { k = 1 } ^ { K } ( \beta _ { k } \sum _ { n = 0 } ^ { N - 1 } e ^ { j 2 \pi n T _ { 0 } f _ { d, k } } e ^ { - j 2 \pi \frac { n d } { N } } \\ & + \frac { q ^ { \prime } ( m, n ) } { v ( m, n ) } e ^ { - j 2 \pi \frac { n d } { N } } ) |,$$

where d = 0 , ..., N -1 . We perform FFT across all rows of the g m,n ( ) and then average the results as shown below:

$$\overline { \nu } ( d ) = \frac { 1 } { M } \sum _ { m = 0 } ^ { M - 1 } v _ { m } ( d ). \quad \quad ( 1 4 ) \quad \ \ 2 0$$

Afterward, we determine the index of the maximum value, denoted as ˆ d k for target k in (14), and the Doppler frequency can be estimated as [13]

$$\hat { \Phi } _ { \hat { d }, k } = \frac { \hat { d } _ { k } } { T _ { s } N }. \quad \quad \ \ ( 1 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The Doppler frequency resolution is

$$\hat { \Phi } _ { \Phi } \hat { f } = \frac { 1 } { T _ { s } N }.$$

and the maximum detectable Doppler frequency estimation is

$$\hat { f } _ { m a x } = \frac { 1 } { T _ { s } }. \quad \quad \ \ ( 1 7 ) \quad \ \ F i \varrho.$$

We can improve the accuracy of Doppler frequency estimation by increasing the FFT point similar to range estimation.

## A. Doppler Ambiguity

PRS was initially designed for positioning in 5G NR and not for sensing. In positioning, time-frequency uncertainties are mitigated by the time-frequency servo loop, which is not the case for sensing. PRS features a periodic structure that can be adjusted by the comb size, as demonstrated in Fig. 3, where an example with a comb size, i.e., K comb = 6 is presented. The presence of empty REs within the PRS structure leads to the emergence of fake peaks in Doppler estimation using the method described in section IV. These peaks introduce ambiguity in distinguishing between the real and fake Doppler frequency shift resulting in Doppler ambiguity as shown in Figs. 4a and 5a. The cost of the approach proposed in [2] to circumvent this ambiguity is the reduction of the maximum detectable range and Doppler frequency shift by a factor of K comb . In scenarios where precise Doppler resolution is crucial, reducing the subcarrier spacing, as defined in (16), becomes necessary. However, this adjustment with larger comb sizes leads to a notable reduction in the maximum detectable Doppler frequency shift. Hence, an alternative method is necessary to alleviate Doppler frequency ambiguity while ensuring a high maximum detectable Doppler frequency.

## B. Proposed Method

The key idea to eliminate ambiguity in Doppler estimation is to leverage DMRS as a reference signal within the PDSCH for ISAC scenarios where the OFDM frame contains PRS for sensing and PDSCH for communication. To this end, we generate an OFDM frame, and based on the sensing and communication requirements of each application, we perform the time and frequency resource allocation. Fig. 3 depicts an example of such resource allocation in the frequency domain between PRS and PDSCH. In this scenario, DMRS serves not only for communication tasks like channel estimation but also for enhancing sensing capabilities by accurately distinguishing

Fig. 3: Resource allocation between PRS and PDSCH showing a slot and 2 PRBs of the whole resource grid with K comb = 6 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-95d1138e-10f7-41c1-8e73-392fabaaf416.jpg)

between real and fake Doppler frequency shifts. Additionally, this improvement is achieved without compromising accuracy or the maximum detectable Doppler frequency. Additionally, it does not require any modifications to the gNB configuration, thus fully aligning with 3GPP standards.

Initially, we set double symbol for DMRS configuration and one additional DMRS symbol. Based on DMRS configuration type, the total number of subcarriers per PRB that can be allocated to DMRS, i.e., M ′ , is 4 or 6 . We define g DMRS ( m,n ) and g PRS ( m,n ) as the result of point-wise division similar to (5) over DMRS and PRS symbols, respectively. Then, we form the following

$$g _ { P R S, D M R S } ^ { m, m ^ { \prime } } ( n ) = g _ { P R S } ^ { m } ( n ) + g _ { D M R S } ^ { m ^ { \prime } } ( n ), \quad ( 1 8 )$$

where g m RS P ( n ) is the m -th row of g PRS ( m,n ) and g m MRS ′ D ( n ) is the m ′ -th row of g DMRS ( m,n ′ ) . Based on TS 38.211 [10], m ′ = { 0 2 4 6 8 10 , , , , , } for DMRS configuration type 1 and m ′ = { 0 1 6 7 , , , } for DMRS configuration type 2 . Then we calculate FFT over m -th row of g m,m ′ PRS,DMRS ( n ) as follows

$$\text{nc}. & v _ { m } ^ { m ^ { \prime } } ( d ) = | \text{FFT} ( g _ { P R S, D M R S } ^ { m, m ^ { \prime } } ( n ) ) | = | \sum _ { k = 1 } ^ { K } ( 2 \beta _ { k } \sum _ { n = 0 } ^ { N - 1 } e ^ { j 2 \pi n T _ { 0 } f _ { d, k } } \\ \text{timation} \\ \text{PDSCH} \\ \text{ins} \ P R S & \text{end, we} \\ \text{in} \ n & \text{ and}$$

where v PRS and v DMRS are transmitted complex symbols of PRS and DMRS , respectively. Afterward, we perform averaging over m as follows

$$\overline { \text{v} } ^ { m ^ { \prime } } ( d ) = \frac { 1 } { M } \sum _ { m = 0 } ^ { M - 1 } v _ { m } ^ { m ^ { \prime } } ( d ). \text{ \quad \ \ } ( 2 0 )$$

Eventually, we obtain the element wise multiplication of v m ′ ( d ) over m ′ as follows

$$\overline { \nu } ( d ) = \prod _ { m ^ { \prime } = 0 } ^ { M ^ { \prime } - 1 } \overline { \nu } ^ { m ^ { \prime } } ( d ). \quad \ \ ( 2 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

At the end, we use (21) to find the index of the maximum value and estimate Doppler. The Doppler frequency shift resolution and the maximum detectable Doppler frequency can be obtained using (16) and (17), respectively. This approach disrupts the periodic pattern of PRS, effectively reducing the magnitude of spurious peaks. Consequently, distinguishing between real and false Doppler frequencies becomes straightforward. If PRS and PDSCH have different resources allocated in the time domain, i.e., slots, we pad zeros to the reference signal with the smaller slots so g m RS P ( n ) and g m MRS ′ D ( n ) can have similar vector sizes.

## V. PARETO OPTIMALITY IN ISAC RESOURCE ALLOCATION

In this section, we introduce an optimization problem to reach Pareto optimal point in time and frequency resource allocation between communication and sensing without Doppler ambiguity through the proposed algorithm. Towards this goal, we define the range resolution metric as c ∆ f , R 0 as the maximum throughput of the received PDSCH per one slot and one PRB, and 1 T s as the Doppler resolution metric. We present the following resource allocation problem to identify the Pareto optimal point.

$$\max _ { m _ { 0 }, n _ { 0 }, \dots, m _ { K }, n _ { K } } & \frac { \alpha _ { 0 } m _ { 0 } n _ { 0 } } { R _ { \max } } R _ { 0 } - \sum _ { k = 1 } ^ { K } ( \frac { \gamma _ { k, 1 } } { d _ { \max } \ m _ { k } } \frac { c } { \Delta f } \\ & + \frac { \gamma _ { k, 2 } } { \nu _ { \max } \ n _ { k } } \frac { 1 } { T _ { s } } ),$$

$$\int _ { n _ { 0 }, \dots, m _ { K }, n _ { K } } K _ { \max } \quad \underline { k = 1 } ^ { \dagger } a _ { \max } \ m _ { k } \ \Delta J \quad \text{$\omega$} \\ + \frac { \gamma _ { k, 2 } } { \nu _ { \max } \ n _ { k } } \frac { 1 } { T _ { s } } ), \quad ( 2 2 a ) \ \text{condu} \\ s. t. \ m _ { i } \in \{ 1, \dots, M _ { \max } - 1 \}, \ i = 0, \dots, K, \ ( 2 2 b ) \quad \text{the $D$} \\ n _ { i } \in \{ 1, \dots, N _ { \max } - 1 \}, \ i = 0, \dots, K, \ ( 2 2 c ) \quad \text{of $30$} \\ \sum _ { i = 0 } ^ { K } m _ { i } = M _ { \max }, \quad ( 2 2 d ) \ \text{DMR} \\ \sum _ { i = 0 } ^ { K } n _ { i } = N _ { \max }, \quad ( 2 2 e ) \ \text{based} \\ \int _ { i = 0 } ^ { K } \int _ { i } \alpha _ { 0 }, \gamma _ { k, 2 } \ \text{represent the communication weight,} \quad 16 \, \mu _ { i } ^ { \dagger } \,. \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cd.$$

where α 0 , γ k, 1 , and γ k, 2 represent the communication weight, range resolution weight, and Doppler resolution weight of target k , respectively. These values can be selected based on the priorities of the application, ensuring that α 0 + ∑ K k =1 ( γ k, 1 + γ k, 2 ) = 1 . M max and N max are available resources in frequency and time domain, respectively. R max, d max, and ν max denote the maximum achievable throughput, range resolution, and Doppler resolution, respectively, when all available resources in the time i.e., N max, and frequency i.e., M max, domains are allocated for communication data transmission or range and Doppler estimation. Additionally, m 0 and n 0 represent the number of PRBs and slots, respectively, allocated for communication, while m k and n k denote the number of PRBs and slots allocated for range and Doppler estimation of the target k . By solving this optimization problem, we

(a) Doppler estimation with PRS only.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-4580cc23-4f80-4cf9-bcf6-db6137e27dbf.jpg)

(b) Doppler estimation with PRS and DMRS combined.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-acdf94d4-89d9-4bca-846e-a1316b74e2fa.jpg)

Fig. 4: Doppler ambiguity elimination with the presence of 2 targets with K comb = 4 , f c = 4 GHz and SCS = 30 kHz.

can identify the Pareto optimal point for resource allocation between communication and sensing.

## VI. SIMULATION RESULTS

To validate our approach, we exploited MATLAB 5G toolboxes, ensuring compatibility with 3GPP standards. We conducted simulations for a scenario where two targets caused the Doppler frequency shift of 4 kHz and 9 kHz. The central frequency was set to 4 GHz with a subcarrier spacing (SCS) of 30 kHz, and the SNR is 0 dB. PRS number of symbols is 12 , DMRS configuration type is 2 , the number of front-loaded DMRS symbols is 2 with 1 additional DMRS symbol position which make 4 symbols in total, PDSCH mapping type is 'A', Code rate used to calculate transport block sizes is 490 / 1024 based on [9]. Virtual resource block (VRB) bundle size is 4 , VRB to PRB interleaving is disabled, and modulation is 16 quadrature amplitude modulation (16-QAM) for PDSCH and β 1 = β 2 = 1 . In Fig. 4a, we can observe the Doppler ambiguity using PRS only, with the real and false peaks of Doppler frequency shifts. However, in Fig. 4b, the magnitudes of the fake Doppler frequencies are significantly suppressed by using the proposed method, and we can easily select the two highest peaks originating from two targets. It is worth mentioning that the number of the targets can be estimated by model order selection based on information-theoretic criteria [15]. In Fig. 5a, we can see an extreme Doppler ambiguity when K comb = 12 while Fig. 5b shows the performance of the proposed method with extreme ambiguity in Doppler estimation. In this case, we can still easily select the highest two peaks and estimate the Doppler frequency shifts of the targets correctly. Once the Doppler ambiguity is mitigated,

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-9bb2e877-8bbd-49fd-adde-2b82572f6dd2.jpg)

0.2

10

5

(a) Doppler estimation with PRS only.

10

15

20

Doppler (kHz)

25

Estimation

0.8

Target 1

(b) Doppler estimation with PRS and DMRS combined.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-40dd27fd-78b6-4488-9d52-51a15086820d.jpg)

Fig. 5: Doppler ambiguity elimination with the presence of 2 targets with K comb = 12 , f c = 4 GHz and SCS = 30 kHz.

the Pareto optimal point, as illustrated in Fig. 6 in time and frequency resource allocation can be found. F in Fig. 6 is the objective function of the optimization problem in (22a). In this simulation, we set M max = 48 , N max = 20 , K = 1 , α 0 = 2 , and γ 1 1 , = γ 1 2 , = 1 . As illustrated in Fig. 6, the Pareto optimal point is attained by allocating 3 slots and 5 PRBs for PRS, with the remaining resources allocated to PDSCH.

## VII. CONCLUSION

In this work, we addressed the Doppler ambiguity appeared using PRS by exploiting DMRS within PDSCH for unconnected targets in ISAC framework. We proposed a novel method to alleviate Doppler ambiguity, which is fully compliant with 3GPP standards and can be integrated into existing 5G networks. Furthermore, our approach enables sensing without imposing additional overhead, requiring modifications to gNB configurations, or compromising Doppler estimation accuracy, while also enhancing the maximum detectable Doppler shift. Finally, we introduced a resource allocation problem between PDSCH and PRS to determine the Pareto optimal point between communication and sensing while effectively eliminating Doppler ambiguity.

## REFERENCES

- [1] J. A. Zhang, M. L. Rahman, K. Wu, X. Huang, Y. J. Guo, S. Chen, and J. Yuan, 'Enabling joint communication and radar sensing in mobile networks-a survey,' IEEE Communications Surveys &amp; Tutorials , vol. 24, no. 1, pp. 306-345, 2021.
- [2] Z. Wei, Y. Wang, L. Ma, S. Yang, Z. Feng, C. Pan, Q. Zhang, Y. Wang, H. Wu, and P. Zhang, '5G PRS-based sensing: A sensing reference signal approach for joint sensing and communication system,' IEEE Transactions on Vehicular Technology , vol. 72, no. 3, pp. 3250-3263, 2022.
- [3] C. D. Ozkaptan, E. Ekici, O. Altintas, and C.-H. Wang, 'OFDM pilotbased radar for joint vehicular communication and radar systems,' in IEEE Vehicular Networking Conference (VNC) , 2018, pp. 1-8.
- [4] Z. Huang, K. Wang, A. Liu, Y. Cai, R. Du, and T. X. Han, 'Joint pilot optimization, target detection and channel estimation for integrated sensing and communication systems,' IEEE Transactions on Wireless Communications , vol. 21, no. 12, pp. 10 351-10 365, 2022.
- [5] L. Ma, C. Pan, Q. Wang, M. Lou, Y. Wang, and T. Jiang, 'A downlink pilot based signal processing method for integrated sensing and communication towards 6G,' in IEEE 95th Vehicular Technology Conference (VTC) , June 2022, pp. 1-5.
- [6] C.-H. Wang, O. Altintas, C. D. Ozkaptan, and E. Ekici, 'Multi-range joint automotive radar and communication using pilot-based OFDM radar,' in IEEE Vehicular Networking Conference (VNC) , 2020, pp. 1-4.
- [7] K. Khosroshahi, P. Sehier, and S. Mekki, 'Leveraging PRS and PDSCH for integrated sensing and communication systems,' 2024, accepted in GLOBECOM 2024 conference. [Online]. Available: https://arxiv.org/abs/2408.00667
- [8] Z. Wei, F. Li, H. Liu, X. Chen, H. Wu, K. Han, and Z. Feng, 'Multiple reference signals collaborative sensing for integrated sensing and communication system towards 5G-A and 6G,' arXiv preprint arXiv:2312.02170 , 2023.
- [9] 3GPP, 'NR; physical layer procedures for data,' 3rd Generation Partnership Project (3GPP), Technical Specification (TS) 38.214, 2024, v18.2.0.
- [10] --, 'NR; physical channels and modulation,' 3rd Generation Partnership Project (3GPP), Technical Specification (TS) 38.211, 2024, v18.2.0.
- [11] S. Buzzi, C. D'Andrea, and M. Lops, 'Using massive mimo arrays for joint communication and sensing,' in 53rd Asilomar Conference on Signals, Systems, and Computers . IEEE, 2019, pp. 5-9.
- [12] K. M. Braun, 'OFDM radar algorithms in mobile communication networks,' Ph.D. dissertation, Karlsruhe, Karlsruher Institut f¨r Technologie u (KIT), 2014.
- [13] L. Pucci, E. Matricardi, E. Paolini, W. Xu, and A. Giorgetti, 'Performance analysis of joint sensing and communication based on 5G new radio,' in IEEE Globecom Workshops , 2021, pp. 1-6.
- [14] J. Li, S. An, J. An, H. Zirath, and Z. S. He, 'OFDM radar range accuracy enhancement using fractional fourier transformation and phase analysis techniques,' IEEE Sensors Journal , vol. 20, no. 2, pp. 1011-1018, 2019.
- [15] A. Mariani, A. Giorgetti, and M. Chiani, 'Model order selection based on information theoretic criteria: Design of the penalty,' IEEE Transactions on Signal Processing , vol. 63, no. 11, pp. 2779-2789, 2015.

Fig. 6: Pareto optimal point in time and frequency resource allocation between communication and sensing.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00679v2-figure-6bfe075f-b168-44f9-a618-be2526a65ce5.jpg)
Keivan Khosroshahi Philippe Sehier Sami Mekki
Doppler Ambiguity Elimination Using 5G Signals in Integrated Sensing and Communication
The industrial point of view towards integrated sensing and communication
(ISAC), the preference is to leverage existing resources and fifth-generation
(5G) infrastructure to minimize deployment costs and complexity. In this
context, we explore the utilization of current 5G new radio (NR) signals
aligned with 3rd generation partnership project (3GPP) standards. Positioning
reference signals (PRS) for sensing and physical downlink shared channel
(PDSCH) for communication have been chosen to form an ISAC framework. However,
PRS-based sensing suffers from Doppler ambiguity when the Doppler frequency
shift is severe. To address this challenge, we introduce a novel method within
the ISAC system that leverages the demodulation reference signal (DMRS) present
in PDSCH to eliminate Doppler ambiguity. Furthermore, we formulate a resource
allocation problem between PRS and PDSCH to achieve a Pareto optimal point
between communication and sensing without Doppler ambiguity. Through
simulations and analysis, we demonstrate the effectiveness of our proposed
method on joint DMRS-PRS exploitation in mitigating Doppler ambiguity and the
efficiency of the resource allocation scheme in achieving Pareto optimality for
ISAC within a 5G NR framework.