## A Rotation-based Method for Precoding in Gaussian MIMOME Channels

Xinliang Zhang, Yue Qi, Student Member, IEEE , and Mojtaba Vaezi, Senior Member, IEEE

Abstract -The problem of maximizing secrecy rate of multipleinput multiple-output multiple-eavesdropper (MIMOME) channels with arbitrary numbers of antennas at each node is studied in this paper. First, the optimization problem corresponding to the secrecy capacity of the MIMOME channel is converted to an equivalent optimization based on Givens rotations and eigenvalue decomposition of the covariance matrix. In this new formulation, precoder is a rotation matrix which results in a positive semi-definite (PSD) covariance matrix by construction. This removes the PSD matrix constraint and makes the problem easier to tackle. Next, a Broyden-Fletcher-Goldfarb-Shanno (BFGS)-based algorithm is developed to find the rotation and power allocation parameters. Further, the generalized singular value decomposition (GSVD)-based precoding is used to initialize this algorithm. The proposed rotation-BFGS method provides an efficient approach to find a near-optimal transmit strategy for the MIMOME channel and outperforms various state-of-the-art analytical and numerical methods. In particular, the rotationBFGS precoding achieves higher secrecy rates than the celebrated GSVD precoding, with a reasonably higher computational complexity. Extensive numerical results elaborate on the effectiveness of the rotation-BFGS precoding. The new framework developed in this paper can be applied to a variety of similar problems in the context of multi-antenna channels with and without secrecy.

Index Terms -Physical layer security, MIMO wiretap channel, secrecy capacity, beamforming, precoding, covariance, rotation.

## I. INTRODUCTION

As a complement to higher-layer security measures, physical layer security has emerged as a significant technique for security in the lowest layer of communication, i.e., the physical layer. Founded on information-theoretic security, which is built on classical Shannon's notion of perfect secrecy, physical layer security can offer unbreakable security, unlike conventional secret-key-based cryptosystems. Physical layer security was laid in the 1970s in Wyner's seminal work on the wiretap channel [2] where the idea of secure communication based on the communication channel itself without using encryption keys was first introduced. In this work, Wyner proved that in a wiretap channel (a channel in which a transmitter conveys information to a legitimate receiver in the presence of an eavesdropper) communication can be both robust to transmission errors ( reliable ) and confidential ( secure ), to a certain degree, provided that the legitimate user's channel is

This paper was partially presented in the IEEE International Symposium on Personal, Indoor, and Mobile Radio Communications, September 2020 [1].

The authors are with the Department of Electrical and Computer Engineering, Villanova University, Villanova, PA, USA (e-mail: { xzhang4, yqi, mvaezi } @villanova.edu).

better than the eavesdropper's channel 1 . He established the capacity of the degraded wiretap channel. Later, Csiszar and Korner [4] generalized this result to arbitrary, not necessarily degraded, wiretap channels.

In the past decades, physical layer security has been applied to enhance the classical wiretap channels (e.g., by including more realistic assumptions) and to study advanced wiretap channels (e.g., quantum communication [5], [6]). Particularly, as multiple-input multiple-output (MIMO) networks continue to flourish worldwide, a significant effort has been made to study the MIMO wiretap channel which allows for the exploitation of space/time/user dimensions of wireless channels for secure communications. Specifically, secrecy capacity of Gaussian multiple-input multiple-output multipleeavesdropper (MIMOME) channels under an average total power constraint was established independently in [7]-[9]. The capacity result is abstracted as an optimization problem over input covariance matrix. This problem is non-convex and its optimal solution is known only for limited settings [10]-[14].

Among notable sub-optimal solutions that can be applied to the MIMOME channel is the generalized singular value decomposition (GSVD)-based precoding [15]. GSVD-based precoding decomposes transmitted channel matrices into several parallel subchannels and confidential information is transmitted over subchannels where the legitimate user is stronger than the eavesdropper. This method gives a closed-form solution for achievable secrecy rate which is relatively fast and is asymptotically optimal at high signal-to-noise ratios (SNRs). However, its performance is not good at certain settings, e.g., when the eavesdropper has a single antenna while other nodes have multiple antennas [14]. Another important suboptimal solution is Li et al. 's alternating optimization and water filling (AOWF) algorithm [16] which alternates the original optimization problem to a convex problem and finds the corresponding Lagrange multipliers in an iterative manner. AOWF is more computationally expensive than GSVD-based precoding but it can provide a better secrecy rate in some settings. The performance of this method also varies depending on the number of antennas at different nodes. For example, its performance is not as good as the GSVD-based precoding when the number of antennas at the eavesdropper is greater than that of the transmitter. There are also other numerical solutions for this optimization problem [17]-[19]. Specifically, in [17], a barrier method based iterative algorithm is tailored to obtain global optimal with guaranteed convergence. However,

1 Later in the 1990s, Maurer proved that secret key generation through public communication over an insecure yet authenticated channel is possible even when a legitimate user has a worse channel than an eavesdropper [3].

it is not straightforward to obtain the barrier parameter and this involves a challenging optimization problem. In addition, this solution is devised based on an upper bound custom-made for the MIMOME channel and its extension to other related problems is not straightforward.

Recently, based on a trigonometric parameterization of the covariance matrix, a closed-form solution for optimal precoding and power allocation of the MIMOME channel with two transmit antennas was obtained in [14], [20]. This approach in finding the optimal covariance matrix is completely different from existing linear beamforming methods. It does not require degradedness condition of [12] and [21], and thus provides the optimal solution for both full-rank and rank-deficient cases in one shot. The above beamforming and power allocation schemes are, however, limited to two transmit antenna cases, and the optimal transmit covariance matrix is still open in general. Givens rotation has been previously investigated for quantizer design in MIMO broadcast channels [22]. Later, Givens rotation is used as a beamformer to maximize sumrate under average signal-to-noise ratio (SNR) with limited feedback in multiuser multiple-input single-output (MISO) channels [23]. However, Givens rotation theory is the first time to be applied on MIMO wiretap channel to the authors' knowledge.

In this paper, we use Givens rotation matrices to generalize the approach of [20] to arbitrary numbers of antennas at each node, and introduce a new method for precoding and power allocation in the MIMOME channel. In this approach, without loss of generality, the precoding matrix is formed using a rotation matrix. Then, the covariance matrix can be seen as an operator to appropriately stretch (by power allocation) and rotate (by precoding) the input symbols to form a transmit signal that best fits the channels of the legitimate user and eavesdropper. The capacity expression is then transformed into optimizing rotation angles and power allocation parameters. This problem is very different from the original problem, which requires exploring new techniques to solve it. One advantage of the new problem formulation is that it converts the matrix covariance constraint (symmetric and positive semidefiniteness (PSD)) into a set of linear constraints and thus simplifies the optimization problem. We then provide a numerical solution for the new optimization problem and show that it outperforms the existing methods in various antenna settings. The method in [14] can be seen as using two-dimension Givens rotation for precoding.

## A. Contributions

The main contributions of this paper are listed below:

- · We use a rotation modeling method for the parameterization of the covariance matrix. This parameterization gives a new representation of the capacity expression for the MIMOME channel for which more efficient solutions can be exploited. Particularly, the PSD constraint on the covariance matrix is removed because in the proposed method the covariance matrix is PSD by construction.
- · For n t = 1 and n t = 2 , the proposed scheme reduces to those in [11] and [14], respectively, for which closed-

form solutions are known. For n t ≥ 3 , finding a closedform solution is still challenging. In such a case, we introduce a Broyden-Fletcher-Goldfarb-Shanno (BFGS)based method to iteratively solve the rotation and power allocation parameters. We name this approach rotationBFGS method in which a rectifier is designed to remove the constraints. Numerical results in different antenna settings confirm that the proposed scheme works better than the well-known GSVD and AOWF. Specifically, the proposed approach outperforms GSVD when n e &lt; n t , and AOWF approach when n e ≥ n t , where n e is the number of antennas at the eavesdropper. Particularly, the gap between the proposed and GSVD-based methods is remarkably high when the eavesdropper has a single antenna.

- · To improve the computational complexity of the proposed algorithm (by reducing the number of iterations), we develop an algorithm to exploit GSVD as an initialization for our rotation-BFGS method. This initialization improves the results and reduces the computational complexity as it reduces the number of iterations in the optimization problem.

## B. Other Related Works

An interesting aspect of the proposed approach is its generality and its great potential for extension to other related problems. The MIMOME channel has turned out as a fundamental tool for the study of physical-layer security in many other related problems throughout the past decade. Many solutions developed for the MIMOME has appeared to be instrumental in designing transmit strategies that maximize the secrecy rate of extensions of this basic channel model to MIMO channels with multiple eavesdroppers [16], secure relaying [24], [25], ergodic secrecy rate [26], finite alphabet signaling [27]-[29], artificial noise [30] and cooperative jammer [31], among others. The rotation method is also applicable to any problem that can be cast as an optimization problem over a covariance matrix. Therefore, it is worth studying the optimal covariance matrix of the MIMOME channel as a general tool for physical-layer security in various MIMO settings.

## C. Organizations and Notations

The remainder of this paper is organized as follows. Section II describes the system model and related works. Section III introduces and elaborates on the Givens rotation and reformulates the secrecy capacity for the MIMOME channel. Section IV, details a rotation-BFGS based algorithm to optimize the achievable secrecy rate. In Section V, numerous simulation results are carried out to demonstrate the effectiveness of the proposed method. Section VI draws the conclusion.

