
directions, Θi, of the Mincoming signals. The Matrix Pencil
algorithm is similar to ESPRIT, however it is a non-statistical
method.
The advantages of the MUSIC algorithm compared to
classical DOA estimation methods include general sensor con-
figurations, ultra-slow sampling, and small array dimensions.
MUSIC can be applied to both narrowband and wideband sig-
nals without prior knowledge of the signals [11]. The MUSIC
algorithm, which is closely related to Prony’s method [12], has
been evaluated as having superior high-resolution performance
at the cost of computation and storage [13]. Using peaks of a
spatial pseudo-spectrum, MUSIC is able to identify the angle-
of-arrival of signals from −90 to +90 degrees. For this reason,
we chose to evaluate our DF disruption technique against an
adversarial node equipped with the MUSIC algorithm.
This paper is organized as follows. Section II defines the
RIS wave propagation model. How to program the surface to
maximize the Signal-to-Noise Ratio (SNR) at the adversarial
node is explained in Section III. In Section IV, we evaluate
the impact of correlated multipath on MUSIC direction-finding
performance, independent on the source of the multipath. In
Section V the RIS propagation model is applied and perfor-
mance is evaluated as a function of RIS and node geometry.
We end the paper with conclusions in Section VI.
II. RIS WAVE PROPAGATION MODEL
A RIS is a type of electromagnetic metasurface [1]. Typi-
cally, it is a sheet of inexpensive and adaptive thin composite
material, which can cover walls, buildings, ceilings, etc. It
consists of individually tuned passive reflectors, making the
phase response of the incoming wave tunable. The individ-
ually controlled unit cell of a RIS incorporates low power
electronic circuits with components such as positive-intrinsic-
negative (PIN) diodes and varactors, where the bias voltage
of the varactor can be tuned. The RIS can be controlled and
programmed using a simple microcontroller such as Raspberry
Pi. The RIS is reconfigurable and can be programmed to
control and modify the incident radio waves by elementary
electromagnetic functions such as reflection, refraction, ab-
sorption, focusing/beamforming, polarization, splitting, analog
processing and collimation [1]. A RIS can also be used for
joint transmit and passive beamforming design [14].
For the direct Line-of-Sight (LOS) path between the RIS
and the adversarial node, the attenuation and the delay of the
received baseband signal are [15],
A0=λcqGAdv
T x GT x
Adv
4πkpAdv −pT xk, τ0=kpAdv −pT xk
c0
,(1)
respectively, where λcis the wavelength of the carrier wave,
GAdv
T x is the antenna gain of the transmit antenna in the
direction of the adversarial node, GT x
Adv is the antenna gain of
the adversarial node in the direction of the transmitter, pAdv
is the 3-dimensional (3D) position of the adversarial node,
pT x is the 3D position of the transmitter, c0is the speed of
light through air, and k·k the l2-norm. For the RIS path, the
attenuation and delay are,
Ai,j =√µλ2
c
16π2qGAdv
i,j Gi,j
Adv Gi,j
T x GT x
i,j
kpAdv −pi,j kkpi,j −pT xk,(2)
τi,j =kpi,j −pT xk+kpAdv −pi,j k
c0
,(3)
respectively, where i, j are the element index, GAdv
i,j is the
antenna gain of the RIS element in the direction of the
adversarial node, Gi,j
Adv is the antenna gain of the adversarial
node in the direction of the RIS element, Gi,j
T x is the antenna
gain of the transmitter in the direction of the RIS element,
GT x
i,j is the antenna gain of the RIS element in the direction
of the transmitter, µ∈[0,1] is the fraction of the incident
energy that is scattered, and pi,j is the RIS element location
with x, y, z-coordinates,
x= 0
y=idy−0.5dy((Q+ 1) mod 2) ∀i= [1, Q]
z=jdz−0.5dz((P+ 1) mod 2) ∀j= [1, P ],
(4)
where dyis the element spacing in the y-direction, dzis
the element spacing in the z-direction, Q is the number of
element columns of the surface and P is the number of rows
of elements. The complex baseband signal received at the
adversarial node, y(t), can be expressed as,
y(t) = A0e−j2πfcτ0x(t−τ0)
+X
i,j
Ai,j e−j2πfcτi,j −jφi,j xt−τi,j −φi,j
2πfc
+n(t)
(5)
where τ0is the time delay from the transmitter to the adver-
sarial node, τi,j is the time delay from the transmitter to the
adversarial node via the RIS elements, φi,j
2πfcis a tunable delay
that enables the anomalous scattering of the incident radio
waves and n(t)is additive white Gaussian noise (AWGN).
III. RIS ARRAY PROCESSING
The phases φi,j are programmed by setting the bias voltage
of the varactor of each surface element. Maximizing the
received power of the baseband signal (5) requires that all
terms have the same phase. In order to maximize the SNR at
the ULA, φi,j is programmed as,
φi,j = 2πfc(τ0−τi,j ) mod 2π(6)
We assume that the location of the adversarial node is known.
IV. SIMULATION RESULTS
Without using the complex RIS model shown in Section
II, we first analyzed the behavior of the MUSIC algorithm
using orthogonal frequency division multiplexing (OFDM)
signals generated as they arrive on the adversarial node’s ULA.
We looked at the ability of the adversarial node to detect
signals in the MUSIC pseduo-spectrum when two correlated
or two uncorrelated signals arrive at the ULA from different