I Introduction
The wireless propagation environment is traditionally viewed as a detrimental medium of communication that distorts the transmitted signals in an uncontrollable manner. Wireless systems have improved consistently over the years through the development of increasingly sophisticated transmission and reception techniques that counteract the undesirable effects of propagation. An alternative to this conventional reactive approach is currently gaining interest, as several works propose to achieve a level of control over the propagation of electromagnetic waves by placing in the environment ‘smart’ passive devices that manipulate the impinging waves according to a desired functionality, such as reflection in an intended direction, full absorption, polarization, or filtering [1, 2, 3]. Furthermore, by endowing such devices with interconnection capability, they can be introduced in control loops that enable adaptation to the changing conditions of the wireless network, thus creating a smart radio environment. Examples of reconfigurable devices that can manipulate EM waves include tunable reflectarrays [4, 1], frequencyselective surfaces [5], or metasurfaces consisting of a thin layer of metamaterials [6]. The latter have particularly attractive properties and are expected to be a key enabler of smart wireless environments [2, 3].
The focus of this paper is on the concept of a large reflective surface (LRS) which, irrespective of the building technology, can be abstracted into an array of passive reflector elements that induce adjustable phase shifts on the reflected signals. An LRS placed somewhere between a source and a destination can greatly enhance the communication between the two when the phase shifts are appropriately configured such that the reflected waves combine coherently at the destination [1, 7, 8, 9]. Most of the studies assume the phase shifts induced by the composite propagation channel through the LRS are perfectly known and the reflector phases are set to the ideal values, with no errors. However, perfect phase estimation or highprecision configuration of the reflection phases is unfeasible. In [10, 11], the LRSbased system is optimized over lowresolution phase shifts.
In this paper, we study the communication through an LRS whose induced phase shifts deviate from the ideal values according to a generic probability distribution. This setup could, for example, represent imperfect phase estimation, quantized reflection phases, or both. We show that the transmission through a large number
of imperfect reflectors is equivalent to a pointtopoint communication over a Nakagami fading channel with the following characteristics: the diversity order grows with , the average SNR scales with relative to the average singlereflector SNR, and both parameters are attenuated by the phase uncertainty. Numerical results for a limited number of reflectors validate the theoretical prediction and show that accurate knowledge or representation of the ideal phase shifts is not necessary to closely approach the ideal performance.Ii System Model
The considered system comprises a singleantenna source , a singleantenna destination and an LRS with reflector elements , which assist the communication between and . Since we want to study the effect of the phase errors at the LRS, we restrict our focus on the transmission through the LRS only and assume there is no direct link between and . Assuming slow and flat fading, let and be the complex fading coefficients of the to and, respectively, to channels. The reflectors are sufficiently spaced apart such that all the coefficients are mutually independent. Furthermore, they have unit power and average magnitudes
(1) 
for all . In the following, we use
(2) 
Denote by the phase shift induced by . The signal received by the destination over the equivalent baseband channel is^{1}^{1}1
The signal is normalized by the standard deviation of the receiver noise.
(3) 
where is the transmitted symbol with mean zero and , models the normalized receiver noise and is the average SNR when only one reflector is present and accounts for path loss, shadowing and reflection loss.
Ideally, given the phases and , the phase is set so as to cancel the overall phase shift , which maximizes the SNR at the receiver [9]. Assuming the phase shifts induced by the channels are not estimated perfectly and/or the desired phases cannot be set precisely (e.g., when only a discrete set of phases is technologically possible), we model the deviation of from the ideal setting by the phase noise , which is randomly distributed on according to a certain circular distribution [12].
Now, we express the received signal as
(4) 
where
(5) 
All the variables in the r.h.s. of (4) and (5) are mutually independent. We assume ,
, are i.i.d. with common characteristic function expressed as a sequence of complex numbers
,(6) 
which are also called trigonometric (or circular) moments
[12]. Note that , for all ; forequality holds if, and only if, the distribution is degenerate (Dirac delta). A distribution that is broad on the circle (signifying high uncertainty) gives small magnitudes of the trigonometric moments (for the uniform distribution,
for ). It is reasonable to assume that the distribution of the error has mean direction zero () and its pdf is symmetric around zero. Due to the symmetry, the characteristic function takes real values, i.e., , for all .The receiver performs ideal coherent detection based on the sufficient statistic
(7) 
where .
Iii Equivalent Scalar Fading Channel
In this section, we find an equivalent representation of the transmission through an LRS as a pointtopoint fading channel. First, we provide the following result.
Lemma 1.
For large , the coefficient given by (5
) has a (non circularly symmetric) complex normal distribution. Its real and imaginary parts,
and , are independent, and and with(8)  
(9)  
(10) 
Proof:
See Appendix A. ∎
The ideal case where the phase shifts induced by the LRSaugmented channel are perfectly estimated and the reflector phases can be set precisely is instantiated as follows.
Corollary 1.
When there are no phase shift errors, i.e., , for all , the coefficient in (5) is real and distributed as .
Proof:
Since the phase errors are zero, we have . This implies that in Lemma 1 and the result follows. ∎
The ideal case of no phase errors and Rayleigh fading (i.e., ) is studied in [9], which Corollary 1 recovers.
Corollary 2.
When the phase shift errors are uniformly distributed on , which corresponds to a complete lack of knowledge about the phases of and , , the coefficient in (5
) has a circularlysymmetric complex normal distribution with mean zero and variance
, which means the equivalent channel resembles Rayleigh fading.Proof:
For the uniform circular distribution, and the result follows from Lemma 1. ∎
When , even though and are Gaussians with nonzero and, respectively, zero means, the magnitude does not have a Rice distribution because in general. We obtain the following result instead.
Theorem 1.
For large and , the magnitude of in (5) has a Nakagami distribution with pdf
(11) 
where the fading parameter is given by
(12) 
Proof:
See Appendix B. ∎
Thus, the transmission through the LRS is statistically equivalent to a direct transmission over a Nakagami fading channel.^{2}^{2}2The magnitude of the fading coefficient has a Nakagami distribution, although its phase is not uniformly distributed as is sometimes assumed for Nakagami fading. However, this is irrelevant, since we consider phase coherent detection. The equivalent channel has larger receive power relative to the transmission through one reflector. The average singlereflector SNR can be very low because of the attenuation of the constituent channels and reflection. Next, we make a few remarks.
Remark 1.
The full distributions of the coefficients , , , are irrelevant to the distribution of , as this depends only on the average magnitudes (1). The individual coefficients may even have different distributions, as long as they are mutually independent and have the respective means.
Furthermore, the phase errors manifest in (11) through the first two trigonometric moments (6) of their common distribution, which are subunitary reals for the considered zeromean symmetric distributions.
Remark 2.
The average SNR, , is attenuated by , such that the broader the phase distribution (larger errors), the lower .
The average SNR does not grow indefinitely with , because it cannot be larger than the transmit signal power divided by the receiver noise.^{3}^{3}3The upper bound would correspond to the contrived case when the emitted waves all arrive at the receiver without any absorption.
Remark 3.
The quantity is w.r.t. .
Iv Performance Analysis
We investigate how the performance of the system is impacted by the system size and phase errors. We consider the following types of errors:

