Waveform Design for the Integrated Sensing, Communication, and Simultaneous Wireless Information and Power Transfer System

Next-generation communication systems demand the 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 (ISACSs) 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 parameter 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 an iterative 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, which is solved using through the SDR technique. The simulation results demonstrate that the ISWSCPTS exhibits stronger parameter estimation performance in high-speed scenarios compared to orthogonal frequency division multiplexing (OFDM), the waveform designed by CP-OTFS AFS demonstrates superior interference resilience, and the beamforming designed by SDR-BD strikes a balance in the overall performance of the ISWSCPTS.


Introduction
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 (ISACSs), 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 ISACSs 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 ISACSs.One such example is a battlefield unmanned aerial vehicle (UAV).In combat environments, UAVs operating under constrained resources need to simultaneously accomplish target perception, transmit information, and enhance endurance.One approach to addressing this challenge is the efficient and rational allocation of power in ISACSs to improve energy efficiency.Some existing works have studied power allocation problems in ISACSs.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 energyharvesting 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.Waveforms based on orthogonal frequency division multiplexing (OFDM), extensively used in fourth-generation (4G) and emerging fifth-generation (5G) cellular systems as well as WiFi networks, are a promising choice due to their simple processing framework and high processing gain.However, OFDM encounters challenges in highmobility scenarios.The OFDM waveform suffers from severe inter-carrier interference (ICI), exacerbated by varying normalized Doppler shifts between the highest and lowest subcarriers.Consequently, achieving synchronization becomes challenging.Similar to communication systems, 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 is 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 that are sparse in the DD domain.Previous studies have shown that OTFS offers significant communication performance improvements 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, the 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, and the beamforming is designed with communication and power transfer quality of service (QoS) constraints in mind.We analyze the shortcomings of existing technologies, as depicted in Table 1, to clarify the motivation, innovation, and significance of this work.

ISACS Lacking expansion on PT. OFDM
Introducing severe ICI In high-speed scenarios.OTFS-based radar in [9] Inaccuracy in estimating TSC.
System in [15] Lacking consideration for high-speed scenarios.
In this work, we propose a framework for a cyclic prefix (CP)-OTFS-based integrated simultaneous wireless sensing, communication, and power transfer system (ISWSCPTS), which incorporates both a search mode (SM) and joint work mode (JWM).This framework integrates the advantages of ISACS and SWIPT, enabling the 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 parameter estimation (MF-TDaPE) algorithm is proposed.This algorithm exhibits superior performance to OFDM in high-speed scenarios.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.

