Preprint
Article

Waveform Design For The Integrated Sensing, Communication and Simultaneous Wireless Information And Power Transfer System

This version is not peer-reviewed.

Submitted:

14 May 2024

Posted:

14 May 2024

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
The next-generation communication systems demand integration of sensing, communication, and power transfer (PT) capabilities, requiring high spectral efficiency, energy efficiency, and low cost, while also necessitating robustness in high-speed scenarios. Integrated sensing and communication systems (ISACS) exhibit the ability to simultaneously perform communication and sensing tasks using a single RF signal, while simultaneous wireless information and power transfer (SWIPT) systems can handle simultaneous information and energy transmission, and orthogonal time frequency space (OTFS) signals are adept at handling high Doppler scenarios. Combining the advantages of these three technologies, a novel cyclic prefix (CP) OTFS based integrated simultaneous wireless sensing, communication, and power transfer system (ISWSCPTS) framework is proposed in this work. Within the ISWSCPTS, the CP-OTFS matched filter (MF) based target detection and parameters estimation (MF-TDaPE) algorithm is proposed to endow the system with sensing capabilities. To enhance the system’s sensing capability, a waveform design algorithm based on CP-OTFS ambiguity function shaping (AFS) is proposed, which is solved by iteratively method. Furthermore, to maximize the system’s sensing performance under communication and PT Quality of Service (QoS) constraints, a semidefinite relaxation (SDR) beamforming design (SDR-BD) algorithm is proposed, solved using through SDR technique. Simulation results demonstrate that the ISWSCPTS exhibits stronger parameter estimation performance in high-speed scenarios compared to orthogonal frequency division multiplexing (OFDM), waveform designed by CP-OTFS AFS demonstrates superior interference resilience, and beamforming designed by SDR-BD strikes a balance in the overall performance of the ISWSCPTS.
Keywords: 
;  ;  ;  ;  ;  ;  

1. Introduction

The next-generation communication systems require the integration of additional functionalities to enhance system spectral efficiency, energy efficiency, and reduce system costs. Integrated sensing and communication systems (ISACS), which converge sensing and communication functionalities within a unified framework, have gained significant traction in recent years [1,2,3,4]. By co-locating sensing and communication tasks and sharing common resources such as antennas and signal processing algorithms, ISACS promise enhanced performance, reduced latency, and improved resource utilization compared to traditional separate implementations. However, the functionalities of ISACS systems still do not fully leverage the potential of wireless radio frequency signals, such as power transfer (PT).
Certain modern wireless applications impose energy transfer requirements on ISACS systems. One such example is battlefield unmanned aerial vehicle (UAV). In combat environments, UAVs operating under constrained resources need to simultaneously accomplish target perception, information transmission, and enhance endurance. One approach to addressing this challenge is the efficient and rational allocation of power in ISACS systems to improve energy efficiency. Some existing works have studied power allocation problems in ISACS systems. Another approach is equipping ISACS with energy harvesting devices to collect energy from received radio frequency (RF) signals, thereby extending the system’s operational lifespan. Regarding the latter, simultaneous wireless information and power transfer (SWIPT), first proposed in [5], has garnered significant attention in recent years [6,7,8]. The core of SWIPT is to use a single signal as a carrier for both information transmission and energy transfer. Communication nodes extract information from the received signal, while energy harvesting nodes collect energy from it. Clearly, inspired by ISACS, SWIPT still has the potential to expand perception capabilities and further improve system efficiency.
In ISACS, the transmitted waveform plays a crucial role in determining system performance. The waveforms based on orthogonal frequency-division multiplexing (OFDM) is a promising choice due to its simple processing framework and high processing gain. However, the performance of OFDM-based waveforms sharply deteriorates in high-speed mobile scenarios [9]. Next-generation communication systems also demand robustness in high-Doppler time-varying scenarios, as required for ISACS. Orthogonal Time Frequency Space (OTFS), a modulation scheme that modulates/demodulates information in the Doppler-delay (DD) domain, has been proposed in recent years [10,11,12]. OTFS represents high-Doppler time-varying channels sparsely in the DD domain. Previous studies have shown that OTFS offers significant communication performance improvement compared to OFDM in high-speed time-varying scenarios [11,12]. OTFS is also utilized in sensing systems and demonstrates excellent performance compared to OFDM [9,13,14]. Due to the aforementioned advantages, cyclic prefix OTFS (CP-OTFS) is adopted as the fundamental waveform in this work.
Currently, there is relatively limited research on the integration of ISACS with SWIPT. In [15], the Cramer-Rao lower bound (CRLB) of the estimation of the targets’ angle is employed as the optimization objective, the beamforming is designed with communication and power transfer quality of service (QoS) constraints in mind. In this work, we propose a framework for cyclic prefix (CP)-OTFS based integrated simultaneous wireless sensing, communication, and power transfer system (ISWSCPTS) that incorporates both search mode (SM) and joint work mode (JWM). This framework integrates the advantages of ISACS and SWIPT, enabling simultaneous perception of targets, information transmission, and energy transfer capabilities. The main contributions of this work are as follows:
  • The framework of ISWSCPTS is proposed. This framework comprises two operational modes: SM and JWM. To enhance sensing capabilities for both SM and JWM, The CP-OTFS matched filter (MF) based target detection and parameters estimation (MF-TDaPE) algorithm is proposed. This algorithm exhibits superior performance to OFDM in high-speed scenarios.
  • An CP-OTFS ambiguity function (AF) shaping (AFS) waveform design algorithm is proposed. Firstly, a novel DD domain AF for CP-OTFS is proposed. Subsequently, aiming to minimize the integrated sidelobe level (ISL) of the proposed AF while adhering to the QoS for communication and PT and constant modulus constraints, a non-convex CP-OTFS waveform design optimization problem is formulated. This problem is then solved through an iterative algorithm to obtain waveforms with superior AF characteristics and interference resilience.
  • The semidefinite relaxation (SDR) beamforming design (SDR-BD) algorithm tailored for ISWSCPTS is proposed. Initially, a non-convex optimization problem is formulated with sensing QoS as the optimization objective and communication and PT QoS as constraints. Subsequently, the SDR technique is employed to solve this problem. The designed waveform optimizes perceptual capability while ensuring communication and PT performance.
