A Bistatic Sensing System in Space-Air-Ground Integrated Networks

Xiangyu Li1,2, Bodong Shang2*, Qingqing Wu1 1Department of Electronic Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
2Eastern Institute for Advanced Study, Eastern Institute of Technology, Ningbo 315200, China
e-mail: xiangyuli@sjtu.edu.cn, bdshang@eitech.edu.cn, qingqingwu@sjtu.edu.cn
Abstract

Sensing is anticipated to have wider extensions in communication systems with the boom of non-terrestrial networks (NTNs) during the past years. In this paper, we study a bistatic sensing system by maximizing the signal-to-interference-plus-noise ration (SINR) from the target aircraft in the space-air-ground integrated network (SAGIN). We formulate a joint optimization problem for the transmit beamforming of low-earth orbit (LEO) satellite and the receive filtering of ground base station. To tackle this problem, we decompose the original problem into two sub-problems and use the alternating optimization to solve them iteratively. Using techniques of fractional programming and generalized Rayleigh quotient, the closed-form solution for each sub-problem is returned. Simulation results show that the proposed algorithm has good convergence performance. Moreover, the optimization of receive filtering dominates the optimality, especially when the satellite altitude becomes higher, which provides valuable network design insights.

Index Terms:
Bistatic sensing, beamforming, LEO satellite, space-air-ground integrated networks, alternating optimization.

I Introduction

I-A Motivations

Sensing has been envisioned to become a key technology in the fifth generation (5G) and beyond wireless communication systems for obtaining emerging applications such as navigation and monitoring [1]. By applying sensing techniques, the quality of service (QoS) in terms of mutual communications as well as continuous detection and tracking of targets can be enhanced [2]. Due to the expansion of global networks and the advancement of air communication and transportation, vital significance has been brought to the task of sensing aircraft for aviation safety, traffic efficiency, and aerial security purposes [3]. By tracking and monitoring the movements of aircraft, not only the airspace traffic can be better controlled, but potential threats or unauthorized airspace entries can be detected promptly as well. Therefore, the enhancement of sensing quality for aircraft remains as a crucial problem.

Traditionally, the implementation of target sensing is based on terrestrial networks (TNs) [4]. However, the direct extension of sensing to wider areas is facing challenges. First, considering the large range and limited number of sensing equipment, it can be difficult to simultaneously monitor a large number of devices and timely focus on a certain target, especially in aerial regions. Second, long distances to the sensed targets will bring high propagation delay and path loss when echo signals are reflected back to the transceiving radar, as well as large investments in specific sensing equipment.

Satellite communications are recognized as a solution to satisfy the growing demand for seamless wireless connectivity and enhance the service availability in a global range [5, 6]. Compared with geostationary-earth orbit (GEO) and medium-earth orbit (MEO) satellites, low-earth orbit (LEO) satellites impose much less stringent requirements, including a lower latency, higher data rate and high visibility because of their lower altitudes [7]. Thus, the LEO satellite can serve as a promising candidate for wide-range ubiquitous sensing, while satisfying the requirements of low-altitude economy.

For sensing aircraft, LEO satellite-assisted sensing systems may have more advantages over terrestrial sensing systems. On the one hand, an LEO satellite could simultaneously sense and detect target aircraft which are tens of kilometers from each other due to its wide coverage. The terrestrial radar base stations (BSs), however, mainly detect targets below their antenna deployment altitude, and those within their separate range of sector. On the other hand, powered by solar energy, satellites can have more energy for continuous sensing tasks, while the process of scanning and target searching for terrestrial radar BSs can be of huge energy consumption.

I-B Related Work

Sensing techniques have been explored in many scenarios of wireless networks during the past decades, especially in TNs. The authors in [8] explore sensing in vehicle transportation settings and confirm the feasibility and necessity of using the functional integrated sensing and communication platform for different sensing modes. To monitor the operations of machines and the environment in the industrial Internet of Things (IIoT), Zhang et al. use crowd sensing technology and propose a Particle Swarm Optimization-Elman Neural Network algorithm for position prediction [9]. Remote sensing enabled by satellite or airplane radar systems can be integrated with machine learning for the detection of airborne particulates and dust sources [10], and be widely used in environmental monitoring [11]. Additionally, WiFi sensing is utilized indoor to recognize dynamic human activities such as breathing and heart beat by applying modeling or learning-based algorithms for over-the-air signal analysis [12].

Integrated sensing and communication (ISAC) systems have recently been explored for TNs. The work of [13] studies intelligent omni surfaces (IOSs)-aided target sensing and user communication for a blocked urban region, where its bistatic setting for the radar system can effectively reduce the impacts of path loss. Under the condition of imperfect channel state information (CSI), an intelligent reflecting surface (IRS)-aided sensing system is explored in the terrestrial setting [14] while no receiving operations are employed. Furthermore, as satellite-assisted communications are becoming a new trend, [15] investigates massive multiple-input multiple-output in ISAC LEO satellite systems using beam squint-aware techniques. The authors in [16] study non-orthogonal multiple access (NOMA)-aided precoding for ISAC by LEO satellites to ensure security. It is worth mentioning that previous works on sensing systems are mainly targeted at TNs or non-terrestrial networks (NTNs) with direct space-ground channels. There has not been a research on the detection of aerial targets in the space-air-ground integrated network (SAGIN).

I-C Contributions

Motivated by the aforementioned works, we aim to use an LEO satellite as the transmitter and improve the sensing performance of a specific aircraft under the interfering sensing signals. Considering the huge path loss of echo signals reflected back from aircraft to the LEO satellite, we utilize a BS as the bistatic radar on the ground to receive the echo signals. The main contributions are summarized as follows.

  • We establish a bistatic sensing model for the SAGIN, which includes the satellite-to-aircraft channel, the aircraft-to-ground channel, and the sensing signal transmission model.

  • We jointly design the transmit beamforming and receive filtering to maxmize the SINR from the target aircraft under the total transmit power constraint and the receive filtering requirement.

  • We derive a closed-form solution for each sub-problem based on fractional programming and generalized Rayleigh quotient, respectively, to solve the proposed alternating problem efficiently.

  • The simulation results show that the receive filtering optimization dominates the overall performance, and higher transmit power and more antenna numbers contributes positively to the system improvement.