Notations: Bold lowercase letters denote column vectors and bold uppercase letters denote matrices. | x | and log 2 ( x ) denote the absolute value and the binary logarithm of the scalar x . A ( i, j ) denotes the entry ( i, j ) of matrix A . Besides, ( A ) T , tr( A ) , | A | , and diag( A ) are the transpose, trace, determinate, and diagonal of the matrix A . [ A ] + replaces the negative elements with zeros. I n is a n × n identity matrix. And E {·} is the expectation of random variables.

Fig. 1: The MIMOME channel with n t , n r , and n e antennas at the transmitter, legitimate receiver, and eavesdropper.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-7f9865b0-3d9e-4090-ae14-8d53681ad3f8.jpg)

## II. SYSTEM MODEL AND RELATED WORKS

## A. System Model

We consider a MIMOME channel with n t antennas at the transmitter, n r antennas at the receiver, and n e antennas at the eavesdropper, as depicted in Fig. 1. The transmitter knows the perfect channel state information (CSI) of both users 2 . The received signals at the legitimate receiver and the eavesdropper can be, respectively, expressed as

$$y _ { r } = H x + w _ { r }, \quad \quad \ \ ( 1 a ) ^ { \ \lim i }$$

$$y _ { e } = G x + w _ { e }, \quad \ \ \ ( 1 b ) \quad ^ { \bullet }$$

in which H ∈ R n r × n t and G ∈ R n e × n t are the channels corresponding to the receiver and eavesdropper, x ∈ R n t is the transmitted signal, and w r ∈ R n r and w e ∈ R n e are independent and identically distributed (i.i.d) Gaussian noises with zero means and identity covariance matrices. A representation of secrecy capacity is given by [9]

$$( P 1 ) \ \mathcal { C } _ { s } = \max _ { \substack { Q \succeq 0, \text{tr} ( Q ) \leq P _ { t } } } \frac { 1 } { 2 } \log _ { 2 } \frac { | I _ { n _ { r } } + \text{HQH} ^ { T } | } { | I _ { n _ { e } } + \text{GQG} ^ { T } | }. \quad ( 2 ) \quad \begin{array} { c } \ w | \\ \ r a \\ \bullet \ n _ { t } \end{array} } \quad \bullet \ n _ { t }$$

Based on Slyvester's determinant theorem, i.e., det( + I XY ) = det( I + YX ) , (2) can also be rewritten as

$$\mathcal { C } _ { s } & = \max _ { Q \succeq 0, \text{tr} ( Q ) \leq P _ { t } } \frac { 1 } { 2 } \log _ { 2 } \frac { | \mathbf I _ { n _ { t } } + \mathbf H ^ { T } \mathbf H Q | } { | \mathbf I _ { n _ { t } } + G ^ { T } G Q | }, \quad ( 3 ) \quad \frac { \bullet \, C } { \mathbf a } \\. & \quad \. & \. \quad \. \quad \. \quad \. \quad \.$$

where Q = E { xx T } ∈ R n t × n t is the covariance matrix of the channel input x and P t is the total transmit power. Q is symmetric and positive semi-definite (PSD) by definition.

The architecture of the linear precoding and power allocation is depicted in Fig. 2. In this figure, s , . . . , s 1 n t are input symbols which are independent and identically distributed Gaussian random variables with zero means and unit variances, λ i s are power allocation coefficients, V is the precoding matrix, and x = [ x , . . . , x 1 n t ] T is the transmit vector whose covariance is Q .

2 A perfect CSI is assumed, as we are deriving the theoretical limits. This may provide an upper bound in terms of achievable secrecy rates. The method we are developing in this paper is, however, applicable to the case with imperfect CSI. In our future works, we will relax this idealized assumption and consider the practical scenarios with imperfect CSI.

Fig. 2: The structure of linear precoding and power allocation. s , . . . , s 1 n t are the independent input symbols, λ , . . . , λ 1 n t are their corresponding allocated powers, and x , . . . , x 1 n t are the transmitted signals.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-34a021bc-47c5-413e-98a5-4a70a8020e31.jpg)

## B. Existing Results

The optimal transmission over the MIMOME channel is still an open problem in general. However, there are a number of notable analytical results for special numbers of antennas as well as numerical results as listed below.

- 1) Analytical Solutions: An analytical capacity-achieving covariance matrix is known only for special cases. These are limited to:
- · n t = 1 : this is the single-input multiple-output (SIMO) case in which Q is a scalar and the optimal solution is either P t or 0 [11].
- · n r = 1 : the so-called multiple-input single-output multiple-eavesdropper (MISOME) channel in which generalized eigenvalue decomposition of H and G achieves the capacity [32].
- · n t = 2 , n r = 2 , and n e = 1 : the optimization problem is shown to be the Rayleigh quotient and optimal signaling, which is the maximum eigenvalue of this problem, is unitrank [10].
- · n t = 2 : in which the secrecy capacity is obtained by modeling the covariance matrix as a 2 × 2 rotation matrix [14], [20].
- · Q is full-rank (which implies H H T -G G T glyph[follows] 0 ) and also P t is greater than a certain threshold [12], [21]: in this case the problem is convex and Karush-Kuhn-Tucker (KKT) conditions are used to find the optimal Q . It is worth noting that, as n t grows, few channel realizations satisfy the above conditions. For n t = 3 , n r = 3 , and n e = 1 , for example, the probability of having a fullrank solution is less than 18 2% . 3 . This value decreases when n e goes up. Therefore, it can be said that an analytical solution for the MIMOME channel is still an open problem in many practical cases.

2) Suboptimal Analytical Solution: For a general MIMOME channel, a sub-optimal solution can be obtained using GSVD-based beamforming [15]. By applying GSVD on H and G , the optimization problem (3) is simplified to a set of parallel non-interfering channels whose optimal power

3 This is obtained by Monte Carlo experiments with 10 6 trails where H and G have the same distributions.

allocation can be obtained using Lagrange multiplier and KKT conditions. Since parallelization using GSVD does not necessarily convert this problem into an equivalent one, GSVDbased beamforming is not the optimal solution in general. It is not even close enough to the capacity in some cases. For example, GSVD-based beamforming achieves less than 70% of the capacity when n t = 3 , n r = 2 , and n e = 1 .

3) Numerical Solutions: There are still important cases of the MIMOME for which optimal Q is unknown. Due to the intractability of the problem in an analytical form, numerical solutions have been developed to tackle this problem. The AOWF [16], which is computationally efficient to implement, is one of them. Despite its effectiveness in many cases, AOWF experiences problems when n e is greater than n t , for example, which is caused by a failure in finding an optimal Lagrange multiplier. We modify this issue in this paper, as we will see later in Section V. The price is a higher time consumption in the modified approach.

In the next section, we apply a rotation-based model for the covariance matrix Q which is a generalization of the solution in [14], [16], from n t = 2 to any arbitrary n t . This model is then used to find transmit signaling that can be used to achieve secrecy capacity of the MIMOME channel regardless of the number of antennas at different nodes.

## III. A ROTATION MODELING OF THE PROBLEM

In the following subsections, we further model matrix V using the rotation matrix related method by reviewing n t = 2 [20] first and then generalizing it to an arbitrary n t with proof.

The covariance matrix Q can be eigendecomposed as

$$Q = V \Lambda V ^ { T },$$

in which Λ is a diagonal matrix and its diagonal elements are the eigenvalues of Q , which are real and non-negative, i.e.,

$$\mathbf \Lambda = \text{diag} ( \lambda _ { i } ), \, \lambda _ { i } \geq 0, \, i = 1, 2, \dots, n _ { t }, \quad ( 5 )$$

and the total power constraint tr( Q ) ≤ P t will be equivalent to

$$\sum _ { i = 1 } ^ { n _ { t } } \lambda _ { i } & \leq P _ { t }. & ( 6 ) \\ \text{is the matrix consumed of $n$ corresponding}$$

Also, V ∈ R n t × n t is the matrix composed of n t corresponding eigenvectors of Q in R n t vector space. Since Q is symmetric, the matrix V is orthonormal .

With this, an immediate implication of the above decomposition is that linear precoding can achieve the capacity of the MIMOME channel.

## A. Rotation Modeling for n t = 2

The capacity region of MIMOME channel with two transmit antennas has been established in [14], [20], using a rotation matrix in two-dimensional (2D) space. In this case, the eigenvalue matrix Λ can be written as

$$\begin{matrix} \text{can be written as } & \quad \text{$\bullet$} \\ \Lambda = \begin{bmatrix} \lambda _ { 1 } & 0 \\ 0 & \lambda _ { 2 } \end{bmatrix}, & \quad & ( 7 ) & \text{and} \\ & \text{order} \end{matrix}$$

and, without loss of generality, the eigenvectors can be written as the following orthonormal (rotation) matrix

The rotation angle θ 12 in rotation matrix V 12 corresponds to the rotation from the direction of the standard basis (unit vector) e 1 = (1 0) , T to the standard basis e 2 = (0 1) , T in R 2 vector space. This rotation is achieved on the plane defined by e 1 and e 2 . Then, the covariance matrix Q can be built using three parameters: the two non-negative eigenvalues λ 1 and λ 2 , and the rotation angle θ 12 . These parameters are complete to represent any arbitrary 2 × 2 covariance matrix as in (4). Then the original optimization problem (3) can be equivalently converted to

$$e \text{ following orthonormal (rotation) matrix} \\ V = V _ { 1 2 } \stackrel { \triangle } { = } \begin{bmatrix} \cos \theta _ { 1 2 } & - \sin \theta _ { 1 2 } \\ \sin \theta _ { 1 2 } & \cos \theta _ { 1 2 } \end{bmatrix}. \quad \text{ (8)} \\ \text{rotation angle} \, \theta _ { 1 2 } \text{ in rotation matrix} \, V _ { 1 2 } \text{ corresponds}$$

