## Flexible Beam Coverage Optimization for Movable-Antenna Array

Dong Wang, Weidong Mei, Boyu Ning, and Zhi Chen

National Key Laboratory of Wireless Communications

University of Electronic Science and Technology of China (UESTC), Chengdu 611731, China

Emails: DongwangUESTC@outlook.com; wmei@uestc.edu.cn; boydning@outlook.com; chenzhi@uestc.edu.cn

Abstract -Fluid antennas (FAs) and movable antennas (MAs) have attracted increasing attention in wireless communications recently. As compared to the conventional fixed-position antennas (FPAs), their geometry can be dynamically reconfigured, such that more flexible beamforming can be achieved for signal coverage and/or interference nulling. In this paper, we investigate the use of MAs to achieve uniform coverage for multiple regions with arbitrary number and width in the spatial domain. In particular, we aim to jointly optimize the MAs' weights and positions within a linear array to maximize the minimum beam gain over the desired spatial regions. However, the resulting problem is non-convex and difficult to be optimally solved. To tackle this difficulty, we propose an alternating optimization (AO) algorithm to obtain a high-quality suboptimal solution, where the MAs' weights and positions are alternately optimized by applying successive convex approximation (SCA) technique. Numerical results show that our proposed MAbased beam coverage scheme can achieve much better performance than conventional FPAs.

## I. INTRODUCTION

Beamforming is a pivotal signal-processing technology for multiple-antenna systems, which is able to enhance/reduce the signal strength over desired/undesired directions by optimizing the amplitude and/or phase of the antenna weight vector (AWV) [1]. However, due to the fixed geometry of conventional fixedposition antenna (FPA) arrays, their steering vectors have inherent correlation/orthogonality among different steering angles, thus resulting in difficulty in achieving flexible beamforming for interference nulling and/or beam coverage.

Recently, fluid antennas (FAs) and movable antennas (MAs) have emerged as a new technology that enables local antenna movement within a confined region at the transmitter/receiver, thereby providing additional spatial degrees of freedom (DoFs) to enhance the wireless communication performance [2]-[6]. Inspired by its compelling benefits, there have been some prior studies focusing on its use for point-to-point communications, multi-user systems, secure communicationes, among others [7][13]. In particular, the authors in [14] and [15] have shown the capability of FAs/MAs to achieve flexible interference nulling, which is able to amplify signals over desired directions and mitigate interference over undesired directions at the same time.

Different from the above papers, we investigate in this paper the use of MAs for flexible beam coverage, aiming to achieve uniform beam gain over a single or multiple regions with

This work was supported by the National Natural Science Foundation of China under Grant 62027806.

Fig. 1. MA-assisted beam coverage for multiple subregions

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00434v2-figure-a521cdd4-71b9-481a-b16f-5a6e73def6bc.jpg)

different width in the spatial domain, as shown in Fig. 1. Note that such a beam coverage problem can find broad applications in wireless communications. For example, due to the mobility of users, the beamwidth can be dynamically adjusted to cover an area based on the distribution of the users within it [16]. Additionally, the number of covered areas may fluctuate over time due to varying network traffic demands. Consequently, more efficient beam coverage approaches are required compared to beamforming design with conventional FPAs only to accommodate diverse requirements for beamwidth and beam number, which motivates this work. In particular, we aim to jointly optimize the MAs' weights and positions within a linear array to maximize the minimum beam gain over the desired spatial regions. However, the resulting problem is non-convex and difficult to be optimally solved. To tackle this challenging problem, we propose an alternating optimization (AO) algorithm to obtain a high-quality suboptimal solution to it, where the MAs' weights and positions are alternately updated by applying the successive convex approximation (SCA) technique. Numerical results show that our proposed MA-based beam coverage scheme can achieve much better coverage performance than conventional FPAs.

Notations: R and C represent the sets of real and complex number, respectively. For a complex number a , Re { a } and | a | denote its real part and modulus, respectively. For a complexvalued vector x x , T , x H , ∥ x ∥ 2 and x ( n ) denote its transpose, conjugate transpose, l 2 -norm, n -th element, respectively. For a matrix A , Tr ( A ) represents its trace. x ∼ CN ( µ R , ) represents that x follows circularly symmetric complex Gaussian (CSCG) distribution with mean µ and covariance matrix R .

## II. SYSTEM MODEL

As depicted in Fig. 1, we consider a transmitter equipped with a linear array with N MAs, whose positions can be flexibly adjusted within a one-dimensional (1D) region denoted as C t . The length of the MA array is set to D . For convenience, we assume a reference origin point x 0 = 0 in the transmit region C t , as shown in Fig. 1. Let x n denote the coordinate of the n -th MA's position, n ∈ N ≜ { 1 2 , , · · · , N } . Accordingly, the antenna position vector (APV) of all N MAs is denoted by x ≜ [ x , x 1 2 , · · · , x N ] ∈ R N × 1 . Let f c and λ denote the operating carrier frequency of the considered system and carrier wavelength, respectively. For a given steering angle θ , the array response vector of the MA array can be determined as a function of the APV x and θ , i.e.,

