Temporal-Assisted Dynamic Beampattern Optimization in Integrated Sensing and Communication Systems

Shengcai Zhou, Student Member, IEEE Luping Xiang, Member, IEEE, and Kun Yang, Fellow, IEEE The authors would like to thank the financial support of Quzhou Government (Grant No.: 2023D005), Natural Science Foundation of China (Grant No. 62132004 and Grant No. 62301122), MOST Major Research and Development Project (Grant No.: 2021YFB2900204), Sichuan Science and Technology Program (Grant No.: 2023NSFSC1375). (Corresponding author: Luping Xiang.)Shengcai Zhou is with the Yangtze Delta Region Institute (Quzhou), University of Electronic Science and Technology of China, Quzhou 324003, China, and also the School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China, email: [email protected] Xiang and Kun Yang are with the State Key Laboratory of Novel Software Technology, Nanjing University, Nanjing 210008, China, and School of Intelligent Software and Engineering, Nanjing University (Suzhou Campus), Suzhou, 215163, China, email: [email protected], [email protected].
Abstract

In this paper, an integrated sensing and communication (ISAC) system is investigated. Initially, we introduce a design criterion wherein sensing data acquired from the preceding time slot is employed for instantaneous optimal beamforming in the succeeding time slot, aiming to enhance the communication rate. Subsequently, the development of optimal beamforming is addressed, and a high-caliber suboptimal resolution is derived utilizing successive convex approximation (SCA) techniques combined with the iterative rank minimization (IRM) methodology. Our evaluations, grounded on numerical analyses, reveal that the communication rate of the introduced beamforming strategy surpasses that of conventional omnidirectional sensing and pilot based approaches.

Index Terms:
ISAC, beamforming, sensing-assisted, temporal-assisted.

I Introduction

The integration of sensing and communication (ISAC) emerges as a focal approach [14, 17], since the radio propagation is able to convey information generated from the transmitter and extract information based on scattered echoes [13]. The ISAC system allows a co-designed waveform to transmit information to the communication user while sensing the radar target. Sensing-assisted communication technology, which enables communication and sensing to cooperate to omit special pilot design , has attracted a lot of attention in recent years[1].

The study of [4] proposed a shared deployment scheme that directly utilizes communication signals as dual-function waveforms, addressing the additional power loss issue caused by separate transmission of radar and communication signals [3]. In addition to the simple beamforming gain metrics [4], [16] also investigated the trade-offs between minimizing the Cramer-Rao Bound (CRB) and trajectory optimization. A further approach involves directly activating different numbers of transmitting antennas at the transmitter to achieve dynamic beampattern with varying beamwidths [15]. Simultaneously, researchers in the Vehicle-to-Infrastructure (V2I) domain studied the sensing-assisted beam tracking rooted in radar echo data, bypassing the traditionally high-overhead pilot-based channel estimation [5, 6]. Initial V2I surveys used beam training to extract precise channel details, yet they still depended on pilot signals and failed to integrate radar-communication combined signals into a cohesive system. Furthermore, [20] describes some popular framework protocol designs for UAVs, discussing performance variations based on factors such as quality of service requirements, target localization, and maneuverability.

This study delves into an optimal beamforming protocol supported by the ISAC framework. Contrary to prior research, our proposed A2G model employs an integrated waveform for persistent environmental sensing, subsequently employing this data to enhance communication. Additionally, we examine a optimization problem about the ramifications of imprecise sensing data on communication rates, manifested in beamforming through beam coverage breadth. An overview of this study’s merits is encapsulated in Table I. Key contributions from this research include:

TABLE I: Contrasting Our Contributions To The State-Of-Art
Contributions this work [3] [4] [5, 6, 19] [2]
Continuous wave
Integrated signal
Beam coverage
Sensing-assisted
  • A sensing-assisted communication protocol tailored for the ISAC system is introduced. This protocol designs the forthcoming transmit beampattern employing the channel data acquired from the preceding time slot to enhance the communication rate.

  • We employ the successive convex approximation (SCA) method [11] to reframe the presented non-convex optimization challenge as a convex optimization problem. Subsequently, the iterative rank minimization (IRM) algorithm [12] is applied to derive high-caliber solutions.

  • Our simulation outcomes highlight the superiority of the devised protocol. When contrasted against the traditional omnidirectional sensing-based method and the pilot-based approach, there’s an observable enhancement in throughput by 5.7%percent\%% and 8.7%percent\%%, respectively.

The rest of this paper is organised as follows. Section II introduces the system model, Section III formulates the beam optimization problem. Section IV provides numerical results, and Section V summarizes this paper.

II System Model

Refer to caption


Figure 1: System model.

As illustrated in Fig. 1, an ISAC system is considered, comprising an aerial base station (ABS) equipped with a transmit uniform linear array (ULA) encompassing Ntsubscript𝑁𝑡N_{t}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT antennas and a distinct receive ULA consisting of Nrsubscript𝑁𝑟N_{r}italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT antennas. This ABS is adept at sensing proximate single-antenna IoT devices and facilitating downlink communications to them. Without loss of generality, it is postulated that the ABS senses and communicates with the IoT devices exclusively over the line of sight (LoS) channels.

For the described system, it is assumed that the protocol endures for N𝑁Nitalic_N = 10 time intervals, each interval length is ΔTΔ𝑇\Delta Troman_Δ italic_T. Initially, ABS has imprecise information about K𝐾Kitalic_K devices in its vicinity, alongside an unknown number of devices without channel information, so it configures a multi-beam ISAC signal 𝐬(t)𝐬𝑡\mathbf{s}(t)bold_s ( italic_t ), having K𝐾Kitalic_K dimensions, expressed as

𝐬n(t)=[s1,n(t),,sK,n(t)]T,subscript𝐬𝑛𝑡superscriptsubscript𝑠1𝑛𝑡subscript𝑠𝐾𝑛𝑡𝑇\mathbf{s}_{n}(t)=[s_{1,n}(t),...,s_{K,n}(t)]^{T},bold_s start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) = [ italic_s start_POSTSUBSCRIPT 1 , italic_n end_POSTSUBSCRIPT ( italic_t ) , … , italic_s start_POSTSUBSCRIPT italic_K , italic_n end_POSTSUBSCRIPT ( italic_t ) ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT , (1)

where the k𝑘kitalic_k-th signal sk,n(t)subscript𝑠𝑘𝑛𝑡s_{k,n}(t)italic_s start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) includes the communication information of the k𝑘kitalic_k-th device at time slot n𝑛nitalic_n, with independent circularly symmetric complex Gaussian (CSCG) random variables with zero mean and unit variance, represented by sk,n(t)𝒞𝒩(0,1)similar-tosubscript𝑠𝑘𝑛𝑡𝒞𝒩01s_{k,n}(t)\sim\mathcal{CN}(0,1)italic_s start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) ∼ caligraphic_C caligraphic_N ( 0 , 1 ).

