sensors

Article
RIS-Assisted Hybrid Beamforming and Connected User Vehicle
Localization for Millimeter Wave MIMO Systems
Md. Abdul Latif Sarker 1 , Woosung Son 2 and Dong Seog Han 3, *

1 Center for ICT & Automotive Convergence, Kyungpook National University, Daegu 41566, Republic of Korea
2 Graduate School of Electronic and Electrical Engineering, Kyungpook National University,
Daegu 41566, Republic of Korea
3 School of Electronics Engineering, Kyungpook National University, Daegu 41566, Republic of Korea
* Correspondence: dshan@knu.ac.kr; Tel.: +82-53-950-6609

Abstract: A reconfigurable intelligent surface (RIS) is a type of metasurface that can dynamically
control the reflection and transmission of electromagnetic waves, such as radio waves, by chang-
ing its physical properties. Recently, RISs have played an important role in intelligently reshaping
wireless propagation environments to improve the received signal gain as well as spectral efficiency
performance. In this paper, we consider a millimeter wave (mmWave) vehicle-to-vehicle (V2V)
multiple-input multiple-output (MIMO) system in which, an RIS is deployed to aid downlink V2V
data transmission. In particular, the line-of-sight path of the mmWave system is affected by blockages,
resulting in higher signaling overhead. Thus, the system performance may suffer due to interruptions
caused by static or mobile blockers, such as buildings, trees, vehicles, and pedestrians. In this paper,
we propose an RIS-assisted hybrid beamforming scheme for blockage-aware mmWave V2V MIMO
systems to increase communication service coverage. First, we propose a conjugate gradient and
location-based hybrid beamforming (CG-LHB) algorithm to solve the user sub-rate maximization
problem. We then propose a double-step iterative algorithm that utilizes an error covariance matrix
splitting method to minimize the effect of location error on the passive beamforming. The proposed
algorithms perform quite well when the channel uncertainty is smaller than 10%. Finally, the simula-
tion results validate the proposed CG-LHB algorithm in terms of the RIS-assisted equivalent channel
for mmWave V2V MIMO communications.
Citation: Sarker, M.A.L.; Son, W.;
Han, D.S. RIS-Assisted Hybrid Keywords: millimeter-wave MIMO; RIS; vehicle-to-vehicle communications; connected autonomous
Beamforming and Connected User vehicle; spectral efficiency; location error
Vehicle Localization for Millimeter
Wave MIMO Systems. Sensors 2023,
23, 3713. https://doi.org/10.3390/
s23073713
1. Introduction
Academic Editor: Changchuan Yin Millimeter wave (mmWave) communication is a promising technology that plays a ma-
Received: 28 February 2023
jor role in fifth-generation (5G) wireless communications owing to its wide bandwidth [1],
Revised: 28 March 2023 sufficient spectral-efficiency [2–5] and high quality of service (QoS) [6–10]. In highly di-
Accepted: 31 March 2023 rective wireless communications, beamforming is a signal processing tool used to steer,
Published: 3 April 2023 shape, and control multiple antenna directions efficiently, which minimizes transmission
losses [11–16]. The use of fully digital beamforming has been widely adopted in traditional
communication systems, but it requires a radio frequency (RF) chain per antenna element,
which increases power consumption, hardware complexity, and costs. As a solution to this
Copyright: © 2023 by the authors. problem, hybrid beamforming has been proposed as an alternative approach for mmWave
Licensee MDPI, Basel, Switzerland. massive multiple-input and multiple-output (MIMO) systems [17,18]. To further develop
This article is an open access article traditional hybrid beamforming, the proposed conjugate gradient-based hybrid beamform-
distributed under the terms and
ing algorithm can deal with the user sub-rate maximization problem and offer the benefit
conditions of the Creative Commons
of significantly improving spectral efficiency.
Attribution (CC BY) license (https://
The authors developed an adaptive algorithm in [2] to estimate the mmWave channel
creativecommons.org/licenses/by/
parameters that exploit the poor scattering nature of this channel. In [2], focus on both
4.0/).

Sensors 2023, 23, 3713. https://doi.org/10.3390/s23073713 https://www.mdpi.com/journal/sensors

Sensors 2023, 23, 3713 2 of 16

the single and multi-path channels to validate the adaptive channel estimation algorithm.
A generalized hybrid architecture was proposed with a small number of RF chains and a
finite number of ADC bits and the authors derived achievable rates with channel inver-
sion and singular value decomposition-based transmission methods. A content-centric
communications scenario was studied and a beamforming-based rate-splitting method
applied to maximize the sum rates of all users in [11]. In addition, joint beamforming and
power allocation was demonstrated for a non-orthogonal multiple access (NOMA)-based
satellite-terrestrial integrated network in [12]. In [13], a receiver prototype was presented
that enables reconfiguration between analog and digital beamforming modes to better deal
with different operating scenarios. The authors designed a true-time-delay analog and
digital beamformer in [13]. Elsewhere, a network was proposed featuring a multi-antenna
array base station (BS) and a reconfigurable intelligent surface (RIS) to deliver information
to information users and power to energy users [14]. Conjugate and zero-forcing beam-
forming was proposed to develop a path-following algorithm to solve the user throughput
maximization problem subject to a realistic constraint on the quality-of-energy-service in
terms of the harvested energy thresholds in [14]. In [15], the authors investigated a secure
transmission design for an intelligent reflecting surface (IRS)-assisted unmanned aerial
vehicles network and solved active and passive beamforming and trajectory optimiza-
tions. In addition, a neural networks-based beam codebook was designed in [16]. There
are three hybrid beamforming architectures investigated in the literature: fully connected,
fixed subarray and dynamic subarray. In [17], an iterative hybrid beamforming design
algorithm was proposed with a fully connected architecture to maximize the weighted
sum-rate performance of mmWave MIMO systems. The authors considered an RIS-aided
mmWave-NOMA downlink system with a hybrid beamforming structure in [18]. They
developed an alternating optimization algorithm to solve the non-convex optimization is-
sue. Finally, geometric mean decomposition-based beamforming for RIS-assisted mmWave
hybrid MIMO systems was proposed in [19]. The authors designed phase shifters for RIS
by maximizing the array gain for the line-of-sight (LoS) channel.
Recently, mmWave communication has included a vehicle-to-everything (V2X) paradigm
to facilitate connected autonomous driving, video entertainment, vehicle platooning, and
dynamic map update [20]. All of these services require high-speed connections to ex-
change safety information (e.g., vehicle location, direction and speed) among adjacent
vehicles since recent V2X standards such as cellular V2X in [21] , do not have advanced
QoS and cannot provide consistent high-speed data transmission. Owing to the sparsity of
the mmWave channel [1], the MIMO system typically requires an LoS connection, which
maintains the received power level. Indeed, the LoS channel between the transmitter and
receiver may be blocked by static or moving objects, such as buildings, trees, vehicles, and
humans [22]. To satisfy the service quality of V2X communication under LoS blockage, RISs
have also been represented as large intelligent surfaces (LISs), which have been investigated
as an emerging technology in [6,7,23]. An RIS controller with beamforming is typically
offered with high-accuracy indoor or outdoor radio localization [24]. In mmWave MIMO
communication, the RIS can easily establish robust connections subject to blockage or
mobility awareness. RISs are mainly used as antenna arrays for increasing coverage area
and reducing energy consumption [25]. Moreover, RISs are applied as reflective relays [26].
The functionalities of RIS controllers increase the quality of the received signal when the
adjustable passive elements can individually steer an incident electromagnetic wave in
any direction.
Based on the RIS element arrays, beamforming and steering concepts have been
investigated in [23–25,27]. Due to the notable benefits of RIS technology, outdoor cov-
erage was analyzed using RISs in [25] and IRS-aided beamforming concepts discussed
in [27]. A channel estimation algorithm was proposed in [28,29], an RIS-aided manifold
optimization algorithm was developed in [30], and an RIS-aided real-time autonomous
vehicle model was investigated in [31]. Elsewhere, researchers proposed ultra-reliable
low-latency communication in a factory automation (FA) scenario to deploy an IRS to
Sensors 2023, 23, 3713 3 of 16