$$a ( x, \theta ) = [ e ^ { j \frac { 2 \pi } { x } x _ { 1 } \cos \theta }, \cdots, e ^ { j \frac { 2 \pi } { x } x _ { N } \cos \theta } ] ^ { T } \in \mathbb { C } ^ { N \times 1 }. \quad ( 1 ) \quad \text{where $1$}$$

To ease the practical implementation, we assume analog beamforming with a constant magnitude at the transmitter, which is given by

$$\omega = \frac { 1 } { \sqrt { N } } [ e ^ { j \phi _ { 1 } }, e ^ { j \phi _ { 2 } }, \cdots, e ^ { j \phi _ { N } } ] ^ { T } \in \mathbb { C } ^ { N \times 1 }, \quad ( 2 )$$

where ϕ n denotes the phase shift at the n -th MA.

Based on the above, the beam gain over the steering angle θ can be expressed as

$$G ( \omega, x, \theta ) = \left | \omega ^ { H } a ( x, \theta ) \right | ^ { 2 }. \text{ \quad \ \ } ( 3 ) \text{ \ \ } \text{algor}.$$

As shown in Fig. 1, in this paper, we aim to achieve uniform beam coverage in the angular domain by jointly optimizing the APV x and the analog beamformer ω . In particular, we consider that the desired coverage region is composed of K disjoint subregions, i.e., R = R ∪ R 1 2 ∪ · · · ∪ R K , with R k = { θ ∈ [ θ ( k ) min , θ ( k ) max ] } , where θ ( k ) min and θ ( k ) max denote the boundary angles of the k -th subregion, with θ ( k ) min &lt; θ ( k ) max . It follows that by changing the number of the subregions (i.e., K ) and the width of each subregion k (i.e., θ ( k ) max -θ ( k ) min ), flexible beam coverage can be achieved.

To achieve uniform beam coverage over the K subregions, we aim to maximize the minimum beam gain over them, i.e.,