Subsequent to the beamforming process, the resultant transmitted signal is

𝐱n(t)=𝐖n𝐬n(t),subscript𝐱𝑛𝑡subscript𝐖𝑛subscript𝐬𝑛𝑡\mathbf{x}_{n}(t)=\mathbf{W}_{n}\mathbf{s}_{n}(t),bold_x start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) = bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_s start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) , (2)

where 𝐖nNt×Ksubscript𝐖𝑛superscriptsubscript𝑁𝑡𝐾\mathbf{W}_{n}\in\mathbb{C}^{N_{t}\times K}bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ∈ blackboard_C start_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT × italic_K end_POSTSUPERSCRIPT denotes the beamforming matrix (alternatively termed as the linear precoder).

II-A Omnidirectional Phase

Initially, the ABS lacks knowledge regarding the channel information of the devices. Thus, the proposed protocol utilizes an omnidirectional beam for both sensing and communication with these devices. In this phase, we assign n=1𝑛1n=1italic_n = 1.

II-A1 Communication Process

According to literature [3], the channel vector from the ABS to the k𝑘kitalic_k-th device is represented as

𝐡k,n=α0dk𝐚(θk),subscript𝐡𝑘𝑛subscript𝛼0subscript𝑑𝑘𝐚subscript𝜃𝑘\mathbf{h}_{k,n}=\frac{\alpha_{0}}{d_{k}}\mathbf{a}(\theta_{k}),bold_h start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG bold_a ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) , (3)

where the term α0subscript𝛼0\alpha_{0}italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT signifies the path loss at a reference distance d0=1subscript𝑑01d_{0}=1italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 m. The distance between the ABS and the k𝑘kitalic_k-th device is dksubscript𝑑𝑘d_{k}italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. The steering vector aimed at device k𝑘kitalic_k is expressed as

𝐚(θk)=[1,ej2πdλsinθk,,ej2πdλ(Nt1)sinθk]T,𝐚subscript𝜃𝑘superscript1superscript𝑒𝑗2𝜋𝑑𝜆subscript𝜃𝑘superscript𝑒𝑗2𝜋𝑑𝜆subscript𝑁𝑡1subscript𝜃𝑘𝑇\mathbf{a}(\theta_{k})=[1,e^{j2{\pi}\frac{d}{\lambda}\sin{\theta_{k}}},...,e^{% j2{\pi}\frac{d}{\lambda}(N_{t}-1)\sin{\theta_{k}}}]^{T},bold_a ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = [ 1 , italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_d end_ARG start_ARG italic_λ end_ARG roman_sin italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , … , italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_d end_ARG start_ARG italic_λ end_ARG ( italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 1 ) roman_sin italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT , (4)

where d𝑑ditalic_d represents the gap between two neighboring antennas and θk(π,π)subscript𝜃𝑘𝜋𝜋\theta_{k}\in(-\pi,\pi)italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ∈ ( - italic_π , italic_π ) is the device’s angular position.

Consequently, the signal received by the k𝑘kitalic_k-th device is

yk,n(t)=𝐡k,nH𝐱n(t)+zC(t),subscript𝑦𝑘𝑛𝑡superscriptsubscript𝐡𝑘𝑛𝐻subscript𝐱𝑛𝑡subscript𝑧𝐶𝑡y_{k,n}(t)=\mathbf{h}_{k,n}^{H}\mathbf{x}_{n}(t)+z_{C}(t),italic_y start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) = bold_h start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_x start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) + italic_z start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT ( italic_t ) , (5)

with zC(t)𝒞𝒩(0,σC2)similar-tosubscript𝑧𝐶𝑡𝒞𝒩0superscriptsubscript𝜎𝐶2z_{C}(t)\sim\mathcal{CN}(0,\sigma_{C}^{2})italic_z start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT ( italic_t ) ∼ caligraphic_C caligraphic_N ( 0 , italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) symbolizing the additive white Gaussian noise (AWGN) at the receiver.

Given a unit total transmit power of PTsubscript𝑃𝑇P_{T}italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and based on the principles of omnidirectional beampattern [7], the orthogonal nature of the signal 𝐱n(t)subscript𝐱𝑛𝑡\mathbf{x}_{n}(t)bold_x start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) is ensured. For scenarios where KNt𝐾subscript𝑁𝑡K\geq N_{t}italic_K ≥ italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, the spatial covariance matrix of the omnidirectional beam is

𝐂=𝔼[𝐱n(t)𝐱nH(t)]=𝐖n𝐖nH=PTNt𝐈Nt,𝐂𝔼delimited-[]subscript𝐱𝑛𝑡subscriptsuperscript𝐱𝐻𝑛𝑡subscript𝐖𝑛subscriptsuperscript𝐖𝐻𝑛subscript𝑃𝑇subscript𝑁𝑡subscript𝐈subscript𝑁𝑡\mathbf{C}=\mathbb{E}[\mathbf{x}_{n}(t)\mathbf{x}^{H}_{n}(t)]=\mathbf{W}_{n}% \mathbf{W}^{H}_{n}=\frac{P_{T}}{N_{t}}\mathbf{I}_{N_{t}},bold_C = blackboard_E [ bold_x start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) bold_x start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) ] = bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_W start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = divide start_ARG italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG start_ARG italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG bold_I start_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUBSCRIPT , (6)

where 𝐈Ntsubscript𝐈subscript𝑁𝑡\mathbf{I}_{N_{t}}bold_I start_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUBSCRIPT denotes the Nt×Ntsubscript𝑁𝑡subscript𝑁𝑡N_{t}\times N_{t}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT × italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT identity matrix, the received SINR of the k𝑘kitalic_k-th device is represented as

γk,n=1K𝐡kH𝐂𝐡kK1K𝐡kH𝐂𝐡k+σC2=α02PTdk2Kα02PT(K1)dk2K+σC2.subscript𝛾𝑘𝑛1𝐾superscriptsubscript𝐡𝑘𝐻subscript𝐂𝐡𝑘𝐾1𝐾superscriptsubscript𝐡𝑘𝐻subscript𝐂𝐡𝑘superscriptsubscript𝜎𝐶2superscriptsubscript𝛼02subscript𝑃𝑇superscriptsubscript𝑑𝑘2𝐾superscriptsubscript𝛼02subscript𝑃𝑇𝐾1superscriptsubscript𝑑𝑘2𝐾superscriptsubscript𝜎𝐶2\gamma_{k,n}=\frac{\frac{1}{K}\mathbf{h}_{k}^{H}\mathbf{C}\mathbf{h}_{k}}{% \frac{K-1}{K}\mathbf{h}_{k}^{H}\mathbf{C}\mathbf{h}_{k}+\sigma_{C}^{2}}=\frac{% \frac{\alpha_{0}^{2}P_{T}}{d_{k}^{2}K}}{\frac{\alpha_{0}^{2}P_{T}(K-1)}{d_{k}^% {2}K}+\sigma_{C}^{2}}.italic_γ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG divide start_ARG 1 end_ARG start_ARG italic_K end_ARG bold_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Ch start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG start_ARG divide start_ARG italic_K - 1 end_ARG start_ARG italic_K end_ARG bold_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Ch start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG divide start_ARG italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_K end_ARG end_ARG start_ARG divide start_ARG italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_K - 1 ) end_ARG start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_K end_ARG + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (7)