create an alternative transmission link, which can enhance transmission reliability [32]. The
performance of single and multiple RIS-assisted systems has also been considered without
a direct path between the transmitter and the receiver in indoor and outdoor propagation
environments [33]. In addition, researchers proposed a geometry-based stochastic channel
model that supports the movements of transceivers and clusters, and the evolution of
clusters was considered in the space domain, where a reflecting coefficients design was
based on the minimum path loss [34]. Elsewhere, researchers provided a comprehensive
overview of state-of-the-art research on emerging and promising RIS/IRS-aided wireless
systems, with an emphasis on signal processing techniques for solving various radio lo-
calization, transmission design, and channel estimation issues [35]. A tutorial overview of
single and multi-IRS-aided wireless networks has also been provided, with an emphasis on
addressing the new and more challenging issues in IRS reflection optimization and channel
acquisition design, in [36]. In [37], a multi-IRS-aided system was studied in which, the IRS
and base station (BS) are managed by a central processing unit to coordinate data transmis-
sion and maximize the weighted sum rate of all cell-edge users by jointly optimizing the
BS’s transmit beamforming and each IRS’s phase shifts, subject to the BS transmit power
budget. The authors solved the non-convex and unit modulus constraint optimization
using an efficient iterative power allocation algorithm. The achievable secrecy performance
of mmWave MIMO systems was studied using RISs in [38], where the authors assumed a
smart environment in which an RIS is placed between the source and the legitimate user
to enhance the main link. The authors, motivated by the aforementioned observations,
studied beamforming optimization for an IRS-aided multi-antenna communication system
by incorporating signal distortions caused by hardware impairments [39]. The authors
optimized both source’s transmit active and passive beamforming to maximize the SNR
received at the destination. Elsewhere, the application of an active refracting RIS-enabled
transmitter for a secure internet of things network was investigated to enhance secure
communication in the considered network and develop an alternating optimization algo-
rithm to optimize the sum secrecy rate by jointly designing the power allocation, transmit
beamforming, and the phase shifts of an RIS in [40].
Due to the high beam training costs and the effect of user location errors in vehicular
MIMO communication systems, we propose an RIS-assisted location-based hybrid beam-
forming algorithm as well as a double-step iterative algorithm for minimizing the effect
of user vehicle location error, which has not yet been investigated in the literature. The
contributions of this study are listed as follows:
• Due to blockage awareness, we first develop an RIS-assisted V2V MIMO channel
model and demonstrate the geometrical relation between the RIS controller and
user vehicle position. By using the distance between the transmitting and receiving
vehicle in relation to the distance of an RIS controller, we estimate the path amplitude
and phase.
• We then design an RIS-assisted low-dimensional equivalent effective channel. In re-
sponse to a sub-design problem of an original beamformer, we propose a conjugate
gradient and location-based hybrid beamforming algorithm. We apply a Karush-
Kuhn-Tucker condition and a Lagrangian method to solve an original problem into a
sub-optimal problem for developing the conventional hybrid beamforming algorithm,
which reduces the user sub-rate maximization problem. The proposed algorithm
attains significant spectral efficiency performance.
• We next consider a covariance splitting method to reform the error covariance matrix.
Since the channel is geometrically contained in location information, we propose a
double-step iterative algorithm that minimizes the effect of location error in user
vehicles. Finally, the simulation results demonstrate the superiority of the proposed
algorithms over their counterparts.
Sensors 2023, 23, 3713 4 of 16

The rest of the paper is organized as follows. The signal and channel model and the
geometrical relation between the RIS controller and connected autonomous vehicle position
are discussed in Section 2. The proposed hybrid beamforming and user vehicle localization
is presented in Section 3. Simulation results are provided in Section 4. The conclusions are
offered in Section 5.

2. Signal and Channel Model

2.1. Signal Model
Consider an access point that serves a point-to-point (P2P) vehicular MIMO network
using mmWave technology, as shown in Figure 1, where a transmit vehicle is operating at
mmWave frequency bands with Nt transmit antennas. The user vehicle is equipped with
Nu receive antennas. Let the transmit and user vehicles be autonomous, which are assumed
to be fully connected with different modules and sensors as the power module is used as
a battery to supply power to the other module. The perception module is recognized for
driving environments and detected objects using sensors such as light detection and ranging
(LiDAR), radars, cameras, a global positioning system (GPS), and inertial measurement
unit (IMU), and other sensors as shown in Figure 2. A WiFi-based mmWave transceiver is
also intelligently connected to the transmit and receive vehicles in Figure 2. We consider an
IEEE 802.11ad Wi-Fi chips with the number of 8 antennas to communicate at the 60 GHz
frequency band.
A blockage scenario is considered in the system model where the direct link between
the transmit and user vehicles is blocked owing to the critical propagation conditions.
Hence, the RIS is deployed in the transmission model to assist and expand the service
coverage of the P2P vehicular MIMO communication. The received signal vector r ∈ C Nu ×1
can be modeled as
r = He x + z, (1)
where x ∈ C Nt ×1 is the transmitted signal vector, which satisfies E[xx H ] = I, He is the
Nu × Nt effective channel from the transmitting vehicle to the user vehicle, z ∼ CN (0, σz2 I)
is the Nu × 1 additive white Gaussian noise vector and σz2 denotes the noise variance.
After applying a hybrid precoding, x is given by

