Skip to main content
SearchLoginLogin or Signup

Joint Wireless Communication and Sensing: A Comprehensive Outlook

The author analyzed the conflict between communications and sensing in terms of bandwidth and power. It is found that the bandwidth is approximately partitioned between communication and sensing, thus making an almost zero-sum game, while the conflict in the power is marginal.

Published onAug 31, 2024
Joint Wireless Communication and Sensing: A Comprehensive Outlook
·

Abstract

Joint communications and sensing (JCS) can improve the efficiency of power, bandwidth and hardware usage. The conflict between communications and sensing is analyzed in terms of bandwidth and power. It is found that the bandwidth is approximately partitioned between communication and sensing, thus making an almost zero-sum game, while the conflict in the power is marginal. The same conclusion holds when the receiver and target are separated. The feasible region and tradeoff boundary of JCS are obtained, based on which the adaptive signaling in JCS is studied.

Introduction

In recent years, more focus has been paid to the joint communications and sensing (JCS), namely the integration of traditional communication and radar systems [1][2]. A major motivation is the development of various cyber physical systems (CPSs), such as autonomous driving and unmanned vehicular (UAV) networks. In such CPSs, both communications and sensing are needed for the decision of control actions that need the information of the environment. Such information can either be exchanged by the CPS nodes using communications, such as broadcasting the position and velocity information that are obtained from local measurement (e.g., using GPS), or be obtained from sensing, such as radar ranging and Doppler based velocity estimation.

In traditional systems, communications and radar are developed, designed and implemented almost independently, although they share some common theoretical foundation and hardware components. However, for most cases, they both use electromagnetic (EM) wave as the medium. Therefore, they can share the same waveform and complete the tasks of communications and sensing in one round of EM wave emission. An example is shown in Fig. 1, in which vehicle A sends out an EM pulse, which consists of a message to vehicle B using modulations. When the message is delivered to vehicle B, the EM wave may be reflected by vehicle C (which could coincide with vehicle B). The reflected EM wave is intercepted by vehicle A and is then used to infer the information of vehicle C, which is essentially the sensing procedure. Such a JCS procedure can be understood in the following two different ways:

  • The communication message is embedded in a sensing signal via modulation.

  • A communication signal is used to illuminate an target for sensing, similarly to a radar beam.

Moreover, the JCS can be considered as a two-way communication channel: in the forward channel, the transmitter ‘pushes’ the information to a receiver, while in the reverse channel, it ‘pulls’ the information from an target.

Figure 1

An example of JCS in vehicular networks

Although there have been plenty of studies on JCS, which will be detailed in the next section, there lacks systematic study on the fundamental aspect of JCS. A key question is: how do the purposes of communications and sensing conflict with each other? In the frequency division multiplexing (FDM) scheme of JCS, in which the communication and sensing use different frequency bands in the signal, the conflict is in the allocations of bandwidth and power. However, it is no longer obvious when the radar and communication use the same waveform. Understanding the conflict between communications and sensing is of key importance for assessing the tradeoff between the two tasks and the corresponding system design, similarly to many other tradeoff problems in communications, such as multiplexing-diversity tradeoff in multi-antenna systems [3].

In this paper, we study the fundamental tradeoff between communications and sensing, based on frequency modulation waveforms that are widely used in wideband radar systems. Through mathematical analysis, we conclude that (a) the signal bandwidth is approximately partitioned between communications and sensing; hence, the shared waveform scheme is similar to FDM or time division multiplexing (TDM) scheme of JCS; the same conclusion holds when the receiver and sensed targeted are separated; (b) communication and sensing have marginal conflict in the power when the receiver and sensed target are separated, thus making the shared waveform scheme more efficient than the FDM and TDM schemes. The feasible region and optimal tradeoff for JCS will be obtained, under the guidance of which we will study the adaptive signaling in JCS.

The remainder of this paper is organized as follows. The existing studies related to this paper will be detailed in Section 2. The signal model, as well as the modulation scheme, will be discussed in Section 3. Then, the performance metrics of JCS will be analyzed in Section 5, based on which the tradeoff of communications and sensing will be obtained and an adaptive signaling scheme will be proposed. The numerical results and conclusions will be given in Sections 6 and 7, respectively.

Related Work

Joint Communications and Sensing:

In this case, the same EM waveform is shared by communications and sensing, namely the signal carries communication messages to the receiver and the reflected signal carries the information of target for sensing. Comprehensive surveys can be found in [1][2]. The integration of communications and radar has been studied for spread spectrum and multicarrier signals [4], for stepped frequency waveforms [5], for OFDM signals [4], for WiMax signals [6], for frequency hopping MIMO signals [7], and for shift keying signals [8]. It has been exploited for traffic safety in [9]. The designs of waveform and beam pattern have been studied in [10] and [11], respectively. In the mmWave band, the joint design has been analyzed from the signal processing perspective in [12]. In [13], the sharing of 79GHz band for radar and communications in vehicles is studied in terms of time and frequency duplexing. A testbed combining radar and communications in the mmWave band has been built in UT-Austin [14][15], mainly for the purpose of automotive applications. In [16], a joint communication and radar scheme is proposed based on duplexing in different beams. It differs from the proposed research by its dedicated hardware and waveform.