$$\mathcal { C } _ { s } = \max _ { \lambda _ { 1 }, \lambda _ { 2 }, \theta _ { 1 2 } } \frac { 1 } { 2 } \log _ { 2 } \frac { | I _ { n _ { t } } + \mathbf H ^ { T } \mathbf H \mathbf V \mathbf A \mathbf V ^ { T } | } { | I _ { n _ { t } } + G ^ { T } \mathbf G \mathbf V \mathbf A \mathbf V ^ { T } | }, \\ \text{s.t.} \quad \lambda _ { 1 } \geq 0, \ \lambda _ { 2 } \geq 0, \ \lambda _ { 1 } + \lambda _ { 2 } \leq P _ { t }. \quad ( 9 b )$$

$$\ s. t. \ \ \lambda _ { 1 } \geq 0, \ \lambda _ {$$

In light of this modeling, an analytical solution for optimal precoding matrix and power allocation scheme are obtained in [20] by finding θ 12 , λ 1 and λ 2 . In the next subsection, we extend this method to the cases for an arbitrary n t .

## B. Generalization to an Arbitrary n t

To generalize the rotation modeling method to an arbitrary n t × n t covariance matrix, Λ ∈ R n t × n t is a diagonal matrix with non-negative elements Λ ( i, i ) glyph[defines] λ i . V is a rotation matrix in R n t × n t vector space which can be obtained by

$$V = \prod _ { i = 1 } ^ { n _ { t } - 1 } \pro$$

in which the basic rotation matrix V ij is the Givens matrix [33] defined as and

$$\begin{array} { l l$$

$$\begin{bmatrix} v _ { i i } & v _ { i j } \\ v _$$

It should be noted that the order of multiplication in (10) is not unique, and a different order will lead to different rotation angles θ ij . In this paper, without loss of generality, we use the order definition in (10).

V ij represents a rotation from the i th standard basis to the j th standard basis in R n t vector space with a rotation angle θ ij . That is, we show that an arbitrary orthogonal matrix V can be represented by (10). Further, an arbitrary covariance matrix Q ∈ R n t × n t can be represented by n t non-negative eigenvalues and 1 2 n t ( n t -1) rotation angles.

Lemma 1. To reach the secrecy capacity of the MIMOME with n t ≥ 2 , it is sufficient to use a PSD diagonal Λ and the rotation matrix V is given in (10) to generate the input covariance matrix Q = VΛV T .

Proof. First, we prove that Q = VΛV T is a covariance matrix. It is straightforward to check that V in (10) is an orthonormal matrix, i.e., VV T = I . Since diagonal elements of Λ are non-negative, Q is symmetric and PSD, i.e., a covariance matrix.

Next, we prove that an arbitrary covariance matrix Q can be written as Q = VΛV T while V is defined as (10). It suffices to find θ ij such that (13) holds for a given orthonormal V

$$\begin{smallmatrix} i j \, \text{such that} \, ( 1 5 )$$

This process is carried out in Algorithm 1. In fact, (13) applies a series of Givens rotation on V . They keep the relative orthogonality and rotates V to identity matrix I .

To better appreciate this, we note that if we expand the product of matrices on the left side of V in (13), we see that V is initially multiplied with V T 12 . Then, the corresponding rotation angle θ 12 can be chosen to set the entry (2 1) , of ( V V 12 ) to zero. Next, V T 13 is multiplied to the new V and θ 13 can be chosen to set the entry (3 1) , to zero. This process continues until the last element under the main diagonal of V , i.e., the entry ( n , n t t -1) , becomes zero. We note that, since V is orthonormal, the upper triangle also will be zero throughout this process. That is, the left side of (13) becomes an identity matrix. This proves Lemma 1. The rotation angles θ ij can be obtained by Algorithm 1 which is a generalization of Algorithm 5.1.3 of [33] for vectors in R n t , where atan2( · , · ) is the four-quadrant inverse tangent denoted bt atan2 in MATLAB.

In certain cases, the eigenvalue decomposition of Q may give an improper rotation matrix [34] whose determinate is -1 and cannot be converted to an identical matrix by rotation. As a result, V ij for (13) does not exist. To deal with this issue, the improper rotation matrix can be converted to a proper rotation matrix by exchanging arbitrary two eigenvectors of V and the corresponding eigenvalues of Λ . As an example, we can define V ′ = VI ′ and Λ ′ = I ′ T ΛI ′ , where

This specific I exchanges the last two eigenvectors and eigenvalues of Q . Since I I ′ ′ T = I , it is clear that

$$\stackrel { \text{p-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-b-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-H-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-l-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-g-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-g-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-h-s-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-h-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-s-h-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-h-s-h-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-h-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-g-h-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-k-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-g-h-$$

$$& \text{values} \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cap \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \downarrow \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \upcup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \, \cup \,.$$

In such a case, V ′ is a rotation matrix and the Q will remain the same. This completes the proof of Lemma 1.

Lemma 1 shows that any covariance matrix can be formed using a rotation matrix and a diagonal power allocation matrix.

## Algorithm 1 Rotation Angles Solution

- 1: Initialize [ V Λ , ] = eig( Q ) , i = 1 ;
- 2: if det( V ) = -1 then
- 3: Exchange first two columns of V ;
- 4: Exchange first two values on diagonal of Λ ;

$$\text{ance} 3 \text{ - Exchange first two columns of } V ; \text{ance} 4 \text{ - Exchange first two values on diagonal of } \Lambda ; \text{$2$} \\ \text{$5$} \text{ end if} \\ \text{$6$} \text{ while } i \leq ( n _ { t } - 1 ) \text{ do} \\ \text{$7$} \text{ - } j = i + 1 ; \text{$8$} \\ \text{$8$} \text{ while } j \leq n _ { t } \text{ do} \\ \text{$9$} \text{ - } \theta _ { i j } = - \text{atan2(-V(j,i),V(i,i));V(i,i));$} \\ \text{$10$} \text{ - } \text{$v_{rot}=I_{n_{t}}$} ; \text{$1$} \\ \text{$1$} \text{ - } \text{$v_{rot}(i,i)=V_{rot}(j,j)=cos\theta_{ij}$} ; \text{$12$} \text{ - } \text{$v_{rot}(j,i)=-V_{rot}(i,j)=sin\theta_{ij}$} ; \text{$13$} \text{ - } \text{$V=V_{rot}V$} ; \text{$14$} \\ \text{$(13)} \text{ 14$} \text{ - } \text{$j=j+1$} ; \text{$15$} \text{ - } \text{end while} \\ \text{$16$} \text{ - } i = i + 1 ; \text{$1$} \\ \text{plies} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{                                                                                                                                                                                                        
    \text{                                                                                                                                                                                                    
   \text{                                                                                                                                                                                                  
   \text{\text{ative} } \text{                                                                                                                                                                                                   
                                                                                                                                                                                                        }                                                                                                                                                                                                        \text{                                                                                                                                                                                                 
      \text{                                                                                                                                                                                               
        \text{$$

- 17: end while
- 18: Output θ ij , ∀ 1 ≤ i &lt; j ≤ n t .

This is useful in many optimization problems, including that of the MIMOME channel in (3), as it removes the PSD constraint of the covariance matrix (i.e., Q glyph[followsequal] 0 ). Instead, we will have a set of linear constraints to make sure that the diagonal elements of the power allocation matrix are non-negative and their sum is not greater than P t .

Specifically, for any n t , the optimization problem (3) is equivalently reformulated as

$$\widetilde { \text{es} } _ { \text{on} \\ \cdot ) } \left ( P 2 \right ) \, \mathcal { C } _ { s } = \max _ { \lambda, \theta } \frac { 1 } { 2 } \log _ { 2 } \frac { | \mathbf I _ { n _ { t } } + \mathbf H ^ { T } \mathbf H \mathbf V \mathbf A \mathbf V ^ { T } | } { | \mathbf I _ { n _ { t } } + G ^ { T } G \mathbf V \mathbf A \mathbf V ^ { T } | }, \quad ( 1 6 a )$$

$$\lambda _ { i } \geq 0, i \in \{ 1, \dots, n _ { t } \},$$

$$\lambda, \vartheta \ 2 & \quad | \mathfrak { m } _ { n _ { t } } + \mathfrak { m } ^ { * } \ \mathfrak { m } v \ \mathfrak { m } v ^ { * } | \\ \text{s.t.} \quad & \sum _ { i = 1 } ^ { n _ { t } } \lambda _ { i } \leq P _ { t }, \quad & \quad ( 1 6 b ) \\ & \quad \lambda _ { i } \geq 0, i \in \{ 1, \dots, n _ { t } \}, \quad & \quad ( 1 6 c )$$

in which we have defined

$$\lambda & \stackrel { \triangle } { = } \{ \lambda _ { i } \}, \, 1 \leq i \leq n _ { t }, \, \quad \quad \quad ( 1 7 \tt a ) \\ \tt a & \stackrel { \triangle } { = } \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gam, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma\, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma \, \Gamma$$

$$\theta & \stackrel { \triangle } { = } \{ \theta _ { i j } \}, \, 1 \leq i < j \leq n _ { t }. \quad \quad ( 1 7 b ) \\. & \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon$$

as the compact form of the parameters of (10)-(12) . 4

Different from (3), this new representation replaces the constraint that Q is symmetric and PSD by linear constraints on λ only while rotation angles can take any number, i.e., θ ∈ R . Then, numerical methods can be applied to optimize the parameters θ and λ to obtain the optimal secrecy rate. For Q with dimension n t , the required number of eigenvalues is n t , whilst this number is 1 2 n t ( n t -1) for the rotation angle. The total number of parameters is 1 2 n t ( n t +1) , which is equal to the number of the elements of the upper triangular of Q .

4 It is worth mentioning that Givens rotation method applied in many problems can finally boil down to finding the optimal rotation parameters. However, the application of the rotation method is different. In [23], only the optimal unitary beamforming without power allocation is evaluated because it constrained the precoding matrices to be unitary for limited feedback. On the other hand, the optimization approach is in different manners. In their iterative algorithm, the rotation angles of each covariance matrix are obtained one by one. In our paper, these parameters can be updated vector-wised once.

Theoretically, the rotation modeling method can provide a systematic approach to traverse Q by traversing λ and θ in finite regions. The feasible region of λ is given in (16b)-(16c) and for each rotation angel the regions can be [0 , 2 π ) 5 . The rotation method can be applied to many other problems, some listed in Section I-B.

As mentioned earlier, for n t = 2 a closed-form solution for (P2) is known in [14]. However, finding closed-form solutions of λ and θ is challenging for n t ≥ 3 . Next, we introduce a structure to solve this problem effectively.

## IV. SOLVING THE NEW OPTIMIZATION PROBLEM

In this section, we develop a novel method to solve (P2). Specifically, we convert (P2) to an unconstrained problem and resort to BFGS [35] to solve it. We should highlight that various iterative optimization methods, such as gradient descent, Newton's method, and quasi-Newton algorithms with interior-point method (that is, MATLAB convex optimization tool fmincon ), may be used to optimize (P2). Among them, Newton's method is the fastest, but it requires the local Hessian matrix and its inverse, which may not exist at some points. Quasi-Newton methods keep the convergence advantage of Newton's method without requiring the local Hessian matrix and its inverse. BFGS algorithm is an outstanding variety of quasi-Newton methods that can converge faster for nonconvex problems [36], [37]. The proposed method is called rotation-BFGS. The block diagram of this method is given in Fig. 3 in which the inputs are channel matrices H and G and average power P t and the outputs are the parameters of the rotation model (i.e., λ ∗ and θ ∗ ) as well as the corresponding secrecy rate R ∗ . As shown in Fig. 3, there are four blocks in the BFGS-based optimizer: GSVD for initialization, BFGS for optimization, eigenvalues rectifier, and objective function (16a) evaluator. In the following subsections, we introduce the functionality of each block and the optimization algorithm.

Fig. 3: The system design of rotation-BFGS method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-984e2777-6c9d-4ef0-9302-f0ad628d9d8a.jpg)

## A. Functionality of Each Block

1) Initialization using GSVD: While initial values of λ and θ can be chosen randomly, efficient initial values can save time by reducing the number of iterations. For this reason, and

5 It can be proved that in some cases the optimal θ ij is in [0 , π ) . See, for example, the case for n t = 2 [14].

knowing that GSVD-based beamforming is a good solution for this problem, we use GSVD to find the initial values ( λ (0) and θ (0) ) for our rotation-BFGS algorithm. The solution given by GSVD-based beamforming provides a precoding matrix E which satisfies:

$$\text{H} E = \Psi _ { r } C, \text{ \quad \ \ } ( 1 8 a )$$

$$G E = \Psi _ { e } D, \text{ \quad \ \ } ( 1 8 b )$$

$$C ^ { T } C + D ^ { T } D = \mathbf I, \text{ \quad \ \ } \text{(1&c)}$$

where E ∈ R n t × q , q = min( n , n t r + n e ) , C C T = diag( c i ) and D D T = diag( d i ) , i ∈ { 1 , . . . , q } , are diagonal matrices, and Ψ r ∈ R n r × n r and Ψ e ∈ R n e × n e are orthonormal matrices. Besides, the power allocation matrix P = diag( p i ) is determined by [15]

$$\text{gmgu} \quad \text{is odd} \quad \text{$vy$ } \text{$1\nu$} \\ \text{radient} \quad \text{as with} \quad \text{$p_{i}=\begin{cases} \max (0, \frac{2 ( c _ { i } - d _ { i } ) / ( \mu e _ { i } ) - 2 } { 1 + \sqrt { 1 - 4 c _ { i } d _ { i } + 4 ( c _ { i } - d _ { i } ) c _ { i } d _ { i } / ( \mu e _ { i } ) } } ), & \text{if $c_{i}>d_{i},$} \\ 0, & \text{otherwise,} \\ \text{lessian} \end{cases}$$

in which p i and e i are the i th diagonal element of P and E E T respectively, and µ is the Lagrange multiplier to ensure

$$\text{tr} ( \text{EPE} ^ { T } ) = P _ { t }. \text{ \quad \ \ } ( 2 0 )$$

GSVD-based Q is then EPE T . Thus, we can determine V (0) and Λ (0) from eigenvalue decomposition of EPE T , i.e., from

$$\mathbf Q ^ { ( 0 ) } = \text{EP} \, \text{$T=V^{(0)} \Lambda^{(0)} \mathbf V^{(0) T}$.}$$

Then, λ (0) = diag( Λ (0) ) and θ (0) can be obtained from V (0) using Algorithm 1.

2) BFGS Optimizer: Different iterative optimization methods, such as gradient descent, Newton's method, and quasiNewton algorithms may be used to optimize (P2). Among the above methods, Newton's method is the fastest method, but it requires a local Hessian matrix and its inverse, which may not exist or hard to obtain. Quasi-Newton methods keep the convergence advantage of Newton's method without requiring the local Hessian matrix and its inverse. An outstanding variety of quasi-Newton methods is the BFGS algorithm [35]. We give a brief introduction to this method in the following.

Given an unconstrained optimization problem

$$\arg \min _ { x } f ( { \mathbf x } ), \quad \quad \quad ( 2 2 )$$

with argument vector x and objective function f ( x ) , BFGS algorithm updates the vector x iteratively according to

$$x ^ { ( k + 1 ) } = x ^ { ( k ) } - \alpha ^ { ( k ) } \mathbf M ^ { ( k ) } \mathbf g ^ { ( k ) }, \ k \geq 0, \quad ( 2 3 )$$

in which

- · α ( k ) is the step size to minimize f ( x ( k +1) ) , which can be obtained via a line search.
- · g ( k ) is the gradient of f ( ) · at x ( k ) .
- · M ( k ) is an approximation of the inverse of the Hessian matrix.

The matrix M is initialized by a unity matrix, i.e., M (0) = I , and is then updated as [35]

$$\text{and is then updated as } [ 3 5 ] & & \text{namely,} \\ M ^ { ( k + 1 ) } = & \left ( \text{I } - \frac { \delta _ { x } ^ { ( k ) } \delta _ { g } ^ { ( k ) T } } { \delta _ { g } ^ { ( k ) T } \delta _ { x } ^ { ( k ) } } \right ) M ^ { ( k ) } \left ( \text{I } - \frac { \delta _ { g } ^ { ( k ) } \delta _ { x } ^ { ( k ) T } } { \delta _ { g } ^ { ( k ) T } \delta _ { x } ^ { ( k ) } } \right ) & & \text{This as an e} \\ & + \frac { \delta _ { x } ^ { ( k ) } \delta _ { x } ^ { ( k ) T } } { \delta _ { g } ^ { ( k ) T } \delta _ { x } ^ { ( k ) } }, & & \text{region } \iota \\ \text{in which}$$

in which

$$\delta _ { x } ^ { ( k ) } & \stackrel { \triangle } { = } x ^ { ( k + 1 ) } - x ^ { ( k ) }, \quad \quad \ \ ( 2 5 a ) \quad \text{the } f \\ \tilde { x } ^ { ( k ) } & \stackrel { \triangle } { = } \tilde { \pi } ^ { ( k + 1 ) } \quad \tilde { \pi } ^ { ( k ) } \quad \ \text{$\varkappa \times \ \text{sec} \alpha$}$$

$$\delta _ { g } ^ { ( k ) } \stackrel { \triangle } { = } g ^ { ( k + 1 ) } - g ^ { ( k ) }, \quad \quad \ \ ( 2 5 b ) \quad \sec _ { \text{folio} } ^ { \cdot }$$

respectively, represent the difference between the arguments and the gradients in two successive iterations. The gradient can be obtained analytically or numerically.

Due to its fast convergence speed and self-correcting property [36], BFGS is widely used in unconstrained optimization problems. However, the problem in this paper is a constrained optimization problem due to (16b) and (16c). To overcome this limitation, we relax the constraints of (P2) when using the BFGS optimizer, but we add a new block called a rectifier, as shown in Fig. 3. The rectifying block is used to ensure that the constraint on λ i s are satisfied as elaborated on in what follows.

3) Rectifying Eigenvalues: To use BFGS, the optimization problem should be unconstrained. However, the eigenvalues in (P2) need to be constrained. To overcome this issue, we rectify the eigenvalues to ensure that the constraints in (16b)-(16c) are satisfied.

Before talking about the rectification process, we highlight that we only optimize the first n t -1 eigenvalues because the last eigenvalue will be obtained by ∑ n t i =1 λ i = P t . That is, we use P t -∑ n t -1 i =1 λ i as the value of the last eigenvalue ( λ n t ) since in the MIMOME channel it is known that optimal solution uses the total power [7], [13], [14].

$$\lambda = r ( \tilde { \lambda }, P _ { t } ), \quad \ \ \ ( 2 6 ) \ \ \ r o t a$$

Suppose ˜ λ ∈ R n t -1 is the vector of first n t -1 unconstrained eigenvalues obtained from the BFGS algorithm. We obtain λ ∈ R n t (the rectified eigenvalue vector) from where the rectifying function r( · , · ) is defined by the following successive processes:

$$\lambda ^ { + } = [ \tilde { \lambda } ] ^ { + },$$

$$\lambda = \left [ \bar { \lambda }, P _ { t } - \sum _ { i = 1 } ^ { n _ { t } - 1 } \bar { \lambda } \right ]. \quad \begin{array} { c c c } \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \alpha & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta &\beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta & \beta \\ \ln \left ( 2 7 a \right ). \left [ \cdot \right ] ^ { + } & \text{is an element-wise operation which forces all} & \quad \text{the rel}$$

$$& \lambda ^ { + } = [ \tilde { \lambda } ] ^ { + }, \\ & \bar { \lambda } = \begin{cases} \lambda ^ { + } \cdot \frac { P _ { t } } { \sum _ { i = 1 } ^ { n _ { t } - 1 } \lambda _ { i } ^ { + } }, & \sum _ { i = 1 } ^ { n _ { t } - 1 } \lambda _ { i } ^ { + } > P _ { t }, & \text{otherwise}, & \text{the $arg$} \\ \Gamma & n _ { t } - 1 & \Gamma & \text{otherwise}, & \text{the cur} \\ \end{cases}$$

In (27a), [ ] · + is an element-wise operation which forces all negative elements of ˜ λ to 0 while keeping non-negative values unaltered. This will take care of the constraints in (16c). But, the elements of the new vector λ + may not still satisfy (16b). In such a case, in (27b), we scale the new vector such that the sum of the eigenvalues does not exceed P t . Finally, in the (27c), the last eigenvalue is added. Thus, the problem (P2)

can be solved using an unconstrained optimization method, namely, the BFGS.

This process illustrated in Fig. 4 for n t = 3 and P t = 5 , as an example. In this figure, the blue (shaded) area denotes the feasible region of the eigenvalues, red points outside the region are unconstrained eigenvalues ˜ = [ λ λ , λ 1 2 , · · · , λ n t -1 ] (the output of the BFGS block), and green squares denote their corresponding rectified eigenvalues. With the proposed rectification in (27a)-(27b), unrestricted inputs will be forced to new points on the boundary of the feasible region. For example, for the point ( -1 4) , we have ˜ = [ λ -1 4] , , ¯ = λ λ + = [0 4] , (the second case in (27b)) and λ = [ λ , λ 1 2 , P -λ 1 -λ 2 ] = [0 4 1] , , , following (27a)-(27b). Similarly, the red point ( -1 , -1) will result in λ = [0 0 5] , , .

Fig. 4: Illustration of eigenvalues rectification for n t = 3 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-8be3f36a-e5bd-49fa-94a3-b8989b879986.jpg)

4) The Objective Function: In this paper, our goal is to maximize the secrecy rate R defined by

$$R ( \lambda, \theta ) & \stackrel { \triangle } { = } \frac { 1 } { 2 } \log _ { 2 } \frac { | I _ { n _ { t } } + V ^ { T } \mathbf H ^ { T } \mathbf H V \Lambda | } { | I _ { n _ { t } } + V ^ { T } \mathbf G ^ { T } \mathbf G V \Lambda | }, \quad ( 2 8 ) \\ \text{which is a function of $A$ and $\Lambda$ as shown in $P$} \times \text{Noting}$$

which is a function of θ and λ as shown in (P2). Noting that our problem is a maximization rather than a minimization problem, in comparison to (22), we define the objective function f ( x ) = -R . With this, we can link BFGS and rotation method together.

## B. Rotation-BFGS Algorithm

As shown in Fig. 3, GSVD provides an initial value for the argument x in the BFGS block. The rectifier and the objective function can be considered together which require the current value of x and provide corresponding achievable rate back to the BFGS. Then, BFGS will update x and check the termination condition.

To link the BFGS to the rotation model, we need to define the relation between x , the arguments inside BFGS, and the parameters λ and θ . This is given by

$$x \stackrel { \triangle } { = } [ \tilde { \lambda }, \theta ],$$

where ˜ λ is the first n t -1 elements of λ and θ is the same as we defined in (17b). In this way, we reduce one argument (eigenvalue) for efficiency. Then, in the k th iteration ( k ≥ 0 )

of the rotation-BFGS algorithm, the relation between argument x ( k ) in the BFGS and the eigenvalues λ ( k ) in precoding is

$$x ^ { ( k ) } \stackrel { \triangle } { = } [ \tilde { \lambda } ^ { ( k ) }, \theta ^ { ( k ) } ],$$

$$\lambda ^ { ( k ) } \stackrel { \triangle } { = } r ( \tilde { \lambda } ^ { ( k ) }, P _ { t } ), \quad \quad \quad ( 3 0 b )$$

in which ˜ λ ( k ) denotes the first n t -1 unconstrained eigenvalues and λ ( k ) is the constrained (rectified) eigenvalues which are obtained by (26)-(27c). Besides, the value of the function is denoted as

$$f ^ { ( k ) } \stackrel { \triangle } { = } f ( { \mathbf x } ^ { ( k ) } ) = - R ( { \mathbf r } ( \tilde { \mathbf \lambda } ^ { ( k ) }, P _ { t } ), \theta ^ { ( k ) } ), \quad ( 3 1 ) \quad \Big | _ { \quad } \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

in which λ ( k ) and θ ( k ) can be obtained from x ( k ) according to (30a)-(30b), and f ( ) · denotes the mapping from x ( k ) to f ( k ) . Finally, the gradient vector g ( k ) with respect to x ( k ) is obtained numerically. Specifically, the i th element of g in the k th iteration is given by

$$g _ { i } ^ { ( k ) } = \frac { f ( \mathbf x ^ { ( k ) } + \epsilon _ { i } ) - f ^ { ( k ) } } { | \epsilon _ { i } | }, \quad ( 3 2 ) \quad \Big | \\ \text{has the same length as $x$ and its all commonent are}$$

where glyph[epsilon1] i has the same length as x and its all component are zero except for the i th element which is a constant glyph[epsilon1] 1 .

## Algorithm 2 Rotation-BFGS Method

- 1: Initialize: glyph[epsilon1] 1 = glyph[epsilon1] 2 = 10 -4 , k = 0 , and M (0) = I ;
- 2: Find V (0) and Λ (0) using GSVD in (21);
- 3: Find λ (0) which is the diagonal of Λ (0) ;
- 5: Find θ (0) using Algorithm 1;
- 4: Find ˜ λ (0) which is the first ( n t -1) elements in λ (0) ;
- 6: Find R (0) using (28);
- 7: Find x (0) , g (0) , and f (0) using (30a)-(32);
- 8: while 1 do
- 9: Find α ( k ) = α ∗ by line search using Algorithm 3;
- 10: Find x ( k +1) using (23);
- 11: Find λ ( k +1) and θ ( k +1) from x ( k ) using (30a)-(30b);
- 12: Find R ( λ ( k +1) , θ ( k +1) ) using (28);
- 13: Find g ( k +1) and f ( k +1) using (30a)-(32);
- 14: Find M ( k +1) using BFGS by (24)-(25b);

$$\begin{smallmatrix} \iota \colon & \text{$\Piu$ is $\alpha$ using $Du$} \cup \nolimits ( \sigma \nu ( \sigma \nu ), & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \text{$\Delta \nu \nolimits ( \sigma \nu \nolimits ) }, & \begin{smallmatrix} \iota \colon & \text{$\Piu$ is $\alpha$} \\ \iota \colon & \text{$\Piu$ is $\alpha$} \end{smallmatrix} \\ \iota \colon & \text{$\Piu$ is $\alpha$} \end{smallmatrix}$$

18:

end if

- 19: k = k +1 ;

20: end while

- 21: Output: R ∗ , λ ∗ , and θ ∗ .

In our proposed rotation-BFGS method, as shown in Fig. 3, GSVD-based beamforming is first applied to find initial values of λ (0) , θ (0) which are obtained by GSVD decomposition as described in Section IV-A1. Thus x (0) = [ ˜ λ (0) , θ (0) ] . We should also highlight that ˜ λ (0) contains the first n t -1 elements of λ (0) .

Algorithm 2 illustrates the details of the proposed optimization process. Within this algorithm, we require a line search. The details of the line search are given in Algorithm 3 in the Appendix. Algorithm 2 will terminate if the secrecy rate at two successive iterations are very close, i.e., when their difference is smaller than a tolerance ( glyph[epsilon1] 2 ). Algorithm 2 guarantees a

TABLE I: Achievable rate (in bps/Hz) of each method for n t = 3 and P t = 30 W .

| Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   |
|-----------------|-----------------|-----------------|-----------------|-----------------|-----------------|-----------------|-----------------|
| 2               | 2               | 3               | 3               | 3               | 3               | 3               | 3               |
| 4               | 4               | 5               | 6               | n r             | 2.58            | n t = 3 1 1.99  | 1.14            |
|                 | 0.39            | 0.23            | 1 2             | 3.91            | 2.99            | 1.75            | 0.63 1.09       |
|                 | 0.71            | 0.43            | 3               | 4.84            | 3.59            | 2.24            | 1.47            |
|                 | 1.01            | 0.70            | 4               | 5.48            | 4.16            | 2.70            | 1.73            |
|                 | 1.24            | 0.89            | 5               | 5.97            | 4.52            | 3.01            | 2.03            |
|                 | 1.45            | 1.09            | 6               | 6.46            | 4.88            | 3.30            | 2.35            |
|                 | 1.76            | 1.28            | Rotation-Fmin   | 2               | n e 3           | n t = 3 1 4     | 5               |
| 6               | 6               | 6               | 6               | 6               | 6               | 6               | 6               |
| 1               | 1               | 1.99            | 1.99            | 1.99            | 1.99            | 1.99            | 1.99            |
| 1.14            | 1.14            | 0.63            | 0.39            | 0.23            |                 | 2.58 2 2.99     | 3.91            |
|                 |                 | 1.75            | 1.09            | 0.71            | 0.43            | 4.84            | 3               |
|                 | 3.59            | 2.24            | 1.47            | 1.01            | 0.70            | n r 5.48        | 4               |
|                 | 4.16            | 2.70            | 1.73            | 1.24            | 0.89            | 5.97            | 5               |
|                 | 4.52            | 3.01            | 2.03            | 1.45            | 1.09            | 6.46 3.30       | 6               |
|                 | 4.88            |                 | 2.35            | 1.76 3          | 1.28 n e        | GSVD [15] = 3 1 | 2               |
| n t             | 4               | 5               | 6               | n r             | 1               | 2.56            | 1.82            |
| 1.12            | 1.12            | 1.12            | 1.12            | 1.12            | 1.12            | 1.12            | 1.12            |
| 0.63            | 0.63            | 0.39            | 0.39            | 0.39            | 0.39            | 0.39            | 0.39            |
| 0.23            | 0.23            | 1.74            | 2               | 2.97            | 2.87 0.70       |                 | 1.09            |
|                 | 0.43            |                 | 3               | 4.41            | 3.50 1.01       | 2.23            | 1.46            |
|                 | 0.70            |                 | 4               | 5.17            | 4.10 1.24       | 2.69            | 1.73            |
|                 | 0.89            |                 | 5               | 5.76            | 4.46 1.45       | 3.00            | 2.02            |
|                 | 1.09            | 3.29            | 6               | 6.28            | 4.83 1.76       |                 | 2.35            |
|                 | 1.28            | AOWF [16]       | 1               | 2               | n e 3 6         | = 3 4           | 5               |
| n t             | n r             | 1               | 2.58            | 1.99            | 1.13            | 0.63            | 0.38            |
| 0.23            | 0.23            | 1.74 2.23       | 1.74 2.23       | 1.74 2.23       | 1.74 2.23       | 1.74 2.23       | 1.74 2.23       |
| 2               | 2               |                 | 2.98            |                 | 1.08            | 3.92 0.69       | 0.42            |
|                 |                 | 4.86            | 3.58            |                 | 1.45            | 3 0.99          | 0.68            |
|                 |                 | 5.49            | 4.15            | 2.69            | 1.70            | 4 1.20          | 0.85            |
|                 |                 | 5.99            | 4.51            | 2.98            | 1.97            | 5 1.39          | 1.02            |
|                 |                 | 6.47            | 4.87            | 3.26            | 2.28            | 6 1.67          | 1.18            |

global convergence based on the global convergence theorem [35, Chapter 7, pp. 196-204], because the optimized function (31) (the objective function together with the rectifier in Fig. 3) is continuous, and the iteration in (23) can lead to a decreasing sequence, i.e., f ( x ( k +1) ) ≤ f ( x ( k ) ) .

We analyze the complexity of the different algorithms here. The computation of matrix multiplications and matrix inverse yields the complexity of O ( L 3 ) where L = max( n , n t r , n e ) . GSVD-based precoding [15], used as the initial point generator, has the complexity of O ( L 3 + L log(1 /glyph[epsilon1] )) [38] in which glyph[epsilon1] is the convergence tolerance of the algorithm, while bisection search requires O (log(1 /glyph[epsilon1] )) iterations [39]. Besides, the BFGS algorithm has the complexity of O ( n 2 ) [37], where n is the size of input variables which is the total number of optimized rotation parameters, i.e., 1 2 n t ( n t +1) . Thus, the overall complexity of Algorithm 2 is O ( n 4 t + L 3 + L log(1 /glyph[epsilon1] )) . AOWF [16, Algorithm 1] yields O ( L 3 glyph[epsilon1] log(1 /glyph[epsilon1] )) , in which O ( 1 glyph[epsilon1] ) and O (log(1 /glyph[epsilon1] )) are the outer layer loop and inner bisection search, respectively. It should also be mentioned that rotationFmin has the same asymptotic complexity as rotation-BFGS. However, for small values of n t , rotation-BFGS has about an order of magnitude smaller complexity than rotation-Fmin, as we show in Section V.

TABLE II: Achievable rate (in bps/Hz) of each method for n t = 4 and P t = 30 W .

| Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   | Rotation-BFGS   |
|-----------------|-----------------|-----------------|-----------------|-----------------|-----------------|-----------------|
|                 |                 | n e             | n e             | n e             | n e             | n e             |
| n t = 4         | n t = 4         | 1               | 2               | 3               | 4               | 5               |
|                 | 1               | 3.04            | 2.63            | 2.06            | 1.25            | 0.77            |
|                 | 2               | 4.72            | 4.01            | 3.12            | 2.02            | 1.32            |
|                 | 3               | 5.91            | 4.97            | 3.78            | 2.60            | 1.72            |
|                 | 4               | 6.81            | 5.68            | 4.40            | 3.05            | 2.08            |
|                 | 5               | 7.63            | 6.28            | 4.82            | 3.46            | 2.50            |
|                 | 6               | 8.24            | 6.76            | 5.22            | 3.78            | 2.74            |
| Rotation-Fmin   | Rotation-Fmin   | Rotation-Fmin   | Rotation-Fmin   | Rotation-Fmin   | Rotation-Fmin   | Rotation-Fmin   |
|                 |                 | n e             | n e             | n e             | n e             | n e             |
| n t = 4         | n t = 4         | 1               | 2               | 3               | 4               | 5               |
|                 | 1               | 3.04            | 2.63            | 2.06            | 1.25            | 0.77            |
|                 | 2               | 4.73            | 4.02            | 3.12            | 2.02            | 1.32            |
|                 | 3               | 5.91            | 4.98            | 3.78            | 2.60            | 1.72            |
|                 | 4               | 6.81            | 5.69            | 4.40            | 3.05            | 2.08            |
|                 | 5               | 7.63            | 6.29            | 4.82            | 3.46            | 2.50            |
|                 | 6               | 8.24            | 6.77            | 5.22            | 3.78            | 2.74            |
| GSVD [15]       | GSVD [15]       | GSVD [15]       | GSVD [15]       | GSVD [15]       | GSVD [15]       | GSVD [15]       |
|                 |                 | n e             | n e             | n e             | n e             | n e             |
| n t = 4         | n t = 4         | 1               | 2               | 3               | 4               | 5               |
|                 | 1               | 3.03            | 2.59            | 1.82            | 1.22            | 0.76            |
|                 | 2               | 4.27            | 3.02            | 2.97            | 1.99            | 1.31            |
|                 | 3               | 3.87            | 4.42            | 3.67            | 2.57            | 1.72            |
|                 | 4               | 5.80            | 5.33            | 4.31            | 3.03            | 2.08            |
|                 | 5               | 6.91            | 6.03            | 4.75            | 3.44            | 2.50            |
|                 | 6               | 7.70            | 6.56            | 5.16            | 3.76            | 2.73            |
| AOWF [16]       | AOWF [16]       | AOWF [16]       | AOWF [16]       | AOWF [16]       | AOWF [16]       | AOWF [16]       |
|                 |                 | n e             | n e             | n e             | n e             | n e             |
|                 | 4               | 1               | 2               | 3               | 4               | 5               |
|                 | 1               | 3.04            | 2.63            | 2.06            | 1.24            | 0.76            |
|                 | 2               | 4.74            | 4.04            | 3.12            | 2.01            | 1.31            |
|                 | 3               | 5.94            | 4.98            | 3.77            | 2.58            | 1.70            |
|                 | 4               | 6.82            | 5.68            | 4.39            | 3.04            | 2.06            |
|                 | 5               | 7.63            | 6.29            | 4.81            | 3.44            | 2.46            |
|                 | 6               | 8.24            | 6.75            | 5.21            | 3.74            | 2.69            |

## V. NUMERICAL RESULT

In this section, extensive numerical results are provided to illustrate the performance of the proposed rotation-BFGS method. Four methods are taken into account:

- · Rotation-BFGS : the proposed rotation-BFGS parameterization solved by the BFGS method.
- · Rotation-Fmin : the proposed rotation modeling solved using MATLAB convex optimization tool fmincon .
- · GSVD : GSVD-based beamforming with optimal power allocation [15], as described in Section IV-A1.
- · AOWF : alternating optimization and water-filling [16].

All results are based on averaging over 1000 realizations of independent H and G . These entries of H and G are generated based on the standard Gaussian distribution, i.e., N (0 1) , . Two performance metrics are considered. In the first part, we focus on the achievable secrecy rate. In the second part, the time consumption is analyzed for different methods in various antenna configurations.

## A. Achievable Secrecy Rate

We evaluate the performance of the Rotation-BFGS, Rotation-Fmin, GSVD, and AOWF approaches with respect to the variation of the number of antennas. The average secrecy rates are listed in Tables I and II respectively for n t = 3 and n t = 4 . As expected, increasing the number of antennas

P t =30 (W), n =3 t

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-f119def2-abd6-44a4-af0f-8b13e144d8b7.jpg)

(a) Relative rate improvement ( % ) of the proposed method to GSVD.

## P t =30 (W), n =3 t

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-51199dd8-9381-4b00-b5af-06bac4abe116.jpg)

(b) Relative rate improvement ( % ) of the proposed method to AOWF.

Fig. 5: Comparisons between the secrecy transmission rates of the proposed method with GSVD and AOWF.

at the eavesdropper decreases secrecy rate. On the contrary, the secrecy rate will increase when the legitimate receiver or transmitter has higher number of antennas. To better appreciate the improvement due to our proposed numerical method, relative secrecy rate improvement between Rotation-BFGS and GSVD and Rotation-BFGS and AOWF is investigated in the following and illustrated in Fig. 5. Let us define relative rate improvement with respect to GSVD and AOWF, respectively, as

$$\eta _ { g } = \frac { R _ { r } - R _ { g } } { R _ { g } } \times 1 0 0 \%, \quad \quad ( 3 3 \tt a )$$

$$\eta _ { a } = \frac { R _ { r } - R _ { a } } { R _ { a } } \times 1 0 0 \%, \text{ \quad \ \ } ( 3 3 b )$$

(a) Achievable secrecy rate with relatively small n e s.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-13867dd8-a534-4e1a-8306-74797e436876.jpg)

(b) Achievable secrecy rate with relatively large n e s.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-380272c4-8d55-48bc-baa2-21beb8c2e8bd.jpg)