The structure of the paper is outlined as follows: Section 2.1 introduces the principles of CP-OTFS. Section 2.2 elaborates on the framework principles, waveform design, and beamforming design. Section 3 presents simulation results and corresponding discussions. Finally, Section 4 concludes with a summary and outlines future work prospects.

2. Materials and Methods

Considering an CP-OTFS based ISWSCPTS equipped with a uniform linear array (ULA) with N t transmitting elements and a single receiving antenna to sense a single target, deliver information to a single communication node (CN), and transfer energy to a single energy receiving node (ERN). As shown in Figure 1. The ISWSCPTS transmits a single information-carrying CP-OTFS signal. Through this multifunctional signal, the ISWSCPTS receives the target echo to complete perception, the CN receives the signal to accomplish information collection, and the ERN receives the signal to complete energy harvesting. The ISWSCPTS is also referred to as sensing node (SN) due to its function of sensing. We details the system model of ISWSCPTS starting from the CP-OTFS model.

2.1. The CP-OTFS Model

2.1.1. Basic Concepts of CP-OTFS

The time–frequency (TF) plane is discretized to a N × M grid as following:
Δ = ( n T , m Δ f ) , n = 0 , . . . , N 1 , m = 0 , . . . , M 1
where T and Δ f = 1 / T are sampling intervals of time and frequency axes, respectively.
The modulated TF samples χ [ n , m ] are transmitted over an OTFS frame with duration T ^ = N T and occupy a bandwidth B = M Δ f . The sampling frequency f s = B .
The DD plane is discretized to an N × M lattice as following:
Ξ = ( k N T , l M Δ f ) , k = 0 , . . . , N 1 , l = 0 , . . . , M 1 .

2.1.2. Input-Output Relationship of CP-OTFS Based ISWSCPTS

   A. The transmitter model
The CP-OTFS based system is depicted in Figure 2. The symbols χ [ k , l ] in DD domain is first transformed into symbols χ ^ [ n , m ] in TF domain by the inverse symplectic finite Fourier transform (ISFFT). Then, the Heisenberg transform is applied to χ ^ to generate time domain discrete signal x ^ [ i ] , i = 0 , . . . , N M 1 . The expression of x ^ [ i ] is
x ^ [ i ] = n N 1 m M 1 k N 1 l M 1 χ [ k , l ] exp { j 2 π ( k n N l m M + m i M ) } g t x [ i n M ] ,
where the g t x [ ι ] is the transmitting rectangular pulse and its expression is
g t x [ ι ] = 1 , 0 ι M 1 0 , o t h e r w i s e
Next, the discrete CP-OTFS transmitting signal x [ i ] is constructed by adding a CP of length M to x [ i ] . Finally, x [ i ] is converted to radio frequency (RF) signal x ( t ) by RF transmitter module.
   B. The channel model
In the application scenarios of ISWSCPTS, the perception target, communication channel, and energy transmission channel can all be modeled using a unified framework based on the DD channel (DDC) model. Assuming there are P DDCs which can be expressed as:
ρ ( τ , f d ) = p = 1 P α p δ ( τ τ p ) δ ( f d f d p ) ,
where α p is the channel gain coefficient, f d p is the Doppler, and τ p is the delay. Let α = [ α 1 , . . . , α p ] T .
In this work, we make the following assumptions:
τ p = l p / M Δ f , l p = 0 , 1 , . . . , M 1 f d p = ( k p ) N / N T , k p = 0 , 1 , . . . , N 1
where l p and k p are integers.
Thus, the discrete form of ρ ( τ , f d ) can be expressed as:
ρ [ k , l ] = p = 1 P α p δ [ k ( k p ) N ] δ [ l l p ]
where
( k p ) N = k p , k p N / 2 N k p , o t h e r w i s e
   C. The receiver model