Sensing Waveform for Communications

The typical radar waveforms, particularly the FMCW (chirp) waveform, have been leveraged to convey information [17][18][19][20][21]. The main approach is to consider the FMCW waveform as the carrier, and the information is modulated similarly to the quadrature amplitude modulation (QAM). The author has carried out experiments using chirp radar signals in the mmWave band for communications [19], in which the communication rate is in the order of hundreds of kbps.

Signal Model

In this section, we introduce the signal model of JCS, using traditional radar waveforms.

System Model

We assume that a transmitter sends out signal s(t)s(t), which is received by a communication receiver and reflected by a target as well. For simplicity, we assume that the receiver is a portion of the target; e.g., a vehicle sends a message to another vehicle while sensing the information of the same vehicle. For the cases in which the communication receiver is in a location different from the target, the major concern is the beamforming which directs different portions of power to the receiver and target.

Carrier Signal

We assume that the traditional radar waveform is used as the carrier, over which the communication information is modulated. With one radar pulse, the transmitted signal is given by

s(t)=A(t)exp(2πtf(τ)dτ+θ0),t[0,Tp],(1)\begin{aligned} s(t)=A(t)\exp\left(2\pi \int_{-\infty}^t f(\tau)d\tau+\theta_0\right), t\in [0,T_p], \end{aligned} \tag{1}

where TpT_p is the pulse duration, AA, ff and θ0\theta_0 are the amplitude, frequency and initial phase, respectively.

The following two types of radar waveforms are considered in this paper:

  • Frequency Modulation Continuous Waveform (FMCW): In the FMCW case, the instantaneous frequency increases linearly with time, namely

    f(t)=St+f0,t[0,Tp],(2)\begin{aligned} f(t)=St+f_0,\qquad t\in [0,T_p], \end{aligned} \tag{2}

    where SS is the frequency increasing slope and f0f_0 is the initial frequency.

  • Step-Frequency (SF): In the SF scheme, the frequency is piecewisely constant and increases after each interval of constant frequency, namely

    f(t)=Δfk+f0,t[(k1)Δt,kΔt],(3)\begin{aligned} f(t)=\Delta f k+f_0, t\in [(k-1)\Delta t,k\Delta t], \end{aligned} \tag{3}

    for k=1,...,Kk=1,...,K, where Δf\Delta f is the increase of frequency between successive intervals, and KK is the total number of intervals, such that KΔt=TpK\Delta t=T_p.

We assume that both the FMCW and SF schemes use the same bandwidth, namely STp=ΔfKST_p=\Delta fK.

Wideband JCS Signal: Modulation Scheme

Based on the carrier signal model in (1), we consider the following generic signal model with information modulation:

s(t)=n=0Ns1An(t(n1)Ts)×exp(2πfn(t(n1)Ts)+θn),(4)\begin{aligned} s(t)&=&\sum_{n=0}^{N_s-1} A_n(t-(n-1) T_s) \\ &\times&\exp\left(2\pi f_n(t-(n-1) T_s)+\theta_n\right), \end{aligned} \tag{4}

where NsN_s is the total number of information symbols in one radar pulse, and TsT_s is the symbol period given by Ts=TpNsT_s=\frac{T_p}{N_s}. Similarly to modern communication systems, the information can be modulated in the following parameters of sinusoidal signals:

  • Complex Amplitude: The complex amplitude AA (including the signal magnitude and phase) can be used to convey information, similarly to the amplitude modulation (PAM) in digital communications. However, the variation of amplitude will decrease the average SNR of reflected sensing signal, thus impairing the performance of sensing.

  • Frequency: The frequency can also be varied for conveying information, in both the FMCW and SF schemes, similarly to frequency shift keying (FSK) in digital communications. It adds randomness to the frequency, thus may affect the performance of sensing.

In this paper, we propose two modulation schemes, based on the wideband radar signal waveforms, by following the same principle of typical digital modulations. There are totally four possible combinations of modulation (amplitude and frequency) and carrier waveform (FMCW and SF). For simplicity, we will consider only the schemes of QAM-FMCW and FSK-SF.

QAM-FMCW

We consider modulating the amplitude and phase, similarly to the quadratic amplitude modulation (QAM) in traditional digital communication systems. In a contrast, with a time-varying frequency the carrier in the QAM-FMCW scheme is not a pure sinusoidal function, due to the requirement of wideband radar sensing. Therefore the transmitted signal is given by

s(t)=n=0Ns1An(tnTs)exp(j2πθ(t)),(5)\begin{aligned} s(t)&=&\sum_{n=0}^{N_s-1} A_{n}(t-nT_s)\exp\left(j2\pi \theta(t)\right), \end{aligned} \tag{5}

where AnA_n is the complex symbol similar to QAM, and θ(t)\theta(t) is the deterministic (thus carrying no information) phase determined by the FMCW scheme, which is given by

θ(t)=0tf(s)ds=12St2+f0t.(6)\begin{aligned} \theta(t)=\int_0^t f(s)ds=\frac{1}{2}St^2+f_0t. \end{aligned} \tag{6}

FSK-SF

In the FSK-SF scheme, the amplitude AA is a constant. For assuring the performance of communications, we assume KNsK\gg N_s and K/NsK/N_s is an integer. We can allow some randomness in the frequency of each symbol period, namely the frequency at time tt is given by