If K<Nt𝐾subscript𝑁𝑡K<N_{t}italic_K < italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, the linear precoder is redundant; hence, K𝐾Kitalic_K antennas can be directly manipulated to transmit K𝐾Kitalic_K ISAC signals individually. The received SINR is same as (7).

In conclusion, the sum rate achievable by the IoT device k𝑘kitalic_k during the n𝑛nitalic_n-th time slot is denoted as

Rk,n=log2(1+γk,n).subscript𝑅𝑘𝑛subscript21subscript𝛾𝑘𝑛R_{k,n}=\log_{2}{(1+\gamma_{k,n})}.italic_R start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 1 + italic_γ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) . (8)

II-A2 Radar Process

According to [18], the echos of MIMO orthogonal waveform are still orthogonal, so radar can recognize the reflection of each device. Consequently, for device k𝑘kitalic_k, the received echo reflection is represented as

𝐫k,n(t)=βk𝐛(θk)𝐚H(θk)𝐱n(tτk,n)+𝐳k,n(t),subscript𝐫𝑘𝑛𝑡subscript𝛽𝑘𝐛subscript𝜃𝑘superscript𝐚𝐻subscript𝜃𝑘subscript𝐱𝑛𝑡subscript𝜏𝑘𝑛subscript𝐳𝑘𝑛𝑡\mathbf{r}_{k,n}(t)=\beta_{k}\mathbf{b}(\theta_{k})\mathbf{a}^{H}(\theta_{k})% \mathbf{x}_{n}(t-\tau_{k,n})+\mathbf{z}_{k,n}(t),bold_r start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) = italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT bold_b ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_x start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t - italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) + bold_z start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) , (9)

where 𝐛(θk)𝐛subscript𝜃𝑘\mathbf{b}(\theta_{k})bold_b ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is receive steering vector and has the same structure as 𝐚(θk)𝐚subscript𝜃𝑘\mathbf{a}(\theta_{k})bold_a ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ), denoted as:

𝐛(θk)=[1,ej2πdλsinθk,,ej2πdλ(Nr1)sinθk]T,𝐛subscript𝜃𝑘superscript1superscript𝑒𝑗2𝜋𝑑𝜆subscript𝜃𝑘superscript𝑒𝑗2𝜋𝑑𝜆subscript𝑁𝑟1subscript𝜃𝑘𝑇\mathbf{b}(\theta_{k})=[1,e^{j2{\pi}\frac{d}{\lambda}\sin{\theta_{k}}},...,e^{% j2{\pi}\frac{d}{\lambda}(N_{r}-1)\sin{\theta_{k}}}]^{T},bold_b ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = [ 1 , italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_d end_ARG start_ARG italic_λ end_ARG roman_sin italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , … , italic_e start_POSTSUPERSCRIPT italic_j 2 italic_π divide start_ARG italic_d end_ARG start_ARG italic_λ end_ARG ( italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - 1 ) roman_sin italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT , (10)

and 𝐳k,n(t)subscript𝐳𝑘𝑛𝑡\mathbf{z}_{k,n}(t)bold_z start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) represents the zero-mean complex additive white Gaussian noise with variance σ2superscript𝜎2\sigma^{2}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Both βk,nsubscript𝛽𝑘𝑛\beta_{k,n}italic_β start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT and τk,nsubscript𝜏𝑘𝑛\tau_{k,n}italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT signify the reflection coefficient and time lag for device k𝑘kitalic_k during time slot n𝑛nitalic_n. The coefficient of reflection is expressed as

βk=εdk2,subscript𝛽𝑘𝜀superscriptsubscript𝑑𝑘2\beta_{k}=\frac{\varepsilon}{d_{k}^{2}},italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = divide start_ARG italic_ε end_ARG start_ARG italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (11)

where ε𝜀\varepsilonitalic_ε denotes the complex radar cross-section (RCS).

Using a conventional matched-filtering methodology [8], signal delays can be estimated and the refined output vectors are

𝐑~k,n=GmβkPTNt𝐛(θk)𝐚H(θk)+𝐙r,subscript~𝐑𝑘𝑛subscript𝐺𝑚subscript𝛽𝑘subscript𝑃𝑇subscript𝑁𝑡𝐛subscript𝜃𝑘superscript𝐚𝐻subscript𝜃𝑘subscript𝐙𝑟\tilde{\mathbf{R}}_{k,n}=G_{m}\beta_{k}\frac{P_{T}}{N_{t}}\mathbf{b}(\theta_{k% })\mathbf{a}^{H}(\theta_{k})+\mathbf{Z}_{r},over~ start_ARG bold_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = italic_G start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT divide start_ARG italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG start_ARG italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG bold_b ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) + bold_Z start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT , (12)

where Gmsubscript𝐺𝑚G_{m}italic_G start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT is signal processing gain, brought by matched-filtering. The noise matrix, 𝐙rsubscript𝐙𝑟\mathbf{Z}_{r}bold_Z start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, consists of independent, zero-mean, and complex Gaussian elements with variance σr2superscriptsubscript𝜎𝑟2\sigma_{r}^{2}italic_σ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Additionally, the model measuring the distance is

τ^k,n=2dkc+zτk,n.subscript^𝜏𝑘𝑛2subscript𝑑𝑘𝑐subscript𝑧subscript𝜏𝑘𝑛\hat{\tau}_{k,n}=\frac{2d_{k}}{c}+z_{\tau_{k,n}}.over^ start_ARG italic_τ end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG 2 italic_d start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG start_ARG italic_c end_ARG + italic_z start_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT . (13)

With delay, the distance d^k,nsubscript^𝑑𝑘𝑛\hat{d}_{k,n}over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT can be simply derived. Employing either the Capon method or the generalized likelihood ratio test (GLRT) [9], the angle between the ABS and device k𝑘kitalic_k is deduced as

θ^k,n=θk,n+zθk,n,subscript^𝜃𝑘𝑛subscript𝜃𝑘𝑛subscript𝑧subscript𝜃𝑘𝑛\hat{\theta}_{k,n}=\theta_{k,n}+z_{\theta_{k,n}},over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT + italic_z start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT , (14)

where θk,nsubscript𝜃𝑘𝑛\theta_{k,n}italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT is equivalent to θksubscript𝜃𝑘\theta_{k}italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. The measurement models of (13) and (14) both take the Gaussian distribution model.

