## Illumination Design for Joint Imaging and Wireless Power Transfer Systems

Qianyu Yang, Student Member, IEEE , Haiyang Zhang, Member, IEEE , Chunguo Li, Senior Member, IEEE , Ruiqi Liu, Member, IEEE , Baoyun Wang, Senior Member, IEEE

Abstract -This paper presents a novel concept termed Integrated Imaging and Wireless Power Transfer (IWPT), wherein the integration of imaging and wireless power transfer functionalities is achieved on a unified hardware platform. IWPT leverages a transmitting array to efficiently illuminate a specific Region of Interest (ROI), enabling the extraction of ROI's scattering coefficients while concurrently providing wireless power to nearby users. The integration of IWPT offers compelling advantages, including notable reductions in power consumption and spectrum utilization, pivotal for the optimization of future 6G wireless networks. As an initial investigation, we explore two antenna architectures: a fully digital array and a digital/analog hybrid array. Our goal is to characterize the fundamental tradeoff between imaging and wireless power transfer by optimizing the illumination signal. With imaging operating in the near-field, we formulate the illumination signal design as an optimization problem that minimizes the condition number of the equivalent channel. To address this optimization problem, we propose an semi-definite relaxation-based approach for the fully digital array and an alternating optimization algorithm for the hybrid array. Finally, numerical results verify the effectiveness of our proposed solutions and demonstrate the trade-off between imaging and wireless power transfer.

Index termsImaging, wireless power transfer, near field, integrated imaging and wireless power transfer (IWPT).

a designated Region of Interest (ROI) [9], [10]. However, in the context of forthcoming wireless systems, the allure of nearfield imaging, particularly through the utilization of extremely large Massive MIMO (XL-MIMO) and high-frequency bands, has garnered attention. Initial investigations have delved into the beamforming challenges intrinsic to this approach [11], hinting at a future where imaging could seamlessly integrate into the IoE sensor network without necessitating a dedicated imaging infrastructure. On the other hand, the concept of radio wave-powered devices has been a focal point in pioneering discussions [12], [13]. Notably, within the realm of wireless communication, there is a growing interest in amalgamating WPT with existing systems, such as simultaneous wireless information and power transfer (SWIPT) [14] and integrated sensing and wireless power transfer (ISWPT) [15]. These technologies aim to streamline communication and power delivery at the transmitter's end to curtail device expenses. Despite these advancements, the integration of imaging and WPT remains unexplored due to the traditionally distinct nature of these fields and their operational paradigms.

## I. INTRODUCTION

Leveraging high-frequency bands and extensive antenna arrays is anticipated to empower 6G systems to predominantly function within the radiation near-field region, offering unparalleled communication and sensing capabilities, flexibility, and resolution [1], [2]. Within this framework, notable research contributions have underscored the potential of holographic imaging [3], a technique widely utilized in various environmental sensing applications such as security inspections and medical diagnoses. Typically integrated into internet of things (IoE) devices for decision-making processes, these applications face limitations imposed by device size and hardware costs, often constrained by battery capacities [4], [5]. Radio frequency-based wireless power transfer (WPT) emerges as a promising technology for wirelessly charging such devices, presenting a viable solution to address energy scarcity concerns [6], [7].

Traditionally, imaging and WPT have been explored as separate domains. Imaging tasks have conventionally relied on technologies like microwaves or visible light, including lidars [8], with limited exploration on repurposing wireless communication systems for imaging unfamiliar objects within

To exploit higher array gain in such scenarios, precoding and combining techniques play a crucial role. The most adaptable solution for an array of radiating elements is the fully digital architecture, where each antenna element connects to a dedicated Radio Frequency (RF) chain, enabling precise beam control. However, the scalability of this approach poses challenges due to the substantial costs associated with expensive RF components and high power consumption [16]. To address the complexity of fully digital architectures, hybrid solutions have emerged, blending traditional antenna technologies with the integration of multiple antennas or panels to an RF chain through dedicated analog circuitry [17]. These hybrid architectures optimize resources by employing fewer RF chains relative to antenna elements. They achieve this by combining low-dimensional digital processing with high-dimensional analog precoding techniques, such as phaseshifter networks [18]. Such methodologies find applications in various large-scale antenna technologies, including near-field communications [19] and wideband communications [20].

In this paper, we propose a novel concept that merges imaging and WPT into a unified hardware platform, terming the envisioned scenario as integrated imaging and wireless power transfer (IWPT). The proposed IWPT system is envisioned for deployment in common IoE settings that necessitate imaging functionalities, like security inspections and nondestructive testing, enabling the simultaneous operation of imaging and charging sensors. The integration of imaging and WPT within a single system presents several advantages. Firstly, it streamlines device costs and optimizes spectrum resources compared to conventional separate designs, as WPT can be achieved by utilizing the imaging beams, obviating the necessity for additional power beam configurations. Moreover, the continuous coverage provided by the illuminating beam over the ROI ensures a stable energy supply to devices within this area. Additionally, harvesting energy from the illumination signals enhances the overall energy efficiency of the platform. Furthermore, the proximity of charge sensors to the ROI typically situates them within the radiation near-field range of the transmitter, which amplifies the efficiency of WPT due to the near-field propagation mechanism, thereby offering additional benefits to the system's energy transfer capabilities.

In the envisioned IWPT scenario, both imaging and WPT functionalities are consolidated within a singular hardware platform, operating on a unified waveform. As these functions utilize the same antenna and transmission power resources, a crucial performance trade-off emerges between imaging and WPT. This trade-off is fundamentally characterized by optimizing the beamforming vector of the transmitted signal within the IWPT system. The primary contributions of this communication can be succinctly summarized as follows:

- · To our knowledge, we are the first to propose the concept of IWPT, which aims to enhance cost-effectiveness and improve spectrum and energy efficiency by leveraging a shared hardware platform, spectrum, and power resources for two functions. This aligns with the anticipated development trajectory of future 6G networks. In our investigation, we explore two array architectures for this novel concept: the fully digital array and the hybrid array.
- · Additionally, we employ the condition number of the equivalent channel as a metric to evaluate the quality of imaging, providing an effective measure for performance assessment. We analyze the trade-off between imaging and WPT by optimizing the transmit beamforming vectors. Although the beamforming design problem is nonconvex, we tackle it iteratively using the semidefinite relaxation approach and successive convex approximation method. Furthermore, for the challenge of coupling digital and analog precoding in the hybrid array, we develop an alternating optimization algorithm to address this issue.
- · We provide numerical results to demonstrate the fundamental performance trade-off between imaging and WPT. Numerical results also verify the effectiveness of the proposed beamforming designs for the novel IWPT system.