phase estimation errors: is modelled as a zeromean von Mises variable whose concentration parameter captures the accuracy of the estimation. In this case, the characteristic function is , where is the modified Bessel function of the first kind and order [12].

quantization errors: when only a discrete set of phases can be configured, , the error is assumed to be uniformly distributed over . For this distribution, we obtain and .
We study these types of errors separately, although the case with both estimation and quantization errors can be easily treated because, assuming they are independent, the resulting and are given by the product of the respective trigonometric moments.
Iva Distribution of the SNR
The performance of the system is fundamentally determined by the distribution of the instantaneous SNR, which equals .
Corollary 3.
For large , the instantaneous SNR,
, is gamma distributed. Its pdf can be expressed as
(13) 
where is given by (12) and is the average SNR.
Proof:
As an example, the pdf (13) is displayed in Fig. 1 for the following scenario: the phase errors have a zeromean von Mises distribution with concentration ; the to channels exhibit Rician fading with unit power and a Rice factor of , which gives , while the fading over the to channels is Rayleigh with unit power, for which . We observe a very good agreement between the theoretical approximation and the Monte Carlo estimate even for a moderate value of . Also, when becomes large, the pdf concentrates around the mean and approaches a Gaussian pdf.
IvB Average error probability
We evaluate through simulations the average error probability, , of the LRS system for BPSK transmission, and the same fading models as in the previous example. In Fig. 2, the results for estimation and quantization errors are plotted, together with the theoretical performance [13] of the Nakagami fading channel given by Theorem 1. We observe a good agreement between the simulations and the theoretical prediction even for a limited . Other results not included here show that the match is closer for , whereas for the theoretical approximation is less accurate, although fairly good. As illustrated in Fig. (a)a, even a vague knowledge of the ideal phase shifts can be beneficial: e.g., when (for which the distribution is broad) the gap to ideal performance is of about dB, while for (for which the distribution is still not very concentrated) is within dB from the ideal case. It is also remarkable to observe in Fig. (b)b that phase quantization of just one bit is already at about dB from the ideal , while with two bits the error probability becomes close to the ideal.
The average error probability for the Nakagami channel at high average SNR is [14]
(14) 
Using , we express , such that the diversity and coding gains w.r.t. the average singlereflector SNR are and . These relationships may be useful to find the number of reflector elements required to achieve a certain or , given the distribution of the phase errors (e.g., for estimation errors). While can be easily established in the case of the diversity gain based on (12), the relationship is more complicated in the case of the coding gain.
V Conclusion
We showed that the composite channel through an LRS with phase errors is equivalent to a pointtopoint channel with Nakagami fading whose parameters are affected by the phase uncertainty through the first two trigonometric moments of the phase error distribution. Despite the phase shift errors, the average SNR still grows with (w.r.t. the singlereflector SNR) and the diversity order grows linearly with ; however, both are attenuated by the phase uncertainty. Also, the distributions of the fading coefficients of the two constituent channels manifest in the equivalent channel only though their average magnitudes. Numerical results for a limited number of reflectors showed good agreement with the theoretical prediction and that the performance is remarkably robust against the phase errors.
Appendix A Proof of Lemma 1
The complex variables , , are i.i.d. with common mean , variance and pseudovariance (see [15] for the secondorder characterization of complex variables). According to the assumptions on the distribution of , the trigonometric moments and, consequently,
. By invoking the central limit theorem, we approximate the distribution of
in (5) for large by a complex normal distribution, . Let and . We have [15](15) 
Being jointly Gaussian with zero covariance, it follows that and are independent. Furthermore, and with ,
(16)  
(17) 
Appendix B Proof of Theorem 1
We first consider the squared magnitude . The normalized variable
has a noncentral chisquared distribution, while
is gamma distributed with shape and scale . Thus, is the sum of a scaled noncentral chisquared variable and a gamma variable. We characterize the distribution of through its cumulant generating function, . Since and are independent (Lemma 1),(18) 
We now approximate for large the first term in the r.h.s. of (18). We develop its Maclaurin series expansion as
where
Given that , we have , such that (18) is well approximated by^{4}^{4}4Alternatively, the Maclaurin expansion of (18
) can be truncated to secondorder, which corresponds to a Gaussian distribution and incurs an error of
as well. However, for moderate values of , the Gaussian approximation is less accurate than the one we developed.(19) 
Expression (19) corresponds to the cumulant generating function of the sum of three independent gamma variables with (shape, scale) parameters , and respectively . While the resulting distribution can be characterized based on [16], we simplify by making a further approximation,
(20) 
which has error. Thus, for large , has a gamma distribution with shape and scale . It follows through variable transformation that has a Nakagami distribution with fading parameter (shape)
and spread .
References
 [1] X. Tan, Z. Sun, J. M. Jornet, and D. Pados, “Increasing indoor spectrum sharing capacity using smart reflectarray,” in 2016 IEEE International Conference on Communications (ICC), May 2016, pp. 1–6.
 [2] C. Liaskos, S. Nie, A. Tsioliaridou, A. Pitsillides, S. Ioannidis, and I. Akyildiz, “A new wireless communication paradigm through softwarecontrolled metasurfaces,” IEEE Communications Magazine, vol. 56, no. 9, pp. 162–169, Sep. 2018.
 [3] M. D. Renzo, M. Debbah, D.T. PhanHuy, A. Zappone, M.S. Alouini, C. Yuen, V. Sciancalepore, G. C. Alexandropoulos, J. Hoydis, H. Gacanin, J. d. Rosny, A. Bounceur, G. Lerosey, and M. Fink, “Smart radio environments empowered by reconfigurable AI metasurfaces: an idea whose time has come,” EURASIP Journal on Wireless Communications and Networking, vol. 2019, no. 1, p. 129, May 2019.
 [4] S. V. Hum, M. Okoniewski, and R. J. Davies, “Modeling and design of electronically tunable reflectarrays,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 8, pp. 2200–2210, Aug 2007.
 [5] L. Subrt and P. Pechac, “Controlling propagation environments using intelligent walls,” in 2012 6th European Conference on Antennas and Propagation (EUCAP), March 2012, pp. 1–5.
 [6] T. Jun Cui, M. Q. Qi, X. Wan, J. Zhao, and Q. Cheng, “Coding metamaterials, digital metamaterials and programmable metamaterials,” Light: Science & Applications, vol. 3, p. e218, 10 2014.
 [7] C. Huang, A. Zappone, M. Debbah, and C. Yuen, “Achievable rate maximization by passive intelligent mirrors,” in 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), April 2018, pp. 3714–3718.
 [8] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network: Joint active and passive beamforming design,” in 2018 IEEE Global Communications Conference (GLOBECOM), Dec 2018, pp. 1–6.
 [9] E. Basar, “Transmission Through Large Intelligent Surfaces: A New Frontier in Wireless Communications,” arXiv eprints, p. arXiv:1902.08463, Feb 2019.
 [10] C. Huang, G. C. Alexandropoulos, A. Zappone, M. Debbah, and C. Yuen, “Energy efficient multiuser MISO communication using low resolution large intelligent surfaces,” in 2018 IEEE Globecom Workshops (GC Wkshps), Dec 2018, pp. 1–6.
 [11] Q. Wu and R. Zhang, “Beamforming optimization for intelligent reflecting surface with discrete phase shifts,” in ICASSP 2019  2019 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), May 2019, pp. 7830–7833.
 [12] K. V. Mardia and P. E. Jupp, Directional Statistics. John Wiley & Sons, 2000.
 [13] M. K. Simon and M.S. Alouini, Digital communication over fading channels. John Wiley & Sons, 2000.
 [14] Zhengdao Wang and G. B. Giannakis, “What determines average and outage performance in fading channels?” in Global Telecommunications Conference, 2002. GLOBECOM ’02. IEEE, vol. 2, Nov 2002, pp. 1192–1196 vol.2.

[15]
B. Picinbono, “Secondorder complex random vectors and normal distributions,”
IEEE Transactions on Signal Processing, vol. 44, no. 10, pp. 2637–2640, Oct 1996. 
[16]
P. G. Moschopoulos, “The distribution of the sum of independent gamma random variables,”
Annals of the Institute of Statistical Mathematics, vol. 37, no. 3, pp. 541–544, Dec 1985.
Comments
There are no comments yet.