Accordingly, zτk,nsubscript𝑧subscript𝜏𝑘𝑛z_{\tau_{k,n}}italic_z start_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT and zθk,nsubscript𝑧subscript𝜃𝑘𝑛z_{\theta_{k,n}}italic_z start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT both represent the Gaussian noise in the measurement with zero mean and variances στk,n2subscriptsuperscript𝜎2subscript𝜏𝑘𝑛\sigma^{2}_{\tau_{k,n}}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT and σθk,n2subscriptsuperscript𝜎2subscript𝜃𝑘𝑛\sigma^{2}_{\theta_{k,n}}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT, respectively. Given the challenge in acquiring variances for τ^k,nsubscript^𝜏𝑘𝑛\hat{\tau}_{k,n}over^ start_ARG italic_τ end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT and θ^k,nsubscript^𝜃𝑘𝑛\hat{\theta}_{k,n}over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT, the Cramér-Rao Bound (CRB) is utilized as the sensing metric, due to its unbiased estimation and ability to provide a lower MSE bound [10]:

στk,n2subscriptsuperscript𝜎2subscript𝜏𝑘𝑛\displaystyle\sigma^{2}_{\tau_{k,n}}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT =CRBτ=α1SNRk,nNtNrκ2,absent𝐶𝑅subscript𝐵𝜏subscript𝛼1𝑆𝑁subscript𝑅𝑘𝑛subscript𝑁𝑡subscript𝑁𝑟superscript𝜅2\displaystyle=CRB_{\tau}=\frac{\alpha_{1}}{SNR_{k,n}N_{t}N_{r}\kappa^{2}},= italic_C italic_R italic_B start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT = divide start_ARG italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_S italic_N italic_R start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (15)
σθk,n2subscriptsuperscript𝜎2subscript𝜃𝑘𝑛\displaystyle\sigma^{2}_{\theta_{k,n}}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT =CRBθ=α2SNRk,nNtNrξ2,absent𝐶𝑅subscript𝐵𝜃subscript𝛼2𝑆𝑁subscript𝑅𝑘𝑛subscript𝑁𝑡subscript𝑁𝑟superscript𝜉2\displaystyle=CRB_{\theta}=\frac{\alpha_{2}}{SNR_{k,n}N_{t}N_{r}\xi^{2}},= italic_C italic_R italic_B start_POSTSUBSCRIPT italic_θ end_POSTSUBSCRIPT = divide start_ARG italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_S italic_N italic_R start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (16)

where pre-set constants are α1subscript𝛼1\alpha_{1}italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and α2subscript𝛼2\alpha_{2}italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, related to system configurations. The radar echo signal-noise ratio (SNR) is denoted as SNRk,n=PTGm|βk|2σ2𝑆𝑁subscript𝑅𝑘𝑛subscript𝑃𝑇subscript𝐺𝑚superscriptsubscript𝛽𝑘2superscript𝜎2SNR_{k,n}=\frac{P_{T}G_{m}|\beta_{k}|^{2}}{\sigma^{2}}italic_S italic_N italic_R start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_G start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT | italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG. As per [10], the squared effective bandwidth is κ2superscript𝜅2\kappa^{2}italic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and the root mean square aperture width of the beampattern is

ξ2=π2dk2cos2θk(Nt21)3λ2.superscript𝜉2superscript𝜋2subscriptsuperscript𝑑2𝑘superscript2subscript𝜃𝑘superscriptsubscript𝑁𝑡213superscript𝜆2\xi^{2}=\frac{\pi^{2}d^{2}_{k}\cos^{2}{\theta_{k}}(N_{t}^{2}-1)}{3\lambda^{2}}.italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) end_ARG start_ARG 3 italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (17)

With (13) and (14), we can derive estimated channel h^k,nsubscript^h𝑘𝑛\hat{\textbf{h}}_{k,n}over^ start_ARG h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT, which is part of the optimization problem in Section III.

II-B Directional Phase

Within this phase, the beampattern is tailored based on acquired parameters, enabling the transmission of a directional beampattern and the setting is established as n=2,,N𝑛2𝑁n=2,...,Nitalic_n = 2 , … , italic_N.

II-B1 Communication Process

The received signal for device k𝑘kitalic_k is represented as

yk,n(t)=𝐡k,nH𝐖n𝐬n(t)+zC(t).subscript𝑦𝑘𝑛𝑡superscriptsubscript𝐡𝑘𝑛𝐻subscript𝐖𝑛subscript𝐬𝑛𝑡subscript𝑧𝐶𝑡y_{k,n}(t)=\mathbf{h}_{k,n}^{H}\mathbf{W}_{n}\mathbf{s}_{n}(t)+z_{C}(t).italic_y start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) = bold_h start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_s start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_t ) + italic_z start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT ( italic_t ) . (18)

In concise terms, the SINR for device k𝑘kitalic_k can be articulated as

γk,n=|𝐡kH𝐰k,n|2i=1,ikK|𝐡kH𝐰i,n|2+σC2,subscript𝛾𝑘𝑛superscriptsuperscriptsubscript𝐡𝑘𝐻subscript𝐰𝑘𝑛2superscriptsubscriptformulae-sequence𝑖1𝑖𝑘𝐾superscriptsuperscriptsubscript𝐡𝑘𝐻subscript𝐰𝑖𝑛2superscriptsubscript𝜎𝐶2\displaystyle\gamma_{k,n}=\frac{|\mathbf{h}_{k}^{H}\mathbf{w}_{k,n}|^{2}}{\sum% _{i=1,i\neq k}^{K}|\mathbf{h}_{k}^{H}\mathbf{w}_{i,n}|^{2}+\sigma_{C}^{2}},italic_γ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG | bold_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_w start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_w start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (19)

where the beamforming vector, denoted as 𝐰k,nsubscript𝐰𝑘𝑛\mathbf{w}_{k,n}bold_w start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT, stands as the k𝑘kitalic_k-th column vector extracted from 𝐖nsubscript𝐖𝑛\mathbf{W}_{n}bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. For time slot n𝑛nitalic_n, the achievable sum rate pertaining to IoT device k𝑘kitalic_k aligns with expression (8)8(\ref{sumrate})( ).

II-B2 Radar Process

According to the assumption of [5], the steering vectors from different angles are asymptotically orthogonal under the massive MIMO regime, so the reflected echoes do not interfere with each other. Similarly, the received echo reflection is characterized as

𝐫k,n(t)=βk𝐛(θk)𝐚H(θk)𝐖n𝐬(tτk,n)+𝐳k,n(t).subscript𝐫𝑘𝑛𝑡subscript𝛽𝑘𝐛subscript𝜃𝑘superscript𝐚𝐻subscript𝜃𝑘subscript𝐖𝑛𝐬𝑡subscript𝜏𝑘𝑛subscript𝐳𝑘𝑛𝑡\mathbf{r}_{k,n}(t)=\beta_{k}\mathbf{b}(\theta_{k})\mathbf{a}^{H}(\theta_{k})% \mathbf{W}_{n}\mathbf{s}(t-\tau_{k,n})+\mathbf{z}_{k,n}(t).bold_r start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) = italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT bold_b ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_W start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT bold_s ( italic_t - italic_τ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) + bold_z start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ( italic_t ) . (20)