Notation: Scalar variables, vectors and matrices are represented with lower letters, lower bold letters, and capital bold letters, respectively (e.g., x , x , and X , respectively). The term C N × N denotes a complex space of dimension N × N . ( ) · H denotes the Hermittian operators. For a matrix R R , ⪰ 0 means that R is positive semidefinite.

The rest of this paper is organized as follows: Sec. II introduces the considered array architecture, the signal model and the performance metrics of two functions. In Sec. III,

Fig. 1. Illustration of a joint imaging and WPT scenario.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-e68d1b4a-2085-42fd-8a42-acf59ce99005.jpg)

the design of beamforming design scheme is discussed, while Sec. IV details numerical simulations. Finally, Sec. V provides concluding remarks.

## II. SYSTEM MODEL

In this section, we provide a detailed description of the propsed ISWPT system. We begin by outlining the signal and channel model in Section II-A. Following that, in Sec. II-B, we introduce the two array architectures utilized throughout this paper. Subsequently, we introduce the imaging and WPT metrics in II-C and formulate the beamforming design problem for ISWPT systems in II-D.

## A. Transmitted Signal and Channel Model

We consider a joint imaging and wireless power transfer system as illustrated in Fig 1, wherein the transmitting and receiving discrete antenna arrays (TX/RX) are colocated and equipped with N antennas, the n -th antenna is located at a n = [ x , y a n a n , z a n ] , n = 1 2 , , · · · , N . The TX illuminates the ROI with the illumination signal and serves M energy receivers simultaneously, where the m -th energy receiver is located at p m = [ x p m , y p m , z p m ] , m = 1 2 , , · · · , M . The RX observes the reflected signal of the illumination signal from the ROI. The ROI is divided into a grid of K square cells of size ∆ , whose positions are r k = [ x , y r k r k , z r k ] , k = 1 2 , , · · · , K . The k -th cell is characterized by a scattering coefficient defined as γ k , related to the radar cross section (RCS) of the scatterer included in the cell (if any), which have RCS k = | γ k | 2 2 λ/ π 4 . If the k -th cell is empty, then γ k = 0 . The magnitude of the scattering coefficient | γ k | is upper bounded by the maximum RCS from a scatterer of area ∆ 2 , which corresponds to the RCS of a perfect electric conductor (PEC) having the same area given by RCSmax = ∆ 2 2 λ/ π 4 . Therefore, we have | γ k | ≤ ∆ .

Given that energy receivers can harvest energy from the illumination signal, there is no necessity for designing dedicated WPT signals. In this context, we contemplate employing solely the illumination signal to facilitate both imaging and WPT functionalities. Consequently, the transmitted baseband signal can be represented as x = [ x , x 1 2 , · · · , x N ] T ∈ C N × 1 . To optimize WPT efficiency, we make the assumption that the transmitter (TX) consistently transmits signals at the maximum power level, such that ∥ x ∥ 2 = P t , where P t signifies the total transmit power budget.

Assuming both imaging and WPT occur in the radiating near-field, i.e., the distance, d , between the TX/RX and the energy receiver or ROI should ensure that d &lt; 2 D 2 λ , where D and λ representing the array diameter and the signal wavelength, respectively. Define H T ∈ C K × N , H R ∈ C N × K and G ∈ C N × M as the imaging transmission channel matrix, the imaging receive channel matrix, and the power transmission channel matrix, respectively. Since the TX and RX is colocated, we have H R = H T T . Then, under the line of sight (LOS) channel condition, the element of the imaging transmission could be expressed as

$$\mathbf H _ { T } [ k, n ] = \mathbf H _ { R } [ n, k ] = F ( \Theta _ { k, n } ) \frac { \lambda } { 4 \pi d _ { k, n } } e ^ { - j \frac { 2 \pi } { \lambda } d _ { k, n } }, \quad ( 1 )$$

and the element of the power transmission channel matrix could be expressed as

$$\mathbf G [ m, n ] = F ( \Theta _ { m, n } ) \frac { \lambda } { 4 \pi d _ { m, n } } e ^ { - j \frac { 2 \pi } { \lambda } d _ { m, n } }, \quad \ ( 2 ) \quad \text{phase} \quad \text{comp} ]$$

where d k,n and d m,n are the distance between the corresponding antenna and the ROI cell or the energy receiver, F (Θ k,n ) and F (Θ m,n ) are the radiation profile, evaluated in the direction of arrival (i.e., Θ k,n and Θ m,n ), which modeled as [21]. For convenience, in the following study, we assume the channel information is known.

## B. Array Architecture

We consider two types of antenna architectures to generate x , depicted in Fig. 2, that are introduced separately in the following subsections. For convenience, the antennas of all architectures are placed in a rectangular arrangement with N d rows and N e columns, with an overall of N = N N d e elements.

a) Fully digital architecture: The implementation of fully digital antennas, where all processing is conducted in the digital domain, is widely preferred for its remarkable flexibility and the extensive freedom it offers for beamforming [22], notwithstanding the potential for increased costs when deploying large-scale arrays. In architectures employing fully digital antennas, each antenna element is linked to a dedicated Radio Frequency (RF) chain, as shown in Fig. 2(a). The transmitted baseband signal, denoted as x ∈ C N × 1 , is directly conveyed to the antenna via the RF chain subsequent to digital precoding.

b) Hybrid architecture: In our investigation, we center our attention on a partially connected hybrid array configuration, as shown in Figure 2(b). Within this arrangement, the quantity of Radio Frequency (RF) chains, designated as N d , is less than the overall count of antenna elements. Each digital signal output links to an RF chain, and subsequently, each RF chain output connects to the transmit antennas through a partially connected analog beamforming network reliant on

(a) Fully digital architecture

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-d9748bda-8a4d-47cf-8204-d7b22a0513a5.jpg)

(b) Partially connected Hybrid architecture

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-12731af5-b7c8-4c1a-b960-ec42ccfed67a.jpg)

Fig. 2. Considered array architectures.

phase shifters. This hybrid array architecture adeptly balances complexity and performance, offering a cost-efficient resolution for expansive antenna arrays. The array is specifically partitioned into N d linear sub-arrays, each consisting of N e elements uniformly arranged. Each linear sub-array is managed by an autonomous phase-shifter network associated with a single RF chain. Consequently, the signal transmitted by the hybrid array is given by

$$x = Q \, \text{w}.$$

