Empirical Validation of the Impedance-Based RIS Channel Model in an Indoor Scattering Environment

Reconfigurable intelligent surfaces (RISs) are widely considered as one of the ongoing revolutions in the network design due to their property to programmatically alter the propagation properties of the radio environment while keeping the overall manufacturing cost affordable. As such, RISs are a candidate technology for the next-generation wireless networks [1]. In this regard, one of the main research topics has been on exploring and further developing accurate channel modelling, which is essential to fully unlock the potentials of this technology.

Recently, the authors in [2] have proposed a mutual coupling and unit cell-aware electromagnetically-consistent channel model based on mutual impedances. Therein, the unit cells of the RIS are modelled as arbitrarily-spaced wire dipoles, whose loads can be tuned to alter the propagation conditions. Such a model has been exploited in [3] to derive the optimal RIS configuration that maximizes the received power in a single-user and single-antenna setup. Moreover, the authors in [4] have generalized this optimization framework for the case of multi-input multiple-output (MIMO) systems. Therein, the multipath components originated by scattering objects in the environment are modelled as an additive statistical component. On the other hand, the authors of [5] exploit the discrete dipole approximation (DDA) [6] to model the scattering objects present in the environment as loaded wire dipoles, thus effectively capturing its interaction with the transmitted signals. Under this setting, the authors of [7] developed a provably convergent and nearly-optimal RIS optimization algorithm based on Gram-Schmidt’s orthogonalization method.

In this paper, we go one-step beyond and provide an empirical validation of the model in [5] by showcasing how its parameters can be tuned to effectively recreate a real-life environment. To this purpose, we exploit a set of channel measurements [8] performed in an office environment wherein a reflective RIS operating at $28$ GHz [9] aids in the communication between a transmitter [10] and a single-antenna user equipment (UE), in the presence of several scattering objects. Additionally, we demonstrate how the considered RIS channel model effectively captures such physical propagation scenario with few parameters. Finally, numerical results demonstrate that performing the RIS optimization using such a tailored channel model can bring significant gains in terms of received power at the UE.

Notation. We use bold font lower and upper case for vectors and matrices, respectively. $(\cdot)^{\mathrm{H}}$ denotes the hermitian transpose operator, while $j=\sqrt{-1}$ is the imaginary number.

We consider the scenario depicted in Fig. 1, wherein a transmitter (TX) equipped with $M=400$ antenna elements is located in the origin of the reference system [10], an RIS equipped with $N=1600$ elements is placed (i.e., its center-point) at coordinates $\begin{bmatrix}-2&-3&0\end{bmatrix}$ m [9], a single-antenna UE is placed at coordinates $\begin{bmatrix}1&-3.6&0\end{bmatrix}$ m, and a number of scattering objects are present in the environment.Moreover, we employ the impedance-based RIS channel model in [5], such that the scattering objects in the environment are modelled as discrete arrays of half-wavelength loaded dipoles. Hence, the $1\times M$ channel vector is given by

where $Z_{\scriptscriptstyle\mathrm{RL}}\in\mbox{$\mathbb{C}$}$ ($\mathbf{Z}_{\scriptscriptstyle\mathrm{TG}}\in\mbox{$\mathbb{C}$}^{M\times M}$) accounts for the load and self-impedance at the UE (transmitter); $\mathbf{z}_{\scriptscriptstyle\mathrm{ROT}}\in\mbox{$\mathbb{C}$}^{1\times M}$ includes the mutual impedances between the transmitter and the UE, both directly and indirectly through the scattering objects; $\mathbf{z}_{\scriptscriptstyle\mathrm{ROS}}\in\mbox{$\mathbb{C}$}^{1\times N}$ ($\mathbf{Z}_{\scriptscriptstyle\mathrm{SOT}}\in\mbox{$\mathbb{C}$}^{N\times M}$) models the mutual impedance between the RIS and the UE (transmitter), including the effect due to the scattering objects; $\mathbf{Z}_{\scriptscriptstyle\mathrm{SS}}\in\mbox{$\mathbb{C}$}^{N\times N}$ and $\mathbf{Z}_{\scriptscriptstyle\mathrm{SOS}}\in\mbox{$\mathbb{C}$}^{N\times N}$ represent the self-impedance among the RIS elements and the self and mutual impedance among the scattering objects and the RIS, respectively; lastly, $\mathbf{Z}_{\scriptscriptstyle\mathrm{RIS}}\in\mbox{$\mathbb{C}$}^{N\times N}$ stands for the diagonal matrix of tunable RIS impedances modelled as $\mathbf{Z}_{\scriptscriptstyle\mathrm{RIS}}=\mathrm{diag}(R_{0}+jx_{\ell}),\quad\ell=1,\ldots,N$, where $R_{0}\geq 0$ is the constant resistance of each RIS load and $x_{\ell}\in\mbox{$\mathbb{R}$}$ is the tunable reactance at the $\ell$-th element. The precoding vector at the transmitter is denoted by $\mathbf{w}\in\mbox{$\mathbb{C}$}^{M\times 1}$. As depicted in Fig. 1, such vector is defined in order to produce a beam pointing at $-35^{\circ}$ in azimuth and $0^{\circ}$ in elevation. Specifically, its beampattern is given in [10]. We remark that the considered measurements collected in [8] involve arrays of patch antennas with a uniform gain of $6$ dBi, while the adopted RIS channel model in [5] assumes dipole antennas. Hence, to compensate for the reduced antenna gain, we design the precoder vector such that $\|\mathbf{w}\|^{2}\leq M\,G$, with $G=6$ dBi. Moreover, in order to reproduce the setting in [8], we fix the RIS configuration in $\mathbf{Z}_{\scriptscriptstyle\mathrm{RIS}}$ by feeding the channel vector in (1) in the absence of both the scattered paths and the direct link, to the SARIS algorithm in [5] (see Section III-B2 for details). The resulting beam at the RIS exhibits a main lobe pointing at $-10^{\circ}$ in azimuth and $0^{\circ}$ in elevation. As previously done for the precoder vector, we then compensate for the patch antenna gain at the RIS by multiplying the channel coefficient corresponding to the path reflected upon the RIS by $\sqrt{G}$. Lastly, the receive signal at the UE is thus given by

