This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Beamforming Design for Joint Target Sensing and Proactive Eavesdroppingthanks: Manuscript received.

Qian Dan, Hongjiang Lei, Ki-Hong Park, Gaofeng Pan and Mohamed-Slim Alouini This work was supported by the National Natural Science Foundation of China under Grant 61971080. (Corresponding author: Hongjiang Lei).Qian Dan is with the School of Communications and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China, and also with the School of Computer Science, Jiangxi University of Chinese Medicine, Nanchang 330004, China (e-mail: danqian@jxutcm.edu.cn).Hongjiang Lei is with the School of Communications and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China (e-mail: leihj@cqupt.edu.cn).Gaofeng Pan is with the School of Cyberspace Science and Technology, Beijing Institute of Technology, Beijing 100081, China (e-mail: gaofeng.pan.cn@ieee.org).Ki-Hong Park and Mohamed-Slim Alouini are with CEMSE Division, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia (e-mail: kihong.park@kaust.edu.sa, slim.alouini@kaust.edu.sa).
Abstract

This work studies the beamforming design in the joint target sensing and proactive eavesdropping (JTSAPE) system. The JTSAPE base station (BS) receives the information transmitted by the illegal transmitter and transmits the waveform for target sensing. The shared waveform also serves as artificial noise to interfere with the illegal receiver, thereby achieving proactive eavesdropping. We firstly optimize the transmitting beam of the BS to maximize the eavesdropping signal-to-interference-plus-noise ratio or minimize the target estimation parameter Cramér-Rao bound, respectively. Then, the joint optimization of proactive eavesdropping and target sensing is investigated, and the normalized weighted optimization problem is formulated. To address the complexity of the original problem, the formulated problem is decomposed into two subproblems: proactive eavesdropping and target sensing, which are solved by the semi-definite relaxation technique. Furthermore, the scenario in which the quality of the eavesdropping channel is stronger than that of the illegal channel is considered. We utilize the sequential rank-one constraint relaxation method and iteration technique to obtain the high-quality suboptimal solution of the beam transmit covariance matrix. Numerical simulation shows the effectiveness of our proposed algorithm.

Index Terms:
Joint target sensing and proactive eavesdropping (JTSAPE), beamforming design, eavesdropping signal-to-noise ratio (SNR), Cramér-Rao bound (CRB).

I Introduction

I-A Background and Related Works

Integrated sensing and communication (ISAC) technology can share spectrum, hardware platforms, and even baseband waveforms and signal processing between communication and sensing to improve the spectral efficiency, energy efficiency, and hardware efficiency of the system for integration gain. Further, the mutual assistance and mutual gain of the two functions can also be used to improve each other’s performance and obtain coordination gain [1]-[4].

Recently, the performance for various ISAC systems was investigated to achieve concurrent sensing and communication, which were classified into the radar-centric systems [5]-[9], the communication-centric systems [10]-[13], and the joint-design systems [14]-[15]. The authors in [5] investigated beamforming designs for the ISAC systems with multiple communication users (CUs) for scenarios with point or extended targets. In particular, the Cramér-Rao bound (CRB) of the azimuth angle or the response matrix of the target was minimized by designing the dual-functional beamforming matrix. The authors in [6] proposed an integrated sensing, communication, and computation over-the-air framework to enable simultaneous communication, target sensing, and over-the-air computation (AirComp). The shared and separated schemes were proposed, and the accuracy of AirComp was optimized by jointly designing the data transmission beamformer and the aggregation beamformer or the data beamformer and the radar sensing beamformer. In [7], a new angular waveform similarity metric was proposed and combined with the integrated main-lobe-to-sidelobe ratio (IMSR) of the transmit beampattern to evaluate the waveform ambiguity properties. The complementary waveform and linear precoder matrix were jointly designed to minimize the weighted sum of IMSR while considering a predefined signal-to-interference-plus-noise ratio (SINR), power, and peak-to-average-power ratio (PAPR) constraints. In [8], the authors investigated the performance of the ISAC systems with multiple targets and CUs. Both the scenarios without and with prior target knowledge at the BS were considered, and the multi-target response matrix and the reflection coefficients/angles were estimated, respectively. The CRB of multi-target estimation was minimized while satisfying the constraints of a minimum multicast communication rate and a maximum transmit power. The authors in [9] proposed a new robust beamforming scheme for ISAC systems with different point targets. The maximum CRB of the estimated direction of arrival (DOA) was minimized by designing the beamforming matrix while ensuring the quality-of-service (QoS) for all the CUs. In [10], the authors proposed a new waveform design based on constructive interference for the ISAC system with multiple CUs. Their results show that the proposed scheme can reduce the transmit power for the given constraint on estimation accuracy or increase the communication signal-to-noise ratio (SNR) with a given power budget. In [11], the joint active and passive beamforming design problem for the reconfigurable intelligent surface (RIS)-assisted full-duplex (FD) ISAC systems with multiple CUs. The sum rate was maximized by jointly optimizing the transmit beamformer, the linear postprocessing filters, the users’ uplink transmission power, and the reflected phase shift of RIS. In [12], two optimization problems were formulated for the FD ISAC systems. The power consumption was minimized and the sum rate was maximized by jointly optimizing the downlink dual-functional transmit signal, the uplink receive beamformers at the BS, and the transmit power at the uplink users. The authors in [13] investigated the SNR-constrained and CRB-constrained joint beamforming and reflection design problems for the RIS-aided ISAC systems. The beamforming matrix and RIS reflection matrix were jointly optimized to maximize the sum of CUs’ communication rate while considering the SNR/CRB, the BS’s transmission power, and the RIS’s unit-modulus constraints. Due to the possibility of mutual interference between the multiplexed signals, there is a trade-off between target sensing and communication when multi-beam transmitted signals were reused for both sensing and communications. In [14], the angle and the reflection coefficient of the target and the CRB of the angle estimation were utilized as the performance metric of the point target sensing scenario. The authors adopted three performance metrics, the trace, the maximum eigenvalue, and the determinant of the CRB matrix, respectively, to estimate the complete target response matrix with the extended target sensing scenario. The Pareto boundary of the achievable CRB-rate region was given. The energy efficiency (EE) of the ISAC with an extended target and a multi-antenna CU was investigated in [15].  The energy consumption at the transmitter with the on-off scheme was minimized by jointly designing the transmit covariance matrix and the time allocation while taking the QoS for the CU and the CRB constraints for the target into account.

Like other wireless communication systems, ISAC systems also face security issues, and the confident signals are accessible to be wiretapped [4]. Since the transmit signals were reused for both sensing and communications, the security issue becomes imperative, especially when the target is a potential eavesdropper (PE) [16]-[21]. In [16], the radar target was assumed to be a PE, and the SNR at the PE was minimized by designing the radar beam pattern while ensuring the SINR requirement at legitimate users (LUs). Both the scenarios with perfect/imperfect target angle and the channel state information (CSI) were considered, respectively. In [17], destructive interference was utilized to suppress the PE in the ISAC systems, and the SINR of the radar signals at the BS was maximized by jointly designing the transmit and receive beamforming matrices. Assuming the estimation error of the azimuth angle of the PE obeys the Gaussian distribution, the authors in [18] proposed a new CSI error model. Both the scenarios with the perfect eavesdropping CSI/a single LU and with imperfect eavesdropping CSI/multiple LUs were considered. The sum secrecy rate of the considered ISAC system was maximized by optimizing a beamforming vector while considering the CRB and QoS constraints. In [19], the scenarios with perfect and imperfect CSI for the ISAC with multiple CUs and PE were considered. The weighted sum of the mean square error (MSE) between the designed and desired beampatterns and the mean-squared cross-correlation pattern was minimized by jointly optimizing the covariance matrix of the transmitted signals. Ref. [20] investigated the joint secure transmit beamforming designs for ISAC systems with a single LU and multiple PEs. The weighted sum of beampattern matching MSE and cross-correlation patterns was minimized by jointly optimizing the information and the sensing covariance matrix. The scenario with bounded CSI errors and Gaussian CSI errors were considered, and the worst-case secrecy rate and the secrecy outage constraints were considered, respectively. In [21], the authors proposed a novel sensing-aided secure scheme for the ISAC systems. The PEs’ direction was estimated using the combined capon and approximate maximum likelihood technique. Then, the dual-functional beamforming and the artificial noise (AN) matrix were jointly optimized to maximize the normalized weighted sum of the secrecy rate and estimated CRB under the estimated angle errors and the system power budget constraints.

In ISAC systems, the coexisted strong radar signals can be utilized as inherent jamming signals to fight the external eavesdroppers and improve the transmission security [4]. The authors in [22] investigated the joint designs of secure transmit beamforming for ISAC systems with multiple eavesdroppers. For the scenario with the eavesdropping CSI, the maximum eavesdropping SINR of multiple eavesdroppers was minimized by jointly designing the transmit beamformers of communication and radar systems while satisfying the QoS of LUs, the SINR constraint of radar target echo, and the transmit power budget. For the scenario without the eavesdropping CSI, an AN-aided transmit beamforming scheme was proposed. The total transmit power of the BS and radar was minimized by jointly designing the transmit beamformer and the covariance matrices of the AN vector at the BS and radar under the SINR constraints of LUs communication and radar target echoes. In [23], an unmanned aerial vehicle (UAV) was utilized as the RIS-aided ISAC base station (BS) and a secure transmission scheme was proposed to maximize the average achievable rate or the EE by jointly designing the transmit power allocation, the scheduling of users and targets, the phase shifts at RIS, as well as the trajectory and velocity of the UAV.

I-B Motivation and Contributions

The above works implement communication and sensing functions within the ISAC system. To the best of the author’s knowledge, no open literature addresses the joint target sensing and proactive eavesdropping (JTSAPE) system. In this work, we try to answer the following questions: How to design the beamforming in the JTSAPE system? How is the effect of joint proactive eavesdropping and target sensing? A trade-off between parameter estimation CRB and eavesdropping SNR is achieved by a weighted design.

For clarity, the contributions of this work are listed as follows.

  1. 1.

    We consider the scenario in which a multi-antenna JTSAPE BS actively wiretaps the illegal links while sensing the illegitimate transmitter. The beam transmitted by the BS to sense the target is also utilized as an AN to interfere with illegal links and enable proactive eavesdropping. We optimize the transmitting beam of the BS to maximize the eavesdropping SNR or minimize the CRB, respectively. The semidefinite relaxation (SDR) method was employed to solve the proactive eavesdropping subproblem and target sensing subproblem after dropping the rank one constraint.

  2. 2.

    To realize the trade-off between proactive eavesdropping and target sensing, the joint optimization of proactive eavesdropping and target sensing is investigated, and the normalized weighted optimization problem is formulated. To address the complexity of the original problem, the formulated problem is decomposed into two subproblems: proactive eavesdropping and target sensing, which are solved by the SDR technique.

  3. 3.

    Relative to [5] and [9] in which the CRB is solely considered as an optimization objective to maximize sensing performance, a trade-off between sensing CRB and eavesdropping SNR is achieved in this work. Relative to [16], [21], and [22] in which the AN is considered, the radar waves are regarded as AN in this work.

I-C Organization