Here w ∈ C N d × 1 is the the digital precoding output for the hybrid array, and the matrix Q ∈ C N d × N represents analog precoding of the hybrid array that satisfies

̸

$$\mathbf Q _ { n, ( i - 1 ) N _ { e } + l } = \begin{cases} q _ { i, l } \in \mathcal { F } & i = n, \\ 0 & i \neq n. \end{cases} \quad \ ( 4 )$$

where q i,l denotes the phase of phase-shifter, which satisfies

$$q _ { i, l } \in \mathcal { F } \stackrel { \triangle } { = } \left \{ e ^ { j \phi _ { i, l } } \ | \ \phi _ { i, l } \in [ 0, 2 \pi ] \right \}, \quad \forall i, l. \quad ( 5 )$$

Note that we did not distinguish the transmission signals of the two array architectures, but unified them as x because they follow the same constraints. The difference is that fully digital antennas can directly design x , while hybrid arrays require a joint design scheme for the digital and analog domains to obtain x as (3).

## C. Imaging and WPT Metric

According to the process of the illumination signal reflection as mentioned, the received signal on RX can be expressed as

$$y = H _ { R } \Gamma H _ { T } x + n,$$

where Γ ∈ C K × K is a diagonal matrix with the diagonal element Γ k,k = γ k , n is is an additive white Gaussian noise with variance σ 2 . Let γ = [ γ , γ 1 2 , · · · , γ K ] T , we have Γ = diag ( γ ) , our aim is to estimate γ from the received signal, and (6) could be rewritten as

$$y = \mathbf H _ { R } \text{diag} \left ( \mathbf H _ { T } \mathbf x \right ) \gamma + \mathbf n.$$

Without losing generality, we assume that N &gt; K , and define H = H R diag ( H x T ) , H can be equivalent regarded as a MIMO channel, then estimate γ from (7) becomes a conventional signal recovery problem. Using the common least squares method, the estimated scattering coefficients can be expressed as

$$\hat { \gamma } = \mathbf H ^ { \dagger } \mathbf y, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where H † is the pseudo-inverse of the equivalent channel matrix H . Once G T and G R are fixed, H would also determined by the illumination signal, so it is possible to improve the estimation accuracy of ˆ γ in imaging by design x .

To enhance the accuracy of estimation in (8), one direct strategy involves scrutinizing the mean square error (MSE), defined as E ∥ γ -ˆ γ ∥ 2 . However, estimating γ can often present challenges due to ill-posedness [11], particularly when H T and H R exhibit rank deficiency in practical near-field scenarios, leading to a rank-deficient composite channel matrix H . The potential ill-posed nature of the problem renders the analysis of MSE overly intricate for signal design purposes. Instead, in this paper, we choose the condition number of the equivalent channel as the metric of imaging.

The condition number provides insight into the efficient characterization of the near field MIMO channel, in general, can be defined as

$$\text{cond} = \frac { \lambda _ { m a x } } { \lambda _ { m i n } }, \text{ \quad \ \ } ( 9 ) \text{ } R ^ { * } \text{ case}$$

where λ max and λ min are the maximum and minimum eigenvalue of H , respectively. When the condition number of a channel approaches 1, the channel is deemed wellconditioned, indicating that its spatial characteristics are in an optimal state. Conversely, a channel tends towards infinite condition numbers when it is rank deficient. Consequently, our objective is to minimize the condition number of the equivalent channel through the design of the illumination signal. This optimization problem can be formulated as