x = F A F D s, (2)

where F A ∈ C Nt × NRF is the analog radio frequency (RF) beamforming matrix,

F D ∈ C NRF × Ns is the digital precoding matrix, ||F A F D ||2F ≤ P, P is the transmit power,
s ∈ C Ns ×1 is the transmitted symbol vector that satisfies E[ss H ] = I, Ns ≤ min( NRF
t , Nu )
RF
when Nt >> NRF t , N >> N u , N t denotes the number of RF chains for the transmit
u RF RF
vehicle side, and NRF u denotes the number of RF chains for the user intelligent vehicle side.

Setting (2) in (1) and the received signal vector ŝ is given after combining

ŝ = H̃e s + z̃, (3)

where H̃e = W D H W H H F F is the N × N effective equivalent channel, z̃ = W H W H z is

A e A D s s D A
the Ns × 1 equivalent noise vector, W D ∈ C NRF × Ns is the baseband digital combiner, and
W A ∈ C Nu × NRF is the analog combiner which follows the similar constraints as F A . If the
transmitted signal follows a Gaussian distribution over the RIS-assisted mmWave V2V
MIMO channel, then using (3), the achievable spectral efficiency is given by

R = log2 det(I + E−1 H̃e H̃eH ), (4)

where E = σz2 W D H W H W W is the noise covariance matrix after combining and the
A A D
operator (·) H denotes the Hermitian.
Sensors 2023, 23, 3713 5 of 16

Figure 1. System model is illustrated by RIS-assisted V2V communications.

Figure 2. Sensors connected autonomous vehicle with WiFi-based mmWave transceiver.

2.2. Channel Model

Let He = GΘC where G is an Nu × NR channel between RIS to the user vehicle as
Lg
G = γg ∑ β̃l au (φul )aRH (θtl , νtl ), (5)
l =1

where β̃ l denotes the complex channel gain at l-th path l

p of G channel, φul represents the
angle of arrival associated with the user vehicle, γg = Nu NR /L g , and θt and νtl denotes
the azimuth and elevation angles of arrival and departure at l-th path, respectively, the
normalized array response vector for the case of the uniform linear array is given by [2–5].

l) l) T
 