Post the implementation of matched-filtering, the models used for measurements remain consistent with (13)13(\ref{estitau})( ) and (14)14(\ref{estitheta})( ). Concludingly, the SNR of the radar’s echoed signal can be delineated as

SNRk,n=Gm|βk|2i=1K|𝐚H(θk)𝐰i,n|2σ2.𝑆𝑁subscript𝑅𝑘𝑛subscript𝐺𝑚superscriptsubscript𝛽𝑘2superscriptsubscript𝑖1𝐾superscriptsuperscript𝐚𝐻subscript𝜃𝑘subscript𝐰𝑖𝑛2superscript𝜎2\displaystyle SNR_{k,n}=\frac{G_{m}|\beta_{k}|^{2}\sum_{i=1}^{K}|\mathbf{a}^{H% }(\theta_{k})\mathbf{w}_{i,n}|^{2}}{\sigma^{2}}.italic_S italic_N italic_R start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG italic_G start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT | italic_β start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) bold_w start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (21)

III Problem Formulation and Solution

In this section, we develop a formulation for the beamforming problem employing sensing-assisted communication to optimize the communication rate for all devices.

III-A Problem Formulation

The channel estimated by an ABS to device k𝑘kitalic_k using position information from the radar echo is given by

𝐡^k,n=α0d^k,n𝐚(θ^k,n).subscript^𝐡𝑘𝑛subscript𝛼0subscript^𝑑𝑘𝑛𝐚subscript^𝜃𝑘𝑛\hat{\mathbf{h}}_{k,n}=\frac{\alpha_{0}}{\hat{d}_{k,n}}\mathbf{a}(\hat{\theta}% _{k,n}).over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = divide start_ARG italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG over^ start_ARG italic_d end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT end_ARG bold_a ( over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) . (22)

The ISAC system postulates that the communication rate with device k𝑘kitalic_k at time slot n𝑛nitalic_n is

R^k,n=log2(1+tr(𝐡^k,n𝐡^k,nH𝐖k,n)i=1,ikKtr(𝐡^k,n𝐡^k,nH𝐖i,n)+σC2)),\displaystyle\hat{R}_{k,n}=\log_{2}{\left(1+\frac{\text{tr}(\hat{\mathbf{h}}_{% k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{k,n})}{\sum_{i=1,i\neq k}^{K}\text{% tr}(\hat{\mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n})+\sigma_{% C}^{2})}\right)},over^ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 1 + divide start_ARG tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG ) , (23)

where the matrix 𝐖k,nsubscript𝐖𝑘𝑛\mathbf{W}_{k,n}bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT is denoted as 𝐰k,n𝐰k,nHsubscript𝐰𝑘𝑛superscriptsubscript𝐰𝑘𝑛𝐻\mathbf{w}_{k,n}\mathbf{w}_{k,n}^{H}bold_w start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT bold_w start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT. Therefore, the optimization problem can be expressed as

max{𝐖k,n}k=1KR^k,nsubscript𝐖𝑘𝑛maxsuperscriptsubscript𝑘1𝐾subscript^𝑅𝑘𝑛\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\underset{\{% \mathbf{W}_{k,n}\}}{\text{max}}\sum\limits_{k=1}^{K}\hat{R}_{k,n}start_UNDERACCENT { bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT } end_UNDERACCENT start_ARG max end_ARG ∑ start_POSTSUBSCRIPT italic_k = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT over^ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT (24a)
s.t.tr(k=1K𝐖k,n)=PTs.t.trsuperscriptsubscript𝑘1𝐾subscript𝐖𝑘𝑛subscript𝑃𝑇\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\text{s.t.}~{}% \text{tr}\left(\sum\limits_{k=1}^{K}\mathbf{W}_{k,n}\right)=P_{T}s.t. tr ( ∑ start_POSTSUBSCRIPT italic_k = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) = italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT (24b)
𝐖k,n0,𝐖k,n=𝐖k,nH,k,formulae-sequencesucceeds-or-equalssubscript𝐖𝑘𝑛0subscript𝐖𝑘𝑛superscriptsubscript𝐖𝑘𝑛𝐻for-all𝑘\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\mathbf{W}_{k,n}% \succeq 0,\mathbf{W}_{k,n}=\mathbf{W}_{k,n}^{H},~{}\forall k,bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ⪰ 0 , bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT = bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT , ∀ italic_k , (24c)
rank(𝐖k,n)=1,k,ranksubscript𝐖𝑘𝑛1for-all𝑘\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\text{rank}(% \mathbf{W}_{k,n})=1,~{}\forall k,rank ( bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) = 1 , ∀ italic_k , (24d)
R^k,nΓk,k,subscript^𝑅𝑘𝑛subscriptΓ𝑘for-all𝑘\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\hat{R}_{k,n}\geq% \Gamma_{k},~{}\forall k,over^ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ≥ roman_Γ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , ∀ italic_k , (24e)
|𝐚H(θ^k,n1)𝐖k,n𝐚(θ^k,n1)𝐚H(θk,ncover)𝐖k,n𝐚(θk,ncover)|superscript𝐚𝐻subscript^𝜃𝑘𝑛1subscript𝐖𝑘𝑛𝐚subscript^𝜃𝑘𝑛1superscript𝐚𝐻superscriptsubscript𝜃𝑘𝑛𝑐𝑜𝑣𝑒𝑟subscript𝐖𝑘𝑛𝐚superscriptsubscript𝜃𝑘𝑛𝑐𝑜𝑣𝑒𝑟\displaystyle|\mathbf{a}^{H}(\hat{\theta}_{k,n-1})\mathbf{W}_{k,n}\mathbf{a}(% \hat{\theta}_{k,n-1})-\mathbf{a}^{H}(\theta_{k,n}^{cover})\mathbf{W}_{k,n}% \mathbf{a}(\theta_{k,n}^{cover})|| bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT ) bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT bold_a ( over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT ) - bold_a start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c italic_o italic_v italic_e italic_r end_POSTSUPERSCRIPT ) bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT bold_a ( italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c italic_o italic_v italic_e italic_r end_POSTSUPERSCRIPT ) |
Bktr(𝐖k,n),k,|θ^k,n1θk,ncover|lσθk,n1,formulae-sequenceabsentsubscript𝐵𝑘trsubscript𝐖𝑘𝑛for-all𝑘for-allsubscript^𝜃𝑘𝑛1superscriptsubscript𝜃𝑘𝑛𝑐𝑜𝑣𝑒𝑟𝑙subscript𝜎subscript𝜃𝑘𝑛1\displaystyle~{}~{}~{}\leq B_{k}\text{tr}(\mathbf{W}_{k,n}),~{}\forall k,~{}% \forall|\hat{\theta}_{k,n-1}-\theta_{k,n}^{cover}|\leq l\sigma_{\theta_{k,n-1}},≤ italic_B start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT tr ( bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) , ∀ italic_k , ∀ | over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT - italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c italic_o italic_v italic_e italic_r end_POSTSUPERSCRIPT | ≤ italic_l italic_σ start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , (24f)