Fig. 6: Secrecy rate of the MIMOME channel versus the transmit power. The proposed algorithms (Rotation-BFGS and Rotation-Fmin) are compared with GSVD and AOWF.

where R r , R g , and R a represent the secrecy rate achieved by the proposed method, GSVD, and AOWF, respectively 6 . As shown in Fig. 5, the proposed approach is capable of achieving the same or better secrecy rate in almost any antenna setting. The darker the cell color, the higher is improvement achieved by the proposed method

There is a clear pattern that when the eavesdropper has a smaller number of antennas than the transmitter, the proposed method outperforms GSVD. The best case is [ n , n r e ] = [2 1] , where the improvement is about 32%. On the other hand, for cases with larger n e , the rotation-BFGS performs better than AOWF. The best case is [ n , n r e ] = [6 6] , where the improvement is 7.5%. This pattern not only exist when n t = 3 but also at higher n t s, at least n t = 4 . Some more detailed comparisons are presented as follows:

6 Rotation-Fmin is similar to rotation-BFGS in terms of the achievable rate (details are in Table I and II), so relative rate improvements of that to GSVD and AOWF are the same as those in rotation-BFGS.

Fig. 7: Time costs of the proposed methods, GSVD and AOWF, for (a) n t = 3 , (b) n t = 4 , with n r and n e from 1 to 6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-996c4e79-4b0e-433b-913b-0b4f0f4d140b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00994v2-figure-c423fd1a-dbd1-4c0c-b05a-ea6b3c9dcb18.jpg)