1 j2πd sin (φu j2π ( Nt −1)d sin (φu
au (φul ) = √ 1, e λ ,...,e λ , (6)
Nt
Sensors 2023, 23, 3713 6 of 16

where the wavelength λ = c/ f c , c is the speed of light, f c is the carrier frequency, and
d = λ/2 is the antenna spacing. The normalized array response vector for uniform planner
array is given by [28]

1 j2πd l l l
a R (θtl , νtl ) = √ [1, e λ ((n1 −1) cos (νt ) sin (θt )+(n2 −1) sin (νt )) ,
Nt (7)
j2πd l l l
...e λ {( NR,y −1) cos ( νt ) sin ( θt )+( NR,z −1) sin ( νt )} T
] ,

where NR = NR,y × NR,z , and C denotes the NR × Nt channel between the transmit vehicle
and the RIS controller, that is
Lc
C = γc ∑ γ̃l aR (θrl , νrl )atH (φtl ), (8)
l =1
√
where γ̃l denotes the complex channel gain at l-th path of C channel, γc = Nt NR /Lc ,
a R (θrl , νrl ) and atH (φtl ) are designed in the same manner of G channel and the NR × NR RIS
element response matrix Θ is given by

1
Θ= √ D, (9)
NR

where D is the diagonal phase matrix as [25],

D = diag{ϑ1 , . . . , ϑ NR }, (10)

where ϑ = [ϑ1 , . . . , ϑ NR ] T , ϑn = e jϕn , n = 1, 2, . . . , NR and ϕn denotes the phase-shift which

is given by
ϕ1,1 · · · ϕ1,NR
 

ϕ =  ... .. .. . (11)

. . 
ϕ NR ,1 ··· ϕ NR ,NR
Thus, the effective channel He can be written as

Lc Lg
He = ∑ ∑ β l γm µl,m au (φRl )atH (θtl ), (12)
m =1 l =1

where µl,m = a RH (θtl , νtl )Θa R (θrm , νrm ) denotes the passive beamforming gain, β l = γg β̃ l ,
γl = γc γ̃l , and µl,m satisfy |γ1 | ≥ |γ2 |, . . . , ≥ |γLc | and | β 1 | ≥ | β 2 |, . . . , ≥ | β L g |.

2.3. Geometrical Relation between the RIS Controller and the Connected Vehicle Position
Let b = (bx , by , bz ), r = (r x , ry , rz ), and u = (u x , uy , uz ) are the position of the transmit-
ting vehicle, the RIS controller and the user vehicle, respectively in the Cartesian coordinate
as shown in Figure 1, and g is the vector of G channel. The g is given by

g = vec(G) = a R (θtaz , νtel ), (13)

where θ az denotes the azimuth of the angle of departure (AoD) and νel denotes the elevation
of AoD at the RIS side. The geometric relation between the RIS controller and the user
vehicle position is measured by
Sensors 2023, 23, 3713 7 of 16

 
 ux − rx 
θtaz = arcsin q , (14)
 ( u x − r x )2 + ( u y − r y )2 

and
uz − rz
 
νtel = arcsin . (15)
||u − r||2
Similarly, we can make a positional relationship between the transmitting vehicle and
the RIS controller as
 
 r x − b x

θraz = arcsin q , (16)
 (r − b )2 + (r − b )2 
x x y y

and
r z − bz
 
νrel = arcsin . (17)
||r − b||2
Using (13), the passive beamforming gain µl,m is computed as

µ = g H Θg = [g g∗ ] H ϑ, (18)

where Θ = diag{ϑ } and the operator (·)∗ denotes the conjugate, and denotes the
element-wise product. The optimal phase-shift ϑ opt is given by

ϑ opt = arg max |[g g∗ ] H ϑ | (19)

ϑ

3. Proposed CG-LHB Algorithm and User Vehicle Localization

3.1. Proposed CG-LHB Algorithm
Let Ns = NRF , F = F A F D ∈ C Nt × Ns , and W = W A W D ∈ C Nu × Ns . The objective func-
tion is to maximize the rate of the user vehicle while satisfying the maximum transmission
power constraint. Using (4) and (12), the sub-rate maximization problem is written as

Ns  
P
∑ log2
H
max f (F, W, ϑ ) = 1+ H̃e H̃e , s.t. ||F||2F ≤ P,
F,W,ϑ i
Ns σz2 (20)
2 2 2
|F(i, j)| = |W(i, j)| = 1, ∀i,j , and |ϑn | = 1, ∀m=1,...,NR ,

where i = 1, 2, . . . , Ns . The sub-rate maximization problem (20) is a sub-problem with

respect to the digital beamformer F D . To maximize the sub-rate, F D should be optimized to
minimize ||F A f D,i ||2 based on the following multi-object optimization problem (20) as

min ||F A f D,i ||2 , i = 1, . . . , Ns ,

f D,i
(21)
s.t. He F A f D,i = I,

where F A , W A and ϑ are fixed. We see from (21), the Ns column of F D can be optimized
separately according to

min ||F A f D,i ||2 , i = 1, . . . , Ns ,

f D,i
(22)
s.t. He F A f D,i = ςi ,
Sensors 2023, 23, 3713 8 of 16

where ςi denotes an all-zero vector except for the i-th entry being unity. Using the stream-
specific error expressions in (22), we apply Karush-Kuhn-Tucker (KKT) conditions and the
Lagrangian method for the convex problem as

∂ f (f D,i , κ ) H
= (f D,i Af D,i ) T + (κ T He F A f D,i ) T , (23)
∂f D,i

where A = F H
A F A denotes the equivalent Ns × Ns analog beamforming matrix, κ is the
Lagrangian multiplier and the corresponding Lagrangian is
H
f (f D,i , κ ) = f D,i Af D,i + κ T (He F A f D,i − ςi ). (24)

If ∂ f f D,i , κ )/∂f D,i = 0, then f D,i and κ is given by

f D,i = −A−1 F H H
A He κ, (25)

and
κ = −[He F A A−1 F H H −1
A H e ] ςi . (26)
Thus, combining (25) and (26), the digital beamforming vector f D,i is given by

f D,i = A−1 F H H −1 H H −1
A H e [ H e F A A F A H e ] ςi . (27)

Since F A is a constant modulus in (27), we optimize the phase of the analog beamformer
for transforming a complex problem into an unconstrained optimization problem. Let the
beamforming design problem can be written by
opt opt
(F A , F D ) = arg min ||Fopt − F A F D || F , s.t. F A ∈ F A , (28)
F A ,F D

where F A is the set of the analog beamformer, Fopt = V1 , which is measured by using a
singular value decomposition method of He = UΓV H in [41]. Consider F A into F A,m as

1
F A,m = √ e jΩm , (29)
Nt

where Ωm denotes the Nt × Ns phase shift matrix, which phases need to be quantized into
the nearest point in (2πm/2b , m = 0, 1, . . . , 2b − 1), if the phase shifters are of infinite or
finite b-bit resolution. Hence, the phase optimization problem can be expressed as
opt
Ωm = arg min f (Ωm ). (30)
Ωm

where Ωm is converged using a conjugate gradient method, which optimizes along the
direction ∆Ωm
Ωm+1 = Ωm + ∆Ωm , (31)
where ∆Ωm = ηm Ξm and ∆Ωm satisfy ∆Ωm ≤ τ, ηm is determined by minimizing f (Ωm+1 ),
Ξm is given by
Ξm = −∇ f (Ωm ) + ξΞm−1 , (32)
and the ξ is given by
∇ f (Ωm ) T ∇ f (Ωm )
ξ= . (33)
∇ f ( Ω m −1 ) T ∇ f ( Ω m −1 )
Sensors 2023, 23, 3713 9 of 16

Using (29) to (33), we summarize the proposed CG-LHB algorithm in Algorithm 1.

Algorithm 1 Proposed CG-LHB Algorithm

1: Input parameters: H, b, τ, J
opt opt
2: Output: F A , F D .
3: Begin: Set m = 0, and Ω0 ∈ R Nt × Ns
4: If m = 0, 1, . . . , ∞,
Ξm = −∇ f (Ωm )
else compute (31)–(33).
5: If ηm is determined by minimizing f (Ωm+1 ) in (31), then
Ω s ( m + 1 ) ← f ( Ω m +1 );
f s ( m + 1 ) = f ( Ω m +1 ).
6: ∆ f = f s ( m + 1) − f s ( m );
∆Ω = Ωs (m + 1) − Ωs (m),
If ∆ f ≤ τ1 ; stop
If ∆Ω T ∆Ω ≤ τ2 ; stop
If m + 1 = J; stop
If ∇ f (Ωm+1 ) T ∇ f (Ωm+1 ) ≤ τ3 ; converged
m ← m+1
Go To Step 4.
7: if b < ∞ then
8: Ωopt = f (Ωopt )
opt
9: Compute F A
opt
10: Compute F D .

3.2. User Vehicle Localization

Let RIS to user channel covariance Σi = E[gi giH ], where i = 1, . . . , Ns } and the channel
vector g = vec(G) is geometrically generated for the different positions as shown in
Section 2.3. Thus, the user vehicle localization problem can be computed as a hypothesis
testing problem as follows:
Hi : g ∼ CN (0, Σi ). (34)
Hence, the optimal likelihood (ML)-based localization problem is given by

î = arg max pi (g), (35)

1≤i ≤ Ns

where pi (g) is the probability density function under G

1 H −1
pi ( g ) = e − g Σi ( α ) g . (36)
π Nu ,Nt |Σi (α)|

We can now evaluate the error-probability performance under Hi , which is given by

Ns
1
Pe =
Ns ∑ Pe,i , Pe,i = P(î 6= i | Hi ). (37)
i =1

where Pe,i is the conditional error probability under Hi .

Due to user mobility and high location errors, (36) contains a high channel state infor-
mation (CSI) error. To overcome this issue, we analyze the location information generated
by GPS. Consider the location of the user vehicle from the GPS represented by

u = û + ∆u, (38)
Sensors 2023, 23, 3713 10 of 16

where u = [u x , uy , uz ] T denotes the actual location vector of the user vehicle, û = [û x , ûy , ûz ] T
is the ego GPS location vector, and ∆u = [∆u x , ∆uy , ∆uz ] T is the localization error vector
and ∆u satisfy ||∆u||2 ≤ e. Based on (13) and (38), we obtain the actual channel vector g
that satisfy
g(u) = ĝ(u) + ∆g(u), (39)
where ∆g ∼ N (0, Σ) and Σ = [σg,x 2 , σ2 , σ2 ] T is approximated as diagonal with σ2 , σ2
g,y g,z g,x g,y
2 being the channel variances along the x, y, and z directions, respectively. Hence,
and σg,z
by using (39), with the location error bound (36), the CSI error minimization problem can
be formulated as