f(t)=(k1)Δf+f0+mk/NsΔfM,(7)\begin{aligned} f(t)=(k-1)\Delta f+f_0+\frac{m_{\lfloor k/N_s\rfloor}\Delta f}{M}, \end{aligned} \tag{7}

for t[(k1)Δt,kΔt]t\in [(k-1)\Delta t,k\Delta t], k=1,...,Kk=1,...,K, where m{0,...,M1}m\in \{0,...,M-1\} is a random variable that carries the information, and each symbol carries log2M\log_2M bits.

Received Signals

We first consider the received signal at the radar transceiver for the purpose of sensing. The elapsed time between the transmission and reception of the JCS signal is given by td=2dct_d=\frac{2d}{c}, where dd is the distance between the transmitter and the target, and cc is the light speed. For simplicity of analysis, we assume that the target is stationary, thus waiving of the Doppler shift. Then, the received reflected signal is given by

r(t)=n=0Ns1An(tnTs)(t(n1)Ts2dc)×exp(j2πθ(t2dc)+θn),(8)\begin{aligned} r(t)&=&\sum_{n=0}^{N_s-1} A_{n}(t-nT_s)\left(t-(n-1) T_s-\frac{2d}{c}\right) \\ &\times&\exp\left(j2\pi \theta\left(t-\frac{2d}{c}\right)+\theta_n\right), \end{aligned} \tag{8}

The signal received by the communication receiver has the same form as (8), except for replacing dd with d2\frac{d}{2}, since we assume that the communication receiver is colocated with the target. We assume that the communication receiver is perfectly time synchronized with the JCS transmitter, such that the traveling time d2c\frac{d}{2c} can canceled at the mixer.

Joint Communications and Sensing

In this section, we introduce the algorithms for communication demodulation and sensing in the JCS scheme.

Demodulation in Communications

We first introduce the demodulation at communications receiver, for the QAM-FMCW and FSK-SF schemes, respectively. As will be seen, the radar carrier can be perfectly removed, provided that the time is perfectly synchronized. Then the baseband signal can be used for detecting the communication message, as in traditional communication receivers.

QAM-FMCW

For the QAM modulation with FMCW carrier, the demodulation is similar to that of traditional linear modulation with sinusoidal carrier. Due to the assumption of perfect time synchronization, the input signal is mixed with local FMCW signal and the output is passed through a low pass filter. The procedure is given by

r(t)=n=0Ns1An(tnTs)exp(jθ(t))exp(jθ(t))=n=0Ns1An(tnTs),(9)\begin{aligned} r(t)&=&\sum_{n=0}^{N_s-1} A_{n}(t-nT_s)\exp\left(j\theta(t)\right)\exp\left(-j\theta(t)\right) \\ &=&\sum_{n=0}^{N_s-1} A_{n}(t-nT_s), \end{aligned} \tag{9}

which is then put in the decision maker of traditional demodulator. Here the perfect cancellation of the carrier phase is due to the assumption of perfect time synchronization.

FSK-SF

The received signal is mixed with a local mixer with the proper center frequency (k1)Δf+f0(k-1)\Delta f+f_0, namely

exp(j2π((k1)Δf+f0+mΔfM)t)×exp(j2π((k1)Δf+f0)t)=exp(j2πmΔfMt),(10)\begin{aligned} &&\exp\left(j2\pi\left((k-1)\Delta f+f_0+\frac{m\Delta f}{M}\right)t\right)\\ &\times&\exp\left(-j2\pi\left((k-1)\Delta f+f_0\right)t\right)\\ &=&\exp\left(j2\pi\frac{m\Delta f}{M}t\right), \end{aligned} \tag{10}

where the output is based on the assumption of perfect time synchronization. Then, the frequency of the mixer output is estimated and the corresponding mm is obtained. Traditional demodulation approaches of FSK can be employed.

Radar Ranging

Now, we study the ranging algorithm using the reflected radar waveform. In contrast to traditional radar ranging algorithms, the major challenge here is how to remove the impact of modulated information.

QAM-FMCW

The received signal at the radar receiver is mixed with local FMCW signal (without modulation) and the output is given by

r(t)=n=0Ns1An(t(n1)Tsτ))exp(j2πτSt),(11)\begin{aligned} r(t)=\sum_{n=0}^{N_s-1} A_n\left(t-(n-1)T_s-\tau)\right)\exp\left(j2\pi\tau St\right), \end{aligned} \tag{11}