where $x\in\mbox{$\mathbb{C}$}$ is the (known) transmit signal, and $n$ is the additive white Gaussian noise coefficient.

In this section, we firstly describe the considered measurements collected in the real-life environment described in Fig. 1. Subsequently, focusing on the channel model in (1), we describe the procedure to fit the associated relevant parameters in order to match the underlying physical propagation environment.

In order to validate the effectiveness of the channel model in (1) and [5], we evaluate the (normalized) power received at the UE-side in dB, which is given by

in the various cases highlighted in Section III-B. The obtained results for the reflected and direct paths are shown in Table II, demonstrating an error of just $1.16\%$ and $0.35\%$ for the two cases, respectively. For the case of the path including the cylindrical reflector, we generate $10^{3}$ different combinations of the reactance $x_{t}$ for the load of each $t$-th dipole in the cylindrical array by drawing values from the uniform distribution $\mathcal{U}[-330,100]$. Furthermore, we vary the total number of dipoles $N_{x}\times N_{y}$ with fixed ratio $N_{y}=2N_{x}$, the inter-element distance $d$, and the real part of the dipole loads $y_{0}$. For each combination, we select the maximum value obtained over the random reactance.

In Fig. 3, we show the receive power at the UE versus the total number of dipoles for the case of $d=\lambda/4$, and for different values of $y_{0}$. We thus conclude that approximately $2450$ dipoles, corresponding to a $35\times 70$ configuration, with $y_{0}=100~{}\Omega$ is a precise enough fitting of the measured data, resulting in an error of $0.14\%$. Moreover, in Table III, we provide the simulated received power for the case of a $35\times 70$ cylindrical array with $y_{0}=100~{}\Omega$ for different values of the inter-element $d$, thus confirming that a dense array accurately represents the physical propagation properties of a real-life object.

While so far we have successfully emulated the measurements in [8] by configuring the RIS as described in Section II, in the following we demonstrate how the algorithm in [5] can lead to superior performance by designing an RIS configuration that explicitly interacts with the scattering environment and creates constructive interference at the UE. In this regard, we set all channel parameters according to the values that provide the best fitting of the physical propagation environment, as described above, and compare the RIS configuration used during the measurements described in Section III-A versus the SARIS algorithm in [5]. As shown in Table IV, the proposed RIS optimization algorithm achieves significant gains in the received power thanks to an accurate modelling and exploitation of the physical properties of the channel.

In this paper, we shed experimentally validate a recently-introduced electromagnetic-compliant channel model for RIS, which takes into account mutual coupling effects and considers the presence of scattering objects. To this end, we have validated its accuracy in representing the physical characteristics of real-world scenarios. Interestingly, we have fine-tuned the parameters of model using a dataset of channel measurements collected from an office environment. By recreating the channel propagation with the DDA, we have demonstrated that the proposed modeling can be leveraged to derive enhanced RIS configurations, leading to substantial improvements in the overall system performance.

This work was supported in part by the EU H2020 RISE-6G project under grant 101017011.