min f (Θ, ∆g) = ||ĝ(u) + ∆g(u)||2 , s.t. ||Θi,i || = 1, ∀i=1,...,NR ,

∆g
(40)
|a R (θ Rl , νRl )i | = 1, ∀i=1,...,NR , and ||∆g||2 ≤ e.

It is noted that (40) still suffers high localization errors. This issue can be addressed
using a double-step iteration method. Let the error covariance matrix Σ possesses and split
two different covariance matrices Σ = Σ1 + Σ2 in (34), where Σ1 = F∗ Λ1 F, Σ2 = F̂∗ Λ2 F̂,
Λ1 and Λ2 are the diagonal matrices holding the eigenvalues of Σ1 and Σ2 , respectively. If
the matrices Σ1 and Σ2 are positive definite and α is a positive constant, then the iterative
covariance matrix Σ(α) is given using by [42]

Σ(α) = M2−1 N2 M1−1 N1 , (41)

where α denotes the positive constant, M2 = (αI + Σ2 ), N2 = (αI − Σ1 ), M1 = (αI + Σ1 )

and N1 = (αI − Σ2 ). Based on (34), we demonstrate the proposed double-step iterative al-
gorithm for the user vehicle localization in Algorithm 2. The double-step at each Σ1 Σ2 Σ1 Σ2
iterations is necessary to obtain an exact solution for the αI + Σ1 and αI + Σ2 matrices, which
strength diagonal properties of the Σ1 and Σ2 matrices. The complexity of the proposed
Algorithm 2 is O( NR log2 NR ) instead of the original complexity O( NR3 ).

Algorithm 2 Proposed Double-Step Iterative Algorithm for User Vehicle Localization

1: Input parameters: k, α
2: Output: ρ(Σ(α)).
3: If x (k) is a double-step sequence, then
2k+1
x ( 2 ) = M1−1 N1 x (k) + M1−1 b, and
2k+1
x (k+1) = M2−1 N2 x ( 2 ) + M2−1 b, where
Σ ∈ C Nv × Nv , Σ = Mq − Nq (q = 1, 2)
and Σx = b for all initial vectors x (0) ∈ C Nu .
4: If k = 0, 1, . . . , then compute
x (k+1) = Σ(α) x (k) + ¸(α)b, where
Σ(α) = M2−1 N2 M1−1 N1 ,
¸(α) = M2−1 (I + N2 M1−1 ), and
ρ(Σ(α)) ≤ σ (α) < 1, ∀α > 0.
5: When α 6= αmax ,
Repeat step 4,
6: When α = αmax ,
Stop convergence, i.e., ρ(Σ(α)) = ρ(I(αmax )).

4. Simulation Results and Discussion

In this section, we compare the proposed RIS-assisted CG-LHB algorithm against the
conventional hybrid beamforming (G-HB) algorithm via computer simulations. Throughout
the simulations, we considered the simulation parameters of Table 1. The signal-to-noise-
Sensors 2023, 23, 3713 11 of 16

ratio (SNR) at the receiver is defined as Γ = P

Ns σz2
||Υ||2F , where Υ = GΘCF A F D , P = 1
watt, σz2= −174 + 10log10 B = −90 dBm and the bandwidth B = 251.188 MHz [28]. The
azimuth angle of departure (AoD) and angle of arrival (AoA) are uniformly distributed
in the (−π, π ) and the elevation AoD and AoA are uniformly distributed in the interval
(−π/2, π/2) [5,41]. Let the coordinates of the transmit intelligent vehicle, the RIS and
the user vehicle are ( xt , yt ) = (0, 0)m, ( x R , y R ) = (110, 20)m, and ( xu , yu ) = (210, 10)m,
respectively. We set the phase-shifters resolution b = 2 bits and the algorithmic parameter
τ = 10−6 , J = 100, ξ = 0.6 and α = 3 [42] in the computer simulations.

Table 1. Simulation Parameters.

Total number of transmit vehicle antennas Nt = 8

Total number of user vehicle antenna Nu = 8
Total number of RIS discrete elements NR = 32( NR,y = 8, NR,z = 4)
Total number of RF chains NRF = 4
Total number of data streams Ns = 4
The number of propagation path 1
Wavelength λ = 1 cm
Light speed c = 3 × 108 m/s
RIS elements distance d = 0.5 cm
Carrier frequency f c = 28 GHz
RIS dimensions 32 × 32
Noise PSD N0 = −174 dBm/Hz
Noise variance σz2 = N0 + 10log10 B

For the line-of-sight channel (l = 0 case), the RIS-assisted channel vector g is given using

NR β̃ 0 $t $r a R (θtaz , νtel ),
p
g= (42)
where $t and $r are the transmit and receive antenna gains, β̃ 0 ∼ CN (0, 10−0.1PL(d2) )
denotes the channel path gain, PL(d2) in the path loss in the RIS to user link and PL(d2) is
defined as [27,43,44]
PL(d2) = 35.6 + 22.2log10 (d2 ) + X , (43)
where d2 denotes the distance between the RIS and user vehicle as shown in Figure 1, and
X ∼ N (0, σX2 ). The parameter σX is set at σX = 5.8 dB, respectively. Similarly, we also can
measure the corresponding channel C and path loss PL(d1 ) according to (42) and (43), and
d1 represents the distance from the transmitting vehicle to the RIS controller.
Figure 3 depicts the received signal power gain versus SNR. We see that the gain
performance of the received signal power is showing the same in the high SNR case. For
the low SNR case, the gain performance of the received signal power is around 0.30 dB as
shown in Figure 3. Figure 4 shows the achievable spectral efficiency versus transmit power
using user location information. The proposed CG-LHB algorithm provides a significant
increase in the spectral efficiency for a 32 × 32 reflecting array compared to the traditional
RIS-assisted hybrid beamforming [25,27] and with and without RIS-assisted hybrid beam-
forming method [3]. The achievable spectral efficiency of the proposed CG-LHB algorithm
is about 1.30 bits/s/Hz, the conventional RIS-assisted hybrid beamforming method is
around 1.01 bits/s/Hz and the hybrid beamforming without RIS is about 0.90 in Figure 4.
For the proposed CG-LHB algorithm case, the achievable spectral efficiency performance
is increased by approximately 0.301 bits/s/Hz/user at 1 watt transmitted power values.
Figure 5 illustrates the achievable spectral efficiency versus NR by evaluating Algorithm 1.
we also plot the spectral efficiency of the conventional hybrid beamforming with the G-HB
algorithm [25,45]. The performance of the spectral efficiency is around 10% at NR = 32.
Figure 6 illustrates the convergence of the proposed algorithm with two gradient
descent methods as in Algorithm 1 at ξ = 0.6. By measurement of ∇ f (Ω) in (32), the
computational complexity of optimizing F A is dominated and it takes O( Nt NRF Ns ) multi-
plications. It is noted that the total algorithm needs O( Nt NRF Ns J ) multiplications, where
Sensors 2023, 23, 3713 12 of 16