$$\min _ { \ } c o n d \quad & \quad \text{the} \\ s. t. \ \| \mathbf x \| ^ { 2 } = P _ { t }. \ \text{ \quad \ \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \text{ \quad \ } \end{ \text{the}$$

For WPT operation, the received power of the m -th user is contingent on the transmitted signal (here is the illumination signal), which can be represented as

$$P _ { m } \left ( { \mathbf R } \right ) = \zeta { \mathbf g } _ { m } { \mathbf R } { \mathbf g } _ { m } ^ { H }, \quad \ \ \ ( 1 1 ) ^ { \quad \ \varkappa _ { P } \mathbf m }$$

where 0 &lt; ζ &lt; 1 is the constant energy conversion efficiency, and R = xx H denotes the covariance matrix of the illumination signal, g m ∈ C 1 × N is the m -th row of G , denotes the wireless channel between the transmitter and the m -th energy receiver.

In general, the goal of WPT is to maximize the sumharvested power, while adhering to the total transmit power constraint, which can be mathematically formulated as the following maximization problem:

$$\underset { x } { \max } \sum _ { m = 1 } ^ { M } P _ { m } \left ( R \right ) = \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \\ s. t. \text{ } R = \text{xx} ^ { H }, \\ \| \text{x} \| ^ { 2 } = P _ { \text{t} }.$$

Following [14], we can obtain a closed-form optimal solution of (12) directly, which expressed as

$$x ^ { * } = \sqrt { P _ { t } } \mathbf u _ { 1 } \left ( G \right ),$$

where u 1 ( G ) denotes the eigenvector corresponding to the maximal eigenvalue of the matrix G .

## D. Problem Formulation

The result derived from (13) is improbable to represent the optimal solution of (10), indicating the presence of an inherent trade-off between imaging and WPT performance in our systems. In this paper, we would like to characterize this fundamental performance trade-off between imaging and WPT by optimizing the illumination signal x . Mathematically, our interested beamforming design problem is formulated as

$$\min \, \text{cond} \, \begin{array} { c } x \\ s. t. & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \| \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) \\ \| x \| ^ { 2 } = P _ { t }. \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRIG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace } \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{arrayarray} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \ end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \leq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c + 1 } \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { \varphi } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \left ( \text{GRG} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end {array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \quad \begin{array} { c } ( 1 4 ) & \zeta \text{trace} \end{array} \end{array} \end{array} \end{array} \end{} \end{} \end{} \end{} \end{} \end{} \end{} \end{}$$

Where E r is the set received power threshold that satisfies E r ∈ [0 , E max ] , here E max is the achievable maximum received power calculated as E max = ζ trace ( GR G ∗ H ) , and R ∗ is the corvirance matrix of x ∗ from (13). The extreme case, e.g., E r = 0 or E r = E max corresponds to the scenario of only imaging or only WPT.

Problem (14) serves as a comprehensive optimization challenge applicable across all antenna architectures. In the forthcoming section, our focus will revolve around investigating the precoding scheme aimed at acquiring x that minimizes (14) within diverse antenna architectures.

## III. JOINT BEAMFORMING DESIGN

In this section, we study beam precoding design to achieve the performance trade-off between imaging and WPT. We first explicitly express the relationship between the objective function of our optimization problem and the optimization variables in Sec.III-A. Subsequently, we delve into resolving the aforementioned issue within the framework of a fully digital array in Section III-B, followed by the derivation of the optimal hybrid phase shifter configuration in Section III-C.

## A. Problem Reformulation

Problem (14) is hard to solve due to the absence of a direct relationship between x and the condition number as defined by (9). Inspired by [23], the problem of minimizing the condition number of H can be equivalent to minimizing the trace of the product of H and its transpose. This equivalence arises from the fact that minimizing the condition number is akin to minimizing the maximum eigenvalue λ max . Specifically, in accordance with the definition of the condition number, we have

$$\text{cond} \leq \left ( \frac { \lambda _ { m a x } } { \lambda _ { m i n } } \right ) ^ { 2 }. \text{ \quad \quad \ \ } ( 1 5 )$$

Thus the minimization problem can be regarded as minimize λ 2 max . Further more, for the semidefinite matrix HH H , we have

$$\lambda _ { m a x } ^ { 2 } < \text{trace} \left ( \text{HH} ^ { H } \right ). \quad \quad \quad \left ( 1 6 \right ) \quad \quad \quad$$

Based on (15) and (16), problem (14) can be equivalent formulated as

$$\min _ { x } \, \text{trace} \, \left ( \text{HH} ^ { H } \right ) \, \text{ \quad \ \ } \left ( 1 7 a \right ) \, \text{ one}$$

$$s. t. \ H = \mathbf H _ { R } \text{diag} \left ( \mathbf H _ { T } \mathbf x \right ) \quad \ \ ( 1 7 b ) \ \stackrel { \dots$$

$$\zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \text{ \quad \ \ } \left ( 1 7 c \right )$$

$$\| \mathbf x \| ^ { 2 } = P _ { \mathbf t }. \quad \ \ \ ( 1 7 d ) \ \ \text{The}$$

In problem (17), the optimization variable x is tied to the objective function through (17b). However, the diagonalization process in (17b) continues to pose challenges in assessing the influence of x on the objective function. To address this issue, we employ a series of matrix operations to eliminate the constraint (17b) in problem (17):

Firstly, the objective function of (17) can be regarded as

$$\text{trace} \left ( \text{HH} ^ { H } \right ) = \text{vec} \left ( \text{H} \right ) ^ { H } \text{vec} \left ( \text{H} \right ), \quad \ \ ( 1 8 ) \quad \text{where}$$

where vec ( ) · denotes the vector stretching operation on matrix, and according to (17b), we have

$$\text{vec} \left ( \mathbf H \right ) = \left ( \mathbf I _ { N } \otimes \mathbf H _ { R } \right ) \text{vec} \left ( \text{diag} \left ( \mathbf H _ { T } x \right ) \right ). \quad \left ( 1 9 \right ) \text{ lower}$$

Substituting (19) into (18), so we have

$$\text{trace} \left ( \text{HH} ^ { H } \right ) & = \text{vec} \left ( \text{diag} \left ( \text{H} _ { T } x \right ) \right ) ^ { H } \left ( \text{I} _ { N } \otimes \text{H} _ { R } \right ) ^ { H } \left ( \text{I} _ { N } \cdots \\ & \quad \otimes \text{H} _ { R } \right ) \text{vec} \left ( \text{diag} \left ( \text{H} _ { T } x \right ) \right )$$

$$\stackrel { ( a ) } { = } | x ^ { H } \mathbf H _ { T } ^ { H } \mathbf H _ { R } ^ { H } \mathbf H _ { R } \mathbf H _ { T } x |, \quad \ \ ( 2 0 ) \sin _ { \infty } \, 1$$

where ( a ) = comes from that the calculation related to zero elements of vec ( diag ( H x T )) in above equation can be ignored.

Therefore, the problem (17) can be rewritten as

$$\min _ { x } | x ^ { H } \mathbf H _ { T } ^ { H } \mathbf H _ { R } \mathbf H _ { T } x | \text{ \quad \ \ } ( 2 1 a )$$

$$s. t. \ \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq$$

$$\| { \mathbf x } \| ^ { 2 } = P _ { \mathbf t }. \quad \quad \quad$$

## B. Optimal Digital Beamforming Design

Problem (21) is non-convex due to the quadratic equality constraint (21b), rendering it intricate to solve directly. To tackle this issue, we use the SDR approach [24] to transform (21). According to the definition of R , (21) can be rewritten as an equivalent quadratic semidefinite programming (QSDP)

problem with rank-1 constraints:

$$\i \, \text{with rank-1 constraints:} \quad \begin{matrix} \min \, \text{trace} \$$

Due to the non-convex rank-one constraint, the problem (22) is still a non-convex optimization problem. Next, we use the difference-of-convex (DC) programming to handle this rankone constraint. Specifically, for the semi-definite matrix R , the rank-1 constraint can be equivalent expressed as

$$\text{trace} ( R ) = \left \| R \right \| _ { 2 }.$$

Therefore, we can transform the rank-1 constraint of R into a penalty term to the objective function of the problem (22), and the transformed problem is formulated as

$$\text{transformed proiem is formulate as} \quad \begin{array} { c } \min \, \text{trace}$$

where ζ ≫ 0 is a penalty factor related to the rank-1 constraint. (24) is still non-convex due to the presence of a 2-norm ∥ R ∥ 2 in the objective function, so we utilize the SCA method to write a first order Taylor expansion of ∥ R ∥ 2 as its lower bound:

$$\| { \mathbf R } \| _ { 2 } & \geq \left \| { \mathbf R } ^ { ( t ) } \right \| _ { 2 } + \text{trace} \left ( { \mathbf u } _ { m a x } ^ { ( t ) } { \mathbf u } _ { m a x } ^ { ( t ) H } ( { \mathbf R } - { \mathbf R } ^ { ( t ) } ) \right ) \\ \dots & \stackrel { \hat { \varrho } } { = } f \left ( { \mathbf R } ^ { ( t ) }, { \mathbf R } \right ),$$

where R ( ) t is the feasible point of the first order Taylor expansion in the t -th iteration, u ( ) t max is the eigenvector corresponding to the largest singular value of the matrix R ( ) t . Hence, by replacing ∥ R ∥ 2 with f ( R ( ) t , R ) in the t -th iteration, (24) can be approximately transformed into

$$\mathfrak { e } \text{ approximately transformed into} \quad \min \, & \text{trace} \left ( H _ { T } ^ { H } H _ { R } ^ { H } H _ { R } H _ { T } R \right ) \\ & \quad + \eta \left ( \text{trace} ( R ) - f ( R ^ { ( t ) }, R ) \right ) \\ s. t. \ \zeta \text{trace} \left ( \text{GRG} ^ { H } \right ) \geq E _ { r } \\ & \quad \text{trace} ( R ) = P _ { t }, \\ & \quad R \succeq 0. \\ \left ( t \right ) \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array} \quad \text{(26)}$$

R ( ) t is iteratively updated. In each iteration, R ( ) t is taken from the solution of problem (26) in the previous iteration, and the resulting R ( ) t would iteratively close to the optimal solution of (24). Meanwhile, we set that R (1) = R ∗ to ensure the initial satisfies all constraints of (26). The resulting problem (26) is a typical QSDP problem, which is convex and can be solved directly using existing convex optimization solvers such as CVX [25].

Once R ( ) t reaches convergence, it can be approximated as a rank 1 matrix, i.e., the optimal solution of problem (22)

would be found. The jointly optimal beam x opt for imaging and WPT could be calculated from the eigendecomposition of R ( ) t directly, which expressed as

$$x ^ { * } = \lambda _ { m a x } ^ { ( t ) } \mathbf u _ { m a x } ^ { ( t ) }, \quad \quad ( 2 7 ) \, \begin{array} { c } \mu \lambda \\ \mathbf w \end{array}$$

where λ ( ) t max is the largest singular value of the matrix R ( ) t , since the transmit power constraint, we have λ ( ) t max = √ P t , thus (27) could be also expressed as x ⋆ = √ P t u ( ) t max .

The worst case complexity to solve the QSDP problem (26) is O ( N 6 5 . log(1 /ϵ ) ) with the primal-dual interior-point algorithm [26], where ϵ is the solution accuracy. We assume that the number of iterations required for the SCA process to reach convergence is K , thus the whole complexity of the beamforming design is expressed as O ( KN 6 5 . log(1 /ϵ ) ) .

## C. Joint Analog and Digital Design for Hybrid array

As outlined in Section II-B, in the scenario of a fully digital array, the optimal beamforming design can be directly acquired by addressing problem (14). Conversely, in the context of a hybrid array, considering (3), the beamforming design predicament necessitates a collective precoding design spanning the digital and analog domains. Analogous to (22), the optimization conundrum (14) under the hybrid array scenario can be reformulated as

$$\L a t e a \text{ as } & \min _ { \substack { | w ^ { H } Q ^ { H } T Q w | \\ Q, w } } & \text{cable} \\ s. t. & \zeta \text{trace} \left ( \text{GRG} \text{$H$} \right ) \geq E _ { r } & \text{cante} \\ & \| \text{Qw} \| ^ { 2 } = P _ { t } & \text{cinem} \\ & ( 4 ), ( 5 ). & \text{const} \\ \Lambda \quad \sim u \sim \Lambda \quad \sim \quad u \sim u$$

## Where T = H H H H H T H R R T , and R = Qww Q H H .

Since Problem (28) involves coupled optimization variables and unit modulus constraints, it is non-convex. To address this, we aim to derive a hybrid solution for Problem (28) that closely approximates the fully-digital solution of another problem, denoted as Problem (22). This approach is quite often employed for optimizing hybrid analog/digital systems. Essentially, our objective is to determine a hybrid precoding mapping that closely resembles the fully-digital precoding in the context of the Frobenius norm.

Specifically, let x ⋆ ∈ C N × 1 be the result of the optimal digital beamforming design by solving problem (26). The resulting surrogate optimization problem is given by

$$\sin _ { Q, w } \| x ^ { * } - Q w \| ^ { 2 } & & \sin _ { 2 } \mu \\ s. t. \ \| w \| ^ { 2 } & = \frac { P _ { t } } { N _ { d } } & & ( 2 9 ) & 3 \colon \\ & ( 4 ), ( 5 ). & & 5 \colon \\ &. \ \cdots &. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

Where the power constraint in (29) comes from the ∥ Qw ∥ 2 = P t in (28), according to (4) and (5), we have ∥ Qw ∥ 2 = trace ( w Q Qw H H ) = N d trace ( w w H ) .

We tackle (29) using the alternating optimization. In particular, for a given Q , the digital precoding result w ⋆ which minimizes (29) is stated in the following least-squares solution:

$$\dots \quad \text{w} ^ { * } = \left ( Q ^ { H } Q \right ) ^ { - 1 } Q ^ { H } x ^ { * } = Q ^ { H } x ^ { * } / N _ { d }, \quad ( 3 0 ) \ \text{proces}$$

the least-squares solution in (30) dose not satisfy the power constraint in (29). Nonetheless, once w and Q are tuned to optimize (29), we can update the digital precoder by multiplying a factor, i.e., w ⋆ = √ P t w ⋆ √ N d ∥ w ⋆ ∥ .

When optimizing Q with a fixed w , we leverage the partially connected structure constraint (4) and the unit-modulus constraint (5) on the matrix Q . Through these constraints, a closed-form solution for Q emerges. Specifically, when w is given, the phase shifter coefficient matrix Q that minimizes (29) can be attained as per the following proposed theorem.

Theorem 1. For a given w , the matrix Q which minimizes (29) is

̸

$$\mathbf Q _ { n, ( i - 1 ) N _ { e } + l } ^ { * } = \begin{cases} q _ { i, l } ^ { * } \in \mathcal { F } & i = n, \\ 0 & i \neq n. \end{cases} \quad ( 3 1 )$$

here q ⋆ i,l is the corresponding phase shifter coefficient which given by

$$q _ { i, l } ^ { * } = e ^ { j \left ( \zeta x _ { ( i - 1 ) N _ { e } + l } ^ { * } - \zeta \mathbf w _ { i } \right ) }.$$

where ∠ · denotes the angle of complex value.

$$P r o o f \colon \text{The proof is given in Appendix A.}$$

The result of Theorem. 1 indicates that the joint precoding result of the partially connected hybrid array can not perfectly match the result of a fully digital scheme. This discrepancy arises due to the insufficient number of RF chains to impeccably regulate the signal amplitudes at the antennas within the context of the considered partially connected architecture. This disparity isn't a consequence of the chosen scheme but rather stems from the inherent distortion induced by the structural constraints of partially connected hybrid arrays. Consequently, this mismatch engenders a suboptimal trade-off in performance for the hybrid array concerning IWPT compared to the fully digital array.

Finally, the overall resulting configuration algorithm for joint digital and analog precoding for the hybrid array is summarized as the following algorithm.

Algorithm 1 Alternating Optimization of Digital and Analog Precoding for Hybrid Array

Initialize: Initialize the matrix Q 0 ∈ C N d × N constrained by the array architecture with random phase set; Iterations limit T max .

- 1: Obtain x ⋆ by solving problem (26).

While t ≤ T max do

- 2: Calculate w t as (30) with the given matrix Q t -1 and update as w t = √ P t w t √ N d ∥ w t ∥ .
- 4: Update the matrix Q t as (31).
- 3: Update the coefficient q t i,l as (32) with the given w t .
- 5: t = t +1

End while

Output: Q = Q T max ; w = w T max .

In Algorithm 1, it can be seen that Q only needs one update to achieve phase matching based on Theorem. 1, while the update of w is only related to Q by (30), so the algorithm only needs two iterations to converge. In addition, since the update process only involves closed form solution calculation, the

Fig. 3. Imaging results in different illumination design schemes: (a) Original pattern; (b) Random; (c) Imaging only; (d) WPT only; (e) IWPT with fully digital array, where the power threshed is E r = 0 15 . E max ; (f) IWPT with hybrid array, where E r = 0 15 . E max .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-07b4048a-e8df-4ebb-8e07-b009115f7a7c.jpg)

complexity of Algorithm 1 mainly comes from the calculation of x ⋆ , which be expressed as O ( KN 6 5 . log(1 /ϵ ) ) .

TABLE I EXPERIMENT PARAMETERS

## IV. NUMERICAL EVALUATIONS

## A. Simulation Parameters

In this section, we provide numerical results to demonstrate the fundamental trade-off between imaging and WPT under a novel IWPT scenario. In our experiments, the TX/RX is set at configured as uniform squared arrays, with the array's reference position situated at (0 , 0 0) , m on the Y Z-plane. The region of ROI forms a rectangle centered at (2 0 0) , , m and is parallel to the Y Z-plane. We direct the TX energy towards three power users concurrently, positioned at (1 5 . , 1 1) , m, (1 , -1 5 . , 0) m, (1 5 . , -1 0) , m, respectively.Given these specifications, the maximum received power for WPT can be computed as E max. Furthermore, the other simulation parameters are given in Table I.

The numerical indicators of imaging performance are represented by root mean square error (RMSE), calculated as

$$R M S E = \sqrt { \frac { 1 } { N _ { t } } \sum _ { t = 1 } ^ { N _ { t } } \| \hat { \gamma } ^ { ( t ) } - \gamma \| ^ { 2 } }, \quad \ \ ( 3 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
R M S E = \sqrt { \frac { 1 } { N _ { t } } \dots \, \dots } \,.$$

where N t is the number of Monte Carlo experiments, which set as N t = 200 in this paper, and ˆ γ ( ) t is the estimate in the t -th experiment.

## B. Numerical Results

To assess the imaging efficiency of the proposed beamforming scheme, we initially establish the scenario where γ k = ∆ when the k -th cell of the Region of Interest (ROI) contains a scatter to be detected; otherwise, γ k = 0 . The image under consideration for estimation is illustrated in Fig. 3(a). Our

| Parameter                    | Value              |
|------------------------------|--------------------|
| Transmission bandwidth       | ∆ f = 120 kHz      |
| Central frequency            | f c = 28 GHz       |
| Wavelength                   | λ = 0 . 01 m       |
| Noise power spectral density | σ 2 = - 170 dBm/Hz |
| Transmitted power            | P t = 30 dBm       |
| Energy conversion efficiency | ζ = 0 . 5          |
| Array antenna row            | N d = 13           |
| Array antenna column         | N e = 13           |
| Antenna Spacing              | 3 λ/ 2             |
| Antenna Configuration        | 13 × 13 antennas   |
| Cell Distribution of ROI     | 10 × 10 cells      |
| Cell Size                    | ∆ = 0 . 1 m        |
| TX/RX Near Field Distance    | d F = 6 . 5 m      |

objective is to estimate the scattering coefficients γ , which represent the original image, under various illumination signal beamforming strategies.

In our results, following cases are considered: 1) Random: In this scenario, the illumination signal is generated as a random signal, serving as a fundamental performance benchmark. The estimated outcome is depicted in Fig. 3(b); 2) Imaging only: The beamforming process focuses solely on optimizing imaging performance. The corresponding illumination is derived by solving (26) with E r = 0 . The resulting estimate is illustrated in Fig. 3(c); 3) WPT only: Here, the illumination signal is directly obtained via (13). The beamforming aims at maximizing the total received power of the power users exclusively. The estimated outcome is shown in Fig. 3(d); 4) IWPT with fully digital array: The optimal signal beam for imaging and Wireless Power Transfer (WPT) is determined by solving (26) from the proposed beamforming design scheme, with the power threshold set at E r = 0 15 . E max. The estimated result is displayed in Fig. 3(e); 5) IWPT with hybrid array: The signal beam is computed using (3), where the digital and analog precoding are obtained from the proposed Algorithm 1.