1) The MIMOME with small n e : GSVD-based beamforming fails to get close to the secrecy capacity. This phenomenon has been verified by previous literature [14]. As can be seen in Fig. 6(a), rotation-BFGS and AOWF can achieve similar results. On the contrary, GSVD clearly has a gap with those methods which results from sub-optimality of GSVD.

2) The MIMOME with large n e : In this setting, AOWF does not perform very well because the Lagrange multiplier of AOWF cannot be obtained properly. As illustrated in Fig. 6(b) and Fig. 5, our proposed method outperforms AOWF in this regime.

From the simulation results, it is seen that the proposed rotation-BFGS is robust in a wide range of practical antennas settings on each node. This is a big advantage as the robustness towards the variation of n e is necessary to guarantee secure communication. We note that the eavesdropper can have any number of antennas, and the solution for wiretap channels should be robust to such variations. This makes the proposed precoding highly competitive the existing solutions of the MIMOME channel.

## B. Time Consumption

Besides computational complexity, the execution time of each algorithm is also important and provides a means to evaluate the complexity. This evaluation is shown in Figs. 7(a)7(b) for n t = 3 and n t = 4 respectively, which are averaged over 1000 channel realizations in each setting. GSVD has, by far, the best time cost since it transforms the problem into finding the best Lagrange multiplier. However, as we illustrated previously, the precoding provided by GSVD may be far from the optimal solution. The proposed rotation-BFGS method clearly outperforms AOWF, whose time consumption is very high, especially when n e ∈ [4 , 5 6] , . The gap becomes larger when n e increases. In such cases, AOWF does not perform well either in terms of achievable rate or time cost. On the other hand, our approach is robust in any antenna setting and more efficient in time cost compared with AOWF. Rotation-Fmin achieves almost the same performance, but its execution time is higher. This is because fmincon includes BFGS and interior-point method, and instead of the interiorpoint method, we use a rectifier which is more efficient in terms of programming.