The remaining organization of this paper is summarized as follows: Sect. II introduces the system model of the joint proactive eavesdropping and target sensing communication. In Sect. III, an efficient tradeoff algorithm is proposed to solve the formulated problem based on SDR and sequential rank-one constraint relaxation (SROCR) techniques. In Sect. IV, the simulation results of the algorithm are analyzed. Finally, Sect. V summarizes the work of the paper.

Notations: Vectors and matrices are represented in lowercase and uppercase boldface letters, respectively. 𝐈\bf I denotes an identity matrix with appropriate dimensions. ()1{\left(\cdot\right)^{-1}} and ()T{\left(\cdot\right)^{T}} represents the inverse and transpose of the matrix, respectively. (){\left(\cdot\right)^{*}} and ()H{\left(\cdot\right)^{H}} represent the conjugate and conjugate transpose of the matrix, respectively. tr()tr\left(\cdot\right) represents the trace of a matrix. ()(i,j)\left(\cdot\right)_{\left(i,j\right)}denoted the (i,j)\left(i,j\right)-th element of the matrix. \odot denotes Hadamard product. Re(){\mathop{\rm Re}\nolimits}\left(\cdot\right) and Im(){\mathop{\rm Im}\nolimits}\left(\cdot\right) signify the real and the imaginary part of a complex-valued matrix, respectively.

II System Model

Refer to caption
Figure 1: System model consisting of an illegal source (S)\left(S\right), a destination (D)\left(D\right), and an ISAC BS (E)\left(E\right).

As shown in Fig. 1, we consider a JTSAPE system consists of a single-antenna illegal transmission source node SS, a single-antenna destination node DD, and a multi-antenna BS EE. There is an illegal link for adaptive transmission between SS and DD [24], [25], [26], and EE transmits signals to detect SS and eavesdrops the information that SS sends to DD. It is assumed that EE is equipped with uniform linear array (ULA) antennas with NtN_{t} transmitting and NrN_{r} receiving antennas and the angle of EE is known a priori based on the results of the detection stage, like many works [18], [20].

The objective of this work is to estimate the target parameter and implement proactive eavesdropping. To evaluate the estimation accuracy, the CRB and the eavesdropping SINR are evaluated as the performance metrics of the radar and proactive eavesdropping, respectively.

II-A The Radar Model

The transmitted signal, 𝐱ENt×1{{\bf x}_{E}\in\mathbb{C}^{N_{t}\times 1}}, at EE is expressed as

𝐱E=𝐰s,\displaystyle{\bf x}_{E}={\bf{w}}{s}, (1)

where 𝐰Nt×1{\bf w}\in\mathbb{C}^{N_{t}\times 1} is the transmit beamformer to be designed, s{s} denotes the symbols with unit power. It assumes that the data symbols are approximately orthogonal to each other, which denotes 𝔼{|s|2}=1{\mathbb{E}}\left\{{{{\left|s\right|}^{2}}}\right\}=1.

Like [12, 22], since the radar host generally has a strong signal separation capability, and a variety of methods can be used for filtering. Then, EE is assumeded that can completely separate the signal and eliminate the noise. When there are KK targets, the echo received by the radar is expressed as

𝐘(𝜽)=𝐀r(𝜽)𝜷𝑨tT(𝜽)𝐱𝐄+𝐙,\displaystyle{\bf{Y}}\left({\bm{\theta}}\right)={{\bf{A}}_{r}}\left({\bm{\theta}}\right){\bm{\beta A}}_{t}^{T}\left({\bm{\theta}}\right){\bf{x}_{E}}+{\bf{Z}}, (2)

where 𝐀r=[𝐚r(θ1),,𝐚r(θK)]{{\bf{A}}_{r}}=\left[{{{\bf{a}}_{r}}\left({{\theta_{1}}}\right),...,{{\bf{a}}_{r}}\left({{\theta_{K}}}\right)}\right], 𝐀t=[𝐚t(θ1),,𝐚t(θK)]{{\bf{A}}_{t}}=\left[{{{\bf{a}}_{t}}\left({{\theta_{1}}}\right),...,{{\bf{a}}_{t}}\left({{\theta_{K}}}\right)}\right], 𝜷=diag([β1,,βK]T){\bm{\beta}}=diag\left({{{\left[{{\beta_{1}},...,{\beta_{K}}}\right]}^{T}}}\right), 𝜽=[θ1,,θK]T{\bm{\theta}}=\left[{{\theta_{1}},...,{\theta_{K}}}\right]^{T}, βi=bRi+jbIi{\beta_{i}}={b_{{R_{i}}}}+j{b_{{I_{i}}}}, and 𝐙CN(0,σ2𝐈){\bf Z}\sim CN\left(0,\sigma^{2}{\bf{I}}\right), 𝐚t(θ)=[1,ej2πdλsin(θ),,ej2πdλ(Nt1)sin(θ)]T{{\bf a}_{t}}\left(\theta\right)=[1,{e^{j\frac{{2\pi d}}{{\lambda}}\sin\left(\theta\right)}},...,{e^{j\frac{{2\pi d}}{{\lambda}}\left({{N_{t}}-1}\right)\sin\left(\theta\right)}}]^{T} and 𝐚r(θ)=[1,ej2πdλsin(θ),,ej2πdλ(Nr1)sin(θ)]T{{\bf a}_{r}}\left(\theta\right)=[1,{e^{j\frac{{2\pi d}}{{\lambda}}\sin\left(\theta\right)}},...,{e^{j\frac{{2\pi d}}{{\lambda}}\left({{N_{r}}-1}\right)\sin\left(\theta\right)}}]^{T} are the steering vector of the transmitting and receiving arrays, dd is the antenna spacing, λ\lambda denotes the wavelength, θ\theta is the angle of the target when the first antenna is utilized as the reference point, σ2\sigma^{2} is noise power, the normalized parameter β=βR+jβI\beta={\beta_{R}}+j{\beta_{I}} depends on the complex amplitude of the target reflection and path loss coefficient of the SDS-D link. According to (2), 𝐘CN(𝐀r(𝜽)𝜷𝐀tH(𝜽)𝐱E,σ2𝐈){\bf Y}\sim CN\left({{{\bf{A}}_{r}}\left(\bm{\theta}\right)\bm{\beta}{\bf{A}}_{t}^{H}\left(\bm{\theta}\right){{\bf x}_{E}},\sigma^{2}{\bf{I}}}\right). The Fisher information matrix (FIM) of a vector (𝝃)\left({\bm{\xi}}\right) is expressed as

𝐅(𝝃)ij=tr[𝐂1(𝝃)𝐂(𝝃)ξi𝐂1(𝝃)𝐂(𝝃)ξj]+2Re[𝝁H(𝝃)ξi𝐂1(𝝃)𝝁(𝝃)ξj],\displaystyle{\bf F}{\left({\bm{\xi}}\right)_{{{ij}}}}={\rm{tr}}\left[{{{\bf{C}}^{-1}}\left({\bm{\xi}}\right)\frac{{\partial{\bf{C}}\left({\bm{\xi}}\right)}}{{\partial{\xi_{i}}}}{{\bf{C}}^{-1}}\left({\bm{\xi}}\right)\frac{{\partial{\bf{C}}\left({\bm{\xi}}\right)}}{{\partial{\xi_{j}}}}}\right]+2{\mathop{\rm Re}\nolimits}\left[{\frac{{\partial{{\bm{\mu}}^{H}}\left({\bm{\xi}}\right)}}{{\partial{\xi_{i}}}}{{\bf{C}}^{-1}}\left({\bm{\xi}}\right)\frac{{\partial{\bm{\mu}}\left({\bm{\xi}}\right)}}{{\partial{\xi_{j}}}}}\right], (3)

where 𝐂(𝝃){{\bf{C}}\left({\bm{\xi}}\right)} and 𝝁(𝝃){{\bm{\mu}}\left({\bm{\xi}}\right)} denote the covariance matrix and mean vector [27]. In this work, the target direction 𝜽\bm{\theta}, the real and imaginary parts of the coefficient 𝜷\bm{\beta} are estimated. Since σ2𝐈\sigma^{2}{\bf{I}} is independent of the estimator [8, 21, 30], therefore, the FIM for 𝜻=[𝜽,𝜷R,𝜷I]{\bm{\zeta}}=\left[{{{\bm{\theta}}},{\bm{\beta}_{R}},{{{\bm{\beta}}}_{I}}}\right] is given in the following Lemma.

Lemma 1.

The FIM for 𝛇{\bm{\zeta}} is expressed as

𝐅=2σ2[Re(𝐅11)Re(𝐅12)Im(𝐅12)ReT(𝐅12)Re(𝐅22)Im(𝐅22)ImT(𝐅12)ImT(𝐅22)Re(𝐅22)],\displaystyle{\bf{F}}=\frac{2}{{{\sigma^{2}}}}\left[{\begin{array}[]{*{20}{c}}{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{11}}}\right)}&{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{12}}}\right)}&{-{\mathop{\rm Im}\nolimits}\left({{{\bf{F}}_{12}}}\right)}\\ {{{{\mathop{\rm Re}\nolimits}}^{T}}\left({{{\bf{F}}_{12}}}\right)}&{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{22}}}\right)}&{-{\mathop{\rm Im}\nolimits}\left({{{\bf{F}}_{22}}}\right)}\\ {-{{{\mathop{\rm Im}\nolimits}}^{T}}\left({{{\bf{F}}_{12}}}\right)}&{-{{{\mathop{\rm Im}\nolimits}}^{T}}\left({{{\bf{F}}_{22}}}\right)}&{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{22}}}\right)}\end{array}}\right], (7)

where F11=(𝐀rH𝐀˙r)(𝛃𝐀˙tH𝐑𝐀t𝛃)+(𝐀rH𝐀r)(𝛃𝐀˙tH𝐑𝐀˙t𝛃)+(𝐀˙rH𝐀˙r)(𝛃𝐀tH𝐑𝐀t𝛃)+(𝐀˙rH𝐀r)(𝛃𝐀tH𝐑𝐀˙t𝛃){F_{11}}=\left({{\bf{A}}_{r}^{H}{{\dot{\bf{A}}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)+\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\dot{\bf{A}}}_{t}}{\bm{\beta}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\dot{\bf{A}}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\dot{\bf{A}}}_{t}}{\bm{\beta}}}\right), F12=(𝐀rH𝐀r)(𝛃𝐀˙tH𝐑𝐀t)+(𝐀˙rH𝐀r)(𝛃𝐀tH𝐑𝐀t){F_{12}}=\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right), F22=(𝐀rH𝐀r)(𝐀tH𝐑𝐀t){F_{22}}=\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right), 𝐀˙r=[𝐚r(θ1)θ1,,𝐚r(θK)θK]{{\dot{\bf{A}}}_{r}}=\left[{\frac{{\partial{{\bf{a}}_{r}}\left({{\theta_{1}}}\right)}}{{\partial{\theta_{1}}}},...,\frac{{\partial{{\bf{a}}_{r}}\left({{\theta_{K}}}\right)}}{{\partial{\theta_{K}}}}}\right], 𝐀˙t=[𝐚t(θ1)θ1,,𝐚t(θK)θK]{{\dot{\bf{A}}}_{t}}=\left[{\frac{{\partial{{\bf{a}}_{t}}\left({{\theta_{1}}}\right)}}{{\partial{\theta_{1}}}},...,\frac{{\partial{{\bf{a}}_{t}}\left({{\theta_{K}}}\right)}}{{\partial{\theta_{K}}}}}\right], and 𝐑=𝐰𝐰H{\bf{R}}={\bf w}{{\bf w}^{H}}.