where τ=2dc\tau=\frac{2d}{c} is the round trip time, which is to be estimated. To obtain the round trip time τ\tau, we can use the following two approaches:

  • Frequency domain estimation: Notice that AnA_n is known to the radar receiver, thus being deterministic, since it is transmitted by the same transmitter. Hence, we can carry out Fourier transform for the mixer output in (11) and obtain

    r^(jw)=n=0Ns1A^n(j(w2πτS))ej(w2πτS)((n1)Ts+τ)=F(j(w2πτS)),(12)\begin{aligned} \hat{r}(jw)&=&\sum_{n=0}^{N_s-1} \hat{A}_n(j(w-2\pi\tau S)) e^{-j(w-2\pi\tau S)((n-1)T_s+\tau)}\\ &=&F(j(w-2\pi \tau S)), \end{aligned} \tag{12}

    where F=A^n(jw)ejw((n1)Ts+τ)F=\hat{A}_n(jw)e^{-jw((n-1)T_s+\tau)}.

    When the mixer output r^\hat{r} is contaminated by Gaussian noise, the maximum likelihood estimation of τ\tau is given by

    τML=argminτr^(jw)F(j(w2πτS))2.(13)\begin{aligned} \tau_{ML}=\arg\min_{\tau}|\hat{r}(jw)-F(j(w-2\pi\tau S))|^2. \end{aligned} \tag{13}
  • Carrier synchronization: We can also consider the signal in (11) as the carrier with frequency τS\tau S modulated by the information symbols AnA_n. For a typical case of d=200md=200m and S=30MHz/usS=30MHz/us, the beat frequency τS=40MHz\tau S=40MHz. If the symbol frequency (e.g., 1MHz) is much smaller than τS\tau S, we can employ the carrier synchronization approach to estimate τ\tau. When NsN_s is sufficiently large, we can calculate the power spectrum density, estimate the center frequency f^c\hat{f}_c, and calculate τ^=f^cS\hat{\tau}=\frac{\hat{f}_c}{S}.

FSK-SF

Since the radar receiver knows the information symbol mm and thus the frequency, it mixes the received signal with local oscillation. The output of the mixer and filter is given by