The received signal r ( t ) is transformed into the discrete-time signal r [ i ] , i = 0 , . . . , N M 1 , by RF receiver module and removing the CP. r [ i ] can be expressed as:
r [ i ] = p = 1 P α p x [ [ i l p ] N M ] exp { j 2 π f d p ( i T M + T ) } + w ^ [ i ]
where [ ] N M denotes modulo N M operation, w ^ [ i ] is the additive Gaussian white noise (AWGN).
r [ i ] is converted to received symbols Υ ^ [ n , m ] in TF domain through the Wigner transform. Υ ^ [ n , m ] can be expressed as:
Υ ^ [ n , m ] = i = 0 N M 1 r [ i ] g r x * [ i n M ] exp ( j 2 π m i M )
where g r x [ ι ] is the receiving rectangular pulse and its expression is
g r x [ ι ] = 1 , 0 ι M 1 0 , o t h e r w i s e
Finally, Υ ^ [ n , m ] is converted into the received symbols Υ [ k , l ] in DD domain through the symplectic finite Fourier transform (SFFT). Υ [ k , l ] can be expressed as:
Υ [ k , l ] = n = 0 N 1 m = 0 M 1 Υ ^ [ n , m ] exp { j 2 π ( k n N l m M ) } + W [ k , l ]
where W is AWGN matrix.
After laborious derivations, (12) can be recast as:
Υ [ k , l ] = k = 0 N 1 l = 0 M 1 χ [ k , l ] H k , l [ k , l ] + W [ k , l ]
where the expression of H k , l [ k , l ] is as follows:
H k , l [ k , l ] = 1 N M 2 ( H I S I + H I C I ) exp ( j 2 π f d T )
In (14), H I S I and H I C I are expressed as follows:
Preprints 106459 i001
Preprints 106459 i002
According to (13), each receiving symbol is the result of summing all input symbols weighted by H k , l [ k , l ] . Clearly, H k , l [ k , l ] is related to variables l, l, k , and l . The magnitude response of H k , l [ k , l ] is as follows:
| H k , l [ k , l ] | = 1 , k = [ k k p ] N , l = [ l l p ] M 0 , o t h e r w i s e
Therefore, the summation operation in (13) can be removed, and it can be simplified as follows:
Υ [ k , l ] = p = 1 P α p χ [ [ k k p ] N , [ l l p ] M ] H k , l [ k p , l p ] + W [ k , l ] ,
where H k , l [ k p , l p ] is defined in (21), the definitions of H 1 and H 2 are (19) and (20), respectively. Note that (19)–(21) are exact expressions without any approximation.
H 1 = N m = 0 M 1 m = 0 M 1 exp { j 2 π l ( m m ) M } i = l p M 1 exp { j 2 π m m + f d p T M }
H 2 = N exp ( j 2 π [ k k p ] N N ) m = 0 M 1 m = 0 M 1 exp { j 2 π l ( m m ) M } i = 0 l p 1 exp { j 2 π m m + f d p T M }
According to (3)–(12), after tedious derivations, Υ [ k , l ] can be expressed as:
H k , l [ k p , l p ] = 1 N M 2 ( H 1 + H 2 ) exp ( j 2 π f d p T )

2.2. Framework of CP-OTFS Based ISWSCPTS

ISWSCPTS encompasses two work modes: SM JWM for perception, communication, and PT. The following provides a detailed description of these two modes.

2.2.1. The SM of ISWSCPTS

In the SM, ISWSCPTS lacks a priori information of DDCs for the target, communication, and power transfer, necessitating estimation of these parameters. In this mode, the ISWSCPTS scans the region of interest in a phased-array manner to acquire states of DDCs from various angles.
The TF-domain received target echo whose angle is the same as the searching direction R S N C N M × 1 can be expressed as:
R S N = P t α a t H ( θ t a r ) η s S M + w S N ,
where P t is the transmitting power, α is the target scattering coefficient (TSC), s S M is the TF-domain transmitting signal in search mode, w S N C N M × 1 is the AWGN matrix with the variance of each entry being σ S N 2 , θ t a r is the angle of the target, a t ( θ t a r ) C N t × 1 is the steering vector and a t ( θ t a r ) can be expressed as:
a t ( θ t a r ) = [ 1 , exp ( j 2 π d sin ( θ t a r ) λ c ) , exp ( j 2 π 2 d sin ( θ t a r ) λ c ) , . . . , exp ( j 2 π ( N t 1 ) d sin ( θ t a r ) λ c ) ] T ,
where d is the antenna spacing, λ c represents the wavelength.
In the phased-array manner, η = a t ( θ t a r ) , thus (22) can be reformulated as:
R S N = P t α N t s S M + w S N ,
According to (18), the vector form of DD-domain echo y S N C N M × 1 can be expressed as:
y S N = P t α N t ( h S N x ^ S N ) + w S N ,
where h S N = vec ( H S N ) , x ^ S N = vec ( X ^ S N ) , vec ( ) vectorizes matrix along the row direction. For a target with delay tap l t a r and Doppler tap k t a r , the ( k , l ) -th entry of H S N can be expressed as:
H S N [ k , l ] = H k , l [ k t a r , l t a r ] .
the ( k , l ) -th entry of X S N can be expressed as:
X ^ S N [ k , l ] = χ [ [ k k t a r ] N , [ l l t a r ] M ] .
The MF algorithm proposed in Reference [9] represents a state-of-the-art approach for OTFS-based target detection and parameter estimation. In this work, we offer a re-engineered version of the MF algorithm, employing an equivalent yet distinct methodology. Note that (25) can be reformulated as:
y S N = P t α N t Λ S N ( k t a r , l t a r ) P S N ( k t a r , l t a r ) x S M + w S N ,
where x S M = vec ( χ ) , Λ S N ( k t a r , l t a r ) = diag ( h S N ) , P S N ( k t a r , l t a r ) is the permutation matrix. The ( i , j ) -th entry of P S N ( k t a r , l t a r ) can be expressed as:
P S N ( k t a r , l t a r ) = 1 , j = [ k k p ] N M + [ l l p ] M 0 , o t h e r w i s e
where i = k M + l . Let Q S N ( k t a r , l t a r ) = Λ S N ( k t a r , l t a r ) P S N ( k t a r , l t a r ) , (28) can be expressed as:
y S N = P t α N t Q S N ( k t a r , l t a r ) x S M + w S N ,
We construct a matrix Ψ C N M × N M . For the i-th row of Ψ where i = k M + l , Ψ ( i , : ) = ( 1 P t N t P S M Q S N ( k , l ) x S M ) H , P S M = x S M H x S M is the power of transmitting signal. Then, the output ξ S N C N M × 1 of proposed MF algorithm can be expressed as:
ξ S N = Ψ y S N .
Peaks will appear in the indices corresponding to k t a r and l t a r in | ξ S N | , hence target presence can be detected through the constant false alarm rate (CFAR) algorithm. Additionally, by utilizing i = k M + l , the delay and Doppler frequency of the target can be estimated. Specifically, if the i ^ -th element of | ξ S N | is a peak, then the corresponding estimated delay tap l ^ t a r , Doppler tap k ^ t a r and TSC α ^ are:
l ^ t a r = [ i ^ ] M ,
l ^ t a r = i ^ M ,
α ^ = ξ S N [ i ^ ] ,
where is the floor operation. Note that Ψ can be precomputed offline. The CP-OTFS MF based target detection and parameters estimation (MF-DaPE) algorithm is summarized in Algorithm 1.
Algorithm 1: The CP-OTFS MF-DaPE algorithm
Input: Ψ , y S N
Output: l ^ t a r , k ^ t a r , α ^
1 Calculate ξ S N through (31)
2 Find peak index in | ξ S N | through CFAR algorithm
3 Calculate l ^ t a r through (32)
4 Calculate k ^ t a r through (33)
5 Calculate α ^ through (34)
6 Return l ^ t a r , k ^ t a r and α ^ .