The random illumination scheme demonstrates a fundamental ability to recover the original image from the received signal overall, albeit with some residual distortion. In contrast, the imaging-only scenario exhibits the most effective image restoration, showcasing how optimizing the equivalent channel conditions can significantly enhance imaging performance, while the WPT-only estimation yields the most distorted results, underscoring the impact of WPT considerations on imaging fidelity. Furthermore, the estimation outcomes of the IWPT scenario with a fully digital array closely resemble those of the imaging-only scenario. This alignment underscores the efficacy of our proposed beamforming optimization design. Conversely, the performance of the IWPT scheme with the hybrid array falls between that of the fully digital setup and the random scheme. This placement corroborates our earlier analysis in Sec. III-C, but it should be emphasized that despite this, the hybrid array outperformed the random scheme while significantly cutting down the required number of RF chains by a factor of 13 compared to the all-digital array.

In Fig.4, we illustrate the trade-off between imaging and Wireless Power Transfer (WPT). For detailed insights, Fig.4(a) and Fig. 4(b) depict the correlation between the condition number and received power, as well as the Root Mean Square Error (RMSE) and received power, respectively. The comparison points 'Random,' 'Imaging only,' and 'WPT only' are derived from averaging 200 Monte Carlo experiments within the mentioned context. It is evident that both the condition number of the equivalent channel and the RMSE of imaging increase as the received power of users rises, indicating a fundamental trade-off between imaging and WPT. Enhancing performance in one aspect typically comes at the expense of the other. Our proposed joint optimal scheme can achieve two extremes-the imaging-only scenario and the WPT-only scenario-and offers a flexible trade-off by adjusting the power threshold. As previously discussed, the performance of the hybrid array setup is suboptimal compared to the fully digital array case, with noticeable fluctuations in the lines as the power threshold increases. This variability occurs as the trade-off gradually shifts towards WPT, suggesting that a certain level of distortion may benefit imaging and vice versa. Nonetheless, this phenomenon does not alter the overarching trend of the trade-off performance.