J is the number of iterations required to converge with a fixed value of parameter τ and
J < 100. Instead of the conventional gradient method (Algorithm 3 in [46]), which re-
quired O( Nt2 NRF
2 J ) multiplications, the proposed algorithm shows the computationally

more efficient.

Receive Signal Power Gain (|| ||2F ) 4500

4000

3500

3000

2500

2000

Ideal case
Conventional RIS-assisted HB algorithm
1500
Proposed RIS-assistd CG-LHB algorithm

-10 -5 0 5 10 15 20 25 30
SNR

Figure 3. Received power gain versus SNR with the proposed CG−LHB algorithm.

4.5
Conventional HB algorithm with RIS
Achievable spectral efficiency (bits/s/Hz)

4 Prposed CG-LHB algorithm with RIS

Conventional HB algorithm without IRS
3.5

2.5

1.5

0.5

0
0 5 10 15 20 25 30
P (dBm)

Figure 4. Spectral efficiency versus transmit power for the proposed RIS−assisted CG−LHB algorithm.

5.5

5
Achievable spectral efficiency (bps/Hz)

4.5

3.5

2.5

1.5 Proposed CG-LHB algorithm with RIS

1
Conventional G-HB algorithm without RIS
Conventional CG-HB algorithm with RIS
0.5
4 8 12 16 20 32
NR

Figure 5. Spectral efficiency versus number of RIS elements (NR ) for the proposed RIS-assisted
CG−LHB algorithm.
Sensors 2023, 23, 3713 13 of 16

-6 Conjugate Gradient-based analog optimization

Analog beamforming optimization

Gradient Descent-based analog optimization
-8

-10

-12

-14

-16

-18

-20
0 20 40 60 80 100 120 140 160 180 200
Iterations number
Figure 6. Convergence of the proposed CG algorithm−based analog beamforming optimization.

In Figure 7, we plot the error probability versus SNR, which relies on the equivalent ef-
fective channel. The effective equivalent channel is geometrically generated for the different
object locations. To achieve a better error-probability performance, we applied a channel
covariance splitting method and proposed a double-step iterative algorithm. In computer
simulation, we evaluate Pe as a function of SNR, where SNR leads the location informa-
tion subject to the channel path loss, which progressively adds extra scatter positions to
previous positions. From Figure 6, we observe that the proposed double-step algorithm
provides higher accuracy than the conventional iterative algorithm [26] in the case of a
multi-antenna user vehicle. Figure 8 compares the average bit-error-rate versus SNR of a
RIS-assisted V2V MIMO system at Ns = 4, Nt = 8 Nu = 8, and NR = 256. We consider
the proposed double-step iterative and conventional iterative algorithms. We observe that
the proposed double-step algorithm significantly outperforms the conventional iterative
algorithm. In Algorithm 2, we used α = 3 to fix the spectral of the covariance matrices. We
also observe that the average bit-error rate is improved by about 1.01 dB at the average
bit-error rate of 0.0001. Figure 9 shows the channel state information (CSI) error bound
in terms of location error e. We used the location error parameter e = 3 m to validate the
effect of the location error on the passive beamforming. If the number of NR is increased,
the CSI error bound e will vary as shown in Figure 9.

100
Proposed double-step iterative algorithm
Conventional iterative algorithm
Error-probability (Pe)

10-1

10-2
0 2 4 6 8 10 12 14 16 18 20
SNR [dB]

Figure 7. The error−probability performance of the V2V communications at 20 dB SNR.

Sensors 2023, 23, 3713 14 of 16

100

10-1

Average bit-error-rate 10-2

10-3
Conventional iterative algorithm
Proposed double-step iterative algorithm
10-4
-40 -38 -36 -34 -32 -30 -28 -26 -24 -22 -20
SNR in [dB]
Figure 8. Average bit−error−rate versus SNR in a RIS system at Ns = 4, Nt = 8 NR = 256.

3
Conventional iterative algorithm when =3 m
Proposed double-step iterative algorithm when =3 m
2.5

2
Error bound,

1.5

0.5

0
1 5 10 15 20 25
Number of RIS elements (N R )

Figure 9. Channel state information error bound versus the number of RIS elements.

5. Conclusions
In this paper, we proposed RIS-assisted CG-LHB and a user vehicle localization
algorithm. Utilizing the CG-LHB algorithm, we were able to significantly improve spec-
tral efficiency in mmWave MIMO systems. In addition, for user vehicle localization, we
considered a channel covariance splitting method and proposed a double-step iterative
algorithm that reduces the effect of location error on the passive beamforming. We validated
the effectiveness of the proposed algorithms using the conventional algorithms. Hence,
the proposed RIS-assisted CG-LHB algorithm can be extended further to the machine
learning-based beam alignment and location error minimization solution of the connected
autonomous vehicles, which will be explored in future studies.

Author Contributions: M.A.L.S. made contributed to the conception and design, acquisition of data,
analysis, and interpretation of data. W.S. contributed to data collection. D.S.H. supervised the work
and helped in the implementation part. All authors have read and agreed to the published version of
the manuscript.
Funding: This work was supported by the Basic Science Research Program through the National
Research Foundation of Korea (NRF) funded by the Ministry of Education (2021R1A6A1A03043144).
Institutional Review Board Statement: Not Applicable.
Informed Consent Statement: Not Applicable.
Sensors 2023, 23, 3713 15 of 16

Data Availability Statement: Not Applicable.

Conflicts of Interest: The authors declare no conflict of interest.