where the transmission power limitation is imposed by equation (24b). The matrices 𝐖k,nsubscript𝐖𝑘𝑛\mathbf{W}_{k,n}bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT are subject to three constraints, which are semidefinite, Hermitian, and rank-one, as specified by equations (24c) and (24d). The term ΓksubscriptΓ𝑘\Gamma_{k}roman_Γ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT denotes the communication rate threshold for the k𝑘kitalic_k-th device, and increasing the value of ΓksubscriptΓ𝑘\Gamma_{k}roman_Γ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT allows more power to be given to remote devices to ensure the minimum communication rate. Drawing from the findings in [9], we integrate sensing and communication paradigms within equation (24f), capturing the influence of sensing accuracy on beamforming strategies. Pertinently, within the context of equation (24f), based on the Gaussian distribution principle, the constant l𝑙litalic_l is generally chosen as 3 to cover the true angle as much as possible. This ensures maximum angle coverage corresponding to θ^k,n1subscript^𝜃𝑘𝑛1\hat{\theta}_{k,n-1}over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT. Consequently, the cumulative angles θk,ncoversuperscriptsubscript𝜃𝑘𝑛𝑐𝑜𝑣𝑒𝑟\theta_{k,n}^{cover}italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c italic_o italic_v italic_e italic_r end_POSTSUPERSCRIPT, fulfilling the condition |θ^k,n1θk,ncover|lσθk,n1subscript^𝜃𝑘𝑛1superscriptsubscript𝜃𝑘𝑛𝑐𝑜𝑣𝑒𝑟𝑙subscript𝜎subscript𝜃𝑘𝑛1|\hat{\theta}_{k,n-1}-\theta_{k,n}^{cover}|\leq l\sigma_{\theta_{k,n-1}}| over^ start_ARG italic_θ end_ARG start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT - italic_θ start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_c italic_o italic_v italic_e italic_r end_POSTSUPERSCRIPT | ≤ italic_l italic_σ start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_k , italic_n - 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, delineate the beam coverage. To achieve a compact transmission beam, the parameter Bksubscript𝐵𝑘B_{k}italic_B start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT is minimized. In short, by constraining the beam gain of the device at possible angles, (24f) makes the optimized beampattern an approximate ideal radar beampattern with width controllable.

III-B Problem Solution

By ignoring rank-one constraint (24d), problem (24) is rewritten as

max{𝐖k,n}k=1Klog2(1+tr(𝐡^k,n𝐡^k,nH𝐖k,n)i=1,ikKtr(𝐡^k,n𝐡^k,nH𝐖i,n)+σC2))\displaystyle\underset{\{\mathbf{W}_{k,n}\}}{\text{max}}\sum\limits_{k=1}^{K}% \log_{2}{\left(1+\frac{\text{tr}(\hat{\mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^% {H}\mathbf{W}_{k,n})}{\sum_{i=1,i\neq k}^{K}\text{tr}(\hat{\mathbf{h}}_{k,n}% \hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n})+\sigma_{C}^{2})}\right)}start_UNDERACCENT { bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT } end_UNDERACCENT start_ARG max end_ARG ∑ start_POSTSUBSCRIPT italic_k = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 1 + divide start_ARG tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT ) end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG )
s.t.(24b),(24c),(24e)and(24f).s.t.24b24c24eand24f\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}\text{s.t.}~{}(\text{\ref{% opti1-b}}),(\text{\ref{opti1-c}}),(\text{\ref{opti1-e}})~{}\text{and}~{}(\text% {\ref{opti1-f}}).s.t. ( ) , ( ) , ( ) and ( ) . (25)

Secondly, we adopt SCA to approximate the non concave objective function in problem (25) to a concave objective function, which is implemented through multiple iterations. We consider that iteration q0𝑞0q\geq 0italic_q ≥ 0, and {𝐖k,n}subscript𝐖𝑘𝑛\{\mathbf{W}_{k,n}\}{ bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT } is depicted as {𝐖k,n(q)}superscriptsubscript𝐖𝑘𝑛𝑞\{\mathbf{W}_{k,n}^{(q)}\}{ bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT } accordingly. Afterwards, we have