Materials and Methods
Consider 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 informationcarrying CP-OTFS signal.Through this multifunctional signal, the ISWSCPTS receives a target echo to complete the perception, the CN receives the signal to accomplish the information collection, and the ERN receives the signal to complete the energy-harvesting.The ISWSCPTS is also referred to as a sensing node (SN) due to its function of sensing.We detail the system model of ISWSCPTS starting from the CP-OTFS model.The time-frequency (TF) plane is discretized to a N × M grid, as follows: 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 = NT and occupy a bandwidth B = M∆ f .The sampling frequency f s = B.
The DD plane is discretized to an N × M lattice as follows [12]: where the g tx [ι] is the transmitting rectangular pulse and its expression is

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.The traditional TD channel model assumes that the channel impulse response (CIR) is time-invariant.However, with increasing mobility and carrier frequency, this time-invariant assumption may no longer hold.Therefore, new channel models need to consider the CIR in both the TD and DD dimensions.As pointed out in [10], DDC exhibits beneficial features such as separability, stability, compactness, and possible sparsity.Hence, this paper adopts DDC.
Assume there are P DDCs which can be expressed as [9] ρ(τ, where α p is the channel gain coefficient, f d p is the Doppler, and τ p is the delay.Let In this work, we make the following assumptions: where l p and k p are integers.Thus, the discrete form of ρ(τ, f d ) can be expressed as where C. The receiver model The received signal r(t) is transformed into the discrete-time signal r[i], i = 0, ..., N M − 1, by the RF receiver module and by removing the CP.r[i] can be expressed as where [•] N M denotes modulo N M operation, and ŵ[i] is the additive Gaussian white noise (AWGN).r[i] is converted to received symbols Υ[n, m] in the TF domain through the Wigner transform.Υ[n, m] can be expressed as where g rx [ι] is the receiving rectangular pulse, and its expression is Finally, Υ[n, m] is converted into the received symbols Υ[k, l] in the DD domain through the symplectic finite Fourier transform (SFFT).Υ[k, l] can be expressed as where W is the AWGN matrix.
After laborious derivations, ( 12) can be recast as where the expression of H k,l [k ′ , l ′ ] is as follows: In ( 14), H ISI and H ICI are expressed as follows: According to (13), each receiving symbol is the result of summing all input symbols weighted by Therefore, the summation operation in (13) can be removed, and it can be simplified as follows: 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.
According to (3)-( 12), after tedious derivations, Υ[k, l] can be expressed as ISWSCPTS encompasses two work modes: the SM and JWM for perception, communication, and PT.The following provides a detailed description of these two modes.

The SM of ISWSCPTS
In the SM, the ISWSCPTS lacks a priori information on DDCs for the target, communication, and power transfer, necessitating the estimation of these parameters.In this mode, the ISWSCPTS scans the region of interest in a phased-array manner to acquire the states of DDCs from various angles.
The TF domain-received target echo whose angle is the same as the searching direction R SN ∈ C N M×1 can be expressed as where P t is the transmitting power, α is the target scattering coefficient (TSC), s SM is the TF domain transmitting signal in the search mode, w SN ∈ C N M×1 is the AWGN matrix with the variance of each entry being σ 2 SN , θ tar is the angle of the target, a t (θ tar ) ∈ C N t ×1 is the steering vector, and a t (θ tar ) can be expressed as where d is the antenna spacing, λ c represents the wavelength.
In the phased-array manner, η = a t (θ tar ); thus, (22) can be reformulated as According to (18), the vector form of the DD domain echo y SN ∈ C N M×1 can be expressed as where h SN = vec(H SN ), xSN = vec( XSN ), vec(•) vectorizes the matrix along the row direction.For a target with delay tap l tar and Doppler tap k tar , the (k, l)-th entry of H SN can be expressed as The (k, l)-th entry of X SN can be expressed as 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 reengineered version of the MF algorithm, employing an equivalent yet distinct methodology.Note that (25) can be reformulated as where x SM = vec(χ), Λ SN (k tar , l tar ) = diag(h SN ), P SN (k tar , l tar ) is the permutation matrix.The (i, j)-th entry of P SN (k tar , l tar ) can be expressed as where i = kM + l.Let Q SN (k tar , l tar ) = Λ SN (k tar , l tar )P SN (k tar , l tar ), (28) can be expressed as We construct a matrix Ψ ∈ C N M×N M .For the i-th row of Ψ where i = kM + l, Ψ(i, :) = ( is the power of transmitting signal.Then, the output ξ SN ∈ C N M×1 of the proposed MF algorithm can be expressed as Peaks will appear in the indices corresponding to k tar and l tar in |ξ SN |; hence, target presence can be detected through the constant false alarm rate (CFAR) algorithm.Additionally, by utilizing i = kM + l, the delay and Doppler frequency of the target can be estimated.Specifically, if the î-th element of |ξ SN | is a peak, then the corresponding estimated delay tap ltar , Doppler tap ktar , and TSC α are where ⌊•⌋ is the floor operation.Note that Ψ can be precomputed offline.

The JWM of ISWSCPTS
In the JWM, the 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 Similarly, the DD domain-received signal in CN can be expressed as follows: where β is the complex channel gain, a t (θ IR ) ∈ C N t ×1 is the steering vector corresponding to angle of CN, Q CN is the communication channel response matrix with delay tap l CN and Doppler tap k CN , w CN ∈ C N M×1 is the AWGN matrix with the variance of each entry being σ 2 CN .The DD domain-received signal in ERN can be expressed as follows: where γ is the power transfer channel gain, a t (θ ERN ) ∈ C N t ×1 is the transmitting steering vector, Q ERN is the communication channel response matrix with delay tap l ERN and Doppler tap k ERN , w ERN ∈ C N M×1 is the AWGN matrix with the variance of each entry being σ 2 ERN .

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 (ISLs).The definition of the traditional radar signal ambiguity function is as follows [16]: 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: where x is the DD domain transmitting signal, and l 0 and k 0 are the delay tap and Doppler tap of interest, respectively.The ISL of Thus, the problem of waveform design for AF shaping is expressed as The optimization objective in (41) represents the sum of the energy of all mismatched delay-Doppler pairs in the AF.The goal of minimizing this objective is to make the output of the MF for detection as small as possible to reduce false alarms while minimizing interference introduced by uninterested delay and Doppler regions.The constraint of (41) is a constant modulus constraint, which is preferred by radar system.Based on proposition In each iteration of CP-OTFS AFS, the required storage space includes variables x (t) JW M , x H JW M , Θ, Γ k 0 ,l 0 (k, l), Γ k 0 ,l 0 (k, l) H , I N M , and v (t) .Therefore, the space complexity of CP-OTFS AFS is O((2(N − 1)(M − 1) + 2)(N M) 2 + 3N M).

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) the 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 the ISWSCPTS beamforming design.Finally, we propose an algorithm to solve the beamforming design optimization problem.

C1. Sensing metric
In the JWM, the ISWSCPTS necessitates the 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 the SNR as the sensing metric.

C2. Communication metric
Due to the impact of communication quality, such as the bit error rate (BER) and channel capacity, being closely related to SNR and in order to maintain consistency with sensing metrics, the SNR is also employed as the communication metric.Similar to derivation of sensing metric, the SNR CN is where

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 where µ ∈ [0, 1] is the energy-harvesting efficiency, G CN = µP t |γ| 2 a r (θ ERN )a H r (θ ERN ).

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 the beamforming design can be expressed as where ρ represents the minimum SNR required for communication, and ϱ denotes the minimum energy requirement for triggering the energy-harvesting process in ERN.
The optimization objective of P 2 is the received SNR of the SN.By optimizing the SNR of the SN, we aim to maximize sensing performance under the QoS constraints of communication and PT, since the detection and parameter estimation performance of the SN improves with an increasing SNR.This optimization objective is related to the target position and η.Since the JWM of the ISWSCPTS has target position information, η can be optimized to maximize the SNR.
Due to the max operation and quadratic constraints, problem P 2 is a non-convex optimization.In this work, P 2 is solved through the SDR technique [18].According to x H Ax = Tr(Axx H ), where Tr(•) denotes the trace operation, P 2 can be transformed into the following equivalent optimization problem: where Γ = ηη H .Then, the rank constraint is dropped to obtain the following relaxed optimization problem: 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 [18], η ⋆ can be obtained as follows: where λ 0 represents the maximum eigenvalue of Γ ⋆ , and q 0 denotes the corresponding eigenvector.Note that the SDR-BD algorithm is a quadratic programming problem.According to [19], the time complexity of quadratic programming problems is O(N 3 ), and the space complexity is O(N 2 ).Therefore, the time complexity of SDR-BD is O((N M) 3 ), and the space complexity is O((N M) 2 ).

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 2.

Simulation of CP-OTFS MF-TDaPE
To begin with, target detection is performed using the CP-OTFS MF-TDaPE algorithm.The target range and velocity are 689.523m and 95.054 m/s, corresponding to l p = 46 and k p = 5, respectively, with a TSC of 0.8247 + 0.4709i.The SNR of this simulation is 10 dB. Figure 3 illustrates the output ξ SN 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 the OFDM waveform, and the estimation method is the fast Fourier transform (FFT) method described in [20].Since both CP-OTFS and OFDM can accurately estimate the range, the main comparison focuses on the velocity estimation performance and TSC estimation performance.The root mean square error (RMSE) is adopted as the performance metric and is calculated as follows: (56 where ŷ is the estimated parameter, y is the true value, and N s is the number of samples.Figures 4 and 5 depict the RMSE results for the velocity and TSC estimation of fastmoving targets with different velocities, respectively.In the experiments for Figures 4 and 5 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.
Interestingly, in Figure 4, when the target's relative velocity is low, such as 38.02 m/s, OFDM performs better at a low SNR.Conversely, at high velocities, OTFS outperforms OFDM.This phenomenon can be explained by OFDM's weak anti-Doppler capability, causing a "misalignment" of the peak used for parameter estimation.At low velocities, this "misalignment" is small.At a low SNR, the noise can "correct" this misalignment.However, at high velocities, the "misalignment" becomes significant, and the noise's "correcting" effect becomes negligible.The phenomenon observed in Figure 5 is also interesting.For OFDM, at low velocities, such as 38.02 m/s, the TSC estimation performance is intuitive.However, at high velocities, such as 114.06 m/s, the results at a high SNR are unexpectedly worse.This is because, at low velocities, the primary source of TSC estimation error is noise, whereas at high velocities, it is Doppler spread.Additionally, the previously mentioned "correcting" effect of noise at a low SNR also plays a role.The explanation for the near-perfect CP-OTFS results in Figures 4 and 5 is as follows: In CP-OTFS systems, the DDC model induces sparsity in the channel and a low correlation between time and frequency domains, resulting in discrete channel impulses in the DD domain.Due to these characteristics, CP-OTFS systems can achieve perfect Doppler estimation performance.In contrast, OFDM experiences strong ICI in high-speed scenarios, leading to significant performance degradation.

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.Figures 6 and 7, respectively, depict the AF of the initial and optimized waveforms.In the experiments for Figures 6 and 7, the initial waveform is set to a random QPSK signal, with the interested delay tap and Doppler tap set to 46 and 5, respectively.Using the initial and optimized signals, the corresponding AFs are obtained according to (39).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 in 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 interference targets in the simulation scenario.In the experiment for Figure 9, there are no interference targets, the target's delay tap and Doppler tap are set to 46 and 5, respectively, and the target's TSC is set to 0.8247 + 0.4709i.The MF-TDaPE algorithm is then used to estimate the target's TSC at different SNRs, and the corresponding RMSE is calculated, completing the comparative experiment.Figure 9 shows the variation in the target TSC with the SNR in the absence of interference, indicating similar performance between the initial QPSK waveform and the optimized waveform.In the experiment for Figure 10, based on the settings of the experiment in Figure 9, several interfering targets are added.The delay taps and Doppler taps of the interfering targets are randomly selected from [0, M − 1] and [0, N − 1] , respectively, and the TSCs are complex Gaussian random variables with a mean of 0 and a variance of 1.The MF-TDaPE algorithm is then used to estimate the targets' TSCs at different SNRs, and the corresponding RMSE is calculated, completing the comparative experiment.Figure 10 illustrates the TSC estimation performance in the presence of interference.It is evident that the performance is degraded.The reason is that the interference target affects TSC estimation through AF sidelobes.At the same time, it can also be observed that the performance of the algorithm deteriorates as the number of interfering targets increases.

Simulation of SDR-BD
This simulation showcases the beam patterns designed by SDR-BD under different system configurations and QoS settings.The specific simulation parameters for this study are detailed in Table 3. 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 an improved angular resolution of the system.

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 param-eter estimation, the CP-OTFS MF-DaPE algorithm is proposed, demonstrating superior performance in high-speed scenarios compared to OFDM.To enhance the robustness of parameter estimation against interference, the CP-OTFS AFS algorithm is employed for waveform design.The waveform designed by CP-OTFS AFS reduces the ISL of the corresponding AF and improves the 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 a single target, CN, and ERN; (2) CP-OTFS AFS incurs a significant computational burden.In future work, we will focus on extending the ISWSCPTS to multiple target, CN, and ERN scenarios and optimizing CP-OTFS AFS to improve its convergence speed.

Figure 2 .
Figure 2. The structure of CP-OTFS.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 an RF transmitter module.

Figure 3 .
Figure 3.The output of MF.

Figure 6 .
Figure 6.AF of initial waveform.The different colors in the picture represent different magnitudes.

Figure 7 .
Figure 7. AF of optimized waveform.The different colors in the picture represent different magnitudes.

Figure 11 .
Figure 11.Designed beamforming of different N t .

Figure 12
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 dB and ϱ = 5 mW, 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 the respective QoS requirements.

Figure 12 .
Figure 12.Designed beamforming of different QoS settings.

Table 1 .
Shortcomings of existing technologies.
• A 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, PT, and constant modulus constraints, a nonconvex 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.

Table 3 .
Specific simulation parameters for SDR-BD.