f(t)={(mkmkk)ΔfM+kΔf,t[t0,t1](mkmkk+1)ΔfM+(k1)Δf,t[t1,t2].(14)\begin{aligned} f(t)=\left\{ \begin{array}{ll} &\frac{(m_k-m_{k-k'})\Delta f}{M}+k'\Delta f, t\in [t_0,t_1]\\ &\frac{(m_k-m_{k-k'+1})\Delta f}{M}+(k'-1)\Delta f, t\in [t_1,t_2] \end{array} \right.. \end{aligned} \tag{14}

where t0=(k1)Δtt_0=(k-1)\Delta t, t1=t0+{τΔt}t_1=t_0+\{\frac{\tau}{\Delta t}\}, t2=kΔtt_2=k\Delta t (where {}\{\cdot\} means the fractional part of the number) and k=τΔtk'=\left\lceil \frac{\tau}{\Delta t}\right\rceil. The timing is illustrated in Fig. 2.

Figure 2

Timing in the mixer output of FSK-SF

The following four steps are used to estimate τ\tau and thus the distance:

  • Step 1. Estimate the frequency as a function of time tt. In the ideal case, we assume that f(t)f(t) is perfectly obtained. In practice, short-time DFT can be used to estimate the frequency by selecting the peak in the spectrum.

  • Step 2. Coarse estimation: In typical setups, the term kΔfk'\Delta f is significantly larger than the term (mkmkk)ΔfM\frac{(m_k-m_{k-k'})\Delta f}{M}, since the latter is only a fine adjustment of the frequency. Therefore, kk' can be estimated as k^=f(t)Δf\hat{k}'=\left\lfloor\frac{f(t)}{\Delta f}\right\rfloor.

  • Step 3. Fine estimation: Search for the optimal t1t_1 (namely the frequency change time of the received signal) such that the MSE of frequency is minimized, namely

    t^1=argmint1t0t1(f^(t)f(tk^,t1))2dt+t1t2(f^(t)f(tk^,t1))2dt,(15)\begin{aligned} \hat{t}_1&=&\arg\min_{t_1} \int_{t_0}^{t_1}(\hat{f}(t)-f(t|\hat{k}',t_1))^2 dt \\ &+&\int_{t_1}^{t_2}(\hat{f}(t)-f(t|\hat{k}',t_1))^2dt, \end{aligned} \tag{15}

    where f^(t)\hat{f}(t) is the estimated frequency using the received signal and f(tk^,t1)f(t|\hat{k}',t_1) is the frequency given in (14) conditioned on k^\hat{k}' and t1t_1. Note that the above integral can be replaced with summations of the frequency samples. When there are not sufficiently many samples for tracking the instantaneous frequency within one symbol period, the fine estimation of frequency cannot be accomplished.

  • Step 4. Fractional part: The estimation of the frequency change time is given by t^1\hat{t}_1. Then, the fraction part of τΔt\frac{\tau}{\Delta t} is given by t^1(k^1)Δt\hat{t}_1-(\hat{k}'-1)\Delta t.

Then, the round trip time τ\tau is estimated as

τ^=k^Δt+t^1(k^1)Δt.(16)\begin{aligned} \hat{\tau}=\hat{k}'\Delta t+\hat{t}_1-(\hat{k}'-1)\Delta t. \end{aligned} \tag{16}

Performance Analysis

In this section, we analyze the performance of the proposed JCS scheme. For simplicity, we focus on only QAM-FMCW.

Communication Performance

The exact analysis of the communication bit error rate of QAM modulation is difficult. We consider the following upper bound for the bit error rate of QAM (Eq. (4.3-30) in [22]):

Pe,MQAM4Q(3log2MM1PrN0),(17)\begin{aligned} P_{e,M-QAM}\leq 4Q\left(\sqrt{\frac{3\log_2 M}{M-1}\frac{P_r}{N_0}}\right), \end{aligned} \tag{17}

where PrP_r is the received signal and N0N_0 is the power spectral density (PSD) of the thermal noise, This upper bound becomes tight when the signal-to-noise ratio PrN0\frac{P_r}{N_0} is sufficiently large.

We leverage the following approximation:

Q(x)112ex22,(18)\begin{aligned} Q(x)\approx \frac{1}{12}e^{-\frac{x^2}{2}}, \end{aligned} \tag{18}

when xx is sufficiently large. Then, the channel capacity of the binary symmetric channel, in terms of bits per channel use is

C1ex223,(19)\begin{aligned} C\approx 1-\frac{e^{-\frac{x^2}{2}}}{3}, \end{aligned} \tag{19}

where x=3log2MM1PtGN0Nsx=\sqrt{\frac{3\log_2 M}{M-1}\frac{P_tG}{N_0N_s}}, PtP_t is the transmit energy of the radar pulse and GG is the channel gain. Therefore, the data throughput is given by

TNsTP(1ex223)log2M,(20)\begin{aligned} T\approx\frac{N_s}{T_P}\left(1-\frac{e^{-\frac{x^2}{2}}}{3}\right)\log_2M, \end{aligned} \tag{20}

since the channel use lasts time TsT_s.

Radar Performance

The performance of FMCW radar has been analyzed in [23][24]. However, the SNR is assumed to be sufficiently high therein; therefore the performance is dependent on only the frequency increasing slope and the pulse duration, thus being independent of the signal power. Moreover, the impact of QAM symbols (which change the amplitude of the FMCW carrier in a random manner) is novel. In the context of JCS, it is difficult to directly analyze the performances of the proposed schemes of frequency domain estimation and carrier synchronization in Section 4.2.

In this paper, we follow the approach in the following lemma on the frequency estimation in sinusoidal signals contaminated by noise, in order to find the lower bound of the estimation error variance.

Lemma 1 (p.57 of [25]). For a sinusoidal function with frequency ff, SNR γ\gamma and NN samples, the Cramer-Rao bound is given by

V[f]12(2π)2γN(N21)12(2π)2γN3(21)\begin{aligned} V[f]\geq \frac{12}{(2\pi)^2 \gamma N(N^2-1)}\approx \frac{12}{(2\pi)^2 \gamma N^3} \end{aligned} \tag{21}

The approach is based on the Fisher Information matrix defined by

Iij=E[(logf(x,θ)θi)(logf(x,θ)θj)],(22)\begin{aligned} I_{ij}=E\left[\left(\frac{\partial \log f(x,\mathbf{\theta})}{\partial \theta_i}\right)\left(\frac{\partial \log f(x,\mathbf{\theta})}{\partial \theta_j}\right)^*\right], \end{aligned} \tag{22}

where ff is the probability density function of the sample.

Proposition 1. When NN samples are taken for the output of the mixer in the QAM-FMCW scheme, the MSE of ranging is lower bounded by

MSEd3c28π2S2γN(N+1(2N+1)).(23)\begin{aligned} MSE_d\geq \frac{3c^2}{8\pi^2S^2\gamma N(N+1(2N+1))}. \end{aligned} \tag{23}

Remark 1. When the number of samples is sufficiently large, we have

MSEd3c216π2S2γN3.(24)\begin{aligned} MSE_d\geq \frac{3c^2}{16\pi^2S^2\gamma N^3}. \end{aligned} \tag{24}

We observe that the lower bound is irrelevant to the number of information symbols NsN_s. This implies that the modulation for communications have not impact on the Cramer-Rao bound of radar sensing performance.

Proof. When noise is taken into account, the mixer output of QAM-FMCW in (11) is given by

r(t)=n=0Ns1Anexp(j2πτSt)+n(t)=n=0Ns1Anexp(jϕ(t))+n(t)(25)\begin{aligned} r(t)&=&\sum_{n=0}^{N_s-1} A_n\exp\left(j2\pi\tau St\right)+n(t) \\ &=&\sum_{n=0}^{N_s-1} A_n\exp\left(j\phi(t)\right)+n(t) \end{aligned} \tag{25}

where τ=2dc\tau=\frac{2d}{c} is the round trip time, n(t)n(t) is the noise, and

ϕ(t)=4πdStc.(26)\begin{aligned} \phi(t)=\frac{4\pi dSt}{c}. \end{aligned} \tag{26}

NN samples are taken for the signal r(t)r(t), and the samples are denoted by {yn}n=1,...,N\{y_n\}_{n=1,...,N}. The sampling period is δt=TPN=TPGNs\delta t=\frac{T_P}{N}=\frac{T_P}{GN_s}.

We calculate the Fisher’s information:

logf(x)d=dn=1N12PNynA(n)exp(jϕ(tn))2=jn=1NA(n)PNn(tn)exp(jϕ(tn))dϕ(tn)=j4πScPNn=1NA(n)nδtexp(jϕ(tn)),(27)\begin{aligned} \frac{\partial \log f(x)}{\partial d}&=&\frac{\partial}{\partial d}\sum_{n=1}^N \frac{1}{2P_N}\left|y_n-A^{(n)} \exp\left(j\phi(t_n) \right)\right|^2\\ &=&j\sum_{n=1}^N \frac{A^{(n)}}{P_N} n(t_n)\exp\left(j\phi(t_n)\right)\frac{\partial }{\partial d}\phi(t_n)\\ &=&\frac{j4\pi S}{cP_N}\sum_{n=1}^N A^{(n)} n\delta t\exp\left(j\phi(t_n)\right), \end{aligned} \tag{27}

where A(n)=AnGA^{(n)}=A_{\left\lceil\frac{n}{G}\right\rceil}.

Therefore, we have

E[logf(x)d(logf(x)d)]=16π2S2c2PNn=1Nn2E[(A(n))2]=8π2S2γN(N+1(2N+1))3c2.(28)\begin{aligned} &&E\left[\frac{\partial \log f(x)}{\partial d}\left(\frac{\partial \log f(x)}{\partial d}\right)^*\right]\\ &=&\frac{16\pi^2S^2}{c^2P_N}\sum_{n=1}^Nn^2E\left[(A^{(n)})^2\right]\\ &=&\frac{8\pi^2S^2\gamma N(N+1(2N+1))}{3c^2}. \end{aligned} \tag{28}

It concludes the proof by applying the Cramer-Rao bound. ◻

Bandwidth and Power

Now, we consider the consumption of bandwidth and power for the tasks of communications and radar sensing. When separate waveforms and frequency bands are allocated to communications and radar (either dynamically or stationarily), the power and bandwidth are split to the two tasks. Our goal here is to study whether they are also split to the two tasks.

Here we consider a generic radar waveform, which is given by

s(t)=m=n=1NsAmngn(t(m1)Tp),(29)\begin{aligned} s(t)=\sum_{m=-\infty}^\infty \sum_{n=1}^{N_s} A_{mn} g_n(t-(m-1)T_p), \end{aligned} \tag{29}

where gng_n is the carrier waveform in the nn-th symbol period within the same pulse, namely

gn(t)=exp(j2π(12St+f0t))Wn(t).(30)\begin{aligned} g_n(t)=\exp\left(j2\pi\left(\frac{1}{2}St+f_0t\right)\right)W_n(t). \end{aligned} \tag{30}

where Wn(t)W_n(t) is the rectangular window function between [nTs,(n+1)Ts][nT_s,(n+1)T_s]

The information sequence {Am,n}\{A_{m,n}\} are i.i.d, which makes the signal a cyclo-stationary random process, due to the repeated carrier waveform in different radar pulses.

Using the same argument as in [26] (Section 3.4-2), the autocorrelation function R(t+τ,t)R(t+\tau,t) is given by

R(t+τ,t)=m=n=1NsPs×gn(t(m1)Tp+τ)gn(t(m1)Tp).(31)\begin{aligned} &&R(t+\tau,t)=\sum_{m=-\infty}^\infty \sum_{n=1}^{N_s}P_s\\ &\times&g_n(t-(m-1)T_p+\tau)g_n(t-(m-1)T_p). \end{aligned} \tag{31}

It is easy to verify

R(t+τ,t)=R(t+τ+Tp,t+Tp),(32)\begin{aligned} R(t+\tau,t)=R(t+\tau+T_p,t+T_p), \end{aligned} \tag{32}

which means that the signal is cyclo-stationary. The term PsP_s is the power of signal. We denote it by PsP_s.

Due to the cyclostationarity, the autocorrelation function can be decomposed to the Fourier series

R(t+τ,t)=n=RTpn(τ)ej2πntTp,(33)\begin{aligned} R(t+\tau,t)=\sum_{n=-\infty}^\infty R^{n}_{T_p}(\tau)e^{ \frac{-j2\pi nt}{T_p}}, \end{aligned} \tag{33}

where the Fourier series are given by

RTpm(τ)=1Tp0TpR(t+τ,t)ej2πmtTpdt=PsTpl=n=0Ns1nTs(n+1)Tsgn(tlTp)gn(t+τlTp)ej2πmtTpdt=PsTpn=0Ns1nTs(n+1)Tsτgn(t)gn(t+τ)ej2πmtTpdt=PsTpn=0Ns1gn(τ)(gn(τ)ej2πmtTp).(34)\begin{aligned} &&R^{m}_{T_p}(\tau)\\ &=&\frac{1}{T_p}\int_{0}^{T_p}R(t+\tau,t)e^{\frac{-j2\pi mt}{T_p}}dt\\ &=&\frac{P_s}{T_p}\sum_{l=-\infty}^\infty\sum_{n=0}^{N_s-1} \\ &&\int_{nT_s}^{(n+1)T_s}g_n(t-lT_p)g^*_n(t+\tau-lT_p)e^{\frac{-j2\pi mt}{T_p}}dt\\ &=&\frac{P_s}{T_p}\sum_{n=0}^{N_s-1} \int_{nT_s}^{(n+1)T_s-\tau}g_n(t)g_n^*(t+\tau)e^{\frac{-j2\pi mt}{T_p}}dt\\ &=&\frac{P_s}{T_p}\sum_{n=0}^{N_s-1} g_n(\tau)\ast \left(g_n(-\tau)e^{\frac{j2\pi mt}{T_p}}\right). \end{aligned} \tag{34}

Notice that the spectrum of gn(t)g_n(t) is given by

gn^(f)=gWn^(f)=sinc(fTs)ej2πnTsfg^(f),(35)\begin{aligned} \widehat{g_n}(f)&=&\widehat{gW_n}(f)\\ &=&\text{sinc}(fT_s)e^{-j2\pi nT_sf}\ast \hat{g}(f), \end{aligned} \tag{35}

where ^\hat{\cdot} means Fourier transform, and the term ej2πnTsfe^{-j2\pi nT_sf} is due to the time shift of the window function WnW_n.

Then, the power spectrum density of the mm-th order is given by

STsm(f)=n=0Ns1gn^(f)gn^(fmTp)=PsTsn=0Ns1sinc(fTs)ej2πnTsfg^(f)×sinc((fmTp)Ts)ej2πnTs(fmTp)g^(fmTp)=n=0Ns1ej2πmnNsPsTssinc(fTs)g^(f)×sinc((fmTs)Ts)g^(fmTs)={NsPsTssinc(fTs)g^(f)2,m=00,m>0.(36)\begin{aligned} \mathcal{S}^{m}_{T_s}(f)&=&\sum_{n=0}^{N_s-1} \widehat{g_n}(f) \widehat{g_n}\left(f-\frac{m}{T_p}\right)\\ &=&\frac{P_s}{T_s}\sum_{n=0}^{N_s-1}\text{sinc}\left(fT_s\right)e^{-j2\pi nT_sf} \ast \hat{g}(f)\\ &\times&\text{sinc}\left(\left(f-\frac{m}{T_p}\right)T_s\right)e^{j2\pi nT_s\left(f-\frac{m}{T_p}\right)} \ast\hat{g}^*\left(f-\frac{m}{T_p}\right)\\ &=&\sum_{n=0}^{N_s-1}e^{-j\frac{2\pi mn}{N_s}}\frac{P_s}{T_s}\text{sinc}\left(fT_s\right) \ast \hat{g}(f)\\ &\times&\text{sinc}\left(\left(f-\frac{m}{T_s}\right)T_s\right)\ast\hat{g}^*\left(f-\frac{m}{T_s}\right)\\ &=&\left\{ \begin{array}{ll} &\frac{N_sP_s}{T_s}\left|\text{sinc}(fT_s)\ast \hat{g}(f)\right|^2,m=0\\ &0,m>0 \end{array} \right.. \end{aligned} \tag{36}

where in the first equality the term fmTsf-\frac{m}{T_s} is due to the term ej2πmtTpe^{\frac{j2\pi mt}{T_p}} in (34) which results in the shift in the frequency spectrum, and the last equality is due to the fact

n=0Ns1ej2πmnNs={Ns,m=00,m>0.(37)\begin{aligned} \sum_{n=0}^{N_s-1}e^{-j\frac{2\pi mn}{N_s}} =\left\{ \begin{array}{ll} &N_s,m=0\\ &0,m>0 \end{array} \right.. \end{aligned} \tag{37}

In particular, when m=0m=0, which means the time average power spectrum density, we have

STs0(f)=PsTssinc(fTs)g^(f)2.(38)\begin{aligned} \mathcal{S}^0_{T_s}(f)=\frac{P_s}{T_s}\left|\text{sinc}\left(fT_s\right) \ast \hat{g}(f)\right|^2. \end{aligned} \tag{38}

An interesting observation on (38) is that the PSD of the modulated radar signal is determined by the convolution of the following two terms:

  • The sinc function, which represents the bandwidth needed by the communications, since the data rate is inversely proportional to the symbol period TsT_s and thus approximately proportional to the bandwidth BcB_c of the sinc function.

  • The radar signal spectrum g^(f)\hat{g}(f), which determines the bandwidth BsB_s needed for radar sensing. Note that the spectrum bandwidth of FMCW signal can be found in [27].

Due to the property of convolution, the total bandwidth needed for the JCS is given by

BtBc+Bs,(39)\begin{aligned} B_t\approx B_c+B_s, \end{aligned} \tag{39}

which means that the total bandwidth is implicitly partitioned to the bandwidths of communications and sensing. Essentially the functions of communications and radar sensing are not sharing the bandwidth. Each degree of freedom is either used for communications, or used for radar sensing. From this viewpoint of bandwidth, the JCS scheme of shared waveform is approximately the same efficient as a FDMA or TDMA separation of communications and sensing.

In summary, the scheme of shared waveform for JCS is similar to that of FDMA and TDMA JCS in terms of bandwidth efficiency. Meanwhile, the shared waveform is more efficient in terms of power, since the FDMA or TDMA does not utilize the reflected power in the time slot / frequency band for communications.

Numerical Simulations

In this section, we provide the numerical simulation results to evaluate the performance of the proposed JCS schemes.

QAM-FMCW

Ranging Performance

We tested both the schemes of frequency domain estimation and carrier synchronization. The default parameters are shown in Table 1, where most of the radar parameters follow that of TI AWR1243BOOST radar in the 77GHz band.

Table 1

distance

100m

chirp duration TpT_p

60us

symbol number NsN_s

8

modulation

16QAM

freq. inc. rate SS

29.98THz/s

pulse number

50,000

Configuration of simulations for QAM-FMCW

Note that the ML estimation in (13) provides the asymptotically performance only when the processing is analog, namely τ\tau is continuous. However, in practice, discrete Fourier transform is needed for the discrete-time samples. Therefore, we change the minimization to ranging all the frequency indices in a discrete manner. For the carrier synchronization scheme, we simply pick the highest power peak in the frequency spectrum of mixer output, and set the corresponding frequency as the estimation f^c\hat{f}_c.

Figure 3

Ranging MSE with different noise powers and distances

We first tested the performance of ranging, in terms of the standard deviation of ranging errors, with respect to various noise power levels (suppose that the signal power is normalized to 1) and target distances. Both the frequency domain approach and carrier synchronization approach are tested. The results are shown in Fig. 3. An interesting observation is that the frequency domain approach performs far worse than the carrier synchronization approach. The underlying reason is that the frequency domain changes rapidly for different delays τ\tau. Since only discrete samples in the frequency domain are used, an estimation τ^\hat{\tau} close to the true value τ\tau does not guarantee a small matching error in (13). In a sharp contrast, the carrier synchronization approach results in much lower ranging errors. Only when the noise power becomes much higher does the ranging error become significant.

Figure 4

Ranging MSE with different noise powers and Ns

In Fig. 4 we compare the performance of ranging with respect to different numbers of symbols in one radar pulse, namely NsN_s. We observe that, as NsN_s increases, the performance of the frequency domain approach is improved, while that of the carrier synchronization is impaired. The detailed reason for this change is still under exploration. Since the performance of the frequency domain approach is unacceptable, we will give up this approach in the subsequent simulations. Then, the higher data rate, caused by the higher number of symbols in one chirp pulse, will result in a negative impact on the performance of ranging, thus yielding a conflict of interest between communications and radar.

Figure 5

Ranging MSE with different noise powers and sampling frequencies

The ranging performance with respect to different sampling frequencies is shown in Fig. 5. We observe that, when the noise power is significant, the higher sampling rate will substantially improve the performance, as being expected. The impact of different modulation schemes is shown in Fig. 6. We observe that the higher modulation order is, the worse the ranging performance is, when the noise power is significant. The QPSK, which has a constant envelop, achieve the optimal performance. It implies that the variation in the signal amplitude cause the performance degradation of ranging.

Figure 6

Ranging MSE with different noise powers and modulation schemes

Communication-Sensing Tradeoff

Figure 7

From the above simulation results, we observe that the order of modulation, which determines the communication data rate, has marginal impact on the sensing performance, unless the noise power is very high. Hence, the major tradeoff between communications and sensing is due to the implicit sharing of bandwidth as shown in Section 5.3. In Fig. 7, we show the tradeoff between communication and sensing when different bandwidths are allocated to the communications while the total bandwidth is kept a constant. The modulation schemes of 16QAM and QPSK are tested. We observe that an increasing channel capacity implies worse performance of sensing (higher ranging MSE). An interesting observation is that the tradeoff curve is better for the QPSK modulation scheme. One possible reason is that the QPSK scheme has a constant envelop, while the phase change does not affect the performance of sensing. When QAM is used, the variance of signal power (thus the SNR) has a negative impact on the ranging performance. Essentially, this is due to the Jensen’s inequality: as shown in Prop. 1, the Cramer-Rao bound MSEcrMSE_{cr}

MSEcr=E[C0γ]C0E[γ],(40)\begin{aligned} MSE_{cr}=E\left[\frac{C_0}{\gamma}\right]\geq \frac{C_0}{E[\gamma]}, \end{aligned} \tag{40}

where C0C_0 is the constant coefficient and γ\gamma is random due to the different power at different symbol periods.

FSK-SF

Table 2

distance

100m

chirp duration TpT_p

60us

symbol number NsN_s

8

modulation

8-FSK

freq. levels

512

samp. freq.

136Msps

Configuration of simulations for FSK-SF

We also tested the performance of FSK-SF, whose parameters are shown in Table 2. In Fig. 8, we have shown the performance of communication and sensing with respect to different values of NsN_s and MM. The upper figure shows the points characterizing the performances. In each cluster of points due to the same NsN_s, different points corresponds to different MM’s. We observed that a higher NsN_s incurs better communication performance but worse sensing performance. An optimal MM can be selected for the optimal tradeoff between communication and sensing. In the middle subplot, we observe that the modulation order MM needs to be an intermediate number in order to maximize the communication performance. It is shown in the lower subplot that the modulation order MM causes marginal performance difference for the ranging.

Figure 8

Performances of FSK-SF with respect to different Ns and M

The performance of communications and sensing is shown in Fig. 9 for different noise powers at the communication and radar receivers. An observation is that the increase of noise power, either at the communication receiver or radar sensor, has marginal impact on the communications and sensing, unless the noise power is sufficiently high.

Figure 9

Performances of FSK-SF with respect to noise powers

Conclusions

In this paper, we have studied the joint communications and sensing by leveraging the radar sensing waveforms. The approaches of QAM-FMCW and FSK-SF are proposed and analyzed. A major conclusion is that the total bandwidth used by the JCS is approximately the sum of the bandwidth used by communications and sensing, thus having no significant advantage over the scheme of separate communications and sensing, in terms of spectral efficiency.

Comments
0
comment
No comments here
Why not start the discussion?