R^k,nsubscript^𝑅𝑘𝑛\displaystyle\hat{R}_{k,n}over^ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT =log2(i=1Ktr(𝐡^k,n𝐡^k,nH𝐖i,n)+σC2)absentsubscript2superscriptsubscript𝑖1𝐾trsubscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻subscript𝐖𝑖𝑛superscriptsubscript𝜎𝐶2\displaystyle=\log_{2}{\left({\textstyle\sum_{i=1}^{K}}\text{tr}(\hat{\mathbf{% h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n})+\sigma_{C}^{2}\right)}= roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
log2(i=1,ikKtr(𝐡^k,n𝐡^k,nH𝐖i,n)+σC2)subscript2superscriptsubscriptformulae-sequence𝑖1𝑖𝑘𝐾trsubscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻subscript𝐖𝑖𝑛superscriptsubscript𝜎𝐶2\displaystyle~{}~{}~{}~{}-\log_{2}{\left({\textstyle\sum_{i=1,i\neq k}^{K}}% \text{tr}(\hat{\mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n})+% \sigma_{C}^{2}\right)}- roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (26)
log2(i=1Ktr(𝐡^k,n𝐡^k,nH𝐖i,n)+σC2)absentsubscript2superscriptsubscript𝑖1𝐾trsubscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻subscript𝐖𝑖𝑛superscriptsubscript𝜎𝐶2\displaystyle\geq\log_{2}{\left({\textstyle\sum_{i=1}^{K}}\text{tr}(\hat{% \mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n})+\sigma_{C}^{2}% \right)}≥ roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
log2(i=1,ikKtr(𝐡^k,n𝐡^k,nH𝐖i,n(q))+σC2)subscript2superscriptsubscriptformulae-sequence𝑖1𝑖𝑘𝐾trsubscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻superscriptsubscript𝐖𝑖𝑛𝑞superscriptsubscript𝜎𝐶2\displaystyle~{}~{}~{}~{}-\log_{2}{\left({\textstyle\sum_{i=1,i\neq k}^{K}}% \text{tr}(\hat{\mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n}^{(q% )})+\sigma_{C}^{2}\right)}- roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
i=1,ikKtr(𝐗k,n(q)(𝐖i,n𝐖i,n(q)))R~k,n(q),superscriptsubscriptformulae-sequence𝑖1𝑖𝑘𝐾trsuperscriptsubscript𝐗𝑘𝑛𝑞subscript𝐖𝑖𝑛superscriptsubscript𝐖𝑖𝑛𝑞superscriptsubscript~𝑅𝑘𝑛𝑞\displaystyle~{}~{}~{}~{}-{\textstyle\sum_{i=1,i\neq k}^{K}}\text{tr}\left(% \mathbf{X}_{k,n}^{(q)}(\mathbf{W}_{i,n}-\mathbf{W}_{i,n}^{(q)})\right)% \triangleq\tilde{R}_{k,n}^{(q)},- ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( bold_X start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT ( bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT - bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT ) ) ≜ over~ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT , (27)

where 𝐗k,n(q)superscriptsubscript𝐗𝑘𝑛𝑞\mathbf{X}_{k,n}^{(q)}bold_X start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT is defined as

𝐗k,n(q)=log2(e)𝐡^k,n𝐡^k,nHi=1,ikKtr(𝐡^k,n𝐡^k,nH𝐖i,n(q))+σC2.superscriptsubscript𝐗𝑘𝑛𝑞subscript2𝑒subscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻superscriptsubscriptformulae-sequence𝑖1𝑖𝑘𝐾trsubscript^𝐡𝑘𝑛superscriptsubscript^𝐡𝑘𝑛𝐻superscriptsubscript𝐖𝑖𝑛𝑞superscriptsubscript𝜎𝐶2\displaystyle\mathbf{X}_{k,n}^{(q)}=\frac{\log_{2}{(e)}\hat{\mathbf{h}}_{k,n}% \hat{\mathbf{h}}_{k,n}^{H}}{{\textstyle\sum_{i=1,i\neq k}^{K}}\text{tr}(\hat{% \mathbf{h}}_{k,n}\hat{\mathbf{h}}_{k,n}^{H}\mathbf{W}_{i,n}^{(q)})+\sigma_{C}^% {2}}.bold_X start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT = divide start_ARG roman_log start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_e ) over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_i = 1 , italic_i ≠ italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT tr ( over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT over^ start_ARG bold_h end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_W start_POSTSUBSCRIPT italic_i , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT ) + italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (28)

In short, we implement first-order Taylor expansion on the second term of (III-B), then take its lower bound, and the obtained objective function in (III-B) will approach convex function. Accordingly, problem (25) is transformed into the following problem (25) in the q𝑞qitalic_q-th iteration.

max{𝐖k,n}subscript𝐖𝑘𝑛max\displaystyle\underset{\{\mathbf{W}_{k,n}\}}{\text{max}}start_UNDERACCENT { bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT } end_UNDERACCENT start_ARG max end_ARG k=1KR~k,n(q)superscriptsubscript𝑘1𝐾superscriptsubscript~𝑅𝑘𝑛𝑞\displaystyle\sum\limits_{k=1}^{K}\tilde{R}_{k,n}^{(q)}∑ start_POSTSUBSCRIPT italic_k = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT over~ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q ) end_POSTSUPERSCRIPT
s.t. (24b),(24c),(24e)and(24f).24b24c24eand24f\displaystyle(\text{\ref{opti1-b}}),(\text{\ref{opti1-c}}),(\text{\ref{opti1-e% }})~{}\text{and}~{}(\text{\ref{opti1-f}}).( ) , ( ) , ( ) and ( ) . (29)

Now the problem (III-B) is convex, which can be optimally solved by convex optimization solvers such as CVX. After multiple iterations, assuming we obtain an approximate iterative solution of . If {𝐖k,n(q)}superscriptsubscript𝐖𝑘𝑛superscript𝑞\{\mathbf{W}_{k,n}^{(q^{\star})}\}{ bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q start_POSTSUPERSCRIPT ⋆ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT } satisfies the rank-one constraint, then it is the solution to the original problem (24). However, in most cases, {𝐖k,n(q)}superscriptsubscript𝐖𝑘𝑛superscript𝑞\{\mathbf{W}_{k,n}^{(q^{\star})}\}{ bold_W start_POSTSUBSCRIPT italic_k , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_q start_POSTSUPERSCRIPT ⋆ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT } is not a rank-one matrix, but a suboptimal solution can be obtained by classic SDR technique. Due to our strict requirements on the beam shape, conventional eigenvalue decomposition and Gaussian randomization [4] cannot meet our needs, so we have decided to adopt IRM.

IV Simulation Results

In this section, we discuss the numerical findings to substantiate the efficacy of the introduced ISAC protocol. With the main reference to [5] and [6], our simulation parameters delineated in Table II. For the sake of consistency, we fixed the antenna spacing to half the wavelength, ensured that the number of transmitting antennas mirrors the number of receiving ones, and maintained an equal noise variance for both communication and radar, represented as d=λ/2𝑑𝜆2d=\lambda/2italic_d = italic_λ / 2, Nt=Nrsubscript𝑁𝑡subscript𝑁𝑟N_{t}=N_{r}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, and σC2=σ2superscriptsubscript𝜎𝐶2superscript𝜎2\sigma_{C}^{2}=\sigma^{2}italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT.

TABLE II: Simulation Parameters
parameter value
Number of transmit antennas Ntsubscript𝑁𝑡N_{t}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 20
Effective bandwidth κ𝜅\kappaitalic_κ (Mhz) 4
Signal processing gain Gmsubscript𝐺𝑚G_{m}italic_G start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT 10
Carrier wavelength λ𝜆\lambdaitalic_λ (m) 0.06
Additive white Gaussian noise σC2superscriptsubscript𝜎𝐶2\sigma_{C}^{2}italic_σ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 1
Interval length ΔTΔ𝑇\Delta Troman_Δ italic_T 0.01
Power loss of reference distance α02superscriptsubscript𝛼02\alpha_{0}^{2}italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (dB) 50

Refer to caption

Figure 2: Proposed ISAC protocol versus the traditional omnidirectional sensing-based method and the pilot-based approach.

Refer to caption

Figure 3: Beampattern at different time slots.

Refer to caption

Figure 4: Communication rate performance of different number of antennas.

Fig. 2 provides an initial assessment, contrasting the communication rate attained through various techniques across time slots. Two devices are considered: Device 1 located at an angle and distance of 60 and 50m, respectively, and Device 2 positioned at 120 and 100m. For reference, ’Sensing-based’ utilizes omnidirectional sensing in the initial phase and then transmits information through communication waveform. For the ’Pilot-based’ method, linear minimum mean squared Error (LMMSE) channel estimation is used before communication. The communication rate for the proposed protocol is initially minimal but elevates as data accuracy improves. By the time slot n𝑛nitalic_n = 4, the communication rate plateaus, evidencing an enhancement of 6.3%percent\%% and 7.9%percent\%% over ’Sensing-based’ and ’Pilot-based’, respectively. Over the duration of N𝑁Nitalic_N time slots, the throughput experiences an elevation of 5.7%percent\%% and 8.7%percent\%% compared to the aforementioned methods, respectively.