In this work, it is assumed that there is one target. The singal target echo sensing signal received by EE is expressed as [19]

𝐘r=𝐚r(θ)β𝐚tT(θ)𝐱E+𝐧r,\displaystyle{\bf Y}_{r}={{\bf a}_{r}}\left(\theta\right)\beta{{\bf a}_{t}^{T}}\left(\theta\right){{\bf x}_{E}}+{\bf{n}}_{r}, (8)

where 𝐧rCN(0,σ2𝐈){\bf{n}}_{r}\sim CN\left(0,\sigma^{2}{\bf{I}}\right) is the additive white Gaussian noise (AWGN) matrix.

Proof.

Please refer to Appendix A. ∎

The CRB matrix is obtained as

𝚽(𝜻)=(𝐅(𝜻))1.\displaystyle{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)=\left({{\bf F}\left({{\bm{\zeta}}}\right)}\right)^{-1}. (9)

II-B The Proactive Eavesdropping Model

It is assumed that a radar signal is used as AN to realize proactive eavesdropping, and the line-of-sight channel between EE and DD is expressed as [18]

𝐡ED=β0dEDα𝐚tT(θD),\displaystyle{{\bf h}_{ED}}=\sqrt{\frac{{{\beta_{0}}}}{{d_{ED}^{\alpha}}}}{{\bf a}_{t}^{T}}\left({{\theta_{D}}}\right), (10)

where β0\beta_{0} represents the channel gain at the reference distance, dEDd_{ED} represents the distance between EE and DD, and α\alpha represents the path loss exponent. The channel between SS and DD is expressed as

hSD=β0dSDα,\displaystyle{h_{SD}}=\sqrt{\frac{{{\beta_{0}}}}{{d_{SD}^{\alpha}}}}, (11)

where dSDd_{SD} represents the distance between SS and DD. The radar signal transmitted from EE also acts as an AN singnals to reduce the illegal link transmission rate. Thus, the received signal at DD is expressed as

yD=hSDxS+𝐡ED𝐱E+nD,\displaystyle{{y}_{D}}={h_{SD}}{x_{S}}+{{\bf h}_{ED}}{\bf x}_{E}+{n_{D}}, (12)

where xSx_{S} denotes the signal sent by SS, and nDN(0,σD2){n_{D}}\sim N\left({0,\sigma_{D}^{2}}\right) is the AWGN. The SNR of received signal at DD is expressed as

γD\displaystyle{\gamma_{D}} =𝔼|hSDxS|2𝔼|𝐡ED𝐱E+nD|2\displaystyle=\frac{\mathbb{E}{{{\left|{{h_{SD}}{x_{S}}}\right|}^{2}}}}{\mathbb{E}{{{\left|{{{\bf h}_{ED}}{{\bf{x}}_{E}}+{n_{D}}}\right|}^{2}}}}
=PS|hSD|2βEdEDα(θD)+σD2,\displaystyle=\frac{{{P_{S}}{{\left|{{h_{SD}}}\right|}^{2}}}}{{{\beta_{E}}d_{ED}^{-\alpha}\mathcal{I}\left({{\theta_{D}}}\right)+\sigma_{D}^{2}}}, (13)

where PSP_{S} is the transmit power of SS, and (θD)=𝐚tH(θD)𝐑𝐚t(θD)\mathcal{I}\left({{\theta_{D}}}\right)={\bf{a}}_{t}^{H}\left({{\theta_{D}}}\right){\bf{R}}{{\bf{a}}_{t}}\left({{\theta_{D}}}\right) signifies the beampattern from EE.

The communication signals received at EE is expressed as

𝐲E=𝐡SExS+𝐧E,\displaystyle{\bf{y}}_{E}={{\bf{h}}_{SE}}{x_{S}}+{\bf{n}}_{E}, (14)

where 𝐧EN(0,σE2𝐈){\bf{n}}_{E}\sim N\left({0,\sigma_{E}^{2}{\bf{I}}}\right), 𝐡SENr×1{{\bf h}_{SE}\in\mathbb{C}^{N_{r}\times 1}} is the SES-E channel. It is assumed that EE adopts the selection combining scheme, the SNR of EE is expressed as

γE=PS|hSE|2σE2,\displaystyle{\gamma_{E}}=\frac{{{P_{S}}{{\left|{{h_{SE}}}\right|}^{2}}}}{{\sigma_{E}^{2}}}, (15)

where hSE=β0dSEαh_{SE}=\sqrt{\frac{{{\beta_{0}}}}{{d_{SE}^{\alpha}}}}, dSEd_{SE} represents the distance between SS and EE.

III Problem Formulation

In this section, we first consider the covariance matrix design for proactive eavesdropping and target sensing, respectively. Then, the joint optimization of proactive eavesdropping and target sensing is considered. A normalized weighted optimization problem is proposed to examine the trade-off between the target parameter estimation and proactive eavesdropping.

III-A Proactive Eavesdropping Only

In this work, it is assumed that there is an adaptive information transmission link between SS and DD, EE can decode the received information without error, and the success of proactive eavesdropping needs to meet γDγE{\gamma_{D}}\leq{\gamma_{E}} [24, 25, 26]. In this subsection, proactive eavesdropping is considered only, which indicates that the transmitting beam is designed to meet the eavesdropping demand. The problem of proactive eavesdropping is formulated as