2.2.2. The JWM of ISWSCPTS

In the JWM, ISWSCPTS accomplishes target tracking, communication, and power transfer. In this work mode, assuming an approximate target location and using the channel state information (CSI) of the CN and ERN as priors is reasonable, as target detection has already been achieved in the SM, and CSI can be obtained by transmitting pilot signals. The ISWSCPTS aims to achieve better performance in the JWM. In this work, we enhance perception performance while ensuring communication and power transfer performance by designing the transmitting CP-OTFS signal and beamforming.
   A. The receiving model of JWM
The DD-domain received echo of the target in TM can be expressed as:
z S N = P t α a t T ( θ t a r ) η Q S N ( k t a r , l t a r ) x J W M + w S N ,
Similarly, the DD-domain received signal in CN can be expressed as follows:
z C N = P t β a t T ( θ C N ) η Q ( k C N , l C N ) x J W M + w C N ,
where β is the complex channel gain, a t ( θ I R ) C N t × 1 is the steering vector corresponding to angle of CN, Q C N is the communication channel response matrix with delay tap l C N and Doppler tap k C N , w C N C N M × 1 is the AWGN matrix with the variance of each entry being σ C N 2 .
The DD-domain received signal in ERN can be expressed as follows:
z E R N = P t γ a t T ( θ E R N ) η Q ( k E R N , l E R N ) x J W M + w E R N ,
where γ is the power transfer channel gain, a t ( θ E R N ) C N t × 1 is the transmitting steering vector, Q E R N is the communication channel response matrix with delay tap l E R N and Doppler tap k E R N , w E R N C N M × 1 is the AWGN matrix with the variance of each entry being σ E R N 2 .
   B. Waveform design by ambiguity function shaping