$$G _ { \min } ( \omega, x ) = \min _ { \theta \in \mathcal { R } } G ( \omega, x, \theta ), \quad \quad ( 4 ) \quad \text{where} \quad \text{rank} ($$

by jointly optimizing the APV x and the analog beamforming vector ω . The associated optimization problem can be formulated as

$$( \text{Pl} 1 ) \quad & \max _ { \omega, x } G _ { \min } ( \omega, x ) \\ \text{s.t.} \quad & 0 \leq x _ { n } \leq D, \quad n \in \mathcal { N }, \quad & ( \text{5a} ) \\ & x _ { n } - x _ { n - 1 } \geq d _ { \min }, \quad n \in \mathcal { N } \langle \{ 1 \}, \quad & ( \text{5b} ) \\ & | \omega ( n ) | = \frac { 1 } { \sqrt { N } }, \quad n \in \mathcal { N }, \quad & ( \text{5c} ) \quad \text{Note t} \quad & \text{is a c}$$

$$| \omega ( n ) | = \frac { 1 } { \sqrt { N } }, \ \ n \in \mathcal { N },$$

where constraint (5a) ensures that the MAs are located within the transmit region C t ; (5b) ensures the minimum distance between adjacent MAs (denoted as d min ) to avoid antenna coupling; constraint (5c) is with respect to (w.r.t.) the nonconvex modulus constraints on the analog beamforming vector.

However, (P1) is generally difficult to be optimally solved due to the continuous coverage region, the coupling of the analog beamformer ω and APV x in objective function, as well as the unit-modulus constraints. Moreover, the objective function appears to be non-convex in both ω and x . Next, we propose an AO algorithm to solve (P1) iteratively.

## III. PROPOSED ALGORITHM

First, due to the continuous nature of the angle within each subregion R k , we first discretize it into L k discrete values, i.e.,

$$\theta _ { l } ^ { ( k ) } = \theta _ { \min } ^ { ( k ) } + \frac { l ^ { ( k ) } - 1 } { L ^ { ( k ) } - 1 } Z _ { k }, l ^ { ( k ) } = 1, \cdots, L ^ { ( k ) }, \quad ( 6 )$$

where L ( k ) represent the number of the sampling points in R k , and Z k = θ ( k ) max -θ ( k ) min represents the width of R k .

Next, we introduce an auxiliary variable t to convert (P1) into its epigraph form, i.e.,

$$( P 2 ) \quad & \max _ { \omega, x, t } t \\ \text{s.t. } \ G ( \omega, x, \theta _ { l } ^ { ( k ) } ) \geq t, \forall k, l \quad \quad \quad ( 7 a ) \\ \text{(5a), (5b), (5c). }$$

However, (P2) is still difficult to solve due to its non-convex constraints w.r.t. both w and x . Next we propose an AO algorithm to optimize them alternately with the SCA technique.

## A. Optimization of ω for Given x

First, we optimize the analog beamforming vector ω in (P2) with a given APV x . Then, (P1) can be simplified as

$$( P 2. 1 ) \quad \max t \quad \omega, t \quad \\ s. t. \quad ( 7 a ), ( 5 c ).$$

To deal with the non-convex constraint (7a) , we first rewrite the beam gain at θ ( k ) l as

$$G ( \omega, x, \theta _ { l } ^ { ( k ) } ) & = \left | \omega ^ { H } a ( x, \theta _ { l } ^ { ( k ) } ) a ^ { H } ( x, \theta _ { l } ^ { ( k ) } ) \omega \right | ^ { 2 } \\ & = \text{Tr} ( R _ { l } ^ { ( k ) } V ),$$

where R ( k ) l = a x ( , θ ( k ) l ) a H ( x , θ ( ) l k ) and V = ωω H , with rank ( V ) = 1 . If we apply the semidefinite relaxation (SDR) to relax the rank-one constraint, (P2.1) can be recast as

$$\begin{matrix} ( P 2. 2 ) & \max t \\ _ { t, V } \\ s. t. & \text{Tr} ( R _ { l } ^ { ( k ) } V ) \geq t, & \forall l, k \\ V ( n, n ) = \frac { 1 } { N }, & n \in \{ 1, \cdots, N \}. \\ V \succeq 0, & ( 9 c ) \\ \intertext { a a n d } \frac { 1 } { 2 } \sum \partial _ { 0 } \sqrt { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac{ 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac. \, \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 3 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 1 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1. \, \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1. \, 0 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1. \to 0 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } \varphi \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 r } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } {2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { z } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { r } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \sum \frac { 1 } { 2 } \$$

Note that as (9a), (9b) and (9c) are all affine functions, (P2.2) is a classical semi-definite program (SDP) problem, which thus can be optimally solved via the interior-point algorithm [13]. However, the optimized matrix V may not be rank-one. Although Gaussian randomization procedures may be applied to yield a rank-one solution, the resulting performance gap with the globally optimal solution may be large. To tackle this issue, we propose a penalty-based algorithm in this paper, which moves the rank-one constraint to the objective function. Specifically, note that the rank-one constraint can be equivalently written as

$$\text{rank} ( V ) = 1 \Leftrightarrow f ( V ) \stackrel { \triangle } { = } \text{Tr} ( V ) - \sigma ( V ) = 0, \quad ( 1 0 ) \quad \text{the cos} \quad.$$

where σ ( V ) represents the maximum singular value of V . As such, we can consider the following objective function by adding a penalty term, i.e.,

$$\max _ { t, V } t - \rho f ( V ), \quad \quad \quad ( 1 1 ) \ \ A s \, \infty _ { \ f u n c }$$

̸

where ρ ≥ 0 represents a penalty parameter ensuring that the objective function is small enough if Tr ( V ) -σ max = 0 . However, such an objective function is still non-convex as the maximum singular value of V , σ ( V ) , is a convex function (instead of being concave) in V . Nonetheless, it enables us to apply the SCA algorithm to obtain a locally optimal solution. In particular, for a given local point V ( ) i in the i -th SCA iteration, we replace f ( V ) as its first-order Taylor expansion, i.e.,

$$f ( V ) \geq \tilde { f } ( V | V ^ { ( i ) } ) & \stackrel { \triangle } { = } \text{Tr} ( V ) - \sigma ( V ^ { ( i ) } ) + \quad \ \ G ( \omega, \\ & \text{Re} \{ \text{Tr} ( s s ^ { H } ) ( V - V ^ { ( i ) } ) \}, \quad ( 1 2 )$$

where s represents the singular vector corresponding to the maximum singular value of V ( ) i . where the equality holds only when V = V ( ) i . Based on the above, the optimization of V in the i -th SCA iteration is given by

$$( P 2. 3 ) \quad & \max _ { t, V } t - \rho \tilde { f } ( V | V ^ { ( i ) } ) \\ & \text{s.t. } ( 9 a ), ( 9 b ), ( 9 c ).$$

It can be seen that the objective function of (P2.3) is currently a linear function of t and V . As such, (P2.3) is an SDP problem which can be optimally solved [17]. Based on our simulation results, by properly setting the value of ρ , the converged solution of V is always rank-one. Next, we proceed to the ( i + 1) -th SCA for V by updating V ( +1) i = V ( ) i . Finally, to extract ω based on the converged V , we can perform the singular value decomposition (SVD) on V as V = U H 1 Λ U 2 , where U 1 , Λ and U 2 are the left eigenvector matrix, the diagonal matrix of the singular values and the right eigenvector matrix of V , respectively. Then we can obtain the analog beamforming solution to (P2.1) as ω = 1 / √ Ne j arg( U H 1 Λ 1 / 2 r ) , where r is the left singular vector corresponding to the largest singular value of V .

## B. Optimization of x for Given ω

Next, we optimize the APV x in (P2) with a given analog beamformer ω . To tackle the non-convex constraint (7a), we first expand G ( ω x , , θ ( k ) l ) into the following form

$$G & ( \omega, x, \theta _ { l } ^ { ( k ) } ) = \left | \omega ^ { H } a ( x, \theta _ { l } ^ { ( k ) } ) \right | ^ { 2 } & \quad \ \ & A _ { l } ^ { ( k ) } \\ & = \sum _ { p = 1 } ^ { N } \sum _ { q = 1 } ^ { N } \frac { 1 } { N } \cos ( \frac { 2 \pi } { \lambda } \cos ( \theta _ { l } ^ { ( k ) } ) ( x _ { p } - x _ { q } ) - ( \phi _ { p } - \phi _ { q } ) ) \quad \ & \text{Furl} \\ & \quad \text{to $C$}$$

$$\underset { \ v e } { e } \underset { p = 1 } { \ e } \sum _ { q = 1 } ^ { N } \sum _ { q = 1 } ^ { N } \frac { 1 } { N } \cos ( u _ { l } ^ { ( k ) } ( x _ { p }, x _ { q } ) ).$$

where cos( u ( k ) l ( x , x p q )) ≜ ( α ( k ) l ) 2 ( x p -x q ) -( ϕ p -ϕ q ) and α ( k ) l = 2 π λ cos( θ ( k ) l ) . The difficulty in dealing with (14) lies in the cosine function therein. To tackle this challenge, we relax (13) by using the SCA again. For a given z 0 ∈ R , its secondorder Taylor expansion can be expressed as

$$\mathrm b y \quad & \cos ( z _ { 0 } ) - \sin ( z _ { 0 } ) ( z - z _ { 0 } ) - \frac { 1 } { 2 } \cos ( z _ { 0 } ) ( z - z _ { 0 } ) ^ { 2 }. \quad ( 1 4 ) \\..$$

As cos( z 0 ) ≤ 1 and ( z -z 0 ) 2 ≥ 0 , we can construct a quadratic function g z z ( | 0 ) to replace cos( z ) as

$$\mathfrak { n } _ { \mathfrak { x } } \neq 0. \quad \cos ( z ) \geq g ( z | z _ { 0 } ) \stackrel { \Delta } { = } \cos ( z _ { 0 } ) - \sin ( z _ { 0 } ) ( z - z _ { 0 } ) - \frac { 1 } { 2 } ( z - z _ { 0 } ) ^ { 2 }. \\ \text{function}$$

Based on the above, for any given x ( ) i in the i th iteration of SCA for the APV x , by substituting (15) into (13), we can construct a surrogate function of G ( ω x , , θ k ) as

$$( 1 2 ) & \quad \ \ G ( \omega, x, \theta _ { l } ^ { ( k ) } ) \geq \sum _ { p = 1 } ^ { N } \sum _ { q = 1 } ^ { N } \frac { 1 } { N } g ( u _ { l } ^ { ( k ) } ( x _ { p }, x _ { q } ) | u _ { l } ^ { ( k ) } ( x _ { p } ^ { ( i ) }, x _ { q } ^ { ( i ) } ) ) \\ & \quad \stackrel { \triangle } { = } x ^ { T } A _ { l } ^ { ( k ) } x + [ b _ { l } ^ { ( k ) } ] ^ { T } x + c _ { l } ^ { ( k ) }, \quad ( 1 6 ) \\ \text{other} \quad & \quad \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where A ( k ) l ∈ R N × N , b ( k ) l ∈ R N × 1 , and c ( k ) l ∈ R are given by

$$A _ { l } ^ { ( k ) } = - ( \alpha _ { l } ^ { ( k ) } ) ^ { 2 } ( I _ { N } - \frac { 1 } { N } \mathbf 1 _ { N } ) \stackrel { \triangle } { = } - ( \alpha _ { l } ^ { ( k ) } ) ^ { 2 } W, \quad ( 1 7 )$$

$$b _ { l } ^ { ( k ) } ( p ) = \frac { 2 \alpha _ { l } ^ { ( k ) } } { N } \sum _ { q = 1 } ^ { N } [ \alpha _ { l } ^ { ( k ) } x _ { p, q } ^ { ( i ) } - \sin ( u _ { l } ^ { ( k ) } ( x _ { p } ^ { ( i ) }, x _ { q } ^ { ( i ) } ) ) ], \ \ ( 1 8 )$$

$$\text{ly a} \\ \text{lem} \\ \text{tion} \\ \text{t} \right ) \text{-th} \\ \text{track} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{\text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text` \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \text{-th} \\ \text{t} \right ) \end{cases}$$

where I N and 1 N represent the identity matrix and all-ones matrix of size N , respectively, x ( ) i p,q ≜ x ( ) i p -x ( ) i q , and W ≜ I -1 N 1 N . As such, constraint (7a) can be relaxed as

$$x ^ { T } A _ { l } ^ { ( k ) } x + [ b _ { k } ^ { ( l ) } ] ^ { T } x + c _ { l } ^ { ( k ) } \geq t, \forall k, l. \quad \ \ ( 2 0 )$$

Next, we show that (21) is a quadratic constraint (QC). To this end, it is noted that both A ( k ) l and W are hermitian matrices for any given k and l and have the same eigenvectors. Thus, A ( k ) l can be diagonalized as

$$\underset { \mathbb { I } \mathbb { I } } { \text{nalog} } \quad A _ { l } ^ { ( k ) } = [ U _ { l } ^ { ( k ) } ] ^ { H } D _ { l } ^ { ( k ) } U _ { l } ^ { ( k ) } \overset { \triangle } { = } - [ U _ { l } ^ { ( k ) } ] ^ { H } ( \alpha _ { l } ^ { ( k ) } ) ^ { 2 } D U _ { l } ^ { ( k ) } \ \ ( 2 1 )$$

where U ( k ) l is an unitary matrix composed of the eigenvectors of A ( k ) l , D ( k ) l and D are diagonal matrices composed of the eigenvalues of A ( k ) l and W , respectively, with D ( k ) l = -( α ( k ) l ) 2 D . Furthermore, any eigenvalue of D is non-negative according to Gerschgorin disk theorem. Hence, the eigenvalue of D ( k ) l is non-positive and A ( k ) l is a negative semi-definite matrix. It follows that constraint (20) is a convex QC w.r.t. x .

Therefore, in the i -th iteration of the SCA for x x , can be optimized by solving the following optimization problem

$$( P 2. 4 ) & \max _ { x, t } t \\ & s. t. & ( 5 a ), ( 5 b ), ( 2 0 ).$$

Since constraints (5a) and (5b) are linear and (20) is a QC w.r.t x , (P2.2) is a classic quadratically constrained quadratic program (QCQP) problem which can be effectively solved via the interior-point algorithm. Next, we proceed to the ( i +1) -th SCA for x by updating x ( +1) i = x ( ) i .

## C. Overall Algorithm and Complexity Analysis

The overall AO algorithm can be executed as follows. Let ω ( j -1) and x ( j -1) denote the values of ω and x at the beginning of the j -th AO iteration. Then, in this AO iteration, we first optimize V via SCA with fixing x = x ( j -1) and V (0) = ω ( j -1) ω H ( j -1) and obtain ω ( j ) by performing SVD on the converged V . Next, we optimize x via SCA with fixing ω = ω ( j ) and x (0) = x ( j -1) and obtain x ( j ) as the converged solution, and the ( j + 1) -th AO iteration follows. The overall procedures of the AO algorithm are summarized in Algorithm 1.

## Algorithm 1 Proposed AO algorithm to solve (P1)

Input: ω (0) and x (0) .

Output: ω x , .

- 1: Initialization: j ← 1 .
- 2: while AO convergence is not reached do
- 3: Initialize i ← 0 and update V (0) = ω ( j -1) ω H ( j -1) , and x = x ( j -1) .
- 4: while SCA convergence for V is not reached do
- 5: Obtain V ( +1) i by solving problem (P2.3).
- 6: Update i ← i +1 .
- 7: end while
- 8: Obtain ω ( j ) based on the SVD on V ( ) i

.

- 9: Initialize i ← 0 and update ω = ω ( j ) , x (0) = x ( j -1) .
- 10: while SCA convergence for x is not reached do
- 11: Obtain x ( +1) i by solving problem (P2.4).
- 12: Update i ← i +1 .
- 13: end while
- 14: Update x ( j ) = x ( ) i .
- 15: j ← j +1 .
- 16: end while
- 17: return x and ω .

Next, we prove that the proposed AO algorithm with SCA is ensured to converge. Specifically, regarding the SCA for V , let v ( ) i ≜ t -ρf ( V ( ) i ) denote the objective value of (11) in the i -th SCA iteration for V . Then, the following inequalities hold, i.e.,

$$( 2 3 )$$

$$v ^ { ( i ) } & \stackrel { \triangle } { = } t - \rho f ( V ^ { ( i ) } ) \stackrel { ( a ) } { = } t - \rho \tilde { f } ( V ^ { ( i ) } | V ^ { ( i ) } ) \quad ( 2 3 ) \quad \text{ this } 1 \\ & \stackrel { ( b ) } { \leq } t - \rho \tilde { f } ( V ^ { ( i + 1 ) } | V ^ { ( i ) } ) \stackrel { ( c ) } { \leq } t - \rho f ( V ^ { ( i + 1 ) } ) \stackrel { \triangle } { = } v ^ { ( i + 1 ) }, \quad \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        } \text{                                                                                                                                                                                                    
   the  \Sigma$$

where the equality ( a ) holds since the first-order Taylor expansion in (12) is tight at V ( ) i ; the inequality ( b ) holds as V ( +1) i is the optimal solution to (P2.3) and thus maximizes the function t -ρf ˜ ( V V | ( ) i ) ; the inequality ( c ) holds due to (12). Based on the above, the sequence { v ( ) i } is non-decreasing and thus ensured to converge. Similarly, it can be shown that the SCA for optimizing the APV x also converges, for which the details are omitted for brevity. It follows that the proposed AO algorithm in Algorithm 1 must converge.

Finally, we analyze the complexity Algorithm 1. It can be shown that the complexity of optimizing the analog beamformer ω and the APV x are both in the order of O ( √ LN N ( 2 + L )) [18], where L = ∑ K k =1 L ( k ) denotes the total number of sampling points.

## D. Initialization

The converged performance of the AO algorithm depends critically on the initialization. In this subsection, we propose an efficient initialization method to detemine x (0) and ω (0) in Algorithm 1. For the initial MA positions, we consider that all MAs are uniformly deployed within the linear array with equal spacing. As such, we set x (0) = [ D N +1 , 2 D N +1 , · · · , ND T N +1 ] .

Next, to determine ω (0) , we consider solving the SDR problem (P2.2) directly with fixing x = x (0) via the interiorpoint algorithm. Denote by V (0) the associated optimal solution to (P2.2). Then, we obtain ω (0) by performing Gaussian randomization based on V (0) , i.e.,

$$\omega ( 0 ) = \frac { 1 } { \sqrt { N } } e ^ { j \arg ( U _ { 0 } ^ { H } \Lambda _ { 0 } ^ { 1 / 2 } r _ { 0 } ) }, \text{ \quad \ \ } ( 2 4 )$$

where U 0 and Λ 0 are unitary matrix and diagonal matrix composed by eigenvectors and eigenvalues of the matrix V (0) , respectively, and vector r 0 ∈ C N × 1 is a random vector following CN (0 , I N ) . By comparing the objective value of (P2.2) under different realizations of r 0 , we set ω (0) as that corresponding to the one yielding the maximum objective value of (P2.2).

## IV. NUMERICAL RESULTS

In this section, we provide numerical results to evaluate the performance of our proposed MA-enabled beam coverage. Unless otherwise stated, the simulation parameters are set as follows. The penalty parameter for updating AWV is set to ρ = 20 . The carrier frequency is f c = 1 GHz. The length of the MA linear array is set to D = 8 λ . Meanwhile, we set the convergence threshold for the AO algorithm, the SCA for V , and the SCA for x as 10 -5 , 0 01 . , and 0 01 . , respectively. For performance comparison, we consider the following two benchmarks:

- 1) FPA : N FPAs are deployed with half-wavelength spacing and the analog beamfroming is optimized based on the SCA algorithm as in Section III-A.
- 2) MAs with fixed analog beamforming (MA-FAB) : In this benchmark, we set the analog beamforming ω = ω (0) as presented in Section III-D and optimize MA positions based on the SCA algorithm as in Section III-B.

FPA Fig. 2. Optimized beam gains with (a) N = 6 and (b) N = 8 , and optimized MA positions with (c) N = 6 and (d) N = 8 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00434v2-figure-dc26be9d-fcf4-418a-979c-8a6cb5f6f8b7.jpg)

0 0.2 0.4 0.6 0.8 1 1.2 MA-FAB Proposed First, by setting K = 1 , θ (1) min = 0 , and θ (1) max = 180 ◦ , i.e., R = [0 180 ] , ◦ , Figs. 2(a) and 2(b) plot the optimized beam gains over R by different schemes for N = 6 and N = 8 , respectively. Note that we have introduced a common shift to the MA positions such that x 1 = 0 without affecting the coverage performance. The optimized MA positions by different schemes are also shown in Figs. 2(c) and 2(d) for N = 6 and N = 8 , respectively. It is observed from Figs. 2(a) and 2(b) that our proposed AO algorithm can achieve a considerably higher max-min beam gain over R as compared to the two benchmarks, thereby resulting in flatter beam gain for both N = 6 and N = 8 . It is also observed that the MA-FAB benchmark outperforms the FPA benchmark, which indicates that antenna positions may play a more significant role than antenna weights in terms of beam coverage. In addition, it is observed from Figs. 2(c) and 2(d) that the optimized MA positions are different in different schemes. In particular, for both N = 6 and N = 8 , the proposed scheme and the MAFAB benchmark result in a larger antenna aperture, i.e., x N -x 1 , than the FPA benchmark, so as to fully exploit the spatial DoFs available. Moreover, the inter-antenna spacing is observed to be non-uniform in the former two schemes employing MAs, unlike the FPA benchmark with an equal antenna spacing. In addition, by comparing Figs. 2(c) and 2(d), it is observed that the antenna aperture under M = 8 is larger than that under M = 6 in the proposed scheme (1.2 versus 1.14), so as to accommodate the movement of more MAs.

The optimized MA positions by the considered schemes are shown in Fig. 3(b). It is observed that in the case of multi-region coverage, the proposed scheme still achieves much higher maxmin beam gain than the other two benchmarks. In particular, the minimum beam gain by our proposed scheme is observed to be even larger than the maximum beam gain by the FPA benchmark in R 1 and R 3 . While in R 2 , although the FPA can achieve a maximum beam gain of 6 dB, its minimum beam gain is considerably lower ( -19 dB) than that by the proposed scheme ( -1 dB). Moreover, the MA-FAB benchmark is observed to achieve a slightly higher minimum beam gain than the proposed scheme in R 2 . This implies that as far as each subregion concerned, the proposed scheme may not achieve the highest minimum beam gain. It is also observed from Fig. 3(b) that similar to Figs. 2(c) and 2(d), the proposed scheme and the MA-FAB benchmark result in a larger antenna aperture than the FPA benchmark. Particularly, x 8 increases from 1.2 in Fig. 2(d) to 1.65 in Fig. 3(b), which suggests that multi-region coverage may demand a larger antenna aperture than its single-region counterpart due to the discontinuous subregions. Interestingly, it is observed that the optimized MA positions by the proposed scheme can be viewed as being composed of three subarrays, with the first and second MAs in the first subarray, the third to fifth MAs in the second subarray, and the sixth to eighth MAs in the third subarray, which helps cover the three subregions.

Next, Fig. 3(a) plots the optimized beam gain over R by setting K = 3 , R 1 = [0 , 30 ] ◦ , R 2 = [70 ◦ , 110 ] ◦ and R 3 = [160 ◦ , 170 ] ◦ . The number of MAs is set to N = 8 .

Finally, by setting K = 1 and θ (1) min = 0 , we plot the max-min beam gain versus θ (1) max in Fig. 4. The number of MAs is set to N = 8 . It is observed that the max-min beam gains by all considered schemes monotonically decrease

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00434v2-figure-710c9ccc-19bf-411e-bc9a-cc16c4f278e8.jpg)

(a)

Fig. 3. (a) Optimized beam gain and (b) optimized MA positions for multiregion coverage.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00434v2-figure-83c7ab1f-0d10-4ad7-bc2a-27884b67b6fc.jpg)

with θ (1) max . This is expected as the total beam gain over the entire spatial domain is a constant. Thus, when θ (1) max or the coverage width increases, the max-min beam gain should reduce accordingly. It is also observed that in the case with a narrowto-moderate coverage width, i.e., θ (1) max ≤ 50 ◦ , the performance gap between the proposed scheme and the MA-FAB benchmark is negligible, while both of them significantly outperform the FPA benchmark. This indicates that adjusting the MA positions suffices to achieve uniform coverage for a narrow region given the proposed initial analog beamforming design. However, with the increase in θ (1) max , the performance gap between the proposed scheme and the MA-FAB benchmark becomes larger. This indicates that the analog beamforming design will play a more significant role for the more challenging wide-region coverage.

## V. CONCLUSION

In this paper, we investigate the beam coverage design with multiple MAs by jointly optimizing their weights and positions to maximize the minimum beam gain over the desired spatial regions. An AO algorithm was proposed to obtain a high-quality suboptimal solution to this problem. Numerical results show that our proposed MA-enabled beam coverage design significantly outperforms the conventional FPAs by more flexibly catering to the geometry of the subregions. It was also shown that the antenna weights may play a more significant role as the width of the region increases.

Fig. 4. Max-min beam gain versus coverage width with K = 1 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00434v2-figure-314d951b-2ca9-4776-a45c-9310655fa479.jpg)

## REFERENCES

- [1] J. Mietzner, R. Schober, L. Lampe, W. H. Gerstacker and P. A. Hoeher, 'Multiple-antenna techniques for wireless communications - a comprehensive literature survey,' IEEE Commun. Surveys Tuts. , vol. 11, no. 2, pp. 87-105, Second Quarter 2009.
- [2] K.-K. Wong, A. Shojaeifard, K.-F. Tong, and Y. Zhang, 'Fluid antenna systems,' IEEE Trans. Wireless Commun. , vol. 20, no. 3, pp. 1950-1962, Mar. 2020.
- [3] K.-K. Wong, K.-F. Tong, and C.-B. Chae, 'Fluid antenna system-part II: Research opportunities,' IEEE Commun. Lett. , vol. 27, no. 8, pp. 19241928, Aug. 2023.
- [4] B. Ning et all., 'Movable Antenna-Enhanced Wireless Communications: General Architectures and Implementation Methods,' https://arxiv.org/pdf/2407.15448 , 2024.
- [5] L. Zhu, W. Ma, and R. Zhang, 'Modeling and performance analysis for movable antenna enabled wireless communications,' IEEE Trans. Wireless Commun. , early access, 2023.
- [6] L. Zhu, W. Ma, and R. Zhang, 'Movable antennas for wireless communication: Opportunities and challenges,' IEEE Commun. Mag. , early access, 2023.
- [7] Y. Wu, D. Xu, D. W. K. Ng, W. Gerstacker, and R. Schober, 'Movable antenna-enhanced multiuser communication: Jointly optimal discrete antenna positioning and beamforming,' in Proc. IEEE Global Commun. Conf. , Kuala Lumpur, Malaysia, Dec. 2023, pp. 7508-7513.
- [8] L. Zhu, W. Ma, B. Ning, and R. Zhang, 'Movable-antenna enhanced multiuser communication via antenna position optimization,' IEEE Trans. Wireless Commun. , 2023, early access.
- [9] W. Mei, X. Wei, B. Ning, Z. Chen and R. Zhang, 'Movable-antenna position optimization: A graph-based approach,' IEEE Wireless Commun. Lett. , vol. 13, no. 7, pp. 1853-1857, July 2024.
- [10] X. Wei, W. Mei, D. Wang, B. Ning and Z. Chen, 'Joint Beamforming and Antenna Position Optimization for Movable Antenna-Assisted Spectrum Sharing,' IEEE Wireless Commun. Lett. , early acess
- [11] G. Hu, Q. Wu, K. Xu, J. Si and N. Al-Dhahir, 'Secure wireless communication via movable-antenna array,' IEEE Signal Processing Letters , vol. 31, pp. 516-520, 2024.
- [12] Z. Cheng, N. Li, J. Zhu, X. She, C. Ouyang and P. Chen, 'Enabling secure wireless communications via movable antennas,' in Proc. ICASSP , Seoul, Korea, 2024, pp. 9186-9190.
- [13] W. Ma, L. Zhu, and R. Zhang, 'MIMO capacity characterization for movable antenna systems,' IEEE Trans. Wireless Commun. , vol. 23, no. 4, pp. 3392-3407, April 2024.
- [14] W. Ma, L. Zhu and R. Zhang, 'Multi-beam forming with movable-antenna array,' IEEE Commun. Lett. , vol. 28, no. 3, pp. 697-701, Mar. 2024.
- [15] L. Zhu, W. Ma and R. Zhang, 'Movable-antenna array enhanced beamforming: Achieving full array gain with null steering,' IEEE Commun. Lett. , vol. 27, no. 12, pp. 3340-3344, Dec. 2023
- [16] B. Ning, T. Wang, C. Huang, Y. Zhang, and Z. Chen, 'Widebeam designs for terahertz massive MIMO: SCA-ATP and S-SARV,' IEEE Internet Things J. , no. 12, pp. 10857-10869, June 2023.
- [17] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge, U.K.: Cambridge Univ. Press, 2009.
- [18] K. Wang, et al. , 'Outage constrained robust transmit optimization for multiuser MISO downlinks: Tractable approximations by conic optimization,' IEEE Trans. Signal Process. , vol. 62, no. 21, pp. 5690-5705, Nov. 2014.
Dong Wang Weidong Mei Boyu Ning Zhi Chen
Flexible Beam Coverage Optimization for Movable-Antenna Array
Fluid antennas (FAs) and movable antennas (MAs) have attracted increasing
attention in wireless communications recently. As compared to the conventional
fixed-position antennas (FPAs), their geometry can be dynamically reconfigured,
such that more flexible beamforming can be achieved for signal coverage and/or
interference nulling. In this paper, we investigate the use of MAs to achieve
uniform coverage for multiple regions with arbitrary number and width in the
spatial domain. In particular, we aim to jointly optimize the MAs weights and
positions within a linear array to maximize the minimum beam gain over the
desired spatial regions. However, the resulting problem is non-convex and
difficult to be optimally solved. To tackle this difficulty, we propose an
alternating optimization (AO) algorithm to obtain a high-quality suboptimal
solution, where the MAs weights and positions are alternately optimized by
applying successive convex approximation (SCA) technique. Numerical results
show that our proposed MAbased beam coverage scheme can achieve much better
performance than conventional FPAs.