References
1. Rangan, S.; Rappaport, T.S.; Erkip, E. Millimeter-Wave Cellular Wireless Networks: Potentials and Challenges. Proc. IEEE 2014,
102, 366–385. [CrossRef]
2. Alkhateeb, A.; Ayach, O.; Leus, G.; Heath, R.W. Channel Estimation and Hybrid Precoding for Millimeter Wave Cellular Systems.
IEEE J. Sel. Top. Signal Process. 2014, 8, 831–846. [CrossRef]
3. Mo, J.; Alkhateeb, A.; Abu-Surra, S.; Heath, R.W. Hybrid Architectures with Few-Bit ADC Receivers: Achievable Rates and
Energy-Rate Tradeoffs. IEEE Trans. Wirel. Commun. 2017, 16, 2274–2287. [CrossRef]
4. Choi, J.; Lee, G.; Alkhateeb, A.; Gatherer, A.; Al-Dhahir, N.; Evans, B.L. Advanced Receiver Architectures for Millimeter-Wave
Communications with Low-Resolution ADCs. IEEE Commun. Mag. 2020, 58, 42–48. [CrossRef]
5. Sarker, M.A.L.; Kader, M.F.; Han, D.S. Rate-Loss Mitigation for a Millimeter-Wave Beamspace MIMO Lens Antenna Array System
Using a Hybrid Beam Selection Scheme. IEEE Syst. J. 2020, 14, 3582–3585. [CrossRef]
6. Chen, Y.; Wang, Y.; Zhang, J.; Renzo, M.D. QoS-Driven Spectrum Sharing for Reconfigurable Intelligent Surfaces (RISs) Aided
Vehicular Networks. IEEE Trans. Wirel. Commun. 2021, 20, 5969–5985. [CrossRef]
7. Jiang, Z.; Rihan, M.; Zhang, P.; Huang, L.; Zhang, J.; Mohamed, E.M. Intelligent Reflecting Surface Aided Dual-Function Radar
and Communication System. IEEE Syst. J. 2022, 16, 475–486. [CrossRef]
8. Li, J.; Huang, Y.; Chen, P.; Qi, X.; Wen, M.; Arslan, H. Composite Multiple-Mode Orthogonal Frequency Division Multiplexing
with Index Modulation. IEEE Trans. Wirel. Commun. 2022. Early Access. . [CrossRef]
9. Li, J.; Dang, S.; Yan, Y.; Peng, Y.; Al-Rubaye, S.; Tsourdos, A. Generalized Quadrature Spatial Modulation and its Application to
Vehicular Networks with NOMA. IEEE Trans. Intell. Transp. Syst. 2021, 22, 4030–4039. [CrossRef]
10. Gao, Y.; Liu, X.; Li, J.; Fang, Z.; Jiang, X.; Huq, K.M.S. LFT-Net: Local Feature Transformer Network for Point Clouds Analysis.
IEEE Trans. Intell. Transp. Syst. 2022, 24, 2158–2168. [CrossRef]
11. Lin, Z.; Lin, M.; de Cola, T.; Wang, J.B.; Zhu, W.P.; Cheng, J. Supporting IoT with Rate-Splitting Multiple Access in Satellite and
Aerial-Integrated Networks. IEEE Internet Things J. 2021, 8, 11123–11134. [CrossRef]
12. Lin, Z.; Lin, M.; Wang, J.B.; de Cola, T.; Wang, J. Joint Beamforming and Power Allocation for Satellite-Terrestrial Integrated
Networks with Non-Orthogonal Multiple Access. IEEE J. Sel. Top. Signal Process. 2019, 13, 657–670. [CrossRef]
13. Spoof, K.; Tenhunen, M.; Unnikrishnan, V.; Stadius, K.; Kosunen, M.; Ryynänen, J. True-Time-Delay Receiver IC with Reconfig-
urable Analog and Digital Beamforming. IEEE Access 2022, 10, 116375–116383. [CrossRef]
14. Yu, H.; Tuan, H.D.; Dutkiewicz, E.; Poor, H.V.; Hanzo, L. RIS-Aided Zero-Forcing and Regularized Zero-Forcing Beamforming in
Integrated Information and Energy Delivery. IEEE Trans. Wirel. Commun. 2022, 21, 5500–5513. [CrossRef]
15. Pang, X.; Zhao, N.; Tang, J.; Wu, C.; Niyato, D.; Wong, K. IRS-Assisted Secure UAV Transmission via Joint Trajectory and
Beamforming Design. IEEE Trans. Commun. 2022, 70, 1140–1152. [CrossRef]
16. Alrabeiah, M.; Zhang, Y.; Alkhateeb, A. Neural Networks Based Beam Codebooks: Learning mmWave Massive MIMO Beams
That Adapt to Deployment and Hardware. IEEE Trans. Commun. 2022, 70, 3818–3833. [CrossRef]
17. Ni, P.; Liu, R.; Li, M.; Liu, Q. User Association and Hybrid Beamforming Designs for Cooperative mmWave MIMO Systems. IEEE
Trans. Signal Inf. Process. Netw. 2022, 8, 641–654. [CrossRef]
18. Xiu, Y.; Zhao, J.; Sun, W.; Renzo, M.D.; Gui, G.; Zhang, Z.; Wei, N. Reconfigurable Intelligent Surfaces Aided mmWave NOMA:
Joint Power Allocation, Phase Shifts, and Hybrid Beamforming Optimization. IEEE Trans. Wirel. Commun. 2021, 20, 8393–8409.
[CrossRef]
19. Ying, K.; Gao, Z.; Lyu, S.; Wu, Y.; Wang, H.; Alouini, M. GMD-Based Hybrid Beamforming for Large Reconfigurable Intelligent
Surface Assisted Millimeter-Wave Massive MIMO. IEEE Access 2020, 8, 19530–19539. [CrossRef]
20. Fu, S.; Zhang, W.; Jiang, Z. A network-level connected autonomous driving evaluation platform implementing C-V2X technology.
China Commun. 2021, 18, 77–88. [CrossRef]
21. Technical Specification Group Service and System Aspects: Enhancement of 3GPP Support for V2X Scenarios; (Release 16); Document TS
22.186, V16.2.0, V14.0.0; ETSI: Sophia Antipolis, France , 2019.
22. Raghavan, V.; Akhoondzadeh-Asl, L.; Podshivalov, V.; Hulten, J.; Tassoudji, M.; Koymen, O.; Sampath, A.; Li, J. Statistical Blockage
Modeling and Robustness of Beamforming in Millimeter-Wave Systems. IEEE Trans. Microw. Theory Tech. 2019, 67, 3010–3024.
[CrossRef]
23. Renzo, M.D.; Zappone, A.; Debbah, M.; Alouini, A.; Yuen, C.; Rosny, J.; Tretyakov, S. Smart Radio Environments Empowered
by Reconfigurable Intelligent Surfaces: How It Works, State of Research, and the Road Ahead. IEEE J. Sel. Areas Commun. 2020,
38, 2450–2525. [CrossRef]
24. Wymeersch, H.; He, J.; Denis, B.; Clemente, A.; Juntti, M. Radio Localization and Mapping with Reconfigurable Intelligent
Surfaces: Challenges, Opportunities, and Research Directions. IEEE Veh. Technol. Mag. 2020, 15, 52–61. [CrossRef]
25. Huang, C.; Zappone, A.; Alexandropoulos, G.C.; Debbah, M.; Yuen, C. Reconfigurable Intelligent Surfaces for Energy Efficiency
in Wireless Communication. IEEE Trans. Wirel. Commun. 2019, 18, 4157–4170. [CrossRef]
Sensors 2023, 23, 3713 16 of 16