The comparison between Fig. 4(a) and Fig. 4(b) reveals a positive correlation between the condition number and the Root Mean Square Error (RMSE) of imaging, particularly noticeable in the hybrid array scenario where the RMSE results exhibit similar fluctuations as the condition number. Despite a non-linear relationship, the RMSE results are indirectly influenced by optimizing the condition number. Given the complexity introduced by the near-field environment, analyzing imaging RMSE directly can be intricate, making the condition number a simpler yet effective optimization target. Furthermore, our joint optimization strategy based on imaging

(a) Accessible condition number versus received power.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-505b9c0e-6a69-48cc-aaed-474c8bf31ed2.jpg)

Fig. 4. Trade-off between iamging and WPT under different power threshold: (a) Condition number - received power; (b) RMSE - received power.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-32336162-10c2-4121-a5c5-d4bc084b0aa5.jpg)

and WPT indicates that even under high power constraints, the imaging performance remains superior to that of random illumination. In Fig. 4(a), the power gating for improved imaging surpasses 12.2 µ W, while in Fig. 4(b), it reaches 14 µ W, showcasing the substantial potential of combining imaging with WPT.

The hybrid array design achieves a reduction in RF chains at the expense of a certain performance trade-off. In Fig.5, we delve into the correlation between the number of RF chains in the hybrid array and the condition number. Keeping other simulation parameters consistent with Sec. IV-A, we vary the number of antenna rows for TX/RX , i.e., N d , and calculate the corresponding condition number. As N d changes, the total antenna count N also changes, allowing us to compare the results with a fully digital array setup having an identical antenna configuration, the later requiring 13 × more RF chains than the hybrid array. We examine two scenarios: imaging only (i.e., the power threshold E r = 0 ) and IWPT with E r = 0 15 . E max . The trend reveals that the condition number decreases with an increase in RF chains (or antenna numbers in the fully digital array), leading to a