In summary, the proposed rotation-BFGS method gives robust precoding and power allocation for the MIMOME channel. By applying BFGS as an optimizer, the rotation method has a better performance compared to other numerical methods like AOWF in secrecy and time complexity. In addition, the framework we provided in this paper can be used to solve many other problems, some listed in Section I-B.

## VI. CONCLUSIONS

In this paper, we have developed a rotation-based method, called rotation-BFGS, for precoding and power allocation of Gaussian MIMOME channels. In this method, the transmit covariance matrix is constructed using Givens rotation matrices. With this construction, the PSD constraint of the transmit covariance matrix is removed and the capacity optimization problem is simplified. The precoding (rotation) matrix and power allocation coefficients have been obtained iteratively using a modified BFGS algorithm. Compared to existing approaches, the proposed method is robust and performs well independent of the number of antennas at each node. The proposed method outperforms the GSVD-based beamforming when n e &lt; n t and AOWF particularly when n e ≥ n t . This approach can also use the results of existing precoding and power allocation methods, such as the GSVD-based approach, as an initial point to expedite finding the solution.

In addition, the proposed rotation-BFGS method has a great potential for finding precoding and power allocation in various other applications, including in MIMO broadcast channel and MIMO channel with and energy harvesting constraints, both with and without secrecy. Future works will focus on further improving the efficiency of solving parameters and extension of this approach to other related problems.