Notations: Boldface, lower-case letters denote column vectors and boldface upper-case letters denote matrices. The superscript ()1superscript1(\cdot)^{-1}( ⋅ ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, ()Tsuperscript𝑇(\cdot)^{T}( ⋅ ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT and ()Hsuperscript𝐻(\cdot)^{H}( ⋅ ) start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT represent the inverse, transpose and Hermitian, respectively. 𝐈Nsubscript𝐈𝑁\mathbf{I}_{N}bold_I start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT refers to the identity matrices of N𝑁Nitalic_N dimensions, and M×Nsuperscript𝑀𝑁\mathbb{C}^{M\times N}blackboard_C start_POSTSUPERSCRIPT italic_M × italic_N end_POSTSUPERSCRIPT denotes the M×N𝑀𝑁M\times Nitalic_M × italic_N complex matrix. For a matrix 𝐀𝐀\mathbf{A}bold_A and a vector 𝐱𝐱\mathbf{x}bold_x, [𝐀]i,jsubscriptdelimited-[]𝐀𝑖𝑗\left[\mathbf{A}\right]_{i,j}[ bold_A ] start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT stands for the (i,j)𝑖𝑗(i,j)( italic_i , italic_j )-th entry, while [𝐱]isubscriptdelimited-[]𝐱𝑖\left[\mathbf{x}\right]_{i}[ bold_x ] start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT represents the i𝑖iitalic_i-th entry.

II System Model

II-A Channel Model

Refer to caption
Fig. 1: An illustration of the space-air-ground sensing system.

As shown in Fig. 1, we consider a space-air-ground sensing system, consisting of a multi-antenna LEO satellite, a cluster of K𝐾Kitalic_K aircraft, a multi-antenna terrestrial BS, and a gateway. Note that aircraft comprises but not limited to airplanes, unmanned aerial vehicles (UAVs), etc. In order to mitigate the path-loss effect in the system, a bistatic radar is considered. The LEO satellite as a sensing signal transmitter is equipped with Mtsubscript𝑀𝑡M_{t}italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT antennas and the BS as a sensing signal receiver is equipped with Mrsubscript𝑀𝑟{{M}_{r}}italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT antennas. Assume that all aircraft are in the circular area of the same plane; the k0subscript𝑘0k_{0}italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT-th aircraft is the sensed target while other K1𝐾1K-1italic_K - 1 aircraft are interference sources111Prior to optimization operations, the search for the existence of all aircraft in the designated area can be achieved by joint target search and communication channel estimation with omni-directional pilot signals [17].. Moreover, the gateway is deployed as a control center for collecting CSI and performing optimization algorithms.

II-A1 Satellite-to-Aircraft (S2A) Channel

The channel between satellite and the k𝑘kitalic_k-th aircraft is denoted as 𝐡S,k1×Mtsubscript𝐡𝑆𝑘superscript1subscript𝑀𝑡{\mathbf{h}_{S,k}}\in{{\mathbb{C}}^{1\times{{M}_{t}}}}bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT ∈ blackboard_C start_POSTSUPERSCRIPT 1 × italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. The mtsubscript𝑚𝑡m_{t}italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT-th element of 𝐡S,ksubscript��𝑆𝑘{\mathbf{h}_{S,k}}bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT is given as [𝐡S,k]mt=β0S2AdS,kαS2Aej(mt1)2πdλsin(ϑS,k)subscriptdelimited-[]subscript𝐡𝑆𝑘subscript𝑚𝑡superscriptsubscript𝛽0𝑆2𝐴superscriptsubscript𝑑𝑆𝑘subscript𝛼𝑆2𝐴superscript𝑒𝑗subscript𝑚𝑡12𝜋𝑑𝜆subscriptitalic-ϑ𝑆𝑘{{\left[{\mathbf{h}_{S,k}}\right]}_{{{m}_{t}}}}=\sqrt{\beta_{0}^{S2A}{{d}_{S,k% }}^{-\alpha_{S2A}}}{{e}^{-j\left({{m}_{t}}-1\right)\frac{2\pi d}{\lambda}\sin% \left({\vartheta}_{S,k}\right)}}[ bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUBSCRIPT = square-root start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_S 2 italic_A end_POSTSUPERSCRIPT italic_d start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - italic_α start_POSTSUBSCRIPT italic_S 2 italic_A end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG italic_e start_POSTSUPERSCRIPT - italic_j ( italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 1 ) divide start_ARG 2 italic_π italic_d end_ARG start_ARG italic_λ end_ARG roman_sin ( italic_ϑ start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT, where β0S2Asuperscriptsubscript𝛽0𝑆2𝐴\beta_{0}^{S2A}italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_S 2 italic_A end_POSTSUPERSCRIPT is the path-loss at a reference distance for the S2A channel, dS,ksubscript𝑑𝑆𝑘{{d}_{S,k}}italic_d start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT stands for the distance between the satellite and the k𝑘kitalic_k-th aircraft, αSAsubscript𝛼𝑆𝐴\alpha_{SA}italic_α start_POSTSUBSCRIPT italic_S italic_A end_POSTSUBSCRIPT is the satellite-to-aircraft path-loss exponent, d𝑑ditalic_d represents the antenna spacing, λ𝜆\lambdaitalic_λ indicates the wave length, and ϑS,ksubscriptitalic-ϑ𝑆𝑘{\vartheta}_{S,k}italic_ϑ start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT denotes the angle of arrival (AoA) from the transmitting antenna to the aircraft. Without loss of generality, we set the wave length as two times the antenna spacing, i.e. λ=2d𝜆2𝑑\lambda=2ditalic_λ = 2 italic_d.

II-A2 Aircraft-to-Ground (A2G) Channel

The uncorrelated Rician channel is used to model the channel responses, and the channel vector from the k𝑘kitalic_k-th aircraft to the BS is given by 𝐡k,B=βk,B𝐠k,Bsubscript𝐡𝑘𝐵subscript𝛽𝑘𝐵subscript𝐠𝑘𝐵\mathbf{h}_{k,{B}}=\sqrt{\beta_{k,B}}\mathbf{g}_{k,{B}}bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT = square-root start_ARG italic_β start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT end_ARG bold_g start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT, where βk,B=β0A2Gdk,BαA2Gsubscript𝛽𝑘𝐵superscriptsubscript𝛽0𝐴2𝐺superscriptsubscript𝑑𝑘𝐵subscript𝛼𝐴2𝐺{\beta}_{k,B}=\beta_{0}^{A2G}{{d}_{k,B}}^{-\alpha_{A2G}}italic_β start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT = italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_A 2 italic_G end_POSTSUPERSCRIPT italic_d start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - italic_α start_POSTSUBSCRIPT italic_A 2 italic_G end_POSTSUBSCRIPT end_POSTSUPERSCRIPT represents the large-scale fading coefficient, β0A2Gsuperscriptsubscript𝛽0𝐴2𝐺\beta_{0}^{A2G}italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_A 2 italic_G end_POSTSUPERSCRIPT is the path-loss at a reference distance for the A2G channel, αAGsubscript𝛼𝐴𝐺\alpha_{AG}italic_α start_POSTSUBSCRIPT italic_A italic_G end_POSTSUBSCRIPT is the aircraft-to-ground path-loss exponent. The small-scale fading vector 𝐠k,BMr×1subscript𝐠𝑘𝐵superscriptsubscript𝑀𝑟1\mathbf{g}_{k,B}\in{\mathbb{C}}^{{{M}_{r}}\times 1}bold_g start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT ∈ blackboard_C start_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT × 1 end_POSTSUPERSCRIPT is denoted as

𝐠k,B=𝐠¯k,B𝒦𝒦+1+𝐠~k,B1𝒦+1,subscript𝐠𝑘𝐵subscript¯𝐠𝑘𝐵𝒦𝒦1subscript~𝐠𝑘𝐵1𝒦1{\mathbf{g}_{k,B}}={\bar{\mathbf{g}}}_{k,B}\sqrt{\frac{\mathcal{K}}{\mathcal{K% }+1}}+{\tilde{\mathbf{g}}}_{k,B}\sqrt{\frac{1}{\mathcal{K}+1}},bold_g start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT = over¯ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT square-root start_ARG divide start_ARG caligraphic_K end_ARG start_ARG caligraphic_K + 1 end_ARG end_ARG + over~ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT square-root start_ARG divide start_ARG 1 end_ARG start_ARG caligraphic_K + 1 end_ARG end_ARG , (1)

where 𝒦𝒦\mathcal{K}caligraphic_K denotes to the Rician factor, 𝐠¯k,Bsubscript¯𝐠𝑘𝐵{\bar{\mathbf{g}}}_{k,B}over¯ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT and 𝐠~k,Bsubscript~𝐠𝑘𝐵{\tilde{\mathbf{g}}}_{k,B}over~ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT refer to the line-of-sight (LoS) and non-LoS (NLoS) components, respectively. Therefore, the A2G channel is written as

𝐡k,B=βk,B(𝐠¯k,B𝒦𝒦+1+𝐠~k,B1𝒦+1).subscript𝐡𝑘𝐵subscript𝛽𝑘𝐵subscript¯𝐠𝑘𝐵𝒦𝒦1subscript~𝐠𝑘𝐵1𝒦1{\mathbf{h}_{k,B}}=\sqrt{{{\beta}_{k,B}}}\left({{\bar{\mathbf{g}}}}_{k,B}\sqrt% {\frac{{\mathcal{K}}}{\mathcal{K}+1}}+{{\tilde{\mathbf{g}}}}_{k,B}\sqrt{\frac{% 1}{\mathcal{K}+1}}\right).bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT = square-root start_ARG italic_β start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT end_ARG ( over¯ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT square-root start_ARG divide start_ARG caligraphic_K end_ARG start_ARG caligraphic_K + 1 end_ARG end_ARG + over~ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT square-root start_ARG divide start_ARG 1 end_ARG start_ARG caligraphic_K + 1 end_ARG end_ARG ) . (2)

We assume that [𝐠¯k,B]mr=ej(mr1)2πdλsin(ϑk,B)subscriptdelimited-[]subscript¯𝐠𝑘𝐵subscript𝑚𝑟superscript𝑒𝑗subscript𝑚𝑟12𝜋𝑑𝜆subscriptitalic-ϑ𝑘𝐵{\left[{{{\bar{\mathbf{g}}}}_{k,B}}\right]}_{{{m}_{r}}}={e}^{-j\left({{m}_{r}}% -1\right)\frac{2\pi d}{\lambda}\sin\left({{\vartheta}_{k,B}}\right)}[ over¯ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_e start_POSTSUPERSCRIPT - italic_j ( italic_m start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - 1 ) divide start_ARG 2 italic_π italic_d end_ARG start_ARG italic_λ end_ARG roman_sin ( italic_ϑ start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT, where ϑk,Bsubscriptitalic-ϑ𝑘𝐵{{\vartheta}_{k,B}}italic_ϑ start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT denotes the AoA from the aircraft to the terrestrial BS, while the Rayleigh distributed term 𝐠~k,Bsubscript~𝐠𝑘𝐵{\tilde{\mathbf{g}}}_{k,B}over~ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT describes the scattered multi-path and 𝐠~k,B𝒞𝒩(𝟎Mr,𝐈Mr)similar-tosubscript~𝐠𝑘𝐵𝒞𝒩subscript0subscript𝑀𝑟subscript𝐈subscript𝑀𝑟{{\tilde{\mathbf{g}}}_{k,B}}\sim\mathcal{C}\mathcal{N}\left({\mathbf{0}_{{{M}_% {r}}}},{\mathbf{I}_{{M}_{r}}}\right)over~ start_ARG bold_g end_ARG start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT ∼ caligraphic_C caligraphic_N ( bold_0 start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT , bold_I start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT ).

II-B Sensing Signal Transmission

Denote 𝐭Mt×1𝐭superscriptsubscript𝑀𝑡1\mathbf{t}\in{{\mathbb{C}}^{{{M}_{t}}\times 1}}bold_t ∈ blackboard_C start_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT × 1 end_POSTSUPERSCRIPT as the transmit beamforming vector and use q𝒞𝒩(0,1)𝑞𝒞𝒩01q\in\mathcal{CN}(0,1)italic_q ∈ caligraphic_C caligraphic_N ( 0 , 1 ) as the information symbol. The transmitted sensing signal 𝐱𝐱\mathbf{x}bold_x from the LEO satellite is given as 𝐱=𝐭q𝐱𝐭𝑞\mathbf{x}=\mathbf{t}qbold_x = bold_t italic_q, and the echo signal yk1×1subscript𝑦𝑘superscript11{{y}_{k}}\in{{\mathbb{C}}^{1\times 1}}italic_y start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ∈ blackboard_C start_POSTSUPERSCRIPT 1 × 1 end_POSTSUPERSCRIPT reflected from the k𝑘kitalic_k-th aircraft is yk=𝐡S,k𝐱=𝐡S,k𝐭qsubscript𝑦𝑘subscript𝐡𝑆𝑘𝐱subscript𝐡𝑆𝑘𝐭𝑞{{y}_{k}}={\mathbf{h}_{S,k}}\mathbf{x}={\mathbf{h}_{S,k}}\mathbf{t}qitalic_y start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT bold_x = bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT bold_t italic_q. Then, the echo signal 𝐲eMr×1subscript𝐲𝑒superscriptsubscript𝑀𝑟1\mathbf{y}_{e}\in{{\mathbb{C}}^{{{M}_{r}}\times 1}}bold_y start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∈ blackboard_C start_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT × 1 end_POSTSUPERSCRIPT received at the BS is expressed as

𝐲esubscript𝐲𝑒\displaystyle\mathbf{y}_{e}bold_y start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT =𝐡k0,B𝐡S,k0𝐭q+kk0K𝐡k,B𝐡S,k𝐭q+𝝎,absentsubscript𝐡subscript𝑘0𝐵subscript𝐡𝑆subscript𝑘0𝐭𝑞superscriptsubscript𝑘subscript𝑘0𝐾subscript𝐡𝑘𝐵subscript𝐡𝑆𝑘𝐭𝑞𝝎\displaystyle={\mathbf{h}_{{{k}_{0}},B}}{\mathbf{h}_{S,{{k}_{0}}}}\mathbf{t}q+% \sum\limits_{k\neq k_{0}}^{K}{\mathbf{h}_{k,B}}{\mathbf{h}_{S,k}}\mathbf{t}q+% \bm{\omega},= bold_h start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT bold_t italic_q + ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT bold_t italic_q + bold_italic_ω , (3)

Define the finite impulse response (FIR) filter at BS as 𝐮𝐮\mathbf{u}bold_u, and 𝐮Mr×1𝐮superscriptsubscript𝑀𝑟1\mathbf{u}\in{{\mathbb{C}}^{{{M}_{r}}\times 1}}bold_u ∈ blackboard_C start_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT × 1 end_POSTSUPERSCRIPT. Then, the filter output at the BS is given by

yo=𝐮H𝐲e.subscript𝑦𝑜superscript𝐮𝐻subscript𝐲𝑒y_{o}=\mathbf{u}^{H}\mathbf{y}_{e}.italic_y start_POSTSUBSCRIPT italic_o end_POSTSUBSCRIPT = bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_y start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT . (4)

Finally, the SINR of the target aircraft is

SINRk0=|𝐮H𝐡k0,B𝐡S,k0𝐭|2𝐮H(kk0K|𝐡k,B𝐡S,k𝐭|2+|𝝎|2)𝐮.subscriptSINRsubscript𝑘0superscriptsuperscript𝐮𝐻subscript𝐡subscript𝑘0𝐵subscript𝐡𝑆subscript𝑘0𝐭2superscript𝐮𝐻superscriptsubscript𝑘subscript𝑘0𝐾superscriptsubscript𝐡𝑘𝐵subscript𝐡𝑆𝑘𝐭2superscript𝝎2𝐮{\text{SINR}_{{{k}_{0}}}}=\frac{{\left|{\mathbf{u}^{H}}{\mathbf{h}_{{{k}_{0}},% {B}}}{\mathbf{h}_{S,{{k}_{0}}}}\mathbf{t}\right|}^{2}}{{\mathbf{u}^{H}}\left(% \sum\limits_{k\neq{{k}_{0}}}^{K}{{{\left|{\mathbf{h}_{k,{B}}}{\mathbf{h}_{S,k}% }\mathbf{t}\right|}^{2}}}+{{\left|\bm{\omega}\right|}^{2}}\right)\mathbf{u}}.SINR start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG | bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_h start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT bold_t | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT bold_t | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | bold_italic_ω | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) bold_u end_ARG . (5)

III Problem Formulation and Solution

In this paper, we assume the target k0subscript𝑘0k_{0}italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT-th aircraft is in emergent conditions and higher QoS of sensing signal is required for better detection and localization. Thus, we aim to maximize the SINRk0subscriptSINRsubscript𝑘0\text{SINR}_{k_{0}}SINR start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT of the echo signal from the target aircraft by jointly designing the transmit beamforming vector 𝐭𝐭\mathbf{t}bold_t and the receive filtering vector 𝐮𝐮\mathbf{u}bold_u. The problem of interest can be written as

𝒫1:max𝐭,𝐮f1(𝐭,𝐮)=SINRk0:𝒫1subscriptmax𝐭𝐮subscript𝑓1𝐭𝐮subscriptSINRsubscript𝑘0\displaystyle\mathcal{P}1:\mathop{\mathrm{max}}\limits_{\bm{\mathrm{t},\mathrm% {u}}}\quad f_{1}\left(\mathbf{t},\mathbf{u}\right)=\text{SINR}_{k_{0}}caligraphic_P 1 : roman_max start_POSTSUBSCRIPT bold_t bold_, bold_u end_POSTSUBSCRIPT italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_t , bold_u ) = SINR start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT (6a)
s.t.R1:𝐭H𝐭Pt,:s.t.R1superscript𝐭𝐻𝐭subscript𝑃𝑡\text{s.t.}\quad\mathrm{R}1:\mathbf{t}^{H}\mathbf{t}\leq P_{t},s.t. R1 : bold_t start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_t ≤ italic_P start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , (6b)
R2:𝐮H𝐮=1,:R2superscript𝐮𝐻𝐮1\quad\quad\mathrm{R}2:\mathbf{u}^{H}\mathbf{u}=1,R2 : bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u = 1 , (6c)

where R1R1\mathrm{R}1R1 is the sensing power budget requirement, R2R2\mathrm{R}2R2 is the receive filtering requirement. Note that the problem (𝒫1𝒫1\mathcal{P}1caligraphic_P 1) is a non-convex fractional programming expression which requires further operation to get the final optimality. In the following, we will divide the problem (𝒫1𝒫1\mathcal{P}1caligraphic_P 1) into two sub-problems to solve it based on the alternating optimization.

III-A Transmit Beamforming Vector Design

For fixed 𝐮𝐮\mathbf{u}bold_u, we first introduce an auxiliary variable ϖitalic-ϖ\varpiitalic_ϖ to track the closed-form solution of 𝐭𝐭\mathbf{t}bold_t. According to [18, 19], the fractional programming problem maxx{f(x)/g(x)}subscript𝑥𝑓𝑥𝑔𝑥{\max}_{x}\left\{f\left(x\right)/g\left(x\right)\right\}roman_max start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT { italic_f ( italic_x ) / italic_g ( italic_x ) } can be converted into max{x,y}{2yf(x)y2g(x)}subscript𝑥𝑦2𝑦𝑓𝑥superscript𝑦2𝑔𝑥{\max}_{\left\{x,y\right\}}\left\{2y\sqrt{f\left(x\right)}-{{y}^{2}}g\left(x% \right)\right\}roman_max start_POSTSUBSCRIPT { italic_x , italic_y } end_POSTSUBSCRIPT { 2 italic_y square-root start_ARG italic_f ( italic_x ) end_ARG - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g ( italic_x ) }, where y𝑦yitalic_y is a new auxiliary variable.

Denote 𝐪=𝐮H𝐡k0,B𝐡S,k0𝐪superscript𝐮𝐻subscript𝐡subscript𝑘0𝐵subscript𝐡𝑆subscript𝑘0\mathbf{q}=\mathbf{u}^{H}{\mathbf{h}_{{{k}_{0}},{B}}}{\mathbf{h}_{S,{{k}_{0}}}}bold_q = bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_h start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT and 𝐫=𝐮H𝐡k,B𝐡S,k𝐫superscript𝐮𝐻subscript𝐡𝑘𝐵subscript𝐡𝑆𝑘\mathbf{r}=\mathbf{u}^{H}{\mathbf{h}_{k,{B}}}{\mathbf{h}_{S,k}}bold_r = bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT, where 𝐪1×Mt𝐪superscript1subscript𝑀𝑡\mathbf{q}\in{{\mathbb{C}}^{1\times{M}_{t}}}bold_q ∈ blackboard_C start_POSTSUPERSCRIPT 1 × italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUPERSCRIPT and 𝐫1×Mt𝐫superscript1subscript𝑀𝑡\mathbf{r}\in{{\mathbb{C}}^{1\times{M}_{t}}}bold_r ∈ blackboard_C start_POSTSUPERSCRIPT 1 × italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. With an auxiliary variable ϖitalic-ϖ\varpiitalic_ϖ, the problem (𝒫1𝒫1\mathcal{P}1caligraphic_P 1) can be rewritten as

𝒫2:max𝐭,ϖ,𝐮f2(𝐭,ϖ;𝐮):𝒫2subscriptmax𝐭italic-ϖ𝐮subscript𝑓2𝐭italic-ϖ𝐮\displaystyle\mathcal{P}2:\mathop{\mathrm{max}}\limits_{\mathbf{t},\varpi,% \mathbf{u}}f_{2}\left(\mathbf{t},\varpi;\mathbf{u}\right)caligraphic_P 2 : roman_max start_POSTSUBSCRIPT bold_t , italic_ϖ , bold_u end_POSTSUBSCRIPT italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_t , italic_ϖ ; bold_u ) (7a)
s.t.(6b),(6c),s.t.6b6c\text{s.t.}(\mathrm{\ref{R1}}),(\mathrm{\ref{R2}}),s.t. ( ) , ( ) , (7b)

where the new objective function is

f2(𝐭,ϖ;𝐮)subscript𝑓2𝐭italic-ϖ𝐮\displaystyle f_{2}\left(\mathbf{t},\varpi;\mathbf{u}\right)italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_t , italic_ϖ ; bold_u ) (8)
=2{ϖH𝐪𝐭}|ϖ|2(kk0K|𝐫𝐭|2+|𝐮H𝝎|2),absent2superscriptitalic-ϖ𝐻𝐪𝐭superscriptitalic-ϖ2superscriptsubscript𝑘subscript𝑘0𝐾superscript𝐫𝐭2superscriptsuperscript𝐮𝐻𝝎2\displaystyle=2\mathcal{R}\left\{\varpi^{H}\mathbf{q}\mathbf{t}\right\}-\left|% \varpi\right|^{2}\left(\sum\limits_{k\neq k_{0}}^{K}{{{\left|\mathbf{r}\mathbf% {t}\right|}^{2}}}+{{\left|\mathbf{u}^{H}\bm{\omega}\right|}^{2}}\right),= 2 caligraphic_R { italic_ϖ start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_qt } - | italic_ϖ | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_rt | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_italic_ω | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ,

and ϖitalic-ϖ\varpiitalic_ϖ is updated according to [20], which is expressed as

ϖ=𝐪𝐭kk0K|𝐫𝐭|2+|𝐮H𝝎|2.italic-ϖ𝐪𝐭superscriptsubscript𝑘subscript𝑘0𝐾superscript𝐫𝐭2superscriptsuperscript𝐮𝐻𝝎2\varpi=\frac{\mathbf{q}\mathbf{t}}{\sum\limits_{k\neq{{k}_{0}}}^{K}{{{\left|% \mathbf{r}\mathbf{t}\right|}^{2}}}+{{\left|{\mathbf{u}^{H}}\bm{\omega}\right|}% ^{2}}}.italic_ϖ = divide start_ARG bold_qt end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_rt | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_italic_ω | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (9)

Accordingly, the optimal 𝐭superscript𝐭\mathbf{t}^{*}bold_t start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT can be expressed as

𝐭=argmax𝐭,ϖ;𝐮f2(𝐭,ϖ;𝐮).superscript𝐭𝐭italic-ϖ𝐮subscript𝑓2𝐭italic-ϖ𝐮{\mathbf{t}^{*}}=\underset{\mathbf{t},\varpi;\mathbf{u}}{\mathop{\arg\max}}\,{% {f}_{2}}\left(\mathbf{t},\varpi;\mathbf{u}\right).bold_t start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = start_UNDERACCENT bold_t , italic_ϖ ; bold_u end_UNDERACCENT start_ARG roman_arg roman_max end_ARG italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_t , italic_ϖ ; bold_u ) . (10)

The Lagrangian of the associated problem (𝒫2𝒫2\mathcal{P}2caligraphic_P 2) is defined as

(𝐭,ϖ,λ;𝐮)=f2(𝐭,ϖ;𝐮)λ(𝐭H𝐭1).𝐭italic-ϖ𝜆𝐮subscript𝑓2𝐭italic-ϖ𝐮𝜆superscript𝐭𝐻𝐭1\mathcal{L}\left(\mathbf{t},\varpi,\lambda;\mathbf{u}\right)=f_{2}\left(% \mathbf{t},\varpi;\mathbf{u}\right)-\lambda\left({\mathbf{t}^{H}}\mathbf{t}-1% \right).caligraphic_L ( bold_t , italic_ϖ , italic_λ ; bold_u ) = italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_t , italic_ϖ ; bold_u ) - italic_λ ( bold_t start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_t - 1 ) . (11)