Fig. 5. Condition number contrast under different RF chain numbers, in IWPT case, the power threshed is set to E r = 0 15 . E max .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00368v1-figure-5e2cfb39-2374-45d2-b396-b312d45eebfd.jpg)

gradual reduction in the performance gap between the hybrid and fully digital arrays to a stable level. This suggests that the performance of the hybrid array hinges on that of the fully digital array it emulates, with a persistent performance gap due to RF chain reduction. The performance gap between hybrid arrays and fully digital arrays is narrower in the imagingonly scenario compared to the IWPT scenario. This difference is a consequence of the inherent trade-off between imaging and WPT, which exacerbates the mismatch when aligning the imaging beam with the hybrid array.

## V. CONCLUSIONS

In this paper, we introduced the concept of Integrated Imaging and Wireless Power Transfer (IWPT), which combines imaging and wireless power transfer functionalities on a unified hardware platform. IWPT utilizes a transmitting array to efficiently illuminate a specific Region of Interest (ROI), enabling the extraction of the ROI's scattering coefficients while simultaneously providing wireless power to nearby users. We characterized the fundamental trade-off between imaging and wireless power transfer of IWPT system by optimizing the illumination signals, under both the fully digital array and hybrid array architectures. With imaging operating in the near-field, we formulated the illumination signal design as an optimization problem aimed at minimizing the condition number of the equivalent channel. To solve this optimization problem, we proposed a semi-definite relaxationbased approach for the fully digital array and an alternating optimization algorithm for the hybrid array. Our numerical results validate the effectiveness of our proposed solutions and highlight the trade-off between imaging and wireless power transfer. These findings underscore the potential of IWPT systems for practical applications, demonstrating their capability to enhance imaging precision and power transfer efficiency in future wireless networks.

## APPENDIX A

## PROOF OF THEOREM 1

To drop the non-convex structure constraint (4), we define that q = vec ( Q ) and V = w T ⊗ I N . The objective function of (29) can be rewritten as

$$\| x ^ { * } - Q w \| ^ { 2 } = \| x ^ {$$

where q is the modified version of q obtained by removing all the zero elements of q and thus q ∈ C N × 1 , V is the modified version of V obtained by removing the elements having the same index as the zero elements of q . By substituting (A.1) into (29) with given w , (29) can be rewritten as a problem about optimizing q which given by

$$\min _ { \bar { q } } \, & \| x ^ { * } - \$$

where q n denotes the n -th element of q . In addition, according to the definition V , we obverse that V is a diagonal matrix with the diagonal elements vector expressed as [ w 1 l N T e , · · · , w N d l T N e ] T . Therefore, (A.2) can be futher decomposed into N d subproblems, each with an identical structure. In particular, each subproblem individually designs the weighting coefficients for a single RF chains, with the i -th subproblem in scalar form given by

$$\begin{array} {$$

Since q n is only related to the phase, the solution to problem (A.3) is q n = e j ( ∠ x ⋆ n -∠ w i ) . Therefore, according to the mapping relationship between q n and q i,l , i.e., q i,l = q ( i -1) N e + l , the optimal phase shifter coefficient q ⋆ i,l is given by

$$q _ { i, l } ^ { * } = e ^ { j \left ( \not \$$

This proves the theorem.

## REFERENCES

- [1] H. Zhang, N. Shlezinger, F. Guidi, D. Dardari, and Y. C. Eldar, '6G wireless communications: From far-field beam steering to near-field beam focusing,' IEEE Communications Magazine , vol. 61, no. 4, pp. 72-77, 2023.
- [2] R. Liu, H. Lin, H. Lee, F. Chaves, H. Lim, and J. Sk¨ld, 'Beginning of o the Journey Toward 6G: Vision and Framework,' IEEE Communications Magazine , vol. 61, no. 10, pp. 8-9, 2023.
- [3] F. Guidi, A. Guerra, and D. Dardari, 'Personal mobile radars with millimeter-wave massive arrays for indoor mapping,' IEEE Transactions on Mobile Computing , vol. 15, no. 6, pp. 1471-1484, 2016.
- [4] X. Lu, P. Wang, D. Niyato, D. I. Kim, and Z. Han, 'Wireless networks with RF energy harvesting: A contemporary survey,' IEEE Communications Surveys Tutorials , vol. 17, no. 2, pp. 757-789, 2015.
- [5] J. Hu, Q. Wang, and K. Yang, 'Energy self-sustainability in fullspectrum 6G,' IEEE Wireless Communications , vol. 28, no. 1, pp. 104111, 2021.
- [6] B. Clerckx, R. Zhang, R. Schober, D. W. K. Ng, D. I. Kim, and H. V. Poor, 'Fundamentals of wireless information and power transfer: From RF energy harvester models to signal and system designs,' IEEE Journal on Selected Areas in Communications , vol. 37, no. 1, pp. 4-33, 2019.
- [7] S. Li and C. C. Mi, 'Wireless power transfer for electric vehicle applications,' IEEE Journal of Emerging and Selected Topics in Power Electronics , vol. 3, no. 1, pp. 4-17, 2015.
- [8] W. Gong, C. Zhao, H. Yu, M. Chen, W. Xu, and S. Han, 'Threedimensional ghost imaging lidar via sparsity constraint,' Scientific reports , vol. 6, no. 1, p. 26133, 2016.
- [9] J. Gao, Y. Qin, B. Deng, H. Wang, and X. Li, 'Novel efficient 3D short-range imaging algorithms for a scanning 1D-MIMO array,' IEEE Transactions on Image Processing , vol. 27, no. 7, pp. 3631-3643, 2018.
- [10] X. Zhuge and A. G. Yarovoy, 'Three-dimensional near-field MIMO array imaging using range migration techniques,' IEEE Transactions on Image Processing , vol. 21, no. 6, pp. 3026-3033, 2012.
- [11] G. Torcolacci, A. Guerra, H. Zhang, F. Guidi, Q. Yang, Y. C. Eldar, and D. Dardari, 'Holographic imaging with XL-MIMO and RIS: Illumination and reflection design,' IEEE Journal of Selected Topics in Signal Processing , pp. 1-16, 2024.
- [12] H. Zhang, N. Shlezinger, F. Guidi, D. Dardari, M. F. Imani, and Y. C. Eldar, 'Near-field wireless power transfer for 6G internet of everything mobile networks: Opportunities and challenges,' IEEE Communications Magazine , vol. 60, no. 3, pp. 12-18, 2022.
- [13] R. Zhang and C. K. Ho, 'MIMO broadcasting for simultaneous wireless information and power transfer,' IEEE Transactions on Wireless Communications , vol. 12, no. 5, pp. 1989-2001, 2013.
- [14] J. Xu, L. Liu, and R. Zhang, 'Multiuser MISO beamforming for simultaneous wireless information and power transfer,' IEEE Transactions on Signal Processing , vol. 62, no. 18, pp. 4798-4810, 2014.
- [15] Q. Yang, H. Zhang, and B. Wang, 'Beamforming design for integrated sensing and wireless power transfer systems,' IEEE Communications Letters , vol. 27, no. 2, pp. 600-604, 2023.
- [16] T. Zirtiloglu, N. Shlezinger, Y. C. Eldar, and R. T. Yazicigil, 'Powerefficient hybrid MIMO reciever with task-specific beamforming using low-resolution ADCs,' in Proc. IEEE ICASSP , 2022.
- [17] O. Levy and N. Shlezinger, 'Rapid hybrid modular receive beamforming via learned optimization,' in IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP) , 2024, pp. 12 826-12 830.
- [18] S. S. Ioushua and Y. C. Eldar, 'A family of hybrid analog-digital beamforming methods for massive mimo systems,' IEEE Transactions on Signal Processing , vol. 67, no. 12, pp. 3243-3257, 2019.
- [19] H. Zhang, N. Shlezinger, F. Guidi, D. Dardari, M. F. Imani, and Y. C. Eldar, 'Beam focusing for near-field multiuser mimo communications,' IEEE Transactions on Wireless Communications , vol. 21, no. 9, pp. 7476-7490, 2022.
- [20] 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 Journal on Selected Areas in Communications , vol. 40, no. 10, pp. 2994-3009, 2022.
- [21] S. W. Ellingson, 'Path loss in reconfigurable intelligent surface-enabled channels,' in 2021 IEEE 32nd Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC) , 2021, pp. 829-835.
- [22] J. S. Herd and M. D. Conway, 'The evolution to modern phased array architectures,' Proc. IEEE , vol. 104, no. 3, pp. 519-529, 2016.
- [23] W. Huang, B. Lei, S. He, C. Kai, and C. Li, 'Condition number improvement of irs-aided near-field mimo channels,' in 2023 IEEE International Conference on Communications Workshops (ICC Workshops) , 2023, pp. 1210-1215.
- [24] Z.-Q. Luo, W.-K. Ma, A. M.-C. So, Y. Ye, and S. Zhang, 'Semidefinite relaxation of quadratic optimization problems,' IEEE Signal Processing Magazine , vol. 27, no. 3, pp. 20-34, 2010.
- [25] M. Grant and S. Boyd, 'CVX: Matlab software for disciplined convex programming, version 2.1,' http://cvxr.com/cvx, Mar. 2014.
- [26] K.-C. Toh, 'An inexact primal-dual path following algorithm for convex quadratic SDP,' Mathematical programming , vol. 112, no. 1, pp. 221254, 2008.
Qianyu Yang Haiyang Zhang Chunguo Li Ruiqi Liu Baoyun Wang
Illumination Design for Joint Imaging and Wireless Power Transfer Systems
This paper presents a novel concept termed Integrated Imaging and Wireless
Power Transfer (IWPT), wherein the integration of imaging and wireless power
transfer functionalities is achieved on a unified hardware platform. IWPT
leverages a transmitting array to efficiently illuminate a specific Region of
Interest (ROI), enabling the extraction of ROI's scattering coefficients while
concurrently providing wireless power to nearby users. The integration of IWPT
offers compelling advantages, including notable reductions in power consumption
and spectrum utilization, pivotal for the optimization of future 6G wireless
networks. As an initial investigation, we explore two antenna architectures: a
fully digital array and a digital/analog hybrid array. Our goal is to
characterize the fundamental trade-off between imaging and wireless power
transfer by optimizing the illumination signal. With imaging operating in the
near-field, we formulate the illumination signal design as an optimization
problem that minimizes the condition number of the equivalent channel. To
address this optimization problem, we propose an semi-definite relaxation-based
approach for the fully digital array and an alternating optimization algorithm
for the hybrid array. Finally, numerical results verify the effectiveness of
our proposed solutions and demonstrate the trade-off between imaging and
wireless power transfer.