## APPENDIX A

## LINE SEARCH ALGORITHM FOR THE BFGS METHOD

The line search method we use is the Golden section search [35] and summarized in Algorithm 3, where α i and f i for

1 ≤ i ≤ 4 are temporary records for line search steps size and function values.

## Algorithm 3 Line Search (Golden Section [35]) for the BFGS

- 1: Requires from Algorithm 2: glyph[epsilon1] 2 , x ( k ) , M ( k ) , g ( k ) , and f ( k ) ;
- 2: Define: x glyph[defines] x ( k ) , d glyph[defines] M ( k ) g ( k ) ;
- 3: Initialize: glyph[epsilon1] 3 = 5 × 10 -4 , τ 1 = 3 , τ 2 = 0 382 . , and τ 3 = 0 618 . ;
- 4: Initialize: α glyph[defines] [ α , α 1 2 , α 3 , α 4 ] = [0 0 0 0 1] , , , . ;
- 6: while α 4 &lt; 20 and f 1 &lt; f 4 do
- 5: Initialize: f glyph[defines] [ f , f 1 2 , f 3 , f 4 ] = [ f ( k ) , 0 0 , , f ( x -α 4 d )] ;
- 7: Let α 4 = τ 1 α 4 and f 4 = f ( x -α 4 d ) ;
- 8: end while
- 9: Let α 2 = τ 2 α 4 and α 3 = τ 3 α 4 ;
- 10: Let f 2 = f ( x -α 2 d ) and f 3 = f ( x -α 3 d ) ;
- 12: Define: m is the index of the smallest element in f ;
- 11: while α 4 -α 1 &gt; glyph[epsilon1] 3 and max( f ) -min( f ) &gt; glyph[epsilon1] 2 do
- 13: Define: q 1 glyph[defines] max(1 , m -1) and q 2 glyph[defines] min(4 , m +1) ;
- 14: Let α 1 = α q 1 and α 4 = α q 2 ;
- 15: Let α 2 = α 1 + τ 2 ( α 4 -α 1 ) and α 3 = α 1 + τ 3 ( α 4 -α 1 ) ;
- 16: if m = 2 then
- 17: Let f 3 = f 2 and f 2 = f ( x -α 2 d ) ;
- 18: else if m = 3 then
- 19: Let f 2 = f 3 and f 3 = f ( x -α 3 d ) ;
- 20: else
- 21: Let f 2 = f ( x -α 2 d ) and f 3 = f ( x -α 3 d ) ;
- 22: end if
- 23: end while
- 24: Output: α ∗ = α m .

## REFERENCES