26. Renzo, M.D.; Ntontin, K.; Song, J.; Danufane, F. Reconfigurable Intelligent Surfaces vs. Relaying: Differences, Similarities, and
Performance Comparison. IEEE Open J. Commun. Soc. 2020, 1, 798–807. [CrossRef]
27. Özdogan, Ö.; Björnson, E.; Larsson, E.G. Intelligent Reflecting Surfaces: Physics, Propagation, and Pathloss Modeling. IEEE Wirel.
Commun. Lett. 2020, 9, 581–585. [CrossRef]
28. Wei, L.; Huang, C.; Alexandropoulos, G.C.; Yuen, C.; Zhang, Z.; Debbah, M. Channel Estimation for RIS-Empowered Multi-User
MISO Wireless Communications. IEEE Trans. Commun. 2021, 69, 4144–4157. [CrossRef]
29. He, J.; Wymeersch, H.; Juntti, M. Channel Estimation for RIS-Aided mmWave MIMO Systems via Atomic Norm Minimization.
IEEE Trans. Wirel. Commun. 2021, 20, 5786–5797. [CrossRef]
30. Wang, P.; Fang, J.; Dai, L.; Li, H. Joint Transceiver and Large Intelligent Surface Design for Massive MIMO mmWave Systems.
IEEE Trans. Wirel. Commun. 2021, 20, 1052–1064. [CrossRef]
31. Ozcan, Y.U.; Ozdemir, O.; Kurt, G.K. Reconfigurable Intelligent Surfaces for the Connectivity of Autonomous Vehicles. IEEE
Trans. Veh. Technol. 2021, 70, 2508–2513. [CrossRef]
32. Ren, H.; Wang, K.; Pan, C. Intelligent Reflecting Surface-Aided URLLC in a Factory Automation Scenario. IEEE Trans. Commun.
2022, 70, 707–723. [CrossRef]
33. Yildirim, I.; Uyrus, A.; Basar, E. Modeling and Analysis of Reconfigurable Intelligent Surfaces for Indoor and Outdoor Applications
in Future Wireless Networks. IEEE Trans. Commun. 2021, 69, 1290–1301. [CrossRef]
34. Sun, Y.; Wang, C.; Huang, J.; Wang, J. A 3D Non-Stationary Channel Model for 6G Wireless Systems Employing Intelligent
Reflecting Surfaces with Practical Phase Shifts. IEEE Trans. Cogn. Commun. Netw. 2021, 7, 496–510. [CrossRef]
35. Pan, C.; Zhi, K.; Hong, S.; Wu, T.; Pan, Y.; Ren, H.; Renzo, M.D.; Swindlehurst, A.; Zhang, R.; Zhang, A. An Overview of Signal
Processing Techniques for RIS/IRS-Aided Wireless Systems. IEEE J. Sel. Top. Signal Process. 2022, 16, 883–917. [CrossRef]
36. Mei, W.; Zheng, B.; You, C.; Zhang, R. Intelligent Reflecting Surface-Aided Wireless Networks: From Single-Reflection to
Multireflection Design and Optimization. Proc. IEEE 2022, 110, 1380–1400. [CrossRef]
37. Li, Z.; Hua, M.; Wang, Q.; Song, Q. Weighted Sum-Rate Maximization for Multi-IRS Aided Cooperative Transmission. IEEE Wirel.
Commun. Lett. 2020, 9, 1620–1624. [CrossRef]
38. Yang, L.; Yang, J.; Xie, W.; Hasna, M.O.; Tsiftsis, T.; Renzo, M.D. Secrecy Performance Analysis of RIS-Aided Wireless Communica-
tion Systems. IEEE Trans. Veh. Technol. 2020, 69, 12296–12300. [CrossRef]
39. Shen, H.; Xu, W.; Gong, S.; Zhao, C.; Ng, D.W.K. Beamforming Optimization for IRS-Aided Communications with Transceiver
Hardware Impairments. IEEE Trans. Commun. 2021, 69, 1214–1227. [CrossRef]
40. Niu, H.; Lin, Z.; Zhu, Z.; Xiao, P.; Nguyen, H.; Lee, I.; Al-Dhahir, N. Joint Beamforming Design for Secure RIS-Assisted IoT
Networks. IEEE Internet Things J. 2023, 10, 1628–1641. [CrossRef]
41. Ayach, O.E.; Rajagopal, S.; Abu-Surra, S.; Pi, Z.; Heath, R.W. Spatially Sparse Precoding in Millimeter Wave MIMO Systems. IEEE
Trans. Wirel. Commun. 2014, 13, 1499–1513. [CrossRef]
42. Bai, Z.; Golub, G.H.; Ng, M.K. On inexact hermitian and skew-Hermitian splitting methods for non-Hermitian positive definite
linear systems. Linear Algebra Its Appl. 2008, 428, 413–440. [CrossRef]
43. Goldsmith, A. Wireless Communications; Cambridge University Press: Cambridge, UK, 2005.
44. Sun, Y.; An, K.; Luo, J.; Zhu, Y.; Zheng, G.; Chatzinotas, S. Outage Constrained Robust Beamforming Optimization for Multiuser
IRS-Assisted Anti-Jamming Communications with Incomplete Information. IEEE Internet Things J. 2022, 9, 13298–13314. [CrossRef]
45. Tang, L.; Zhu, Y.; Wu, H.; Fu, Q. Fast algorithms for sparse frequency waveform design with sidelobe constraint. Digit. Signal
Process. 2017, 69, 140–153. [CrossRef]
46. Sohrabi, F.; Yu, W. Hybrid Digital and Analog Beamforming Design for Large-Scale Antenna Arrays. IEEE J. Sel. Top. Signal
Process. 2016, 10, 501–513. [CrossRef]

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual
author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to
people or property resulting from any ideas, methods, instructions or products referred to in the content.