The ambiguity function (AF) is crucial metrics for assessing radar signals. To enhance the target tracking performance of the system, we formulate an optimization problem aimed at reshaping the AF to achieve lower integral side-lobe levels (ISL). The definition of the traditional radar signal ambiguity function is as follows:
A ( τ , f d ) = x ( t ) x * ( t τ ) exp ( j 2 π f d t ) d t ,
where x ( t ) is the transmitting signal, τ is the delay and f d is the Doppler frequency. However, the output signal of the OTFS system belongs to the delay-Doppler (DD) domain, and the traditional time-domain approach cannot be used to define the OTFS ambiguity function. For OTFS, a discrete AF in the DD domain has been proposed, with the expression as follows:
A k 0 , l 0 [ k , l ] = | x H Q ( k , l ) Q ( k 0 , l 0 ) x | 2 ,
where x is the DD-domain transmitting signal, l 0 and k 0 are delay tap and Doppler tap of interest, respectively. The ISL of A k 0 , l 0 [ k , l ] is:
ISL = k = 0 , k k 0 N 1 l = 0 , l l 0 M 1 | x H Q ( k , l ) Q ( k 0 , l 0 ) x | 2 ,
Thus, the problem of waveform design for AF shaping is expressed as:
P 0 : min x J W M k = 0 , k k 0 N 1 l = 0 , l l 0 M 1 | x J W M H Q ( k , l ) Q ( k 0 , l 0 ) x J W M | 2 s . t . | x J W M ( i ) | = 1 , i = 1 , 2 , . . . , N M
The constraint of (41) is constant modulus constraint, which is preferred by radar system. Based on the proposition 1 presented in [16], P 0 can be handled by sequentially solving the following approximation problem:
P t + 1 : max x J W M u ( x J W M , x J W M ( t ) ) s . t . | x J W M ( i ) | = 1 , i = 1 , 2 , . . . , N M
where x J W M ( t ) denotes the t-th iteration solution of the proposed iteration algorithm. The expression of u ( x J W M , x J W M ( t ) ) is:
u ( x J W M , x J W M ( t ) ) = x J W M H ( λ ( t ) I N M Θ ( x J W M ( t ) ) ) x J W M ,
where λ ( t ) satisfies that λ ( t ) β λ max ( Θ ( x J W M ( t ) ) ) . Let Γ k 0 , l 0 ( k , l ) = Q ( k , l ) Q ( k 0 , l 0 ) , Θ ( x ) is expressed as:
Θ ( x ) = 1 2 ( N M 1 ) k = 0 , k k 0 N 1 l = 0 , l l 0 M 1 ( Γ k 0 , l 0 ( k , l ) H x x H Γ k 0 , l 0 ( k , l ) + Γ k 0 , l 0 ( k , l ) x x H Γ k 0 , l 0 ( k , l ) H ) .
Furthermore, according to proposition 2 in [16], P t + 1 can be solved through the following problem:
P 1 : max x ^ J W M Re ( x ^ J W M H v ( t ) ) s . t . | x ^ J W M ( i ) | = 1 , i = 1 , 2 , . . . , N M
where v ( t ) can be expressed as:
v ( t ) = ( λ ( t ) I N M Θ ( x J W M ( t ) ) ) x J W M ( t ) .
Then, the closed-form solution for P 1 is:
x ^ J W M = exp ( j arg ( v ( t ) ) ) ,
where exp ( ) and arg ( ) are applied element-wise to the vectors.
Note that the optimal waveform x J W M can also be utilized in SM. This waveform design algorithm is called CP-OTFS AFS, and it is summarized in Algorithm 2.
Algorithm 2: CP-OTFS AFS
Input: Initial x J W M ( 0 ) and stop condition ϵ
Output: The optimized waveform x J W M
1 Let t 0 and x J W M ( t ) x J W M ( 0 )
2 v ( t ) ( λ ( t ) I N M Θ ( x J W M ( t ) ) ) x J W M ( t )
3 x ^ J W M exp ( j arg ( v ( t ) ) )
4 t t + 1
5 x J W M ( t ) x ^ J W M
6 If | x J W M ( t ) x J W M ( t 1 ) | ϵ , return x J W M = x J W M ( t ) . Otherwise, return to step 2.
C. Beamforming Design for ISWSCPTS
In this section, we investigate the beamforming design for ISWSCPTS. Within the beamforming design challenges of ISWSCPTS, two key considerations emerge: 1) optimization of sensing performance; 2) fulfilling the requirements for power transfer in the context of ERN considerations.
We begin by deriving metrics for sensing, communication, and power transfer. Subsequently, we formulate the optimization problem for ISWSCPTS beamforming design. Finally, we propose an algorithm to solve the beamforming design optimization problem.
C1. Sensing metric
In the JWM, the ISWSCPTS system necessitates continuous estimation of target parameters, and the precision of target parameter estimation is closely tied to the signal-to-noise ratio (SNR). Consequently, we employ SNR as the sensing metric.
According to (25), the SNR SN is
SNR SN = E { | | P t α a t H ( θ t a r ) η ( h S N x S N ) | | F 2 } E { | | w S N | | F 2 } .
Based on the properties of | h S N ( j ) | 2 = 1 and | x S N ( j ) | 2 = 1 , through a series of calculations, (48) can be simplified to:
SNR SN = P t | α | 2 η H a r ( θ t a r ) a r H ( θ t a r ) η σ S N 2 = η T G SN η .
where G SN = P t | α | 2 σ S N 2 a r ( θ t a r ) a r H ( θ t a r ) .
C2. Communication metric
Due to the impact of communication quality, such as bit error rate (BER) and channel capacity, being closely related to SNR, and in order to maintain consistency with snesing metrics, SNR is also employed as the communication metric. Similar to derivation of sensing metric, the SNR CN is
SNR CN = E { | | P t β a t H ( θ C N ) η ( h C N x C N ) T | | F 2 } E { | | w C N | | F 2 } = P t | β | 2 η H a t ( θ C N ) a t H ( θ C N ) η σ C N 2 = η H G CN η ,
where G CN = P t | β | 2 σ C N 2 a t ( θ C N ) a t H ( θ C N ) .
C3. Power transfer metric
The ERN collects energy from signals emitted by the SN, thus the harvested energy is employed as the metric of Power transfer. Due to the negligible power of noise compared with transmitting signal, the harvested energy can be expressed as:
E ERN = E { μ | | P t γ a t H ( θ E R N ) η ( h E R N x E R N ) T | | F 2 } = μ P t | γ | 2 η H a r ( θ E R N ) a r H ( θ E R N ) η = η H G ERN η ,
where μ [ 0 , 1 ] is the energy harvesting efficiency, G CN = μ P t | γ | 2 a r ( θ E R N ) a r H ( θ E R N ) .
C4. SDR-BD algorithm
As we want to optimize the sensing performance of ISWSCPTS while meeting the basic requirements for communication and power transfer, the problem of beamforming design can be expressed as:
P 2 : max η η H G SN η s . t . η H G CN η ρ η H G ERN η ϱ η H η 1
where ρ represents the minimum SNR required for communication, ϱ denotes the minimum energy requirement for triggering the energy harvesting process in ERN.
Due to the max operation and quadratic constraints, problem P 2 is non-convex optimization. In this work, P 2 is solved through SDR technique [17]. According to x H A x = Tr ( A x x H ) , where Tr ( ) denotes the trace operation, P 2 can be transformed into the following equivalent optimization problem:
P 3 : max Γ Tr ( G SN Γ ) s . t . Tr ( G CN Γ ) ρ Tr ( G ERN Γ ) ϱ Tr ( Γ ) 1 rank ( Γ ) = 1 Γ ( 0 )
where Γ = η η H .
Then, the rank constraint is dropped to obtain the following relaxed optimization problem:
P 4 : max Γ Tr ( G SN Γ ) s . t . Tr ( G CN Γ ) ρ Tr ( G ERN Γ ) ϱ Tr ( Γ ) 1 Γ ( 0 )
The objective function and constraints in P 4 are all affine, thus making P 4 a convex optimization problem. P 4 can be solved by Matlab using the convex optimization toolbox CVX[]. If the rank of the solution Γ obtained in P 4 is not equal to 1, then feasible solutions η need to be extracted from Γ . Following [1], η can be obtained as follows:
η = λ 0 q 0
where λ 0 represents the maximum eigenvalue of Γ , and q 0 denotes the corresponding eigenvector.