Based on the Lagrangian, four conditions, also known as Karush-Kuhn-Tucker (KKT) conditions, are listed as

𝐭H𝐭1=0,superscript𝐭𝐻𝐭10{\mathbf{t}^{H}}\mathbf{t}-1=0,bold_t start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_t - 1 = 0 , (12a)
λ0,𝜆0\lambda\geq 0,italic_λ ≥ 0 , (12b)
λ(𝐭H𝐭1)=0,𝜆superscript𝐭𝐻𝐭10\lambda\left({\mathbf{t}^{H}}\mathbf{t}-1\right)=0,italic_λ ( bold_t start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_t - 1 ) = 0 , (12c)
t(𝐭,ϖ,λ;𝐮)=2ϖ𝐪H2𝚵𝐭2λ𝐈Mt𝐭=0,subscript𝑡𝐭italic-ϖ𝜆𝐮2italic-ϖsuperscript𝐪𝐻2𝚵𝐭2𝜆subscript𝐈subscript𝑀𝑡𝐭0{{\nabla}_{t}}\mathcal{L}\left(\mathbf{t},\varpi,\lambda;\mathbf{u}\right)=2% \varpi{\mathbf{q}^{H}}-2{\bm{\Xi}}\mathbf{t}-2\lambda{\mathbf{I}_{{{M}_{t}}}}% \mathbf{t}=0,∇ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT caligraphic_L ( bold_t , italic_ϖ , italic_λ ; bold_u ) = 2 italic_ϖ bold_q start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT - 2 bold_Ξ bold_t - 2 italic_λ bold_I start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUBSCRIPT bold_t = 0 , (12d)

where 𝚵=ϖ2kk0K𝐫H𝐫Mt×Mt𝚵superscriptitalic-ϖ2superscriptsubscript𝑘subscript𝑘0𝐾superscript𝐫𝐻𝐫superscriptsubscript𝑀𝑡subscript𝑀𝑡\bm{\Xi}={{\varpi}^{2}}\sum\limits_{k\neq{{k}_{0}}}^{K}{\mathbf{r}^{H}\mathbf{% r}}\in{{\mathbb{C}}^{{{M}_{t}}\times{{M}_{t}}}}bold_Ξ = italic_ϖ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT bold_r start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_r ∈ blackboard_C start_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT × italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. Then, the closed-form solution can be expressed as

𝐭=ϖ(𝚵+λ𝐈Mt)1𝐪H.superscript𝐭italic-ϖsuperscript𝚵𝜆subscript𝐈subscript𝑀𝑡1superscript𝐪𝐻{\mathbf{t}^{*}}=\varpi{{\left(\bm{\Xi}+\lambda{\mathbf{I}_{{{M}_{t}}}}\right)% }^{-1}}{\mathbf{q}^{H}}.bold_t start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = italic_ϖ ( bold_Ξ + italic_λ bold_I start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_q start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT . (13)

The optimal dual variable λ𝜆\lambdaitalic_λ can be searched efficiently via bisection with appropriate upper and lower bound [21]. The calculation of 𝐭𝐭\mathbf{t}bold_t is represented in Algorithm 1 based on the bisection method, where ε1subscript𝜀1\varepsilon_{1}italic_ε start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and ε2subscript𝜀2\varepsilon_{2}italic_ε start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are tolerance values.

Algorithm 1 Transmit Beamforming Vector Update Algorithm
1:Input: λminsubscript𝜆min\lambda_{\text{min}}italic_λ start_POSTSUBSCRIPT min end_POSTSUBSCRIPT, λmaxsubscript𝜆max\lambda_{\text{max}}italic_λ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT, and λ=12(λmin+λmax)𝜆12subscript𝜆minsubscript𝜆max\lambda=\frac{1}{2}\left(\lambda_{\text{min}}+\lambda_{\text{max}}\right)italic_λ = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_λ start_POSTSUBSCRIPT min end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT )
2:repeat
3:     update 𝐭𝐭\mathbf{t}bold_t using equation (13)
4:     if 𝐭>Pdnorm𝐭subscript𝑃𝑑\|\mathbf{t}\|>P_{d}∥ bold_t ∥ > italic_P start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT then
5:         λmin=λsubscript𝜆min𝜆\lambda_{\text{min}}=\lambdaitalic_λ start_POSTSUBSCRIPT min end_POSTSUBSCRIPT = italic_λ
6:     else
7:         λmax=λsubscript𝜆max𝜆\lambda_{\text{max}}=\lambdaitalic_λ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT = italic_λ
8:     end if
9:until |λminλmax|<ε1subscript𝜆minsubscript𝜆maxsubscript𝜀1\left|\lambda_{\text{min}}-\lambda_{\text{max}}\right|<\varepsilon_{1}| italic_λ start_POSTSUBSCRIPT min end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT max end_POSTSUBSCRIPT | < italic_ε start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and |𝐭Pd|<ε2norm𝐭subscript𝑃𝑑subscript𝜀2\left|\|\mathbf{t}\|-P_{d}\right|<\varepsilon_{2}| ∥ bold_t ∥ - italic_P start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT | < italic_ε start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

III-B Receive Filtering Vector Design

When 𝐭𝐭\mathbf{t}bold_t is fixed, we design the solution of 𝐮𝐮\mathbf{u}bold_u for (𝒫1𝒫1\mathcal{P}1caligraphic_P 1) using the generalized Rayleigh quotient method and the stagnation point method [13].

Denote 𝐍=|𝐡k0,B𝐡S,k0𝐭|2𝐍superscriptsubscript𝐡subscript𝑘0𝐵subscript𝐡𝑆subscript𝑘0𝐭2\mathbf{N}={\left|{\mathbf{h}_{{{k}_{0}},{B}}}{\mathbf{h}_{S,{{k}_{0}}}}% \mathbf{t}\right|}^{2}bold_N = | bold_h start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT bold_t | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and 𝐃=kk0K|𝐡k,B𝐡S,k𝐭|2+|𝝎|2𝐃superscriptsubscript𝑘subscript𝑘0𝐾superscriptsubscript𝐡𝑘𝐵subscript𝐡𝑆𝑘𝐭2superscript𝝎2\mathbf{D}=\sum\limits_{k\neq k_{0}}^{K}{{{\left|{\mathbf{h}_{k,{B}}}{\mathbf{% h}_{S,k}}\mathbf{t}\right|}^{2}}}+{{\left|\bm{\omega}\right|}^{2}}bold_D = ∑ start_POSTSUBSCRIPT italic_k ≠ italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_K end_POSTSUPERSCRIPT | bold_h start_POSTSUBSCRIPT italic_k , italic_B end_POSTSUBSCRIPT bold_h start_POSTSUBSCRIPT italic_S , italic_k end_POSTSUBSCRIPT bold_t | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | bold_italic_ω | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, the optimal 𝐮superscript𝐮\mathbf{u}^{*}bold_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT can be written as

𝐮=argmax𝐮𝐮H𝐍𝐮𝐮H𝐃𝐮,superscript𝐮𝐮superscript𝐮𝐻𝐍𝐮superscript𝐮𝐻𝐃𝐮{\mathbf{u}^{*}}=\underset{\mathbf{u}}{\mathop{\arg\max}}\,\frac{\mathbf{u}^{H% }\mathbf{N}\mathbf{u}}{\mathbf{u}^{H}\mathbf{D}\mathbf{u}},bold_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = underbold_u start_ARG roman_arg roman_max end_ARG divide start_ARG bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Nu end_ARG start_ARG bold_u start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Du end_ARG , (14)

where 𝐍𝐍\mathbf{N}bold_N is a Hermitian matrix while 𝐃𝐃\mathbf{D}bold_D is a positive definite Hermitian matrix. Consider a unitary matrix 𝐕𝐕\mathbf{V}bold_V which satisfies 𝐕H=𝐕1superscript𝐕𝐻superscript𝐕1\mathbf{V}^{H}=\mathbf{V}^{-1}bold_V start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT = bold_V start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT and 𝐕H𝐕=𝐈Mrsuperscript𝐕𝐻𝐕subscript𝐈subscript𝑀𝑟\mathbf{V}^{H}\mathbf{V}=\mathbf{I}_{M_{r}}bold_V start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_V = bold_I start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT, we can re-express 𝐃𝐃\mathbf{D}bold_D as 𝐃=𝐕diag(λ1,λ2,,λMr)𝐕H𝐃𝐕diagsubscript𝜆1subscript𝜆2subscript𝜆subscript𝑀𝑟superscript𝐕𝐻\mathbf{D}=\mathbf{V}\text{diag}\left({{\lambda}_{1}},{{\lambda}_{2}},\ldots,{% \lambda}_{{M}_{r}}\right){\mathbf{V}^{H}}bold_D = bold_V diag ( italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , … , italic_λ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) bold_V start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT, and denote 𝐋=𝐕diag(λ1,λ2,,λMr)𝐕H𝐋𝐕diagsubscript𝜆1subscript𝜆2subscript𝜆subscript𝑀𝑟superscript𝐕𝐻\mathbf{L}=\mathbf{V}\text{diag}\left(\sqrt{{\lambda}_{1}},\sqrt{{\lambda}_{2}% },\ldots,\sqrt{{\lambda}_{{M}_{r}}}\right){\mathbf{V}^{H}}bold_L = bold_V diag ( square-root start_ARG italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG , square-root start_ARG italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG , … , square-root start_ARG italic_λ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG ) bold_V start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT, where for i=1,2,,Mr𝑖12subscript𝑀𝑟i=1,2,\ldots,M_{r}italic_i = 1 , 2 , … , italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, λisubscript𝜆𝑖{\lambda}_{i}italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT are positive eigenvalues of 𝐕𝐕\mathbf{V}bold_V. Based on eigenvalue decomposition theory, we have 𝐃=𝐋𝐋𝐃𝐋𝐋\mathbf{D}=\mathbf{L}\cdot\mathbf{L}bold_D = bold_L ⋅ bold_L. Take 𝐩=𝐋𝐮𝐩𝐋𝐮\mathbf{p}=\mathbf{L}\mathbf{u}bold_p = bold_Lu, the original objection function can be expressed as an equivalent form of

𝒫3:max𝐩𝐩H𝐋1𝐍𝐋1𝐩\mathcal{P}3:\quad\mathop{\mathrm{max}}\limits_{\bm{\mathrm{p}}}\quad{\mathbf{% p}^{H}}{\mathbf{L}^{-1}}\mathbf{N}{\mathbf{L}^{-1}}\mathbf{p}caligraphic_P 3 : roman_max start_POSTSUBSCRIPT bold_p end_POSTSUBSCRIPT bold_p start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_NL start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_p (15a)
s.t.𝐩H𝐩=C,s.t.superscript𝐩𝐻𝐩𝐶\text{s.t.}\quad\mathbf{p}^{H}\mathbf{p}=C,s.t. bold_p start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_p = italic_C , (15b)

where C𝐶Citalic_C is a constant which is independent of the final optimality. To attain its optimal solution, the Lagrangian L(𝐩,ν)=𝐩H𝐋1𝐍𝐋1𝐩ν(𝐩H𝐩C)𝐿𝐩𝜈superscript𝐩𝐻superscript𝐋1superscript𝐍𝐋1𝐩𝜈superscript𝐩𝐻𝐩𝐶L\left(\mathbf{p},\nu\right)={\mathbf{p}^{H}}{\mathbf{L}^{-1}}\mathbf{N}{% \mathbf{L}^{-1}}\mathbf{p}-\nu\left({\mathbf{p}^{H}}\mathbf{p}-C\right)italic_L ( bold_p , italic_ν ) = bold_p start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_NL start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_p - italic_ν ( bold_p start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_p - italic_C ) is introduced. By following the optimality condition 𝐩(𝐩,ν)=2𝐋1𝐍𝐋1𝐩2ν𝐩=0𝐋1𝐍𝐋1𝐩=ν𝐩subscript𝐩𝐩𝜈2superscript𝐋1superscript𝐍𝐋1𝐩2𝜈𝐩0superscript𝐋1superscript𝐍𝐋1𝐩𝜈𝐩{{\nabla}_{\mathbf{p}}}\mathcal{L}\left(\mathbf{p},\nu\right)=2{\mathbf{L}^{-1% }}\mathbf{N}{\mathbf{L}^{-1}}\mathbf{p}-2\nu\mathbf{p}=0\Rightarrow{\mathbf{L}% ^{-1}}\mathbf{N}{\mathbf{L}^{-1}}\mathbf{p}=\nu\mathbf{p}∇ start_POSTSUBSCRIPT bold_p end_POSTSUBSCRIPT caligraphic_L ( bold_p , italic_ν ) = 2 bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_NL start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_p - 2 italic_ν bold_p = 0 ⇒ bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_NL start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_p = italic_ν bold_p, the optimal solution 𝐩superscript𝐩\mathbf{p}^{*}bold_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT is the eigenvector corresponding to the maximum eigenvalue of 𝐋1𝐍𝐋1superscript𝐋1superscript𝐍𝐋1{\mathbf{L}^{-1}}\mathbf{N}{\mathbf{L}^{-1}}bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_NL start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. Then, the closed-form analytical solution of the previous equation (14) is