𝒫1:\displaystyle\mathcal{P}_{1}:\; max𝐑γD\displaystyle\mathop{\max}\limits_{\bf{R}}{\gamma_{D}}
s.t.\displaystyle{\mathrm{s.t.}} γDγE,\displaystyle{\gamma_{D}}\leq{\gamma_{E}}, (16a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (16b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (16c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (16d)

where (16a) is the premise of successful eavesdropping, (16b) is the constraint on the total power, where P0P_{0} denotes the power budget of EE, (16c) indicates that the covariance matrix 𝐑\bf{R} is a positive semi-definite matrix, and (16d) indicates that the rank of 𝐑\bf{R} is one.

𝒫1\mathcal{P}_{1} is challengeable to solve because the convexity of (16a) is unknown and (16d) is non-convex. By substituting (13) and (15) into (16a), which is rewritten as

dEDαβE(hSD2σE2hSE2σD2)(θD).\displaystyle\frac{{d_{ED}^{\alpha}}}{{{\beta_{E}}}}\left({\frac{{h_{SD}^{2}\sigma_{E}^{2}}}{{h_{SE}^{2}}}-\sigma_{D}^{2}}\right)\leq\mathcal{I}\left({{\theta_{D}}}\right). (17)

Based on (13), it can be found that maximing γD{\gamma_{D}} equals to minimizing (θD)\mathcal{I}\left({{\theta_{D}}}\right). Then, 𝒫1\mathcal{P}_{1} is equivalent to

𝒫1.1:\displaystyle\mathcal{P}_{1.1}:\; min𝐑(θD)\displaystyle\mathop{\min}\limits_{\bf{R}}\mathcal{I}\left({{\theta_{D}}}\right)
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, dEDαβE(hSD2σE2hSE2σD2)(θD)\displaystyle\frac{{d_{ED}^{\alpha}}}{{{\beta_{E}}}}\left({\frac{{h_{SD}^{2}\sigma_{E}^{2}}}{{h_{SE}^{2}}}-\sigma_{D}^{2}}\right)\leq\mathcal{I}\left({{\theta_{D}}}\right) (18a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (18b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (18c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (18d)

where (θD)\mathcal{I}\left({{\theta_{D}}}\right) is non-negative. With the same method as [16], igoring (18d), 𝒫1.1\mathcal{P}_{1.1} can be solved directly by the convex optimization toolbox.

III-B Target Sensing Only

CRB is the lower bound of the unbiased estimation variance, and therefore minimizing CRB is a strategy for improving the accuracy of the estimation. In this subsection, target sensing is considered only, which indicates that the transmitting beam is designed to minimize the CRB.

III-B1 Pefect Target Direction

In this subsection, it is assumed that the given target angle is perfect. The formulated problem is expressed as

𝒫2.1:\displaystyle\mathcal{P}_{2.1}:\; min𝐑|𝚽(𝜻)|\displaystyle\mathop{\min}\limits_{\bf{R}}{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, (θS)(θn)γs,θnΩ,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega}, (19a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (19b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (19c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (19d)

where θS\theta_{S} denotes the direction of SS, θn\theta_{n} prescribes the sidelobe region, the left-hand side (LHS) of (19a) is the difference between the gain of the target direction and the gain of other directions, Ω{\Omega} represents other angles in the detection range except the target angle θS\theta_{S}, and γs\gamma_{s} is the constant threshold. The purpose of (19a) is to ensure that the gain in the target direction is superior to that in the other directions. The objective function of 𝒫2.1\mathcal{P}_{2.1} is to minimize the determinant of 𝚽(𝜻){{\bm{\Phi}}}\left({{\bm{\zeta}}}\right), which conforms to the convex limit. 𝒫2.1\mathcal{P}_{2.1} is a standard SDR problem after ignoring the rank-1 constraint and can be solved by CVX toolbox.

III-B2 Impefect Target Direction

In this subsection, more practical scenario is considered wherein the given target angle θS{\theta_{S}} is imperfect and there is an uncertain angle Δθ\Delta\theta. Specifically, the angle of the radar’s main detection is changed in Ω1=[θSΔθ,θS+Δθ]{\Omega_{1}}=\left[{{\theta_{S}}-\Delta\theta,{\theta_{S}}+\Delta\theta}\right]. Then, the following optimization problem is formulated

𝒫2.2:\displaystyle\mathcal{P}_{2.2}:\; min𝐑|𝚽(𝜻)|\displaystyle\mathop{\min}\limits_{\bf{R}}{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, (θS)(θn)γs,θnΩ2,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega_{2}}, (20a)
(θm)(1+φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\leq\left({1+\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (20b)
(θm)(1φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\geq\left({1-\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (20c)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (20d)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (20e)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (20f)

where Ω2\Omega_{2} represents other angles except the main lobe direction, φ\varphi is constraint constant, 0φ10\leq\varphi\leq 1. (20a) represents the gain gap threshold constraint between the target direction and side lobes, (20b) and (20c) are both stable mainlobe gain intervals to prevent a large gain difference which affects the detection performance. 𝒫2.2\mathcal{P}_{2.2} under an uncertain target sensing angle can be solved by CVX after ignoring (20f).

III-C Jointly Optimization of Proactive Eavesdropping and Target Sensing

In this section, a normalized weighted optimization problem that reveals the performance tradeoff between proactive eavesdropping and target sensing is proposed. The design of a weighted optimization between the radar CRB and the eavesdropping rate presents the challenge that the two performance metrics have different units and potentially different magnitudes. Inspired by [21], the CRB and the eavesdropping rate are normalized with their respective bound obtained in 𝒫1.1\mathcal{P}_{1.1} and 𝒫2.1\mathcal{P}_{2.1}, respectively.

III-C1 Joint Optimization Proactive Eavesdropping and Target Sensing

To achieve a tradeoff between the eavesdropping performance and the target sensing, while taking the uncertain angle of SS into account, the weighted optimization problem is formulated as

𝒫3:\displaystyle\mathcal{P}_{3}:\; min𝐑ρ(θD)~+(1ρ)|𝚽(𝜻)||𝚽~|\displaystyle\mathop{\min}\limits_{\bf{R}}\rho\frac{{\mathcal{I}\left({{\theta_{D}}}\right)}}{{{{\rm{\tilde{\mathcal{I}}}}}}}+\left({1-\rho}\right)\frac{{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}}{{\left|{{{\tilde{\bm{\Phi}}}}}\right|}}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (21a)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (21b)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (21c)
dEDαβE(hSD2σE2hSE2σD2)(θD),\displaystyle\frac{{d_{ED}^{\alpha}}}{{{\beta_{E}}}}\left({\frac{{h_{SD}^{2}\sigma_{E}^{2}}}{{h_{SE}^{2}}}-\sigma_{D}^{2}}\right)\leq\mathcal{I}\left({{\theta_{D}}}\right), (21d)
(θS)(θn)γs,θnΩ2,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega_{2}}, (21e)
(θm)(1+φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\leq\left({1+\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (21f)
(θm)(1φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\geq\left({1-\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (21g)

where ~\mathcal{\tilde{I}} is the result of the optimization problem 𝒫1.1\mathcal{P}_{1.1}, 𝚽~{{{\tilde{\bm{\Phi}}}}} is the result of the optimization problem 𝒫2.1\mathcal{P}_{2.1} (𝒫2.2)\left(\mathcal{P}_{2.2}\right), and ρ\rho signifies the factor that determines the proactive eavesdropping and target sensing performance weights, 0ρ10\leq\rho\leq 1.

𝒫3\mathcal{P}_{3} is an SDR problem by neglecting (21e), and the CVX toolbox can be utilized to solve the problem. The above-solving process is summarized as Algorithm 1.

1:  Input: wS{w_{S}}, wDw_{D}, wEw_{E}, θS\theta_{S}, θD\theta_{D}, NrN_{r}, NtN_{t}, ηs\eta_{s}, φ\varphi
2:  Output: 𝐑\bf{R}
3:  Obtain ~{{{\rm{\tilde{\mathcal{I}}}}}} by solving (𝒫1\mathcal{P}_{1})
4:  Obtain |𝚽~|{\left|{{\bm{\tilde{\Phi}}}}\right|} by solving the problem (𝒫2.1\mathcal{P}_{2.{1}} or 𝒫2.2\mathcal{P}_{2.{2}})
5:  Obtain the optimal solution by solving (𝒫3\mathcal{P}_{3})
Algorithm 1 An algorithm for non-zero jamming power

III-C2 The interference power is zero

It is unnecessary for EE to send noise to DD when the channel quality of SES-E is better than that of SDS-D, i.e. hSD2σD2hSE2σE2\frac{{h_{SD}^{2}}}{{\sigma_{D}^{2}}}\leq\frac{{h_{SE}^{2}}}{{\sigma_{E}^{2}}} in (21d) [29]. In these scenarios, the eavesdropping rate is equal to the receiving rate of DD without interference and 𝒫3\mathcal{P}_{3} is transformed into a target sensing sub-problem, which is expressed as

𝒫3.1:\displaystyle\mathcal{P}_{3.1}:\; min𝐑|𝚽(𝜻)|\displaystyle\mathop{\min}\limits_{\bf{R}}{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, (θD)=0,\displaystyle\mathcal{I}\left({{\theta_{D}}}\right)=0, (22a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (22b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (22c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (22d)
(θS)(θn)γs,θnΩ2,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega_{2}}, (22e)
(θm)(1+φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\leq\left({1+\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (22f)
(θm)(1φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\geq\left({1-\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (22g)

where (22a) denotes the interference from EE at DD equal to zero.

𝒫3.1\mathcal{P}_{3.1} is a standard convex problem when (22d) is omitted, which can be solved directly by CVX. The solving process is summarized as Algorithm 2.

1:  Input: wS{w_{S}}, wDw_{D}, wEw_{E}, θS\theta_{S}, θD\theta_{D}, NrN_{r}, NtN_{t}, ηs\eta_{s}, φ\varphi
2:  Output: 𝐑\bf{R}
3:  Obtain the optimal solution by solving (𝒫3.1\mathcal{P}_{3.{\mathrm{1}}})
Algorithm 2 An algorithm for zero jamming power

III-C3 The Rank-1 Constraint

The rank-1 constraint is non-convex, which means the optimization problem can not be directly solved by the relevant toolbox. The above problems are solved by omitting the rank-1 constraint and the rank of the obtained solution may be larger than one.

Based on Refs. [16] and [21], the feasible suboptimal solution of the original problem can be obtained by utilizing specific methods based on the SDR technology, and rank-1 constraint is solved by Gaussian randomization or singular value decomposition, etc. It must be noted that these methods may result in a certain degree of distortion, which is tolerable in many cases. However, these methods cannot be applied to our system model since the recovered result may not satisfy γDγE{\gamma_{D}}\leq{\gamma_{E}}, which results in failure eavesdropping.

For the obtained solution with rank(𝐑~)>1{\rm{rank}}\left({\bf{\tilde{R}}}\right)>1, an iterative method is utilized, which is expressed as

𝐑^=𝐃^𝚺^𝐃^𝟏,\displaystyle{\bf{\hat{R}}}=\bf{\hat{D}}{\bf{\hat{\Sigma}}}{\bf{\hat{D}}^{-1}}, (23)

where 𝚺^=diag(0,0,,λ){\bf{\hat{\Sigma}}}={\rm diag}\left({0,0,...,\lambda}\right) denotes the diagonal matrix composed of eigenvalues and 𝐃^{\bf{\hat{D}}} is the eigenvector matrix. Firstly, the eigenvalue matrix 𝚺~=diag(λ1,,λNt){\bf{\tilde{\Sigma}}}={\rm diag}\left({{\lambda_{1}},...,{\lambda_{{N_{t}}}}}\right) is obtained by eigenvalue decomposition. According to the relationship between the matrix eigenvalue and the matrix trace, the maximum eigenvalue of 𝐑~{\bf{\tilde{R}}} must be less than the trace. The eigenvector of maximum eigenvalue of 𝐑~{\bf{\tilde{R}}} is expressed as 𝐮{\bf u} and there must be 𝐮H𝐑~𝐮<tr(𝐑~){{\bf u}^{H}}{\bf{\bf{\tilde{R}}}}{\bf u}<{\rm tr}\left({{\bf{\bf{\tilde{R}}}}}\right). Let w=𝐮H𝐑~𝐮tr(𝐑~){w}=\frac{{{{\bf u}^{H}}{\bf{\bf{\tilde{R}}}}{\bf u}}}{{{\rm tr}\left({{\bf{\bf{\tilde{R}}}}}\right)}}, by using the iterative method, the rank-1 constraint is transformed into

𝐮(j)H𝐑𝐮(j)w(j)tr(𝐑),\displaystyle{{\bf u}^{\left({j}\right)}}^{H}{\bf{R}}{{\bf u}^{\left({j}\right)}}\geq{w^{\left({j}\right)}}{\rm{tr}}\left({\bf{R}}\right), (24)

where ()(j){\left(\right)^{\left(j\right)}} denotes jj-th iteration value. Then 𝒫3\mathcal{P}_{3} and 𝒫3.1\mathcal{P}_{3.1} are reformulated as

𝒫3b:\displaystyle\mathcal{P}_{3{\rm b}}:\; min𝐑ρ(θD)~+(1ρ)|𝚽(𝜻)||𝚽~|\displaystyle\mathop{\min}\limits_{\bf{R}}\rho\frac{{\mathcal{I}\left({{\theta_{D}}}\right)}}{{{{\rm{\tilde{\mathcal{I}}}}}}}+\left({1-\rho}\right)\frac{{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}}{{\left|{{{\tilde{\bm{\Phi}}}}}\right|}}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, 𝐮(j)H𝐑𝐮(j)w(j)tr(𝐑),\displaystyle{{\bf u}^{\left({j}\right)}}^{H}{\bf{R}}{{\bf u}^{\left({j}\right)}}\geq{w^{\left({j}\right)}}{\rm{tr}}\left({\bf{R}}\right), (25a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (25b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (25c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (25d)
dEDαβE(hSD2σE2hSE2σD2)(θD)\displaystyle\frac{{d_{ED}^{\alpha}}}{{{\beta_{E}}}}\left({\frac{{h_{SD}^{2}\sigma_{E}^{2}}}{{h_{SE}^{2}}}-\sigma_{D}^{2}}\right)\leq\mathcal{I}\left({{\theta_{D}}}\right) (25e)
(θS)(θn)γs,θnΩ2,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega_{2}}, (25f)
(θm)(1+φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\leq\left({1+\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (25g)
(θm)(1φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\geq\left({1-\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (25h)

and

𝒫3.1b:\displaystyle\mathcal{P}_{3.{\rm 1b}}:\; min𝐑|𝚽(𝜻)|\displaystyle\mathop{\min}\limits_{\bf{R}}{\left|{{\bm{\Phi}}}\left({{\bm{\zeta}}}\right)\right|}
s.t.\displaystyle{\mathrm{s.t.}}\,\,\, 𝐮(j)H𝐑𝐮(j)w(j)tr(𝐑),\displaystyle{{\bf u}^{\left({j}\right)}}^{H}{\bf{R}}{{\bf u}^{\left({j}\right)}}\geq{w^{\left({j}\right)}}{\rm{tr}}\left({\bf{R}}\right), (26a)
tr(𝐑)P0,\displaystyle{\rm{tr}}({\bf{R}})\leq{P_{0}}, (26b)
𝐑0,\displaystyle{\bf{R}}\succeq 0, (26c)
rank(𝐑)=1,\displaystyle{\rm rank}\left({\bf{R}}\right)=1, (26d)
(θS)(θn)γs,θnΩ2,\displaystyle\mathcal{I}\left({{\theta_{S}}}\right)-\mathcal{I}\left({{\theta_{n}}}\right)\geq{\gamma_{s}},\forall{\theta_{n}}\in{\Omega_{2}}, (26e)
(θm)(1+φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\leq\left({1+\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (26f)
(θm)(1φ)(θS),θmΩ1,\displaystyle\mathcal{I}\left({{\theta_{m}}}\right)\geq\left({1-\varphi}\right)\mathcal{I}\left({{\theta_{S}}}\right),\forall{\theta_{m}}\in{\Omega_{1}}, (26g)
(θD)=0,\displaystyle\mathcal{I}\left({{\theta_{D}}}\right)=0, (26h)

respectively. The initial value of w(0){w^{\left(0\right)}} is obtained by solving 𝒫3\mathcal{P}_{3} (𝒫3.1)(\mathcal{P}_{3.{\rm 1}}), and w(0)=𝐮H𝐑~𝐮tr(𝐑~)[0,1]{{w^{\left(0\right)}}}=\frac{{{{\bf u}^{H}}{\bf{\bf{\tilde{R}}}}{\bf u}}}{{{\rm tr}\left({\bf{\tilde{R}}}\right)}}\in\left[{0,1}\right]. The optimal result of 𝒫3b\mathcal{P}_{3{\rm b}} (𝒫3.1b)(\mathcal{P}_{3.{\rm 1b}}) is denoted as opt(j)opt^{(j)} in jj-th iteration, which are solved by an iterative algorithm summarized as Algorithm 3. After obtain 𝐑\bf{R} by solving 𝒫3.b(𝒫3.1b)\mathcal{P}_{3.{\rm b}}(\mathcal{P}_{3.{\rm 1b}}), the beamforming vector, 𝐰{\bf w}, is obtained through the eigenvalue decomposition.

1:  Input: wS{w_{S}}, wDw_{D}, wEw_{E}, θS\theta_{S}, θD\theta_{D}, NrN_{r}, NtN_{t}, ηs\eta_{s}, φ\varphi, ϑ\vartheta, τ\tau, ρ\rho
2:  Output: 𝐑\bf{R}
3:  Initialization: set initial variables; tolerance ε>0\varepsilon>0 and iteration number j=0j=0;
4:  Executing Algorithm 1 or Algorithm 2 to obtain 𝐑{\bf R}, w(0)=𝐮H𝐑𝐮tr(𝐑){w^{(0)}}=\frac{{{{\bf u}^{H}}\bf R{\bf u}}}{{{\rm tr}\left(\bf R\right)}} and δ(0)[0,1𝐮H𝐑𝐮tr(𝐑)]{\delta^{(0)}}\in\left[{0,1-\frac{{{{\bf u}^{H}}{\bf{R}}{{\bf u}}}}{{{\rm tr}\left({\bf{R}}\right)}}}\right]
5:  Repeat
6:  Solve 𝒫3b\mathcal{P}_{3{\rm b}} or 𝒫3.1b\mathcal{P}_{3.{\rm 1b}}
7:  If solvable, obtain 𝐑(j+1){{\bf{R}}^{\left({j+1}\right)}} and δ(j+1)=δ(j){\delta^{\left({j+1}\right)}}={\delta^{\left({j}\right)}}
8:  else 𝐑(j+1)=𝐑(j){{\bf{R}}^{\left({j+1}\right)}}={{\bf{R}}^{\left(j\right)}} and δ(j+1)=12δ(j){\delta^{\left({j+1}\right)}}=\frac{1}{2}{\delta^{\left({j}\right)}}
9:  w(j+1)=min(𝐮(j+1)H𝐑(j+1)𝐮(j+1)tr(𝐑(j+1))+δ(j+1),1){w^{\left({j+1}\right)}}=\min\left({\frac{{{{\bf u}^{\left({j+1}\right)}}^{H}{{\bf{R}}^{\left({j+1}\right)}}{{\bf u}^{\left({j+1}\right)}}}}{{{\rm tr}\left({{{\bf{R}}^{\left({j+1}\right)}}}\right)}}+{\delta^{\left({j+1}\right)}},1}\right)
10:  j=j+1j=j+1
11:  Until the algorithm converges that opt(j+1)opt(j)ϑ\|opt^{(j+1)}-opt^{(j)}\|\leq\vartheta and w(j+1)τw^{\left({j+1}\right)}\geq\tau or reaches the maximum number of iterations
12:  Obtain the optimal solution to the problem
Algorithm 3 An iterative algorithm for joint optimization

IV Numerical Results and Discussion

In this section, simulation results are presented to prove the convergence and effectiveness of our algorithm. The positions of EE, SS, and DD are expressed as wE=[xE,yE]{w_{E}}=\left[{{x_{E}},{y_{E}}}\right], wS=[xS,yS]{w_{S}}=\left[{{x_{S}},{y_{S}}}\right], and wD=[xD,yD]{w_{D}}=\left[{{x_{D}},{y_{D}}}\right], respectively. The positive direction of the X axis is used as the reference direction, and the azimuth angles of SS and DD relative to EE are expressed as θS{\theta_{S}} and θD{\theta_{D}}, respectively. The direction of the radar detection range is [90,90]\left[{-{{90}^{\circ}},{{90}^{\circ}}}\right], radar resolution of 1 degrees, σ2=1\sigma^{2}=1, and other parameters are listed in Table I. In all the simulation figures, ‘PEO,’ ‘TSO,’ and ‘JPT’ correspond to the results for proactive eavesdropping only, target sensing only, and joint optimization of proactive eavesdropping and target sensing, respectively.

TABLE I: Simulation Parameters
  Parameters Value Parameters Value
wS{w}_{S} [500,0]{\left[500,0\right]} wD{w}_{D} [250,250]{\left[250,250\right]}
wE{w}_{E} [0,0]{\left[0,0\right]} θS\theta_{S} 00^{\circ}
θD\theta_{D} 4545^{\circ} β0\beta_{0} -30 dB
σD2=σE2\sigma_{D}^{2}=\sigma_{E}^{2} -80 dBm φ\varphi 0.05
γs\gamma_{s} 0.01 α\alpha 2.7
PSP_{S} 30 dBm ϑ\vartheta 10410^{-4}
τ\tau 0.999 Nr=NtN_{r}=N_{t} 12
 
Refer to caption
(a) Beampatterns for different schemes.
Refer to caption
(b) Beampatterns versus varying θ\theta.
Refer to caption
(c) CRB with varying Δθ\Delta\theta.
Refer to caption
(d) γD\gamma_{D} for varying Δθ\Delta\theta.
Figure 2: The results versus varying Δθ\Delta\theta with ρ=0.5\rho=0.5 and Nt=12N_{t}=12.
Refer to caption
(a) Beampatterns versus varying NtN_{t}.
Refer to caption
(b) CRB versus varying NtN_{t}.
Refer to caption
(c) γD\gamma_{D} for varying NtN_{t}.
Figure 3: The results versus varying NtN_{t} with ρ=0.5\rho=0.5 and Δθ=5\Delta\theta=5.

Fig. 2(a) shows the beampattern gain with different schemes. It can be observed that, in PEO scheme where only active eavesdropping is carried out, in the direction where DD is located (θD=45)\left({\theta_{D}={45^{\circ}}}\right), the beampattern gain needs to meet the corresponding threshold to reduce the received SINR of DD because the quality of eavesdropping channel is worse than that of SS-DD. In the TSO scheme where only the target sensing is considered, the beampattern gain is highest in the direction of the sensing target (θS=0)\left({\theta_{S}={0^{\circ}}}\right). In JPT, where proactive eavesdropping and target sensing are considered simultaneous, the direction of SS needs to have the highest beampattern gain, and the direction of DD must meet the threshold requirement. Fig. 2(b) plots the beam gain corresponding to the uncertainty with the JPT scheme. We can find that when the target angle is perfect (Δθ=0)\left({\Delta\theta=0}\right), the main lobe of the beam is the narrowest. With the increase of Δθ{\Delta\theta}, the main lobe width becomes more expansive and the gain decreases because the energy is spread over a larger area. However, in the direction of DD, the gain must be not less than the requirement. Fig. 2(c) demonstrates the CRB versus varying Δθ\Delta\theta under different schemes. The CRB of PEO is impendent to Δθ\Delta\theta since only proactive eavesdropping is considered. When target sensing is taken into account, the CRB increases with the increase of the uncertainty angle because the main lobe widens, and the beam gain decreases. As P0P_{0} increases, the CRB of JPT is gradually lower than TSO when the target angle is perfect (Δθ=0)\left({\Delta\theta=0}\right), indicating that the estimation effect of JPT is more accurate. This is because the constraint in the DD’s direction increases the main lobe gain. In the high-Δθ\Delta\theta region, the CRB of TSO is lower than that of JPT because the energy allocated to the DD direction gradually decreases with the increase of Δθ\Delta\theta and the beam gain decreases. In addition, when the transmitting power of the radar increases (P0=35)\left({{P_{0}}=35}\right) dBm, the corresponding CRB decreases because the transmitting power of the radar is high, and the gain in the DD’s direction is also increased. Fig. 2(d) plots the γD\gamma_{D} versus varying Δθ\Delta\theta with different schemes. We can find that, since proactive eavesdropping is not considered in the TSO scheme, the beam gain in DD’s direction may be less than the threshold when P0P_{0} is lower, resulting in eavesdropping failure. On the contrary, with the higher P0P_{0}, the gain in DD’s direction is greater than the threshold, and the eavesdropping is successful, but the eavesdropping SINR is relatively low. The DD’s SINR of the PEO and JPT schemes is a constant, which denotes that successful eavesdropping is achieved and the maximum eavesdropping rate is obtained.

Fig. 3(a) shows the relationship between beam gain and NtN_{t} of the JPT scheme. The results show that the resolution of beam gain increases with increasing NtN_{t}. The main lobe gain also increases because the ability to control the energy distribution increases. Fig. 3(b) demonstrates the relationship between the CRB and NtN_{t} of the JPT scheme. The CRB decreases with the increase of NtN_{t}, which signifies that the sensing accuracy increases. This is because the main lobe gain increases as NtN_{t} increases. Fig. 3(c) plots the SINR at DD for different schemes versus varying NtN_{t}. It is observed that the TSO scheme gradually succeeds in eavesdropping, and eavesdropping SINR decreases with the increase of NtN_{t} since only target sensing is considered.

Refer to caption
(a) Beampatterns versus varying ρ\rho.
Refer to caption
(b) CRB versus varying ρ\rho.
Refer to caption
(c) γD\gamma_{D} for varying ρ\rho.
Figure 4: The results versus varying ρ\rho with Nt=12N_{t}=12 and Δθ=5\Delta\theta=5.

Fig. 4(a) shows the beam gain of the JPT scheme with varying ρ\rho. As shown in the figure, when ρ=1\rho=1, the JPT scheme degenerates into the PEO scheme, but there is a small main lobe in the SS’s detection because JPT introduces constraints on the gain of the main lobe. Fig. 4(b) plots the relationship between CRB and ρ\rho for different schemes. The CRB with ρ=1\rho=1 is the largest because of the same reason of Fig. 4(a). Because the energy of the main lobe is independent of the ρ\rho, the CRB is a constant. The CRB of TSO is better than that of JPT because the gain constraint in the DD direction is absent. Fig. 4(c) plots the γD\gamma_{D} versus varying ρ\rho with different schemes. When ρ=0\rho=0 and P0=35P_{0}=35 dBm, the JPT scheme eavesdropping is successful but the eavesdropping rate is low because the JPT scheme degenerates into the TSO scheme with ρ=0\rho=0 and the gain in DD’s direction is greater than the threshold.

Refer to caption
(a) CRB versus varying P0P_{0}.
Refer to caption
(b) γD\gamma_{D} versus varying P0P_{0}.
Figure 5: The results versus varying P0P_{0} with Nt=12N_{t}=12 and Δθ=5\Delta\theta=5.

Fig. 5(a) demonstrates the relationship between CRB and radar power budget for various schemes. It can be observed that the CRB decreases linearly with the increase of P0P_{0}, which denotes that increasing the power budget can improve the accuracy of sensing. The CRB of the TSO scheme tends to be constant in the high-P0P_{0} region, indicating that when the transmitting power reaches a certain level, simply increasing the transmitting power cannot improve the sensing accuracy. Fig. 5(b) plots the γD\gamma_{D} of various schemes versus the varying P0P_{0}. The γD\gamma_{D} of the TSO scheme decreases with the increase of P0P_{0} because only the target sensing is considered, the gain in the direction of the illegal receiver is not constrained, and no proactive eavesdropping is carried out. When ρ=0\rho=0, the eavesdropping SINR of the JPT scheme also decreases with the increasing P0P_{0} because the JPT scheme degenerates into the TSO scheme. The eavesdropping SINR of PEO and JPT schemes with non-zero ρ\rho is a constant, which signifies proactive eavesdropping also is successfully performed.

Refer to caption
Figure 6: The beampatterns with zero interference.

Fig. 6 plots the beam gain when the interference power is equal to zero. At this time, because the quality of the eavesdropping channel is better than that of the illegal communication channel, it is not necessary to transmit power in DD’s direction.

V Conclusion

This paper studied proactive eavesdropping in the JTSAPE system wherein the multi-antenna BS transmits the signal to realize target sensing which was utilized as AN to realize proactive eavesdropping. Two waveform design schemes were proposed to maximize the eavesdropping SNR or minimize the CRB, respectively. Then, a normalized weighted optimization problem was formulated to optimize proactive eavesdropping and target sensing jointly. After omitting the Rank-1 constraint, the non-convex problem was solved by the SDR technique. The SROCR method is utilized to solve continuously and iteratively to obtain Rank-1 beamforming solution.

Appendix A Proof of Lemma 1

Based on (3), the FIM for 𝜻=[𝜽,𝜷R,𝜷I]{\bm{\zeta}}=\left[{{{\bm{\theta}}},{{\bm{\beta}}_{R}},{{{\bm{\beta}}}_{I}}}\right] is expressed as

𝐅\displaystyle{\bf{F}} =[𝐅(𝜻)ij]\displaystyle=\left[{{\bf{F}}{{\left({\bm{\zeta}}\right)}_{ij}}}\right]
=[2σ2Retr[(𝐀r(𝜽)𝜷𝐀tH(θ)𝐱E)Hζi(𝐀r(𝜽)𝜷𝐀tH(𝜽)𝐱E)ζj]]\displaystyle=\left[{\frac{2}{{{\sigma^{2}}}}{\rm Retr}\left[{\frac{{\partial{{\left({{{\bf{A}}_{r}}\left({\bm{\theta}}\right){\bm{\beta}}{\bf{A}}_{t}^{H}\left(\theta\right){{\bf{x}}_{E}}}\right)}^{H}}}}{{\partial{\zeta_{i}}}}\frac{{\partial\left({{{\bf{A}}_{r}}\left({\bm{\theta}}\right){\bm{\beta}}{\bf{A}}_{t}^{H}\left({\bm{\theta}}\right){{\bf{x}}_{E}}}\right)}}{{\partial{\zeta_{j}}}}}\right]}\right]
=[𝐅(𝜽,𝜽)𝐅(𝜽,𝜷R)𝐅(𝜽,𝜷I)𝐅(𝜷R,𝜽)𝐅(𝜷R,𝜷R)𝐅(𝜷R,𝜷I)𝐅(𝜷I,𝜽)𝐅(𝜷I,𝜷R)𝐅(𝜷I,𝜷I)],\displaystyle=\left[{\begin{array}[]{*{20}{c}}{{\bf{F}}\left({{\bm{\theta}},{\bm{\theta}}}\right)}&{{\bf{F}}\left({{\bm{\theta}},{{\bm{\beta}}_{R}}}\right)}&{{\bf{F}}\left({{\bm{\theta}},{{\bm{\beta}}_{I}}}\right)}\\ {{\bf{F}}\left({{{\bm{\beta}}_{R}},{\bm{\theta}}}\right)}&{{\bf{F}}\left({{{\bm{\beta}}_{R}},{{\bm{\beta}}_{R}}}\right)}&{{\bf{F}}\left({{{\bm{\beta}}_{R}},{{\bm{\beta}}_{I}}}\right)}\\ {{\bf{F}}\left({{{\bm{\beta}}_{I}},{\bm{\theta}}}\right)}&{{\bf{F}}\left({{{\bm{\beta}}_{I}},{{\bm{\beta}}_{R}}}\right)}&{{\bf{F}}\left({{{\bm{\beta}}_{I}},{{\bm{\beta}}_{I}}}\right)}\end{array}}\right], (30)

wherein each of the elements is a KKK*K matrix. In particular, the first element is expressed as

𝐅(𝜽,𝜽)=[F(θ1,θ1)F(θ1,θK)F(θi,θj)F(θK,θ1)F(θK,θK)],\displaystyle{\bf{F}}\left({{\bm{\theta}},{\bm{\theta}}}\right)=\left[{\begin{array}[]{*{20}{c}}{F\left({{\theta_{1}},{\theta_{1}}}\right)}&{...}&{F\left({{\theta_{1}},{\theta_{K}}}\right)}\\ {...}&{F\left({{\theta_{i}},{\theta_{j}}}\right)}&{...}\\ {F\left({{\theta_{K}},{\theta_{1}}}\right)}&{...}&{F\left({{\theta_{K}},{\theta_{K}}}\right)}\end{array}}\right], (34)

where F(θi,θj)F\left({{\theta_{i}},{\theta_{j}}}\right) is expressed as

F(θi,θj)\displaystyle F\left({{\theta_{i}},{\theta_{j}}}\right) =2σ2Retr[(𝐀r𝜷𝑨tT𝐱E)Hθi(𝐀r𝜷𝑨tT𝐱E)θj]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[{\frac{{\partial{{\left({{{\bf{A}}_{r}}{\bm{\beta A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}^{H}}}}{{\partial{\theta_{i}}}}\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{\theta_{j}}}}}\right]
=(a)2σ2Retr[(𝐱EHeieiT𝐀˙t𝜷H𝐀rH+𝐱EH𝐀t𝜷HeieiT𝐀˙rH)×(𝐀˙rejejT𝜷𝐀tT𝐱E+𝐀r𝜷𝐀˙tTejejT𝐱E)]\displaystyle\mathop{=}\limits^{\left(a\right)}\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[\left({{\bf{x}}_{E}^{H}{e_{i}}e_{i}^{T}{\dot{\bf{A}}}_{t}^{*}{{\bm{\beta}}^{H}}{\bf{A}}_{r}^{H}+{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{{\bm{\beta}}^{H}}{e_{i}}e_{i}^{T}{\dot{\bf{A}}}_{r}^{H}}\right)\times\left({{{{\dot{\bf{A}}}}_{r}}{e_{j}}e_{j}^{T}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}+{{\bf{A}}_{r}}{\bm{\beta}}{\dot{\bf{A}}}_{t}^{T}{e_{j}}e_{j}^{T}{{\bf{x}}_{E}}}\right)\right]
=2σ2Retr[𝐱EHeieiT𝐀˙t𝜷H𝐀rH𝐀˙rejejT𝜷𝐀tT𝐱E+𝐱EHeieiT𝐀˙t𝜷H𝐀rH𝐀r𝜷𝐀˙tTejejT𝐱E\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[{{\bf{x}}_{E}^{H}{e_{i}}e_{i}^{T}{\bf{\dot{A}}}_{t}^{*}{{\bm{\beta}}^{H}}{\bf{A}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}{e_{j}}e_{j}^{T}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}+{\bf{x}}_{E}^{H}{e_{i}}e_{i}^{T}{\bf{\dot{A}}}_{t}^{*}{{\bm{\beta}}^{H}}{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}{\bm{\beta}}{\bf{\dot{A}}}_{t}^{T}{e_{j}}e_{j}^{T}{{\bf{x}}_{E}}}\right.
+𝐱EH𝐀t𝜷HeieiT𝐀˙rH𝐀˙rejejT𝜷𝐀tT𝐱E+𝐱EH𝐀t𝜷HeieiT𝐀˙rH𝐀r𝜷𝐀˙tTejejT𝐱E]\displaystyle+\left.{{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{{\bm{\beta}}^{H}}{e_{i}}e_{i}^{T}{\bf{\dot{A}}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}{e_{j}}e_{j}^{T}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}+{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{{\bm{\beta}}^{H}}{e_{i}}e_{i}^{T}{\bf{\dot{A}}}_{r}^{H}{{\bf{A}}_{r}}{\bm{\beta}}{\bf{\dot{A}}}_{t}^{T}{e_{j}}e_{j}^{T}{{\bf{x}}_{E}}}\right]
=(b)2σ2Re[eiT(𝐀rH𝐀˙r)ejejT(𝜷𝐀tT𝐱E𝐱EH𝐀˙t𝜷H)ei+eiT(𝐀rH𝐀r)ejejT(𝜷𝐀˙tT𝐱E𝐱EH𝐀˙t𝜷H)ei\displaystyle\mathop{=}\limits^{\left(b\right)}\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[{e_{i}^{T}\left({{\bf{A}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}}\right){e_{j}}e_{j}^{T}\left({{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}{\bf{x}}_{E}^{H}{\bf{\dot{A}}}_{t}^{*}{{\bm{\beta}}^{H}}}\right){e_{i}}+e_{i}^{T}\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right){e_{j}}e_{j}^{T}\left({{\bm{\beta}}{\bf{\dot{A}}}_{t}^{T}{{\bf{x}}_{E}}{\bf{x}}_{E}^{H}{\bf{\dot{A}}}_{t}^{*}{{\bm{\beta}}^{H}}}\right){e_{i}}}\right.
+eiT(𝐀˙rH𝐀˙r)ejejT(𝜷𝐀tT𝐱E𝐱EH𝐀t𝜷H)ei+eiT(𝐀˙rH𝐀r)ejejT(𝜷𝐀˙tT𝐱E𝐱EH𝐀t𝜷H)ei]\displaystyle\left.{+e_{i}^{T}\left({{\bf{\dot{A}}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}}\right){e_{j}}e_{j}^{T}\left({{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{{\bm{\beta}}^{H}}}\right){e_{i}}+e_{i}^{T}\left({{\bf{\dot{A}}}_{r}^{H}{{\bf{A}}_{r}}}\right){e_{j}}e_{j}^{T}\left({{\bm{\beta}}{\bf{\dot{A}}}_{t}^{T}{{\bf{x}}_{E}}{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{{\bm{\beta}}^{H}}}\right){e_{i}}}\right]
=(c)2σ2Re[L(𝐀rH𝐀˙r)ij(𝜷𝐀˙tH𝐑𝐀t𝜷)ij+L(𝐀rH𝐀r)ij(𝜷𝐀˙tH𝐑𝐀˙t𝜷)ij\displaystyle\mathop{=}\limits^{\left(c\right)}\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[{L{{\left({{\bf{A}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{\dot{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)}_{ij}}+L{{\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{\dot{A}}}_{t}^{H}{{\bf{R}}^{*}}{{{\bf{\dot{A}}}}_{t}}{\bm{\beta}}}\right)}_{ij}}}\right.
+L(𝐀˙rH𝐀˙r)ij(𝜷𝐀tH𝐑𝐀t𝜷)ij+L(𝐀˙rH𝐀r)ij(𝜷𝐀tH𝐑𝐀˙t𝜷)ij]\displaystyle\left.{+L{{\left({{\bf{\dot{A}}}_{r}^{H}{{{\bf{\dot{A}}}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)}_{ij}}+L{{\left({{\bf{\dot{A}}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{{\bf{\dot{A}}}}_{t}}{\bm{\beta}}}\right)}_{ij}}}\right]
=2σ2[Re(𝐅11)]ij,\displaystyle=\frac{2}{{{\sigma^{2}}}}{\left[{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{11}}}\right)}\right]_{ij}}, (35)

where F11=(𝐀rH𝐀˙r)(𝜷𝐀˙tH𝐑𝐀t𝜷)+(𝐀rH𝐀r)(𝜷𝐀˙tH𝐑𝐀˙t𝜷)+(𝐀˙rH𝐀˙r)(𝜷𝐀tH𝐑𝐀t𝜷)+(𝐀˙rH𝐀r)(𝜷𝐀tH𝐑𝐀˙t𝜷){F_{11}}=\left({{\bf{A}}_{r}^{H}{{\dot{\bf{A}}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)+\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\dot{\bf{A}}}_{t}}{\bm{\beta}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\dot{\bf{A}}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\dot{\bf{A}}}_{t}}{\bm{\beta}}}\right), [𝐁]ij{\left[{\bf{{B}}}\right]_{ij}} denotes the (i,j)(i,j) element of the matrix 𝐁\bf{{B}}, eie_{i} denotes the ii-th column of the identity matrix, step (a)\left(a\right) holds by following

(𝐀r𝜷𝐀tT𝐱E)θi=𝐀˙reieiT𝜷𝐀tT𝐱E+𝐀r𝜷𝐀˙tTeieiT𝐱E,\displaystyle\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{\theta_{i}}}}={{\dot{\bf{A}}}_{r}}{e_{i}}e_{i}^{T}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}+{{\bf{A}}_{r}}{\bm{\beta}}{\dot{\bf{A}}}_{t}^{T}{e_{i}}e_{i}^{T}{{\bf{x}}_{E}}, (36)

step (b)\left(b\right) holds by following

tr[(𝐱EHeieiT𝐀˙t𝜷H𝐀rH)(𝐀˙rejejT𝜷𝐀tT𝐱E)]\displaystyle{\mathop{\rm tr}\nolimits}\left[{\left({{\bf{x}}_{E}^{H}{e_{i}}e_{i}^{T}{\dot{\bf{A}}}_{t}^{*}{{\bm{\beta}}^{H}}{\bf{A}}_{r}^{H}}\right)\left({{{{\dot{\bf{A}}}}_{r}}{e_{j}}e_{j}^{T}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}\right] =eiT(𝐀rH𝐀˙r)ejejT(𝜷𝐀tT𝐱E𝐱EH𝐀˙t𝜷H)ei\displaystyle=e_{i}^{T}\left({{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right){e_{j}}e_{j}^{T}\left({{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}{\bf{x}}_{E}^{H}{\dot{\bf{A}}}_{t}^{*}{{\bm{\beta}}^{H}}}\right){e_{i}}
=L(𝐀rH𝐀˙r)ij(𝜷𝐀˙tH𝐑𝐀t𝜷)ij,\displaystyle=L{\left({{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right)_{ij}}{\left({{{\bm{\beta}}^{*}}{\dot{\bf{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)_{ij}}, (37)

and step (c)\left(c\right) holds by following

L(𝐀rH𝐀˙r)ij(𝜷𝐀˙tH𝐑𝐀t𝜷)ij=[𝐀rH𝐀˙r𝜷𝐀˙tH𝐑𝐀t𝜷]ij.\displaystyle L{\left({{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right)_{ij}}{\left({{{\bm{\beta}}^{*}}{\dot{\bf{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)_{ij}}={\left[{{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}\odot{{\bm{\beta}}^{*}}{\dot{\bf{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right]_{ij}}. (38)

With some simple algebraic manipulations, the element of 𝐅(𝜽,𝜷R){\bf{F}}\left({{\bm{\theta}},{{\bm{{\bm{\beta}}}}_{R}}}\right) is expressed as

F(θi,bRj)\displaystyle F\left({{\theta_{i}},{b_{{R_{j}}}}}\right) =2σ2Retr[(𝐀r𝜷𝐀tT𝐱E)Hθi(𝐀r𝜷𝐀tT𝐱E)bRj]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[{\frac{{\partial{{\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}^{H}}}}{{\partial{\theta_{i}}}}\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{b_{{R_{j}}}}}}}\right]
=2σ2Re[L(𝐀rH𝐀r)ij(𝜷𝐀˙tH𝐑𝐀t)ij+L(𝐀˙rH𝐀r)ij(𝜷𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[L{{\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\dot{\bf{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}\right.\left.+L{{\left({{\dot{\bf{A}}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}^{t}}}\right)}_{ij}}\right]
=2σ2[Re(𝐅12)]ij,\displaystyle=\frac{2}{{{\sigma^{2}}}}{\left[{{\mathop{\rm Re}\nolimits}\left({{{\bf{F}}_{12}}}\right)}\right]_{ij}}, (39)

where F12=(𝐀rH𝐀r)(𝜷𝐀˙tH𝐑𝐀t)+(𝐀˙rH𝐀r)(𝜷𝐀tH𝐑𝐀t){F_{12}}=\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}\dot{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)+\left({\dot{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right) and (𝐀r𝜷𝐀tT𝐱E)bRj=𝐀rejejT𝐀tT𝐱E\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{b_{{R_{j}}}}}}={{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}.

The element of 𝐅(𝜽,𝜷I){\bf{F}}\left({{\bm{\theta}},{{\bm{{\bm{\beta}}}}_{I}}}\right) is expressed as

F(θi,bIj)\displaystyle F\left({{\theta_{i}},{b_{{I_{j}}}}}\right) =2σ2Retr[(𝐀r𝜷𝐀tT𝐱E)Hθi(𝐀r𝜷𝐀tT𝐱E)bIj]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[{\frac{{\partial{{\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}^{H}}}}{{\partial{\theta_{i}}}}\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{b_{{I_{j}}}}}}}\right]
=2σ2Re[L(j𝐀rH𝐀r)ij(𝜷𝐀˙tH𝐑𝐀t)ij+L(j𝐀˙rH𝐀r)ij(𝜷𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[L{{\left({j{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\dot{\bf{A}}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}\right.\left.+L{{\left({j{\dot{\bf{A}}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{{\bm{\beta}}^{*}}{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}^{t}}}\right)}_{ij}}\right]
=2σ2[Im(𝐅12)]ij,\displaystyle=-\frac{2}{{{\sigma^{2}}}}{\left[{{\mathop{\rm Im}\nolimits}\left({{{\bf{F}}_{12}}}\right)}\right]_{ij}}, (40)

where F22=(𝐀rH𝐀r)(𝐀tH𝐑𝐀t){F_{22}}=\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)\odot\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right) and (𝐀r𝜷𝐀tT𝐱E)bIj=j𝐀rejejT𝐀tT𝐱E\frac{{\partial\left({{{\bf{A}}_{r}}{\bm{\beta}}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}}{{\partial{b_{{I_{j}}}}}}=j{{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}.

The element of 𝐅(𝜷R,𝜽){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{R}},{\bm{\theta}}}\right)} is expressed as

F(bRi,θj)\displaystyle F\left({{b_{{R_{i}}}},{\theta_{j}}}\right) =2σ2Re[L(𝐀rH𝐀r)ij(𝐀tH𝐑𝐀˙t𝜷)ij+L(𝐀rH𝐀˙r)ij(𝐀tH𝐑𝐀t𝜷)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[L{{\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{{\dot{\bf{A}}}}_{t}}{\bm{\beta}}}\right)}_{ij}}\right.\left.+L{{\left({{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)}_{ij}}\right]
=2σ2[ReT(𝐅12)]ij.\displaystyle=\frac{2}{{{\sigma^{2}}}}{\left[{{{{\mathop{\rm Re}\nolimits}}^{T}}\left({{{\bf{F}}_{12}}}\right)}\right]_{ij}}. (41)

The element of 𝐅(𝜷R,𝜷R){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{R}},{{\bm{{\bm{\beta}}}}_{R}}}\right)} is expressed as

F(bRi,bRj)\displaystyle F\left({{b_{{R_{i}}}},{b_{{R_{j}}}}}\right) =2σ2Retr[(𝐱EH𝐀teieiT𝐀rH)(𝐀rejejT𝐀tT𝐱E)]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\mathop{\rm tr}\nolimits}\left[{\left({{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{e_{i}}e_{i}^{T}{\bf{A}}_{r}^{H}}\right)\left({{{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}\right]
=2σ2Re[L(𝐀rH𝐀r)ij(𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}}\left[{L{{\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}}\right]
=2σ2[Re(𝐅22)]ij.\displaystyle=\frac{2}{{{\sigma^{2}}}}{\left[{{\mathop{\rm Re}}\left({{{\bf{F}}_{22}}}\right)}\right]_{ij}}. (42)

The element of 𝐅(𝜷R,𝜷I){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{R}},{{\bm{{\bm{\beta}}}}_{I}}}\right)} is expressed as

F(bRi,bIj)\displaystyle F\left({{b_{{R_{i}}}},{b_{{I_{j}}}}}\right) =2σ2Retr[(𝐱EHeieiT𝐀t𝐀rH)(j𝐀rejejT𝐀tT𝐱E)]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}}{\mathop{\rm tr}}\left[{\left({{\bf{x}}_{E}^{H}{e_{i}}e_{i}^{T}{\bf{A}}_{t}^{*}{\bf{A}}_{r}^{H}}\right)\left({j{{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}\right]
=2σ2Re[L(j𝐀rH𝐀r)ij(𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}}\left[{L{{\left({j{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}}\right]
=2σ2[ImT(𝐅22)]ij.\displaystyle=-\frac{2}{{{\sigma^{2}}}}{\left[{{{{\mathop{\rm Im}}}^{T}}\left({{{\bf{F}}_{22}}}\right)}\right]_{ij}}. (43)

The element of 𝐅(𝜷I,𝜽){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{I}},{\bm{\theta}}}\right)} is expressed as

F(bIi,θj)\displaystyle F\left({{b_{{I_{i}}}},{\theta_{j}}}\right) =2σ2Re[L(j𝐀rH𝐀˙r)ij(𝐀tH𝐑𝐀t𝜷)ij+L(j𝐀rH𝐀˙r)ij(𝐀tH𝐑𝐀t𝜷)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[L{{\left({-j{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)}_{ij}}\right.\left.+L{{\left({-j{\bf{A}}_{r}^{H}{{{\dot{\bf{A}}}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}{\bm{\beta}}}\right)}_{ij}}\right]
=2σ2[ImT(𝐅12)]ij.\displaystyle=-\frac{2}{{{\sigma^{2}}}}{\left[{{{{\mathop{\rm Im}\nolimits}}^{T}}\left({{{\bf{F}}_{12}}}\right)}\right]_{ij}}. (44)

The element of 𝐅(𝜷I,𝜷R){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{I}},{{\bm{{\bm{\beta}}}}_{R}}}\right)} is expressed as

F(bIi,bRj)\displaystyle F\left({{b_{{I_{i}}}},{b_{{R_{j}}}}}\right) =2σ2Retr[(j𝐱EH𝐀teieiT𝐀rH)(𝐀rejejT𝐀tT𝐱E)]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\rm tr}\left[{\left({-j{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{e_{i}}e_{i}^{T}{\bf{A}}_{r}^{H}}\right)\left({{{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}\right]
=2σ2Re[L(j𝐀rH𝐀r)ij(𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[{{{L\left({-j{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}}\right]
=2σ2[ImT(𝐅22)]ij.\displaystyle=-\frac{2}{{{\sigma^{2}}}}{\left[{{{{\mathop{\rm Im}\nolimits}}^{T}}\left({{{\bf{F}}_{22}}}\right)}\right]_{ij}}. (45)

The element of 𝐅(𝜷I,𝜷I){{\bf{F}}\left({{{\bm{{\bm{\beta}}}}_{I}},{{\bm{{\bm{\beta}}}}_{I}}}\right)} is expressed as

F(bIi,bIj)\displaystyle F\left({{b_{{I_{i}}}},{b_{{I_{j}}}}}\right) =2σ2Retr[(j𝐱EH𝐀teieiT𝐀rH)(j𝐀rejejT𝐀tT𝐱E)]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}{\rm tr}\left[{\left({-j{\bf{x}}_{E}^{H}{\bf{A}}_{t}^{*}{e_{i}}e_{i}^{T}{\bf{A}}_{r}^{H}}\right)\left({j{{\bf{A}}_{r}}{e_{j}}e_{j}^{T}{\bf{A}}_{t}^{T}{{\bf{x}}_{E}}}\right)}\right]
=2σ2Re[L(𝐀rH𝐀r)ij(𝐀tH𝐑𝐀t)ij]\displaystyle=\frac{2}{{{\sigma^{2}}}}{\mathop{\rm Re}\nolimits}\left[L{{{\left({{\bf{A}}_{r}^{H}{{\bf{A}}_{r}}}\right)}_{ij}}{{\left({{\bf{A}}_{t}^{H}{{\bf{R}}^{*}}{{\bf{A}}_{t}}}\right)}_{ij}}}\right]
=2σ2[ReT(𝐅22)]ij.\displaystyle=\frac{2}{{{\sigma^{2}}}}{\left[{{{{\mathop{\rm Re}\nolimits}}^{T}}\left({{{\bf{F}}_{22}}}\right)}\right]_{ij}}. (46)

References

  • [1] F. Liu, Y. Cui, C. Masouros, J. Xu, T. X. Han, Y. C. Eldar, and S. Buzzi, “Integrated sensing and communications: Towards dual-functional wireless networks for 6G and beyond,” IEEE J. Sel. Areas Commun., vol. 40, no. 6, pp. 1728-1767, Jun. 2022.
  • [2] F. Dong, F. Liu, Y. Cui, W. Wang, K. Han, and Z. Wang, “Sensing as a service in 6G perceptive networks: A unified framework for ISAC resource allocation,” IEEE Trans. Wireless Commun., vol. 22, no. 5, pp. 3522-3536, May 2023.
  • [3] Z. Wei, H. Qu, Y. Wang, X. Yuan, H. Wu, Y. Du, K. Han, N. Zhang, and Z. Feng, “Integrated sensing and communication signals toward 5G-A and 6G: A survey,” IEEE Internet Things J., vol. 10, no. 13, pp. 11068-11092, Jul. 2023.
  • [4] Z. Wei, F. Liu, C. Masouros, N. Su, and A. P. Petropulu, “Toward multi-functional 6G wireless networks: Integrating sensing, communication, and security,” IEEE Commun. Mag., vol. 60, no. 4, pp. 65-71, Apr. 2022.
  • [5] F. Liu, Y.-F. Liu, A. Li, C. Masouros, and Y. C. Eldar, “Cramér-Rao bound optimization for joint radar-communication beamforming,” IEEE Trans. Signal Process., vol. 70, pp. 240 - 253, Dec. 2022.
  • [6] X. Li, F. Liu, Z. Zhou, G. Zhu, S. Wang, K. Huang, and Y. Gong, “Integrated sensing, communication, and computation over-the-air: MIMO beamforming design,” IEEE Trans. Wireless Commun., vol. 22, no. 8, pp. 5383-5398, Aug. 2023.
  • [7] C. Wen, Y. Huang, L. Zheng, W. J. Liu, and T. N. Davidson, “Transmit waveform design for dual-function radar-communication systems via hybrid linear-nonlinear precoding,” IEEE Trans. Signal Process., vol. 71, pp. 2130-2145, Jun. 2023.
  • [8] Z. Ren, Y. Peng, X. Song, Y. Fang, L. Qiu, L. Liu, D. W. K. Ng, and J. Xu, “Fundamental CRB-rate tradeoff in multi-antenna ISAC systems with information multicasting and multi-target sensing,” IEEE Trans. Wireless Commun., vol. 23, no. 4, pp. 3870-3885, Apr. 2024.
  • [9] Z. Zhao, L. Zhang, R. Jiang, X.-P. Zhang, X. Tang, and Y. Dong, “Joint beamforming scheme for ISAC systems via robust Cramér–Rao bound optimization,” IEEE Wireless Commun. Lett., vol. 13, no. 3, pp. 889-893, Mar. 2024.
  • [10] T. Liu, Y. Guo, L. Lu, and B. Xia, “Waveform design for integrated sensing and communication systems based on interference exploitation,” J. Commun. Inform. Netw., vol. 7, no. 4, pp. 447-456, Dec. 2022.
  • [11] Y. Guo, Y. Liu, Q. Wu, X. Li, and Q. Shi, “Joint beamforming and power allocation for RIS aided full-duplex integrated sensing and uplink communication system,” IEEE Trans. Wireless Commun., vol. 23, no. 5, pp. 4627-4642, May 2024.
  • [12] Z. He, W. Xu, H. Shen, D. W. K. Ng, Y. C. Eldar, and X. You, “Full-duplex communication for ISAC: Joint beamforming and power optimization,” IEEE J. Sel. Areas Commun., vol. 41, no. 9, pp. 2920-2936, Sep. 2023.
  • [13] R. Liu, M. Li, Q. Liu, and A. L. Swindlehurst, “SNR/CRB-constrained joint beamforming and reflection designs for RIS-ISAC systems,” IEEE Trans. Wireless Commun., doi: 10.1109/twc.2023.3341429, 2023.
  • [14] H. Hua, T. X. Han, and J. Xu, “MIMO integrated sensing and communication: CRB-rate tradeoff,” IEEE Trans. Wireless Commun., vol. 23, no. 4, pp. 2839-2854, Apr. 2024.
  • [15] G. Wu, Y. Fang, J. Xu, Z. Feng, and S. Cui, “Energy-efficient MIMO integrated sensing and communications with on-off non-transmission powers,” IEEE Internet Things J., vol. 11, no. 7, pp. 12177-12191, Apr. 2024.
  • [16] N. Su, F. Liu, and C. Masouros, “Secure radar-communication systems with malicious targets: Integrating radar, communications and jamming functionalities,” IEEE Trans. Wireless Commun., vol. 20, no. 1, pp. 83-95, Jan. 2021.
  • [17] N. Su, F. Liu, Z. Wei, Y.-F. Liu, and C. Masouros, “Secure dual-functional radar-communication transmission: Exploiting interference for resilience against target eavesdropping,” IEEE Trans. Wireless Commun., vol. 21, no. 9, pp. 7238-7252, Sept. 2022.
  • [18] H. Jia, X. Li, and L. Ma, “Physical layer security optimization with Cramér-Rao bound metric in ISAC systems under sensing-specific imperfect CSI model,” IEEE Trans. Veh. Technol., doi: 10.1109/tvt.2023.3347527, Jun. 2023.
  • [19] F. Dong, W. Wang, X. Li, F. Liu, S. Chen, and L. Hanzo, “Joint beamforming design for dual-functional MIMO radar and communication systems guaranteeing physical layer security,” IEEE Trans. Green Commun. Netw., vol. 7, no. 1, pp. 537-549, Mar. 2023.
  • [20] Z. Ren, L. Qiu, J. Xu, and D. W. K. Ng, “Robust transmit beamforming for secure integrated sensing and communication,” IEEE Trans. Commun., vol. 71, no. 9, pp. 5549-5564, Sep. 2023.
  • [21] N. Su, F. Liu, and C. Masouros, “Sensing-assisted eavesdropper estimation: An ISAC breakthrough in physical layer security,” IEEE Trans. Wireless Commun., vol. 23, no. 4, pp. 3162-3174, Apr. 2024.
  • [22] J. Chu, R. Liu, M. Li, Y. Liu, and Q. Liu, “Joint secure transmit beamforming designs for integrated sensing and communication systems,” IEEE Trans. Veh. Technol., vol. 72, no. 4, pp. 4778-4791, Apr. 2023.
  • [23] J. Zhang, J. Xu, W. Lu, N. Zhao, X. Wang, and D. Niyato, “Secure tansmission for IRS-aided UAV-ISAC networks,” IEEE Trans. Wireless Commun., doi: 10.1109/twc.2024.3390169, 2024.
  • [24] F. Feizi, M. Mohammadi, Z. Mobini, and C. Tellambura, “Proactive eavesdropping via jamming in full-duplex multi-antenna systems: Beamforming design and antenna selection,” IEEE Trans. Commun., vol. 68, no. 12, pp. 7563-7577, Dec. 2020.
  • [25] D. Xu and H. Zhu, “Jamming-assisted legitimate eavesdropping and secure communication in multicarrier interference networks,” IEEE Syst. J., vol. 16, no. 1, pp. 954-965, Mar. 2022.
  • [26] G. Hu, Z. Li, J. Si, K. Xu, Y. Cai, D. Xu, and N. Al-Dhahir, “Analysis and optimization of STAR-RIS-assisted proactive eavesdropping with statistical CSI,” IEEE Trans. Veh. Technol., vol. 72, no. 5, pp. 6850-6855, May 2023.
  • [27] S. M. Kay, Fundamentals of Statistical Signal Processing: Estimation Theory. Englewood Cliffs, NJ: Prentice-Hall, Inc., 1993.
  • [28] I. Bekkerman and J. Tabrikian, “Target detection and localization using MIMO radars and sonars,” IEEE Trans. Signal Process., vol. 54, no. 10, pp. 3873-3883, Oct. 2006.
  • [29] X. Jiang, H. Lin, C. Zhong, X. Chen, and Z. Zhang, “Proactive eavesdropping in relaying systems,” IEEE Signal Process. Lett., vol. 24, no. 6, pp. 917-921, Jun. 2017.
  • [30] J. Li, L. Xu, P. Stoica, K. W. Forsythe, and D. W. Bliss, “Range compression and waveform optimization for MIMO radar: A Cramér–Rao bound based study,” IEEE Trans. Signal Process., vol. 56, no. 1, pp. 218-232, Jan. 2008.