3. Results and Discussions

This section conducts simulations to evaluate the performance of CP-OTFS MF-TDaPE , CP-OTFS AFS, and SDR-BD. The basic simulation parameters are summarized in Table 1.

3.1. Simulation of CP-OTFS MF-TDaPE

To begin with, the target detection is performed using CP-OTFS MF-TDaPE algorithm. The target range and velocity are 689.523 m and 95.054 m/s, corresponding to l p = 46 and k p = 5 , respectively, with a TSC of 0.8247 + 0.4709 i . The SNR of this simulation is 10 dB. Figure 3 illustrates the output ξ S N of the MF, where a significant peak appears at positions l = 46 and k = 5 . The target can be correctly detected using the CFAR algorithm.
To validate the parameter estimation performance of the CP-OTFS MF-TDaPE, we conduct parameter estimation for fast-moving targets with different velocities. The benchmark is OFDM waveform, and the estimation method is the fast Fourier transform (FFT) method described in [18]. Since both CP-OTFS and OFDM can accurately estimate the range, the main comparison focuses on velocity estimation performance and TSC estimation performance. The root mean square error (RMSE) is adopted as the performance metric, calculated as follows:
R M S E = ( n = 1 N s | y ^ y | 2 N s ) 1 2
where y ^ is the estimated parameter, y is the true value, N s is the number of samples.
Figure 4 and Figure 5 depict the RMSE results for velocity and TSC estimation of fast-moving targets with different velocities, respectively. It is evident that for fast-moving targets, the performance of velocity estimation sharply deteriorates with increasing velocity in OFDM, while OTFS consistently maintains accurate estimation. Regarding TSC estimation, although OFDM does not exhibit sensitivity to velocity variations, its estimation results still pale in comparison to the nearly error-free estimates provided by CP-OTFS. The results presented demonstrate the robustness of the ISWSCPTS framework in parameter estimation for high-speed motion scenarios.

3.2. Simulation of CP-OTFS AFS

The initial waveform symbols for this simulation are quadrature phase shift keying (QPSK) symbols. We first compare the AF of the waveform optimized by CP-OTFS AFS with the initial waveform. Figure 6 and Figure 7 respectively depict the AF of the initial and optimized waveforms. It can be observed that the optimized waveform exhibits a reduction of approximately 2-3 dB in side lobes in comparison to the initial waveform. Figure 8 illustrates the variation of the ISL with the number of iterations, showing a gradual decrease with successive iterations of the algorithm. However, the convergence speed of CP-OTFS AFS is slow, and there is still potential for improvement.
In fact, the performance of CP-OTFS MF-TDaPE in estimating TSC deteriorates with the emergence of interference. To illustrate this point, we introduced an interference target in the simulation scenario with delay and Doppler taps of 30 and 3, respectively, and the TSC of 1 for this target. Figure 9 shows the variation of the target TSC with SNR in the absence of interference, indicating similar performance between the initial QPSK waveform and the optimized waveform. Figure 10 illustrates the TSC estimation performance in the presence of interference. It is evident that the performance degradation appears. The reason is that the interference target affects TSC estimation through AF sidelobes. At the same time, it can also be observed that starting from an SNR greater than 6 dB, the optimized waveform exhibits stronger resistance to interference.

3.3. Simulation of SDR-BD

This simulation showcases the beam patterns of designed by SDR-BD under different system configurations and QoS settings. The specific simulation parameters for this study are detailed in Table 2.
Figure 11 illustrates the beamforming design results of SDR-BD under different configurations of N t . The simulation result demonstrates that as N t increases, the beamforming becomes narrower in the directions of the target, CN, and ERN, indicating stronger directivity. Additionally, the beamforming exhibits lower sidelobes. Particularly, the highest sidelobe level of the beamforming near the target with N t = 64 is approximately 6.5 dB lower than that with N t = 16 . This phenomenon arises from the increase in system antenna aperture with N t , resulting in improved angular resolution of the system.
Figure 12 compares the beamforming with N t = 16 designed by SDR-BD under different communication and PT QoS settings. The simulation result indicates that increasing the QoS for communication and PT inevitably leads to a deterioration in sensing capability. Specifically, compared to the case where ρ = 15 d B and ϱ = 5 m W , an increase of 5 dB in ρ and 5 mW in ϱ results in a reduction of approximately 4.8 dB in the antenna gain in the direction of the target. This outcome suggests that the beamforming designed by SDR-BD achieves a trade-off among sensing, communication, and PT in ISWSCPTS. Therefore, in practical application scenarios, the overall system performance optimization can be achieved according to respective QoS requirements.