Fig. 3 visually represents how the beampattern fluctuates across time slots. Retaining parameters from Fig. 2, at n𝑛nitalic_n = 2, Device 1 exhibits a reduced beam coverage relative to Device 2, an outcome of its proximal location yielding heightened sensing precision. As time slots progress, the beam coverage visibly contracts, with the beampattern gain consistently amplifying. By the point of n𝑛nitalic_n = 4, the beampattern becomes invariant.

The influence of varying antenna quantities on the beampattern is depicted in Fig. 4. While keeping other parameters static, a reduction in antenna numbers correlates with a notable communication rate decline for every time slot. Relative to our proposed protocol, the communication rates for ’Ntsubscript𝑁𝑡N_{t}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = Nrsubscript𝑁𝑟N_{r}italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = 15’ and ’Ntsubscript𝑁𝑡N_{t}italic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = Nrsubscript𝑁𝑟N_{r}italic_N start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = 10’ decline by 9.8%percent\%% and 29.2%percent\%%, respectively. Such reductions are attributed as follows. Primarily, the antenna quantity influences the variance in angle and distance measurements, which translates to measurement accuracy. Furthermore, a decrement in gain directly compromises the SINR of the radar echo, further influencing accuracy in subsequent time slots.

V Conclusion

In this paper, we have proposed a sensing-assisted communication protocol based on the ISAC system. The central proposition is the enhancement of the downlink communication rate by dynamically refining the transmit beampattern of the current time slot, capitalizing on sensory data procured from the preceding slot. Additionally, the integration of the SCA technique and the IRM algorithm has facilitated the resolution of the proposed optimization challenge. The numerical results validate the advantages of the proposed method. Further, the lower complexity algorithms is worth studying in the future.

References

  • [1] L. Fan, Y. Weijie, Y. Jinhong, Z. J. Andrew, F. Zesong, and Z. Jianming, “Radar-communication spectrum sharing and integration: Overview and prospect,” 2020. [Online]. Available: https://api.semanticscholar.org/CorpusID:235863912
  • [2] Y. Liu, G. Liao, Z. Yang, and J. Xu, “Multiobjective optimal waveform design for OFDM integrated radar and communication systems,” Signal Processing, vol. 141, pp. 331–342, 2017. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0165168417302360
  • [3] Z. Lyu, G. Zhu, and J. Xu, “Joint maneuver and beamforming design for UAV-enabled integrated sensing and communication,” IEEE Transactions on Wireless Communications, vol. 22, no. 4, pp. 2424–2440, 2023.
  • [4] F. Liu, C. Masouros, A. Li, H. Sun, and L. Hanzo, “MU-MIMO communications with MIMO radar: From co-existence to joint transmission,” IEEE Transactions on Wireless Communications, vol. 17, no. 4, pp. 2755–2770, 2018.
  • [5] F. Liu, W. Yuan, C. Masouros, and J. Yuan, “Radar-assisted predictive beamforming for vehicular links: Communication served by sensing,” IEEE Transactions on Wireless Communications, vol. 19, no. 11, pp. 7704–7719, 2020.
  • [6] Z. Du, F. Liu, and Z. Zhang, “Sensing-assisted beam tracking in V2I networks: Extended target case,” in ICASSP 2022 - 2022 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), 2022, pp. 8727–8731.
  • [7] D. R. Fuhrmann and G. San Antonio, “Transmit beamforming for MIMO radar systems using signal cross-correlation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 44, no. 1, pp. 171–186, 2008.
  • [8] B. Friedlander, “On transmit beamforming for MIMO radar,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 4, pp. 3376–3388, 2012.
  • [9] J. Li and P. Stoica, “MIMO radar with colocated antennas,” IEEE Signal Processing Magazine, vol. 24, no. 5, pp. 106–114, 2007.
  • [10] A. Liu, Z. Huang, M. Li, Y. Wan, W. Li, T. X. Han, C. Liu, R. Du, D. K. P. Tan, J. Lu, Y. Shen, F. Colone, and K. Chetty, “A survey on fundamental limits of integrated sensing and communication,” IEEE Communications Surveys & Tutorials, vol. 24, no. 2, pp. 994–1034, 2022.
  • [11] Q. T. Dinh and M. Diehl, “Local convergence of sequential convex programming for nonconvex optimization,” 2010. [Online]. Available: https://api.semanticscholar.org/CorpusID:7036853
  • [12] C. Sun, Y. Liu, R. Dai, and D. J. Grymin, “Two approaches for path planning of unmanned aerial vehicles with avoidance zones,” Journal of Guidance Control and Dynamics, vol. 40, pp. 2076–2083, 2017. [Online]. Available: https://api.semanticscholar.org/CorpusID:125533974
  • [13] Y. Liu and K. Yang, “Communication, sensing, computing and energy harvesting in smart cities,” IET Smart Cities, vol. 4, no. 4, pp. 265–274, 2022. [Online]. Available: https://ietresearch.onlinelibrary.wiley.com/doi/abs/10.1049/smc2.12041
  • [14] L. Xiang, K. Xu, J. Hu, and K. Yang, “Green beamforming design for integrated sensing and communication systems: A practical approach using beam-matching error metrics,” IEEE Transactions on Vehicular Technology, vol. 73, no. 4, pp. 5935–5940, 2024.
  • [15] H. L. V. Trees, “Optimum array processing: Part IV of detection, estimation, and modulation theory,” 2002. [Online]. Available: https://api.semanticscholar.org/CorpusID:197455156
  • [16] X. Jing, F. Liu, C. Masouros, and Y. Zeng, “ISAC from the sky: UAV trajectory design for joint communication and target localization,” 2022.
  • [17] Y. Cui, F. Liu, X. Jing, and J. Mu, “Integrating sensing and communications for ubiquitous IoT: Applications, trends, and challenges,” IEEE Network, vol. 35, no. 5, pp. 158–167, 2021.
  • [18] L. Xu, J. Li, and P. Stoica, “Radar imaging via adaptive MIMO techniques,” in 2006 14th European Signal Processing Conference, 2006, pp. 1–5.
  • [19] L. Xiang, K. Xu, J. Hu, C. Masouros, and K. Yang, “Robust NOMA-assisted OTFS-ISAC network design with 3-D motion prediction topology,” IEEE Internet of Things Journal, vol. 11, no. 9, pp. 15 909–15 918, 2024.
  • [20] K. Meng, Q. Wu, J. Xu, W. Chen, Z. Feng, R. Schober, and A. L. Swindlehurst, “UAV-enabled integrated sensing and communication: Opportunities and challenges,” IEEE Wireless Communications, vol. 31, no. 2, pp. 97–104, 2024.