- [1] X. Zhang, Y. Qi, and M. Vaezi, 'A new precoding for the MIMO Gaussian channel with multi-antenna eavesdroppers,' in Proc. IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC) , pp. 1-6, 2020.
- [2] A. D. Wyner, 'The wire-tap channel,' Bell system technical journal , vol. 54, no. 8, pp. 1355-1387, 1975.
- [3] U. M. Maurer, 'Secret key agreement by public discussion from common information,' IEEE Transactions on Information Theory , vol. 39, no. 3, pp. 733-742, 1993.
- [4] I. Csisz´ ar and J. Korner, 'Broadcast channels with confidential messages,' IEEE Transactions on Information Theory , vol. 24, no. 3, pp. 339-348, 1978.
- [5] N. Cai, A. Winter, and R. W. Yeung, 'Quantum privacy and quantum wiretap channels,' Problems of Information Transmission , vol. 40, no. 4, pp. 318-336, 2004.
- [6] J. Wu, Z. Lin, L. Yin, and G.-L. Long, 'Security of quantum secure direct communication based on Wyner's wiretap channel theory,' Quantum Engineering , vol. 1, no. 4, p. e26, 2019.
- [7] A. Khisti and G. W. Wornell, 'Secure transmission with multiple antennas-part II: The MIMOME wiretap channel,' IEEE Transactions on Information Theory , vol. 11, no. 56, pp. 5515-5532, 2010.
- [8] F. Oggier and B. Hassibi, 'The secrecy capacity of the MIMO wiretap channel,' IEEE Transactions on Information Theory , vol. 57, no. 8, pp. 4961-4972, 2011.
- [9] T. Liu and S. Shamai, 'A note on the secrecy capacity of the multipleantenna wiretap channel,' IEEE Transactions on Information Theory , vol. 55, no. 6, pp. 2547-2553, 2009.
- [10] S. Shafiee, N. Liu, and S. Ulukus, 'Towards the secrecy capacity of the Gaussian MIMO wire-tap channel: The 2-2-1 channel,' IEEE Transactions on Information Theory , vol. 55, no. 9, pp. 4033-4039, 2009.

| [11]   | P. Parada and R. Blahut, 'Secrecy capacity of SIMO and slow fading channels,' in Proc. IEEE International Symposium on Information Theory (ISIT) , pp. 2152-2155, 2005.                                                                                                                    |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [12]   | S. A. A. Fakoorian and A. L. Swindlehurst, 'Full rank solutions for the MIMO Gaussian wiretap channel with an average power constraint,' IEEE Transactions on Signal Processing , vol. 61, no. 10, pp. 2620- 2631, 2013.                                                                   |
| [13]   | S. Loyka and C. D. Charalambous, 'Optimal signaling for secure com- munications over Gaussian MIMOwiretap channels,' IEEE Transactions on Information Theory , vol. 62, no. 12, pp. 7207-7215, 2016.                                                                                       |
| [14]   | M. Vaezi, W. Shin, and H. V. Poor, 'Optimal beamforming for Gaussian MIMO wiretap channels with two transmit antennas,' IEEE Transactions on Wireless Communications , vol. 16, no. 10, pp. 6726-6735, 2017.                                                                               |
| [15]   | S. A. A. Fakoorian and A. L. Swindlehurst, 'Optimal power allocation for GSVD-based beamforming in the MIMO Gaussian wiretap channel,' in Proc. IEEE International Symposium on Information Theory (ISIT) , pp. 2321-2325, 2012.                                                           |
| [16]   | Q. Li, M. Hong, H.-T. Wai, Y.-F. Liu, W.-K. Ma, and Z.-Q. Luo, 'Transmit solutions for MIMO wiretap channels using alternating opti- mization,' IEEE Journal on Selected Areas in Communications , vol. 31, no. 9, pp. 1714-1727, 2013.                                                    |
| [17]   | S. Loyka and C. D. Charalambous, 'An algorithm for global maxi- mization of secrecy rates in Gaussian MIMO wiretap channels,' IEEE Transactions on Communications , vol. 63, no. 6, pp. 2288-2299, 2015.                                                                                   |
| [18]   | J. Steinwandt, S. A. Vorobyov, and M. Haardt, 'Secrecy rate maximiza- tion for MIMO Gaussian wiretap channels with multiple eavesdroppers via alternating matrix POTDC,' in Proc. IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP) , pp. 5686-5690, 2014. |
| [19]   | J. Li and A. Petropulu, 'Transmitter optimization for achieving se- crecy capacity in Gaussian MIMO wiretap channels,' arXiv preprint arXiv:0909.2622 , 2009.                                                                                                                              |
| [20]   | M. Vaezi, W. Shin, H. V. Poor, and J. Lee, 'MIMO Gaussian wiretap channels with two transmit antennas: Optimal precoding and power allocation,' in Proc. IEEE International Symposium on Information Theory , pp. 1708-1712, 2017.                                                         |
| [21]   | S. Loyka and C. D. Charalambous, 'On optimal signaling over secure MIMO channels,' in Proc. IEEE International Symposium on Informa- tion Theory (ISIT) , pp. 443-447, 2012.                                                                                                               |
| [22]   | M. A. Sadrabadi, A. K. Khandani, and F. Lahouti, 'Channel feedback quantization for high data rate MIMO systems,' IEEE Transactions on Wireless Communications , vol. 5, no. 12, pp. 3335-3338, 2006.                                                                                      |
| [23]   | R. De Francisco and D. T. Slock, 'An optimized unitary beamforming technique for MIMO broadcast channels,' IEEE Transactions on Wire- less Communications , vol. 9, no. 3, pp. 990-1000, 2010.                                                                                             |

- [24] Y. Yang, Q. Li, W.-K. Ma, J. Ge, and P. Ching, 'Cooperative secure beamforming for AF relay networks with multiple eavesdroppers,' IEEE Signal Processing Letters , vol. 20, no. 1, pp. 35-38, 2013.
- [25] J. Huang and A. L. Swindlehurst, 'Cooperative jamming for secure communications in MIMO relay networks,' IEEE Transactions on Signal Processing , vol. 59, no. 10, pp. 4871-4884, 2011.
- [26] J. Li and A. P. Petropulu, 'On ergodic secrecy rate for Gaussian MISO wiretap channels,' IEEE Transactions on Wireless Communications , vol. 10, no. 4, pp. 1176-1187, 2011.
- [27] Y. Wu, C. Xiao, Z. Ding, X. Gao, and S. Jin, 'Linear precoding for finitealphabet signaling over MIMOME wiretap channels,' IEEE Transactions on Vehicular Technology , vol. 61, no. 6, pp. 2599-2612, 2012.
- [28] W. Zeng, C. Xiao, M. Wang, and J. Lu, 'Linear precoding for finitealphabet inputs over MIMO fading channels with statistical CSI,' IEEE Transactions on Signal Processing , vol. 60, no. 6, pp. 3134-3148, 2012.
- [29] S. R. Aghdam, A. Nooraiepour, and T. M. Duman, 'An overview of physical layer security with finite-alphabet signaling,' IEEE Communications Surveys and Tutorials , 2018.
- [30] S. R. Aghdam and T. M. Duman, 'Joint precoder and artificial noise design for MIMO wiretap channels with finite-alphabet inputs based on the cut-off rate,' IEEE Transactions on Wireless Communications , vol. 16, no. 6, pp. 3913-3923, 2017.
- [31] X. Yang and A. L. Swindlehurst, 'Limited rate feedback in a MIMO wiretap channel with a cooperative jammer,' IEEE Transactions on Signal Processing , vol. 64, no. 18, pp. 4695-4706, 2016.
- [32] A. Khisti and G. W. Wornell, 'Secure transmission with multiple antennas I: The MISOME wiretap channel,' IEEE Transactions on Information Theory , vol. 56, no. 7, pp. 3088-3104, 2010.
- [33] G. H. Golub and C. F. Van Loan, Matrix Computations . Baltimore, MD, USA: The Johns Hopkins University Press, 2012.
- [34] S. Matthies, 'Orientations and rotations: Computations in crystallographic textures,' Journal of Applied Crystallography , vol. 38, no. 6, pp. 1042-1043, 2005.
- [35] D. G. Luenberger, Y. Ye, et al. , Linear and nonlinear programming , vol. 2. New York: Springer, 1984.
- [36] R. H. Byrd and J. Nocedal, 'A tool for the analysis of quasi-Newton methods with application to unconstrained minimization,' SIAM Journal on Numerical Analysis , vol. 26, no. 3, pp. 727-739, 1989.
- [37] J. Nocedal and S. Wright, Numerical optimization . New York, NY, USA: Springer, 2006.
- [38] D. Park, 'Weighted sum rate maximization of MIMO broadcast and interference channels with confidential messages,' IEEE Transactions on Wireless Communications , vol. 15, no. 3, pp. 1742-1753, 2015.
- [39] S. Boyd, S. P. Boyd, and L. Vandenberghe, Convex optimization . Cambridge, U.K.: Cambridge university press, 2004.
Xinliang Zhang Yue Qi Mojtaba Vaezi
A Rotation-based Method for Precoding in Gaussian MIMOME Channels
The problem of maximizing secrecy rate of multiple-input multiple-output
multiple-eavesdropper (MIMOME) channels with arbitrary numbers of antennas at
each node is studied in this paper. First, the optimization problem
corresponding to the secrecy capacity of the MIMOME channel is converted to an
equivalent optimization based on Givens rotations and eigenvalue decomposition
of the covariance matrix. In this new formulation, precoder is a rotation
matrix which results in a positive semi-definite (PSD) covariance matrix by
construction. This removes the PSD matrix constraint and makes the problem
easier to tackle. Next, a Broyden-Fletcher-Goldfarb-Shanno (BFGS)-based
algorithm is developed to find the rotation and power allocation parameters.
Further, the generalized singular value decomposition (GSVD)-based precoding is
used to initialize this algorithm. The proposed rotation-BFGS method provides
an efficient approach to find a near-optimal transmit strategy for the MIMOME
channel and outperforms various state-of-the-art analytical and numerical
methods. In particular, the rotation BFGS precoding achieves higher secrecy
rates than the celebrated GSVD precoding, with a reasonably higher
computational complexity. Extensive numerical results elaborate on the
effectiveness of the rotation-BFGS precoding. The new framework developed in
this paper can be applied to a variety of similar problems in the context of
multi-antenna channels with and without secrecy.