4. Conclusions

In this work, the ISWSCPTS framework, which integrates ISACS, SWIPT, and OTFS, is proposed. ISWSCPTS comprises two work modes, SM and JWM, which utilize a single CP-OTFS waveform. SM is employed for target detection and parameter estimation, while JWM is used for target tracking, communication, and PT. For target detection and parameters estimation, CP-OTFS MF-DaPE algorithm is proposed, demonstrating superior performance in high-speed scenarios compared to OFDM. To enhance the robustness of parameters estimation against interference, CP-OTFS AFS algorithm is employed for waveform design. The waveform designed by CP-OTFS AFS reduces the ISL of corresponding AF and improves interference resilience. The SDR-BD algorithm is proposed to enhance the overall perceptual capability of the system by designing beamforming to maximize perceptual capability under communication and PT QoS constraints.
However, this work still has the following limitations: (1) The ISWSCPTS scenarios include only single target, CN, and ERN; (2) The CP-OTFS AFS incurs significant computational burden. In future work, we will focus on extending ISWSCPTS to multiple targets, CNs, and ERNs scenarios, and optimizing CP-OTFS AFS to improve its convergence speed.

Author Contributions

Conceptualization, Qilong Miao; methodology, Qilong Miao; formal analysis, Qilong Miao; investigation, simulation and analysis, Qilong Miao; writing original draft preparation, Qilong Miao and Chenfei Xie; writing,review and editing, Weimin Shi, Yong Gao and Lu Chen. Allauthors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article material, further inquiries can be directed to the corresponding author.

Acknowledgments

This work was supported by Sichuan Science and Technology Program under Grant 2023NSFSC0450, and the 111 Project under Grant B17008.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PT Power Transfer
ISACS Integrated Sensing and Communication System
SWIPT Simultaneous Wireless Information and Power Transfer
OTFS Orthogonal Time Frequency Space
RF Radio Frequency
CP Cyclic Prefix
ISWSCPTS Integrated Simultaneous Wireless Sensing, Communication, and Power Transfer System
MF Matched Filter
AF Ambiguity Function
AFS AF Shaping
MF-TDaPE MF based Target Detection and Parameters Estimation
QoS Quality of Service
SDR Semidefinite Relaxation
SDR-BD SDR Beamforming Design
OFDM Orthogonal Frequency Division Multiplexing
UAV Unmanned Aerial Vehicle
SM Search Mode
JWM Joint Work Mode
ISL Integrated Sidelobe Level
ULA Uniform Linear Array
CN Communication Node
ERN Energy Receiving Node
SN Sensing Node
TF Time-Frequency
DD Doppler-Delay
DDC DD Channel
SFFT Symplectic Finite Fourier Transform
ISFFT Inverse SFFT
TSC Target Scattering Coefficient
CFAR Constant False Alarm Rate
DDC DD Channel
CSI Channel State Information
FFT Fast Fourier Transform
RMSE Root Mean Square Error
SNR Signal to Noise Ratio
QPSK Quadrature Phase Shift Keying