𝐮=𝐋1𝐩𝐩,superscript𝐮superscript𝐋1superscript𝐩normsuperscript𝐩\mathbf{u}^{*}=\frac{\mathbf{L}^{-1}\mathbf{p}^{*}}{\|\mathbf{p}^{*}\|},bold_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = divide start_ARG bold_L start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG ∥ bold_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ∥ end_ARG , (16)

with the attained 𝐩superscript𝐩\mathbf{p}^{*}bold_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT. Finally, the maximum SINR of the sensing signal from the target k0subscript𝑘0k_{0}italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT-th aircraft is represented as

SINRk0=(𝐮)H𝐍𝐮(𝐮)H𝐃𝐮.superscriptsubscriptSINRsubscript𝑘0superscriptsuperscript𝐮𝐻superscript𝐍𝐮superscriptsuperscript𝐮𝐻superscript𝐃𝐮\text{SINR}_{{{k}_{0}}}^{*}=\frac{{{\left({\mathbf{u}^{*}}\right)}^{H}}\mathbf% {N}{\mathbf{u}^{*}}}{{{\left({\mathbf{u}^{*}}\right)}^{H}}\mathbf{D}{\mathbf{u% }^{*}}}.SINR start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = divide start_ARG ( bold_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Nu start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG ( bold_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_Du start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG . (17)

III-C Algorithm Summary and Complexity Analysis

Combining the above derivations, the proposed alternating optimization from joint transmit beamforming and receive filtering is summarized in Algorithm 2, where ε3subscript𝜀3\varepsilon_{3}italic_ε start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT is the tolerance value for iterated SINR.

Thereinto, with the matrix inverse operation and the denoted iteration number Iλsubscript𝐼𝜆I_{\lambda}italic_I start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT of the bisection search for λ𝜆\lambdaitalic_λ, the computational complexity to update 𝐭𝐭\mathbf{t}bold_t is approximately calculated as 𝒪(IλMr)𝒪subscript𝐼𝜆subscript𝑀𝑟\mathcal{O}\left(I_{\lambda}{M_{r}}\right)caligraphic_O ( italic_I start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ); while the update of 𝐮𝐮\mathbf{u}bold_u, which comprises an eigenvalue decomposition, requires the complexity of 𝒪(Mr3)𝒪superscriptsubscript𝑀𝑟3\mathcal{O}\left({M_{r}}^{3}\right)caligraphic_O ( italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ). Hence, the complexity of the proposed algorithm at each iteration is 𝒪(IλMr+Mr3)𝒪subscript𝐼𝜆subscript𝑀𝑟superscriptsubscript𝑀𝑟3\mathcal{O}\left(I_{\lambda}{M_{r}}+{M_{r}}^{3}\right)caligraphic_O ( italic_I start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT + italic_M start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ).

Algorithm 2 Alternating Optimization Algorithm for 𝒫1𝒫1\mathcal{P}1caligraphic_P 1
1:Initialize: 𝐭𝐭\mathbf{t}bold_t, 𝐮𝐮\mathbf{u}bold_u, set i=0𝑖0i=0italic_i = 0, SINR(0)=0superscriptSINR00\text{SINR}^{(0)}=0SINR start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = 0
2:repeat
3:     update ϖitalic-ϖ\varpiitalic_ϖ using equation (9)
4:     update 𝐭𝐭\mathbf{t}bold_t using equation (13)
5:     update 𝐮𝐮\mathbf{u}bold_u using equation (16)
6:     i=i+1𝑖𝑖1i=i+1italic_i = italic_i + 1
7:     calculate SINR(i)superscriptSINR𝑖\text{SINR}^{(i)}SINR start_POSTSUPERSCRIPT ( italic_i ) end_POSTSUPERSCRIPT with equation (5)
8:until |SINR(i)SINR(i1)|<ε3superscriptSINR𝑖superscriptSINR𝑖1subscript𝜀3\left|\text{SINR}^{(i)}-\text{SINR}^{(i-1)}\right|<\varepsilon_{3}| SINR start_POSTSUPERSCRIPT ( italic_i ) end_POSTSUPERSCRIPT - SINR start_POSTSUPERSCRIPT ( italic_i - 1 ) end_POSTSUPERSCRIPT | < italic_ε start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT.

IV Simulation Results

In this section, simulation results are provided to validate the proposed algorithm in the studied bistatic space-air-ground sensing system. Unless otherwise noted, we set the S2A channel path-loss αS2A=2.0subscript𝛼𝑆2𝐴2.0\alpha_{S2A}=2.0italic_α start_POSTSUBSCRIPT italic_S 2 italic_A end_POSTSUBSCRIPT = 2.0, the A2G channel path-loss αA2G=2.2subscript𝛼𝐴2𝐺2.2\alpha_{A2G}=2.2italic_α start_POSTSUBSCRIPT italic_A 2 italic_G end_POSTSUBSCRIPT = 2.2, and the noise power is assumed to be 1×10121superscript10121\times 10^{-12}1 × 10 start_POSTSUPERSCRIPT - 12 end_POSTSUPERSCRIPT W. Moreover, the LEO satellite and the BS are set to be located at (10,20,300)1020300\left(-10,20,300\right)( - 10 , 20 , 300 ) km and (0,0,0)000\left(0,0,0\right)( 0 , 0 , 0 ) km, respectively. A total of K=4𝐾4K=4italic_K = 4 aircraft, including one target aircraft and three other interfering aircraft, are uniformly distributed in a circular aerial region with a radius of 10101010 km centering at (20,20,10)202010\left(20,20,10\right)( 20 , 20 , 10 ) km, i.e. all of them are at an altitude of 10101010 km above the ground. In addition, four benchmarks are provided to evaluate the performance of the proposed algorithm: 1) Optimization of receive filtering vector alone; 2) Optimization of the transmit beamforming vector alone; 3) Random receive filtering vectors and transmit beamforming vectors; 4) Transceiving mode where LEO satellite transmits and receives the sensing signal.

Refer to caption
Fig. 2: Convergence validation.
Refer to caption
Fig. 3: Optimal SINR VS satellite altitude.
Refer to caption
Fig. 4: Optimal SINR VS satellite antenna number.

Fig. 4 presents the SINR performance of the received echo signal at the filter output. It is obvious that the SINR increases automatically and converges after only 3 iterations, which validates the efficiency of the proposed algorithm. Additionally, comparing the proposed algorithm with benchmarks, we notice that by optimizing the receive filtering vector alone, a relatively high SINR value can be reached, but the maximum value cannot be obtained. However, optimizing the transmit beamforming vector alone, though it works, has incremental contributions toward the overall performance. This is because the effect of optimization over the transmit beamforming vector can be weakened by the large-scale fading, while the ground receiver’s local optimization over the receive filtering vector is dominant.

In Fig. 4, the optimal SINR under the bistatic mode is compared with that under the transceiving mode at different satellite altitudes, and that by optimizing the receive filtering vector alone. In all cases, the optimal SINR declines as the satellite altitude increases. However, it is noted that the optimization under the bistatic mode can bring a significantly higher SINR than that under the transceiving mode. This has clearly emphasized the effectiveness of using the bistatic mode rather than the transceiving mode for sensing in SAGIN. Moreover, the SINR by optimizing the receive filtering vector alone is lower than that by the joint optimization algorithm, while the gap in between narrows as the satellite altitude increases. This can be due to the fact that the path loss over a longer distance gradually weakens the effect of the transmit beamforming optimization.

Fig. 4 describes the optimal SINR against different numbers of satellite antennas under different transmit power. It is clear that with a fixed antenna number, the increase of the transmit power will brings higher SINR. Moreover, for different transmit power, the growth rate of SINR when increasing the antenna number can be different; the more power there is, the higher growth rate of SINR will be. This is because adding more beams may additionally mitigate interference power, and the rise in desired signal due to increased transmit power is not directly proportional to the increase in interference power. This provides valuable network design insights for the balance of the satellite antenna number and the transmit power when resources are limited.

V Conclusion

In this paper, we studied for the first time the sensing signal performance in the SAGIN, where the LEO satellite and the terrestrial BS were separately used as a transmitter and a radar receiver to construct a bistatic sensing system. Following this setup, we jointly optimized the transmit beamforming vector and the receive filtering vector to maximize the SINR of echo signal from the target, while the closed-form expressions of both vectors were returned. Simulation results validated that the proposed algorithm can effectively improve the sensing performance in the SAGIN; meanwhile, it is manifested that the overall performance can be significantly influenced by the satellite altitude, transmit power and antenna numbers.

VI Acknowledgement

This work was supported in part by the YongRiver Scientific and Technological Innovation Project under Grant 2023A-187-G. The work of Qingqing Wu was supported by the National Natural Science Foundation of China under Grant 62371289 and Grant 62331022.

References

  • [1] A. Yazar, S. Dogan-Tusha, and H. Arslan, “6g vision: An ultra-flexible perspective,” ITU J. Future Evol. Technol., vol. 1, no. 1, pp. 121–140, Jan. 2020.
  • [2] X. Li, S. Guo, and S. Malik, “Sensing-aided peer-to-peer millimeter-wave communication,” Int. J. Inf. Technol., vol. 16, no. 4, pp. 2069–2075, Apr. 2024.
  • [3] R. Shrestha, I. Oh, and S. Kim, “A survey on operation concept, advancements, and challenging issues of urban air traffic management,” Frontiers Future Transp., vol. 2, p. 1, 2021.
  • [4] K. Meng, Q. Wu, C. Masouros, W. Chen, and D. Li, “Intelligent surface empowered integrated sensing and communication: From coexistence to reciprocity,” arXiv:2311.00418, 2023.
  • [5] B. Shang, X. Li, Z. Li, J. Ma, X. Chu, and P. Fan, “Multi-connectivity between terrestrial and non-terrestrial mimo systems,” IEEE Open J. Commun. Soc., vol. 5, pp. 3245–3262, 2024.
  • [6] B. Shang, X. Li, C. Li, and Z. Li, “Coverage in cooperative leo satellite networks,” J. Commun. Inf. Netw., vol. 8, no. 4, pp. 329–340, Dec. 2023.
  • [7] Z. Liu, L. Yin, W. Shin, and B. Clerckx, “Max-min fair energy-efficient beam design for quantized isac leo satellite systems: A rate-splitting approach,” arXiv:2402.09253, 2024.
  • [8] X. Cheng, D. Duan, S. Gao, and L. Yang, “Integrated sensing and communications (isac) for vehicular communication networks (vcn),” IEEE Internet Things J., vol. 9, no. 23, pp. 23 441–23 451, Dec. 2022.
  • [9] J. Zhang, Q. Yan, X. Zhu, and K. Yu, “Smart industrial iot empowered crowd sensing for safety monitoring in coal mine,” Digit. Commun. Netw., vol. 9, no. 2, pp. 296–305, Apr. 2023.
  • [10] D. J. Lary, A. H. Alavi, A. H. Gandomi, and A. L. Walker, “Machine learning in geosciences and remote sensing,” Geosci. Front., vol. 7, no. 1, pp. 3–10, 2016.
  • [11] H. Messer, A. Zinevich, and P. Alpert, “Environmental monitoring by wireless communication networks,” Science, vol. 312, no. 5774, pp. 713–713, 2006.
  • [12] Y. Ma, G. Zhou, and S. Wang, “Wifi sensing with channel state information: A survey,” ACM Comput. Surv., vol. 52, no. 3, pp. 1–36, Jun. 2019.
  • [13] Z. Zhang, W. Chen, Q. Wu, Z. Li, X. Zhu, and J. Yuan, “Intelligent omni surfaces assisted integrated multi-target sensing and multi-user mimo communications,” IEEE Trans. Commun., pp. 1–1, 2024.
  • [14] T. Tian, H. Deng, G. Li, and R. Ma, “Active-passive beamforming with imperfect csi for irs-assisted sensing system,” IEEE Signal Process. Lett., vol. 30, pp. 1052–1056, Aug. 2023.
  • [15] L. You, X. Qiang, C. G. Tsinos, F. Liu, W. Wang, X. Gao, and B. Ottersten, “Beam squint-aware integrated sensing and communications for hybrid massive mimo leo satellite systems,” IEEE J. Sel. Areas Commun., vol. 40, no. 10, pp. 2994–3009, Oct. 2022.
  • [16] M. Huang, F. Gong, G. Li, N. Zhang, and Q.-V. Pham, “Secure precoding for satellite noma-aided integrated sensing and communication,” IEEE Internet Things J., pp. 1–1, 2024.
  • [17] F. Liu, C. Masouros, A. P. Petropulu, H. Griffiths, and L. Hanzo, “Joint radar and communication design: Applications, state-of-the-art, and the road ahead,” IEEE Trans. Commun., vol. 68, no. 6, pp. 3834–3862, Jun. 2020.
  • [18] H. Guo, Y.-C. Liang, J. Chen, and E. G. Larsson, “Weighted sum-rate maximization for reconfigurable intelligent surface aided wireless networks,” IEEE Trans. Wireless Commun., vol. 19, no. 5, pp. 3064–3076, May 2020.
  • [19] D. L. Dampahalage, K. S. Manosha, N. Rajatheva, and M. Latva-Aho, “Weighted-sum-rate maximization for an reconfigurable intelligent surface aided vehicular network,” IEEE Open J. Commun. Soc., vol. 2, pp. 687–703, 2021.
  • [20] K. Shen and W. Yu, “Fractional programming for communication systems—part ii: Uplink scheduling via matching,” IEEE Trans. Signal Process., vol. 66, no. 10, pp. 2631–2644, May 2018.
  • [21] S. P. Boyd and L. Vandenberghe, Convex optimization.   Cambridge, U.K.: Cambridge Univ. Press, 2004.