References

  1. Hassanien, A.; Amin, M.G.; Aboutanios, E.; Himed, B. Dual-Function Radar Communication Systems: A solution to the spectrum congestion problem. IEEE Signal Processing Magazine 2019, 36, 115–126. [Google Scholar] [CrossRef]
  2. Mishra, K.V.; Bhavani Shankar, M.; Koivunen, V.; Ottersten, B.; Vorobyov, S.A. Toward Millimeter-Wave Joint Radar Communications: A Signal Processing Perspective. IEEE Signal Processing Magazine 2019, 36, 100–114. [Google Scholar] [CrossRef]
  3. Liu, F.; Masouros, C.; Petropulu, A.P.; Griffiths, H.; Hanzo, L. Joint Radar and Communication Design: Applications, State-of-the-Art, and the Road Ahead. IEEE Transactions on Communications 2020, 68, 3834–3862. [Google Scholar] [CrossRef]
  4. Liu, F.; Cui, Y.; Masouros, C.; Xu, J.; Han, T.X.; Eldar, Y.C.; Buzzi, S. Integrated Sensing and Communications: Toward Dual-Functional Wireless Networks for 6G and Beyond. IEEE Journal on Selected Areas in Communications 2022, 40, 1728–1767. [Google Scholar] [CrossRef]
  5. Varshney, L.R. Transporting information and energy simultaneously. 2008 IEEE International Symposium on Information Theory, 2008, pp. 1612–1616. [CrossRef]
  6. Amjad, M.; Chughtai, O.; Naeem, M.; Ejaz, W. SWIPT-Assisted Energy Efficiency Optimization in 5G/B5G Cooperative IoT Network. Energies 2021, 14. [Google Scholar] [CrossRef]
  7. Andrawes, A.; Nordin, R.; Abdullah, N.F. Energy-Efficient Downlink for Non-Orthogonal Multiple Access with SWIPT under Constrained Throughput. Energies 2020, 13. [Google Scholar] [CrossRef]
  8. Choi, H.H.; Lee, J.R. Energy-Neutral Operation Based on Simultaneous Wireless Information and Power Transfer for Wireless Powered Sensor Networks. Energies 2019, 12. [Google Scholar] [CrossRef]
  9. Raviteja, P.; Phan, K.T.; Hong, Y.; Viterbo, E. Orthogonal Time Frequency Space (OTFS) Modulation Based Radar System. 2019 IEEE Radar Conference (RadarConf), 2019, pp. 1–6. [CrossRef]
  10. Wei, Z.; Yuan, W.; Li, S.; Yuan, J.; Bharatula, G.; Hadani, R.; Hanzo, L. Orthogonal Time-Frequency Space Modulation: A Promising Next-Generation Waveform. IEEE Wireless Communications 2021, 28, 136–144. [Google Scholar] [CrossRef]
  11. Hadani, R.; Rakib, S.; Tsatsanis, M.; Monk, A.; Goldsmith, A.J.; Molisch, A.F.; Calderbank, R. Orthogonal Time Frequency Space Modulation. 2017 IEEE Wireless Communications and Networking Conference (WCNC), 2017, pp. 1–6. [CrossRef]
  12. Raviteja, P.; Phan, K.T.; Hong, Y.; Viterbo, E. Interference Cancellation and Iterative Detection for Orthogonal Time Frequency Space Modulation. IEEE Transactions on Wireless Communications 2018, 17, 6501–6515. [Google Scholar] [CrossRef]
  13. Gaudio, L.; Kobayashi, M.; Caire, G.; Colavolpe, G. On the Effectiveness of OTFS for Joint Radar Parameter Estimation and Communication. IEEE Transactions on Wireless Communications 2020, 19, 5951–5965. [Google Scholar] [CrossRef]
  14. Gaudio, L.; Kobayashi, M.; Bissinger, B.; Caire, G. Performance Analysis of Joint Radar and Communication using OFDM and OTFS. 2019 IEEE International Conference on Communications Workshops (ICC Workshops), 2019, pp. 1–6. [CrossRef]
  15. Li, X.; Yi, X.; Zhou, Z.; Han, K.; Han, Z.; Gong, Y. Multi-User Beamforming Design for Integrating Sensing, Communications, and Power Transfer. 2023 IEEE Wireless Communications and Networking Conference (WCNC), 2023, pp. 1–6. [CrossRef]
  16. Yang, J.; Cui, G.; Yu, X.; Xiao, Y.; Kong, L. Cognitive Local Ambiguity Function Shaping With Spectral Coexistence. IEEE Access 2018, 6, 50077–50086. [Google Scholar] [CrossRef]
  17. Luo, Z.q.; Ma, W.k.; So, A.M.c.; Ye, Y.; Zhang, S. Semidefinite Relaxation of Quadratic Optimization Problems. IEEE Signal Processing Magazine 2010, 27, 20–34. [Google Scholar] [CrossRef]
  18. Sturm, C.; Wiesbeck, W. Waveform Design and Signal Processing Aspects for Fusion of Wireless Communications and Radar Sensing. Proceedings of the IEEE 2011, 99, 1236–1259. [Google Scholar] [CrossRef]
Figure 1. The structure of ISWSCPTS.
Figure 1. The structure of ISWSCPTS.
Preprints 106459 g001
Figure 2. The structure of CP-OTFS.
Figure 2. The structure of CP-OTFS.
Preprints 106459 g002
Figure 3. The output of MF.
Figure 3. The output of MF.
Preprints 106459 g003
Figure 4. Velocity RMSE vs relative velocity.
Figure 4. Velocity RMSE vs relative velocity.
Preprints 106459 g004
Figure 5. TSC RMSE vs relative velocity.
Figure 5. TSC RMSE vs relative velocity.
Preprints 106459 g005
Figure 6. AF of initial waveform.
Figure 6. AF of initial waveform.
Preprints 106459 g006
Figure 7. AF of optimized waveform.
Figure 7. AF of optimized waveform.
Preprints 106459 g007
Figure 8. ISL vs iterations.
Figure 8. ISL vs iterations.
Preprints 106459 g008
Figure 9. TSC RMSE vs SNR under interference-free conditions.
Figure 9. TSC RMSE vs SNR under interference-free conditions.
Preprints 106459 g009
Figure 10. TSC RMSE vs SNR under interference conditions.
Figure 10. TSC RMSE vs SNR under interference conditions.
Preprints 106459 g010
Figure 11. Designed beamforming of different N t .
Figure 11. Designed beamforming of different N t .
Preprints 106459 g011
Figure 12. Designed beamforming of different N t .
Figure 12. Designed beamforming of different N t .
Preprints 106459 g012
Table 1. Basic simulation parameters.
Table 1. Basic simulation parameters.
Symbol Parameter Value
f c Carrier frequency 77 GHz
N Number of Doppler samples 16
N Number of delay samples 64
B Total bandwidth 10 MHz
Δ f Subcarrier spacing 156.250 kHz
Δ R Range resolution 14.99 m
Δ V Velocity resolution 19.01 m/s
R m a x Unambiguous range 959.336 m
V m a x Unambiguous velocity ±152.086 m/s
Table 2. Specific simulation parameters for SDR-BD.
Table 2. Specific simulation parameters for SDR-BD.
Symbol Parameter Value
P t Transmitting power 1 W
σ C N 2 Variance of noise in CN 0 dBm
σ S N 2 Variance of noise in SN 10 dBm
α TSC 0.8247+0.4709i
β Communication channel gain 0.01
γ PT channel gain 0.02 m
μ Energy harvesting efficiency 0.1
ρ SNR required in CN Configured
ϱ Energy requirement in SN Configured
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

187

Views

68

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated