This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Anchor-Assisted Intelligent Reflecting Surface Channel Estimation for Multiuser Communications

Xinrong Guan12, Qingqing Wu3, Rui Zhang4 1Communications Engineering College, Army Engineering University of PLA, Nanjing, 210007, China
2 Postdoctoral Station, Shenzhen Electric Appliance Company, Shenzhen, 518022, China
3 Faculty of Science and Technology, University of Macau, Macau, 999078, China
4 Department of Electrical and Computer Engineering, National University of Singapore, 117583, Singapore
Email: [email protected], [email protected], [email protected]
Abstract

Due to the passive nature of Intelligent Reflecting Surface (IRS), channel estimation is a fundamental challenge in IRS-aided wireless networks. Particularly, as the number of IRS reflecting elements and/or that of IRS-served users increase, the channel training overhead becomes excessively high. To tackle this challenge, we propose in this paper a new anchor-assisted two-phase channel estimation scheme, where two anchor nodes, namely A1 and A2, are deployed near the IRS for helping the base station (BS) to acquire the cascaded BS-IRS-user channels. Specifically, in the first phase, the partial channel state information (CSI), i.e., the element-wise channel gain square, of the BS-IRS link is obtained by estimating the BS-IRS-A1/A2 channels and the A1-IRS-A2 channel, separately. Then, in the second phase, by leveraging such partial knowledge of the BS-IRS channel that is common to all users, the individual cascaded BS-IRS-user channels are efficiently estimated. Simulation results demonstrate that the proposed anchor-assisted channel estimation scheme is able to achieve comparable mean-squared error (MSE) performance as compared to the conventional scheme, but with significantly reduced channel training time.

I Introduction

Recently, intelligent reflecting surface (IRS) has emerged as a promising technology to achieve high spectral and energy efficiency for future wireless networks [1], [2]. Specifically, IRS is a uniform planar array composed of a large number of low-cost, passive, and tunable reflecting elements. By adaptively varying the reflection coefficient of each element based on the user dynamic channels, IRS can achieve high beamforming and interference suppression gains cost-effectively [3]. As such, IRS has been studied recently in various wireless systems, including non-orthogonal multiple access (NOMA)[4], [5], simultaneous wireless information and power transfer (SWIPT) [6], [7], secrecy communications [8], [9], and so on.

To reap the performance gain of IRS, accurate channel state information (CSI) is required. However, the passive nature of IRS makes channel estimation fundamentally challenging in IRS-aided wireless networks. This is because without radio frequency (RF) chains, IRS can neither transmit nor receive pilot signals, thus the base station (BS)-IRS and IRS-user channels cannot be estimated separately. One alternative approach is to estimate the cascaded BS-IRS-user channel via element-wise on-off operation at each reflecting element [10], [11] or time-varying reflection patterns [12], [13]. However, by applying the above methods to estimate the cascaded channels of multiple users consecutively, the required pilot overhead is the product of the number of IRS reflecting elements and that of users, which is prohibitively large for the case of large IRS serving a high density of users nearby (e.g., in a hot spot scenario).

Refer to caption
Figure 1: IRS-assisted multiuser MISO communication system.

To tackle this problem, we propose in this paper a new anchor-assisted two-phase channel estimation scheme for IRS-aided multiuser communications, which can significantly reduce the channel training overhead by decoupling the estimation of cascaded channels and exploiting multiple antennas at the BS. As shown in Fig. 1, two anchor nodes, namely A1 and A2, are deployed near the IRS to assist its channel estimation. In the first phase, the cascaded BS-IRS-A1/A2 and A1-IRS-A2 channels are estimated separately, based on which the partial CSI of the BS-IRS link, in terms of the square of each IRS element’s channel with the BS (thus, with a +/- sign uncertainty), is obtained. In the second phase, with such partial CSI of the BS-IRS channel that is common to all users, the individual cascaded BS-IRS-user channels are efficiently estimated. Particularly, when the number of antennas at the BS (MM) is no smaller than that (NN) of the IRS reflecting elements, i.e. MNM\!\geq\!N, we show that the proposed scheme can estimate all users’ cascaded channels in the second phase using only one pilot symbol. Besides, even for the case of M<N\!M<N\!, the minimum training overhead (in terms of number of pilot symbols) is shown to be NKM\left\lceil\frac{NK}{M}\right\rceil, where KK is the number of IRS-served users and \left\lceil\cdot\right\rceil denotes the ceiling operation, which is significantly lower than NKNK for the conventional schemes [10, 11, 12, 13]. Moreover, in the case that the IRS-anchor channel is line-of-sight (LoS) by properly deploying the anchor, it is shown that the proposed scheme can be further simplified such that deploying only one anchor is sufficient.

Refer to caption
Figure 2: Chanel estimation and data transmission protocol.

II System Model

As shown in Fig. 1, we consider an IRS-assisted multiple-input-single-output (MISO) communication system, which consists of a BS, an IRS and KK users. The number of antennas at the BS and that of reflecting elements at the IRS are denoted by MM and NN, respectively. The channels from the BS to IRS and User kk are denoted by 𝐇bsM×N{\bf{H}}_{bs}\in{\mathbb{C}^{M\times N}} and 𝐡bukM×1{\bf{h}}_{bu_{k}}\in{\mathbb{C}^{M\times 1}}, respectively, while that from the IRS to User kk is denoted by 𝐡sukN×1{\bf{h}}_{su_{k}}\in{\mathbb{C}^{N\times 1}}. We assume that two single-antenna anchor nodes111In practice, anchors can be idle user terminals and/or dedicated nodes such as adjacent IRS controllers. A1 and A2 are deployed near the IRS to assist in the channel estimation. The channels from the BS to A1 and A2 are denoted by 𝐡ba1M×1{\bf{h}}_{ba_{1}}\!\in\!{\mathbb{C}^{M\times 1}} and 𝐡ba2M×1{\bf{h}}_{ba_{2}}\!\in\!{\mathbb{C}^{M\times 1}}, respectively, those from the IRS to A1 and A2 are denoted by 𝐡sa1N×1{\bf{h}}_{sa_{1}}\!\in\!{\mathbb{C}^{N\times 1}} and 𝐡sa2N×1{\bf{h}}_{sa_{2}}\!\in\!{\mathbb{C}^{N\times 1}}, respectively, and that from A1 to A2 is denoted by ha1a2{{h}}_{a_{1}a_{2}}. As a result, the cascaded BS-IRS-A1/A2/User kk channels are denoted by 𝐇bsa1=𝐇bsdiag(𝐡sa1){\bf H}_{bsa_{1}}\!=\!{\bf H}_{bs}\text{diag}({\bf h}_{sa_{1}}), 𝐇bsa2=𝐇bsdiag(𝐡sa2){\bf H}_{bsa_{2}}\!=\!{\bf H}_{bs}\text{diag}({\bf h}_{sa_{2}}) and 𝐇bsuk=𝐇bsdiag(𝐡suk){\bf H}_{bsu_{k}}={\bf H}_{bs}\text{diag}({\bf h}_{su_{k}}), respectively, while the cascaded A1-IRS-A2 channel is denoted by 𝐡a1sa2=𝐡sa2Tdiag(𝐡sa1){\bf h}_{a_{1}sa_{2}}={\bf h}_{sa_{2}}^{T}\text{diag}({\bf h}_{sa_{1}}). Moreover, the phase-shift matrix of the IRS at time slot ii is denoted by 𝚽i=diag(v1,i,v2,i,.,vN,i){\mathbf{\Phi}}_{i}=\text{diag}\left(v_{1,i},v_{2,i},....,v_{N,i}\right), where vn,iv_{n,i} is the reflection coefficient of the nn-th IRS element at time slot ii, n=1,,Nn=1,...,N. Since IRS is a passive reflecting device, we assume that the channel reciprocity holds for each link between IRS and any other node. The quasi-static flat-fading channel model is assumed for all the links involved.

III Anchor-assisted Channel Estimation

The proposed channel estimation and data transmission protocol is shown in Fig. 2, where TbsT_{bs} and TsuT_{su} are the length of channel coherence block of 𝐇bs{\bf{H}}_{bs} and 𝐡suk{\bf{h}}_{su_{k}}, k\forall k, respectively, and L=KL=K for MNM\geq N while L=KNML=\left\lceil{\frac{KN}{{M}}}\right\rceil for M<NM<N. In practice, TbsT_{bs} is usually much larger than TsuT_{su} due to the fixed locations of the BS and IRS once deployed, while users can move randomly near the IRS. The channel estimation consists of one off-line phase (Phase I) and multiple on-line phases (each termed Phase II). In Phase I, the cascaded BS-IRS-A1, BS-IRS-A2, A1-IRS-A2 channels are estimated separately, which requires 2(N+1)2(N+1) pilot symbols in total. In each Phase II, the BS estimates the cascaded BS-IRS-user channels with K+LK+L pilot symbols. For the estimation of 𝐡bsuk{\bf{h}}_{bsu_{k}} in Phase II, the efficiency is largely improved, especially when MM increases up to MNM\geq N. The detailed channel estimation scheme is elaborated as follows.

III-A Phase I: Off-line Estimation of 𝐇bs𝐇bs{\bf{H}}_{bs}\odot{\bf{H}}_{bs}

III-A1 Step 1

A1 transmits pilot symbol a1,ia_{1,i} with power pp at time slot ii, then the received signals at the BS and A2 are respectively written as

𝐲b(i)=p(𝐡ba1+𝐇bs𝚽i𝐡sa1)a1,i+𝐳b(i),\small{\bf{y}}_{b}^{(i)}=\sqrt{p}\left({\bf{h}}_{ba_{1}}+{\bf H}_{bs}{\bf{\Phi}}_{i}{\bf{h}}_{sa_{1}}\right)a_{1,i}+{\bf{z}}_{b}^{(i)},\vspace{-1mm} (1)
y2(i)=p(ha1a2+𝐡sa2T𝚽i𝐡sa1)a1,i+z2(i).\small{{{y}}_{2}^{(i)}}=\sqrt{p}\left({{h}}_{a_{1}a_{2}}+{\bf{h}}_{sa_{2}}^{T}{\bf{\Phi}}_{i}{\bf{h}}_{sa_{1}}\right)a_{1,i}+{z_{2}^{(i)}}.\vspace{-1mm} (2)

Let 𝐯i=[v1,i,v2,i,.,vN,i]T{\bf{v}}_{i}=[v_{1,i},v_{2,i},....,v_{N,i}]^{T}, 𝐯~i=[1,𝐯iT]T{\bf{\tilde{v}}}_{i}=[1,{\bf{v}}_{i}^{T}]^{T}, 𝐇~ba1=[𝐡ba1,𝐇bsa1]{\bf{\tilde{H}}}_{ba_{1}}=[{\bf{h}}_{ba_{1}},{\bf{H}}_{bsa_{1}}], and 𝐡~a1a2=[ha1a2,𝐡a1sa2]{\bf{\tilde{h}}}_{a_{1}a_{2}}=[h_{a_{1}a_{2}},{\bf{h}}_{a_{1}sa_{2}}]. Then, (1) and (2) are rewritten as

𝐲b(i)=p𝐇~ba1𝐯~ia1,i+𝐳b(i),\small{\bf{y}}_{b}^{(i)}=\sqrt{p}{\bf{\tilde{H}}}_{ba_{1}}{{\bf{\tilde{v}}}_{i}}{a_{1,i}}+{\bf{z}}_{b}^{(i)},\vspace{-1mm} (3)
y2(i)=p𝐡~a1a2𝐯~ia1,i+z2(i).\small{{{y}}_{2}^{(i)}}=\sqrt{p}{\bf{\tilde{h}}}_{a_{1}a_{2}}{{\bf{\tilde{v}}}_{i}}{a_{1,i}}+{z_{2}^{(i)}}.\vspace{-1mm} (4)

There are in total (N+1)M(N+1)M unknowns in 𝐇~ba1{\bf{\tilde{H}}}_{ba_{1}} and the BS has MM observations at each time slot, thus A1 has to transmit N+1N+1 pilot symbols at least for the BS to estimate 𝐇~ba1{\bf{\tilde{H}}}_{ba_{1}}. Similarly, it can be shown that at least N+1N+1 pilot symbols are needed for A2 to estimate 𝐡~a1a2{\bf{\tilde{h}}}_{a_{1}a_{2}}. By setting a1,i=1a_{1,i}=1, i=1,,N+1i=1,...,N+1 and denoting 𝐕~=[𝐯~1,,𝐯~N+1]{{\bf{\tilde{V}}}}=[{{\bf{\tilde{v}}}_{1}},...,{{\bf{\tilde{v}}}_{N+1}}], the overall received signals at the BS and A2 during the N+1N+1 time slots are given by

𝐘b=[𝐲b(1),,𝐲b(N+1)]=p𝐇~ba1𝐕~+𝐙b,\small{{\bf{Y}}_{b}}=[{\bf{y}}_{b}^{(1)},...,{\bf{y}}_{b}^{(N+1)}]=\sqrt{p}{\bf{\tilde{H}}}_{ba_{1}}{{\bf{\tilde{V}}}}+{\bf{Z}}_{b}, (5)
𝐲2=[y2(1),,y2(N+1)]=p𝐡~a1a2𝐕~+𝐳2,\small{{\bf{y}}_{2}}=[{{{y}}_{2}^{(1)}},...,{{{y}}_{2}^{(N+1)}}]=\sqrt{p}{\bf{\tilde{h}}}_{a_{1}a_{2}}{{\bf{\tilde{V}}}}+{{\bf{z}}_{2}}, (6)

where 𝐙b=[𝐳b(1),,𝐳b(N+1)]{\bf{Z}}_{b}=[{\bf{z}}_{b}^{(1)},...,{\bf{z}}_{b}^{(N+1)}] and 𝐳2=[z2(1),,z2(N+1)]{{\bf{z}}_{2}}=[{z_{2}^{(1)}},...,{z_{2}^{(N+1)}}], respectively. By properly constructing 𝐕~{{\bf{\tilde{V}}}} such that rank(𝐕~)=N+1\text{rank}({{\bf{\tilde{V}}}})=N+1, 𝐇~ba1{\bf{\tilde{H}}}_{ba_{1}} and 𝐡~a1a2{\bf{\tilde{h}}}_{a_{1}a_{2}} can be respectively estimated as

𝐇^ba1=[𝐡^ba1,𝐇^bsa1]=1p𝐘b𝐕~1,\small{\bf{\hat{H}}}_{ba_{1}}=[{\bf{\hat{h}}}_{ba_{1}},{\bf{\hat{H}}}_{bsa_{1}}]=\frac{1}{\sqrt{p}}{{\bf{Y}}_{b}}{{\bf{\tilde{V}}}^{-1}},\vspace{-1mm} (7)
𝐡^a1a2=[h^a1a2,𝐡^a1sa2]=1p𝐲2𝐕~1.\small{\bf{\hat{h}}}_{a_{1}a_{2}}=[{\hat{h}}_{a_{1}a_{2}},{\bf{{\hat{h}}}}_{a_{1}sa_{2}}]=\frac{1}{\sqrt{p}}{{\bf{y}}_{2}}{{\bf{\tilde{V}}}^{-1}}.\vspace{-2mm} (8)

Practically, 𝐕~{{\bf{\tilde{V}}}} can be constructed based on the discrete Fourier transform (DFT) matrix [13], i.e.,

𝐕~=[1111ej2πN+1ej2πN+1N1ej2πN+12ej2πN+12N......1ej2πN+1Nej2πN+1N2].\small{{\bf{\tilde{V}}}}=\left[\begin{aligned} &1~{}~{}~{}&~{}1~{}~{}~{}~{}~{}~{}~{}&...~{}~{}&~{}1~{}~{}~{}~{}~{}\\ &1~{}~{}~{}&e^{-j\frac{2\pi}{N+1}}~{}~{}~{}&...~{}~{}&e^{-j\frac{2\pi}{N+1}N}\\ &1~{}~{}~{}&e^{-j\frac{2\pi}{N+1}2}~{}~{}~{}&...~{}~{}&e^{-j\frac{2\pi}{N+1}2N}\\ &.&~{}.~{}~{}~{}~{}~{}~{}~{}&...~{}~{}&~{}.~{}~{}~{}~{}~{}\\ &.&~{}.~{}~{}~{}~{}~{}~{}~{}&...~{}~{}&~{}.~{}~{}~{}~{}~{}\\ &1~{}~{}~{}&e^{-j\frac{2\pi}{N+1}N}~{}~{}~{}&...~{}~{}&e^{-j\frac{2\pi}{N+1}N^{2}}\end{aligned}\right].\vspace{-1mm} (9)

In this case, 𝐕~1{{\bf{\tilde{V}}}^{-1}} can be efficiently computed as by 𝐕~1=1N+1𝐕~H{{\bf{\tilde{V}}}^{-1}}=\frac{1}{N+1}{{\bf{\tilde{V}}}^{H}}.

[𝐲(1)𝐲(2)..𝐲(N)]=p[𝐇bs𝚽1x1,1𝐇bs𝚽1x2,1𝐇bs𝚽1x3,1𝐇bs𝚽1xM,1𝐇bs𝚽2x1,2𝐇bs𝚽2x2,2𝐇bs𝚽2x3,2𝐇bs𝚽2xM,2.......𝐇bs𝚽Nx1,N𝐇bs𝚽Nx2,N𝐇bs𝚽Nx3,N𝐇bs𝚽NxM,N]𝐁MN×MN[𝐡su1𝐡su2..𝐡suM]\small\left[{\begin{array}[]{*{20}{c}}{{{\bf{y}}^{(1)}}}\\ {{{\bf{y}}^{(2)}}}\\ \begin{array}[]{l}.\\ .\end{array}\\ {{{\bf{y}}^{\left(N\right)}}}\end{array}}\right]=\sqrt{p}\underbrace{\left[{\begin{array}[]{*{20}{c}}{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{1,1}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{2,1}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{3,1}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{M,1}}}\\ {{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{1,2}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{2,2}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{3,2}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{M,2}}}\\ .&.&.&{...}&.\\ .&.&{}\hfil&.&{}\hfil\\ {{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N}}{x_{1,N}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N}}{x_{2,N}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N}}{x_{3,N}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N}}{x_{M,N}}}\end{array}}\right]}_{{\bf{B}}\in{\mathbb{C}^{MN\times MN}}}\left[{\begin{array}[]{*{20}{c}}{{{\bf{h}}_{s{u_{1}}}}}\\ {{{\bf{h}}_{s{u_{2}}}}}\\ \begin{array}[]{l}.\\ .\end{array}\\ {{{\bf{h}}_{s{u_{M}}}}}\end{array}}\right]\vspace{-3mm} (21)

III-A2 Step 2

With 𝐇~ba2=[𝐡ba2,𝐇bsa2]{\bf{\tilde{H}}}_{ba_{2}}=[{\bf{h}}_{ba_{2}},{\bf{H}}_{bsa_{2}}], A2 transmits at least N+1N+1 pilot symbols so that the BS can estimate 𝐇~ba2{\bf{\tilde{H}}}_{ba_{2}} as 𝐇^ba2=[𝐡^ba2,𝐇^bsa2]{\bf{\hat{H}}}_{ba_{2}}=[{\bf{\hat{h}}}_{ba_{2}},{\bf{\hat{H}}}_{bsa_{2}}], similarly as in Step 1 and thus the details are omitted. With 𝐡^a1sa2{\bf\hat{h}}_{a_{1}sa_{2}} fed back from A2, BS obtains the estimated BS-IRS-A1, BS-IRS-A2 and A1-IRS-A2 channels, which are given by

𝐇^bsa1\displaystyle{\bf{\hat{H}}}_{bsa_{1}} =𝐇^bsdiag(𝐡^sa1),\displaystyle={\bf{\hat{H}}}_{bs}{\rm{diag}}({\bf{\hat{h}}}_{sa_{1}}), (10)
𝐇^bsa2\displaystyle{\bf{\hat{H}}}_{bsa_{2}} =𝐇^bsdiag(𝐡^sa2),\displaystyle={\bf{\hat{H}}}_{bs}{\rm{diag}}({\bf{\hat{h}}}_{sa_{2}}), (11)
𝐡^a1sa2\displaystyle{\bf{\hat{h}}}_{a_{1}sa_{2}} =𝐡^sa2Tdiag(𝐡^sa1).\displaystyle={\bf{\hat{h}}}_{sa_{2}}^{T}{\rm{diag}}({\bf{\hat{h}}}_{sa_{1}}). (12)

Based on (10)-(12), the BS computes

𝐇^bs𝐇^bs=𝐇^bsa1𝐇^bsa2(diag(𝐡^a1sa2))1,{\bf{\hat{H}}}_{bs}\odot{\bf{\hat{H}}}_{bs}={\bf{\hat{H}}}_{bsa_{1}}\odot{\bf{\hat{H}}}_{bsa_{2}}({\rm{diag}}({\bf{\hat{h}}}_{a_{1}sa_{2}}))^{-1}, (13)

where \odot denotes the Hadamard product. By defining 𝐆=𝐇^bsa1𝐇^bsa2(diag(𝐡^a1sa2)1){\bf{G}}={\bf{\hat{H}}}_{bsa_{1}}\odot{\bf{\hat{H}}}_{bsa_{2}}({\rm{diag}}({\bf{\hat{h}}}_{a_{1}sa_{2}})^{-1}), 𝐇^bs𝐇^bs{\bf{\hat{H}}}_{bs}\odot{\bf{\hat{H}}}_{bs} is rewritten as

𝐇^bs𝐇^bs=[𝐆11𝐆12𝐆1N𝐆21𝐆22𝐆2N......𝐆M1𝐆M2𝐆MN],\small{\bf{\hat{H}}}_{bs}\odot{\bf{\hat{H}}}_{bs}=\left[\begin{aligned} &{\bf{G}}_{11}&{\bf{G}}_{12}~{}~{}~{}&...&{\bf{G}}_{1N}~{}\\ &{\bf{G}}_{21}&{\bf{G}}_{22}~{}~{}~{}&...&{\bf{G}}_{2N}~{}\\ &~{}~{}.&.~{}~{}~{}~{}~{}&...&.~{}~{}~{}\\ &~{}~{}.&.~{}~{}~{}~{}~{}&...&.~{}~{}~{}\\ &{\bf{G}}_{M1}&{\bf{G}}_{M2}~{}~{}&...&{\bf{G}}_{MN}\end{aligned}\right], (14)

where 𝐆mn=𝐇^bs(m,n)2{\bf G}_{mn}\!=\!{\bf\hat{H}}_{bs}(m,n)^{2}. Letting gmn=|𝐆mn|ej𝐆mn2g_{mn}=\sqrt{|{\bf{G}}_{mn}|}e^{j\frac{\angle{\bf{G}}_{mn}}{2}}, then each element in 𝐇^bs{\bf{\hat{H}}}_{bs} can be obtained as

𝐇^bs(m,n)=±gmn,m,n,{\bf{\hat{H}}}_{bs}(m,n)=\pm g_{mn},\forall m,n,\vspace{-1mm} (15)

i.e., we recover each 𝐇^bs(m,n){\bf{\hat{H}}}_{bs}(m,n) but with a +/- sign uncertainty. However, we will show later that such partial CSI estimated is sufficient for resolving the cascaded channels 𝐇bsuk{\bf H}_{bsu_{k}}’s of all users uniquely.

III-B Phase II: On-line Estimation of 𝐡buk{\bf{h}}_{bu_{k}} and 𝐇bsuk{\bf{H}}_{bsu_{k}}

III-B1 Step 1

KK users sequentially transmit one pilot symbol while the BS estimates 𝐡buk{\bf{h}}_{bu_{k}}, k\forall k, respectively, where the details are omitted for brevity.

III-B2 Step 2

Users transmit max(K,KNM)\text{max}\left(K,\left\lceil{\frac{KN}{M}}\right\rceil\right) pilot symbols while the BS estimates 𝐇bsuk{\bf{H}}_{bsu_{k}}, k\forall k, respectively. Denoting the pilot symbol transmitted from User kk at time slot ii by xk,ix_{k,i}, the received signal at the BS is given by

𝐲b(i)=pk=1K(𝐡buk+𝐇bs𝚽i𝐡suk)xk,i+𝐳b(i).\displaystyle{\bf{y}}_{b}^{(i)}=\sqrt{p}\sum_{k=1}^{K}\left({\bf{h}}_{bu_{k}}+{\bf H}_{bs}{\bf{\Phi}}_{i}{\bf{h}}_{su_{k}}\right)x_{k,i}+{\bf{z}}_{b}^{(i)}.\vspace{-1mm} (16)

Let 𝐇bs𝚽i𝐡suk=𝐇bsuk𝐯i{\bf H}_{bs}{\bf{\Phi}}_{i}{\bf{h}}_{su_{k}}={\bf{H}}_{bsu_{k}}{\bf{v}}_{i}, where 𝐇bsuk=𝐇bsdiag(𝐡suk){\bf H}_{bsu_{k}}={\bf H}_{bs}\text{diag}({\bf h}_{su_{k}}) and 𝚽i=diag(𝐯i){\bf{\Phi}}_{i}=\text{diag}({\bf{v}}_{i}). Then, by removing the signal from the direct channel, 𝐲b(i){\bf{y}}_{b}^{(i)} can be re-expressed as

𝐲¯b(i)=pk=1K𝐇^bsuk𝐯ixk,i+𝐳¯b(i),{\bf{\bar{y}}}_{b}^{(i)}=\sqrt{p}\sum_{k=1}^{K}{\bf{\hat{H}}}_{bsu_{k}}{\bf{v}}_{i}x_{k,i}+{\bf{\bar{z}}}_{b}^{(i)},\vspace{-3mm} (17)

where

𝐳¯b(i)=𝐳b(i)+pk=1K(𝐡buk𝐡^buk+(𝐇bsuk𝐇^bsuk)𝐯i)xk,i{\bf{\bar{z}}}_{b}^{(i)}\!=\!{\bf{z}}_{b}^{(i)}+\sqrt{p}\sum_{k=1}^{K}\left({\bf{h}}_{bu_{k}}\!-\!{\bf{\hat{h}}}_{bu_{k}}\!+\!({\bf H}_{bsu_{k}}\!-\!{\bf{\hat{H}}}_{bsu_{k}}){\bf{v}}_{i}\right)x_{k,i}\vspace{-2mm}

is the effective noise, including the channel estimation error in 𝐡^buk{\bf{\hat{h}}}_{bu_{k}} and 𝐇^bsuk{\bf{\hat{H}}}_{bsu_{k}}. In conventional schemes [10, 11, 12, 13], NN pilot symbols are required to estimate each 𝐇bsuk{\bf{H}}_{bsu_{k}}. However, in our proposed scheme, by leveraging the partial CSI of 𝐇^bs{\bf{\hat{H}}}_{bs} obtained in Phase I, the efficiency in estimating 𝐇bsuk{\bf{H}}_{bsu_{k}} can be greatly improved. Specifically, we consider the following two cases.

Case 1: MNM\geq N. In this case, users send pilot symbols consecutively for the BS to estimate 𝐡suk{\bf{h}}_{su_{k}}, independently. Let 𝚽k=𝐈{\bf{\Phi}}_{k}={\bf{I}}, xk,k=1x_{k,k}=1 and xk1,k=0,k1kx_{k_{1},k}=0,k_{1}\neq k. Then, 𝐲¯b(k){\bf{\bar{y}}}_{b}^{(k)} is rewritten as

𝐲¯b(k)=p𝐇^bs𝐡suk+𝐳~b(k),\displaystyle{\bf{\bar{y}}}_{b}^{(k)}=\sqrt{p}{\bf{\hat{H}}}_{bs}{\bf{h}}_{su_{k}}+{\bf{\tilde{z}}}_{b}^{(k)}, (18)

where 𝐳~b(k)=𝐳b(k)+p(𝐡buk𝐡^buk+(𝐇bsuk𝐇^bsuk)𝐯k){\bf{\tilde{z}}}_{b}^{(k)}={\bf{z}}_{b}^{(k)}+\sqrt{p}\left({\bf{h}}_{bu_{k}}-{\bf{\hat{h}}}_{bu_{k}}+({\bf H}_{bsu_{k}}-{\bf{\hat{H}}}_{bsu_{k}}){\bf{v}}_{k}\right) and 𝐯k=[1,1,,1]T{{\bf v}_{k}}=[1,1,...,1]^{T}.

Since 𝐇^bs(m,n){\bf{\hat{H}}}_{bs}(m,n), m,n\forall m,n, has two possible values, 𝐡suk{\bf{h}}_{su_{k}} cannot be uniquely estimated. However, we show that the estimation of the cascaded channel 𝐇bsuk=𝐇bsdiag(𝐡suk){\bf H}_{bsu_{k}}={\bf H}_{bs}\text{diag}({\bf h}_{su_{k}}) can be uniquely obtained by resorting to the following proposition.

Proposition 1.

By setting 𝐇^bs{\bf{\hat{H}}}_{bs} as 𝐖{\bf{W}} where 𝐖(1,n)=g1n{\bf{W}}(1,n)=g_{1n} and 𝐖(m,n)=𝐇^bsa1(m,n)𝐇^bsa1(1,n)g1n{\bf{W}}(m,n)=\frac{{\bf{\hat{H}}}_{bsa_{1}}(m,n)}{{\bf{\hat{H}}}_{bsa_{1}}(1,n)}g_{1n}, m,n\forall m,n, 𝐡suk{\bf{h}}_{su_{k}} can be estimated as

𝐡^suk=1p(𝐖H𝐖)1𝐖H𝐲¯b(k),k=1,,K.\displaystyle{\bf{\hat{h}}}_{su_{k}}=\frac{1}{\sqrt{p}}({\bf{W}}^{H}{\bf{W}})^{-1}{\bf{W}}^{H}{\bf{\bar{y}}}_{b}^{(k)},k=1,...,K. (19)

Then, 𝐇bsuk{\bf{H}}_{bsu_{k}} is uniquely estimated as

𝐇^bsuk=𝐇^bsdiag(𝐡^suk).{\bf{\hat{H}}}_{bsu_{k}}={\bf{\hat{H}}}_{bs}{\rm{diag}}({\bf{\hat{h}}}_{su_{k}}). (20)

Proof: Please see Appendix A. \square

Based on Proposition 1, it takes the BS KK pilot symbols to estimate all BS-IRS-user cascaded channels.

Case 2: M<NM<N. In this case, if all users transmit pilot symbols one by one, it takes the BS at least NM\left\lceil{\frac{N}{M}}\right\rceil pilot symbols to estimate each 𝐡suk{\bf{h}}_{su_{k}} and thus the total overhead is KNMK\left\lceil{\frac{N}{M}}\right\rceil. Thus, we propose a scheme based on orthogonal pilot symbols and orthogonal phase shifts for reducing the total overhead to KNM\left\lceil{\frac{KN}{M}}\right\rceil, which is detailed as follows.

Step (a): Divide the KK users into L1=KML_{1}=\left\lceil{\frac{K}{M}}\right\rceil groups, such that there are MM users in each of the first L11L_{1}-1 groups and M1M_{1} (M1MM_{1}\leq M) users in the last group, i.e., K=(L11)M+M1K=(L_{1}-1)M+M_{1}.

Step (b): For each of the first L11L_{1}-1 groups, NN pilot symbols are used to estimate the BS-IRS-user cascaded channels. Take the first group as an example. The received signals at the BS during NN time slots by removing those from the direct path and neglecting the noise can be written as (21) shown at the bottom of this page, where 𝚽i=diag(𝐯i){\bf{\Phi}}_{i}=\text{diag}({\bf v}_{i}) is the phase-shift matrix in the iith time slot and xk,ix_{k,i} is the pilot symbol transmitted by the kkth user in the iith time slot, kMk\leq M, iNi\leq N. As long as we design Φi\Phi_{i} and xk,ix_{k,i} properly such that 𝐁\bf B given in (21) is full-rank, the IRS-user channels can be estimated as

[𝐡^su1𝐡^su2..𝐡^suM]=1p(𝐁H𝐁)1𝐁H[𝐲¯b(1)𝐲¯b(2)..𝐲¯b(N)].\small\left[{\begin{array}[]{*{20}{c}}{{{\bf{\hat{h}}}_{s{u_{1}}}}}\\ {{{\bf{\hat{h}}}_{s{u_{2}}}}}\\ \begin{array}[]{l}.\\ .\end{array}\\ {{{\bf{\hat{h}}}_{s{u_{M}}}}}\end{array}}\right]=\frac{1}{\sqrt{p}}{\left({{{\bf{B}}^{H}}{\bf{B}}}\right)^{-1}}{{\bf{B}}^{H}}\left[{\begin{array}[]{*{20}{c}}{{{\bf{\bar{y}}}_{b}^{(1)}}}\\ {{{\bf{\bar{y}}}_{b}^{(2)}}}\\ \begin{array}[]{l}.\\ .\end{array}\\ {{{\bf{\bar{y}}}_{b}^{\left(N\right)}}}\end{array}}\right]. (22)

For example, we can design the pilot symbols transmitted by the MM users during the NN time slots as

𝐗=[11111ejθej2θej(M1)θ.......1ej(N1)θej2(N1)θej(M1)(N1)θ],\small{\bf X}\!=\!\left[{\begin{array}[]{*{20}{c}}1&1&1&{...}&1\\ 1&{{e^{-j\theta}}}&{{e^{-j2\theta}}}&{...}&{{e^{-j\left({M-1}\right)\theta}}}\\ .&.&.&{...}&.\\ .&.&{}\hfil&.&{}\hfil\\ 1&{{e^{-j\left({N-1}\right)\theta}}}&{{e^{-j2\left({N-1}\right)\theta}}}&{...}&{{e^{-j{\left({M-1}\right){\left({N-1}\right)}}\theta}}}\end{array}}\right],

where 𝐗ik=xk,i{\bf X}_{ik}=x_{k,i} is the pilot symbol transmitted by User kk in the ii-th time slot, while the phase shifts during the NN time slots are given by

𝐕=[11111ejθej2θej(N1)θ.......1ej(N1)θej2(N1)θej(N1)2θ],\small{\bf V}=\left[{\begin{array}[]{*{20}{c}}1&1&1&{...}&1\\ 1&{{e^{-j\theta}}}&{{e^{-j2\theta}}}&{...}&{{e^{-j\left({N-1}\right)\theta}}}\\ .&.&.&{...}&.\\ .&.&{}\hfil&.&{}\hfil\\ 1&{{e^{-j\left({N-1}\right)\theta}}}&{{e^{-j2\left({N-1}\right)\theta}}}&{...}&{{e^{-j{{\left({N-1}\right)^{2}}}\theta}}}\end{array}}\right],

where 𝐕ni=vn,i{\bf V}_{ni}=v_{n,i} is the phase shift of the nn-th reflecting element in the ii-th time slot.

It should be noted that similar to Case 1, the exact 𝐁\bf B is unknown due to the fact that we only have the estimation of 𝐇bs𝐇bs{\bf{H}}_{bs}\odot{\bf{H}}_{bs}, instead of 𝐇bs{\bf{H}}_{bs}. However, it is shown in Proposition 1 that the cascaded channels can be uniquely estimated as

[𝐇^bsu1,,𝐇^bsuM]=𝐁[diag(𝐡^su1),,diag(𝐡^suM)].[{\bf{\hat{H}}}_{bsu_{1}},...,{\bf{\hat{H}}}_{bsu_{M}}]={\bf{B}}[{\rm{diag}}({\bf{\hat{h}}}_{su_{1}}),...,{\rm{diag}}({\bf{\hat{h}}}_{su_{M}})].

Step (c): For the last group, if M1=MM_{1}=M, the process is the same as that in Step (b) and another NN pilot symbols are needed. In this case, we have L1=KML_{1}=\frac{K}{M} and thus the total overhead is L1N=KNML_{1}N=\frac{KN}{M}.

On the other hand, if M1<MM_{1}<M, another N1=M1NMN_{1}=\left\lceil{\frac{M_{1}N}{M}}\right\rceil pilot symbols are required, M1<N1<NM_{1}<N_{1}<N. Specifically, 𝐗\bf X and 𝐕\bf V are given by

𝐗=[11111ejθej2θej(M11)θ.......1ej(N11)θej2(N11)θej(M11)(N11)θ],\small{\bf X}\!\!=\!\!\left[{\begin{array}[]{*{20}{c}}1&1&1&{...}&1\\ 1&{{e^{-j\theta}}}&{{e^{-j2\theta}}}&{...}&{{e^{-j\left({M_{1}-1}\right)\theta}}}\\ .&.&.&{...}&.\\ .&.&{}\hfil&.&{}\hfil\\ 1&{{e^{-j\left({N_{1}-1}\right)\theta}}}&{{e^{-j2\left({N_{1}-1}\right)\theta}}}&{...}&{{e^{-j{\left({M_{1}\!-\!1}\right){\left({N_{1}\!-\!1}\right)}}\theta}}}\end{array}}\right],
𝐕=[11111ejθej2θej(N11)θ.......1ej(N1)θej2(N1)θej(N11)(N1)θ],\small{\bf V}\!=\!\left[{\begin{array}[]{*{20}{c}}1&1&1&{...}&1\\ 1&{{e^{-j\theta}}}&{{e^{-j2\theta}}}&{...}&{{e^{-j\left({N_{1}-1}\right)\theta}}}\\ .&.&.&{...}&.\\ .&.&{}\hfil&.&{}\hfil\\ 1&{{e^{-j\left({N-1}\right)\theta}}}&{{e^{-j2\left({N-1}\right)\theta}}}&{...}&{{e^{-j{{\left({N_{1}-1}\right)\left({N-1}\right)}}\theta}}}\end{array}}\right],

and thus 𝐁𝐌𝐍𝟏×𝐌𝟏𝐍\bf B\in{\mathbb{C}}^{MN_{1}\times M_{1}N} is re-written as

𝐁=[𝐇bs𝚽1x1,1𝐇bs𝚽1x2,1𝐇bs𝚽1xM1,1𝐇bs𝚽2x1,2𝐇bs𝚽2x2,2𝐇bs𝚽2xM1,2......𝐇bs𝚽N1x1,N1𝐇bs𝚽N1x2,N1𝐇bs𝚽N1xM1,N1].\small{\bf B}\!\!=\!\!\!\left[{\begin{array}[]{*{20}{c}}{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{1,1}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{2,1}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{1}}{x_{M_{1},1}}}\\ {{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{1,2}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{2,2}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{2}}{x_{M_{1},2}}}\\ .&.&{...}&.\\ .&.&.&{}\hfil\\ {{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N_{1}}}{x_{1,N_{1}}}}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N_{1}}}{x_{2,N_{1}}}}&{...}&{{{\bf{H}}_{bs}}{{\bf{\Phi}}_{N_{1}}}{x_{M_{1},N_{1}}}}\end{array}}\right]\!.

It can be verified that MN1>M1NMN_{1}>M_{1}N and 𝐁\bf B is column full-rank, i.e., rank(𝐁)=M1N\text{rank}({\bf B})=M_{1}N, such that the left inverse of 𝐁\bf B exists and thus the IRS-user channels can be estimated similarly as (22). The corresponding overhead is (L11)N+M1NM=M(L11)N+M1NM=KNM(L_{1}-1)N+\left\lceil{\frac{M_{1}N}{M}}\right\rceil=\left\lceil{\frac{M(L_{1}-1)N+M_{1}N}{M}}\right\rceil=\left\lceil{\frac{KN}{M}}\right\rceil.

As a result, the overhead in estimating the cascaded BS-IRS-user channels for the case of M<NM<N is KNM\left\lceil{\frac{KN}{M}}\right\rceil.

Refer to caption
Figure 3: Simulation Setup.

III-C Overall Training Overhead

To summarize, the total training overhead of the proposed scheme is 2(N+1+K)2(N+1+K) for MNM\geq N and 2(N+1)+K+KNM2(N+1)+K+\left\lceil\frac{KN}{M}\right\rceil for M<NM<N, i.e. 2(N+1)+K+max(K,KNM)2(N+1)+K+\text{max}(K,\left\lceil\frac{KN}{M}\right\rceil).

IV Special Case: LoS IRS-Anchor Channel

We assume that the IRS and/or the anchor node “A” can be properly deployed so that the IRS-A channel is LoS. In this case, 𝐡ra{\bf{h}}_{ra} is known a priori based on the knowledge of the positions of IRS and A, thus only one anchor node is sufficient for the proposed scheme and its training overhead can be further reduced.

IV-A Phase I: Off-line Estimation of 𝐇bs{\bf{H}}_{bs}

Let anchor A transmit N+1N+1 pilot symbols while BS estimates 𝐡ba{\bf h}_{ba} and 𝐇bsa=𝐇bsdiag(𝐡ra){\bf H}_{bsa}={\bf H}_{bs}\text{diag}({\bf h}_{ra}), respectively. Since 𝐡ra{\bf{h}}_{ra} is known, 𝐇bs{\bf{H}}_{bs} can be recovered from 𝐇^bsa{\bf{\hat{H}}}_{bsa} as

𝐇^bs=𝐇^bsadiag(𝐡ra)1.{\bf{\hat{H}}}_{bs}={\bf{\hat{H}}}_{bsa}{\rm{diag}}({\bf{h}}_{ra})^{-1}.\vspace{-2mm} (23)

IV-B Phase II: On-line Estimation of 𝐡buk{\bf{h}}_{bu_{k}} and 𝐇bsuk{\bf{H}}_{bsu_{k}}

Users transmit pilot symbols while the BS estimates 𝐡buk{\bf{h}}_{bu_{k}} and 𝐇bsuk{\bf{H}}_{bsu_{k}}, respectively. Specifically, we estimate the channel 𝐡suk{\bf{h}}_{su_{k}} first based on the estimation of 𝐇^bs{\bf{\hat{H}}}_{bs} and then obtain 𝐇bsuk=𝐇^bsdiag(𝐡^suk){\bf{H}}_{bsu_{k}}={\bf\hat{H}}_{bs}\text{diag}({\bf\hat{h}}_{su_{k}}), similar to Proposition 1. The only difference is that we have 𝐇^bs{\bf{\hat{H}}}_{bs} in this case instead of 𝐇^bs𝐇^bs{\bf{\hat{H}}}_{bs}\odot{\bf{\hat{H}}}_{bs}.

Combining Phases I and II, the total training overhead of the proposed scheme in this special case is N+1+2KN+1+2K for MNM\geq N and N+1+K+KNMN+1+K+\left\lceil\frac{KN}{M}\right\rceil for M<NM<N, i.e. N+1+K+max(K,KNM)N+1+K+\text{max}(K,\left\lceil\frac{KN}{M}\right\rceil).

Refer to caption
(a) Training time versus MM, with K=20K=20.
Refer to caption
(b) Normalized training time versus KK.
Refer to caption
(c) MSE versus pp, with K=20K=20.
Figure 4: Performance comparison between the proposed scheme and the benchmark scheme [13].

V Numerical Results

The simulation setup is shown in Fig. 3. It is assumed that BS, IRS (the central point), A1, A2 and A are located at (50, 0, 20), (0, 100, 2), (2, 99, 0), (2, 101, 0) and (2, 100, 0) in meter (m), respectively. We assume that the system operates on a carrier frequency of 750 MHz with the wavelength λc=0.4\lambda_{c}=0.4 m and the path loss at the reference distance d0=1d_{0}\!=\!1 m is given by L0=30L_{0}=30 dB. Suppose that the IRS is equipped with a uniform planar array with 6 rows and 10 columns, and the element spacing is Δd=3λc/8\Delta d\!=\!3\lambda_{c}/8; thus, we have N=60N=60. The noise power is set as σ02=\sigma_{0}^{2}\!=\!-105 dBm. The channel from the BS to A1 is generated by 𝐡ba1=L0dba1cba1𝐠ba1{{\mathbf{h}}_{ba_{1}}}\!=\!\sqrt{{L_{0}}d_{ba_{1}}^{-{c_{ba_{1}}}}}{{\bf{g}}_{ba_{1}}}, where dba1d_{ba_{1}} denotes the distance from the BS to A1 and 𝐠ba1{{\bf{g}}_{ba_{1}}} is the small-scale fading component. The same channel model is adopted for all other channels in general. Particularly, Rayleigh fading is assumed for the channels among the BS, IRS, A1, A2 and each User kk with the path loss exponents set as 3, whereas in the special channel case in Section IV, the channel between IRS and A is assumed to be LoS, with the path loss exponent set as 2. We also apply the scheme proposed in [13] for each of the users to estimate their channels consecutively, which serves as the benchmark.

Fig. 4(a) shows the required training time (in terms of number of pilot symbols) versus the number of antennas, MM, at the BS. It is observed that as MM increases, the proposed scheme significantly reduce the training overhead as compared to that of the benchmark scheme (which is independent of MM). This is because the proposed scheme exploits the multiple antennas at the BS for joint IRS channel estimation, whereas in the benchmark scheme the BS antennas estimate their associated channels independently in parallel. Note that when M=1M=1, the training overhead of the proposed scheme is even larger than that of the benchmark scheme. This is because additional pilot symbols are transmitted by anchors in Phase I, while the training efficiency in Phase II is not improved since max(K,NKM)=NK\text{max}(K,\left\lceil\frac{NK}{M}\right\rceil)=NK in the case of M=1M=1. Moreover, it is observed that the proposed scheme under the special case of LoS IRS-anchor channel is more efficient as compared to the general channel case.

Fig. 4(b) shows the training time of the proposed scheme normalized by that of the benchmark scheme versus KK, with M=10M=10 and 60, respectively. One can observe that the performance gap between the two schemes becomes larger as KK increases. This is because to accommodate one more user, the additional pilot overhead required by the benchmark scheme is N+1N+1, while that by the proposed scheme is 1+max(1,NM)1+\text{max}(1,\left\lceil\frac{N}{M}\right\rceil). As a result, as KK increases, the pilot reduction by using the proposed scheme also increases. Also note that similar to Fig. 4(a), when KK is very small, the proposed scheme is even worse than the benchmark scheme. This reason is that additional 2(N+1)2(N+1)/N+1N+1 pilot symbols are required in Phase I, regardless of KK, while in the benchmark scheme, only N+1N+1 pilot symbols are sufficient when K=1K=1.

Fig. 4(c) shows the normalized mean-squared error (MSE) of the estimations of 𝐡buk{\bf h}_{bu_{k}} and 𝐇bsuk{\bf H}_{bsu_{k}} versus the transmit power of pilot symbols in the on-line phase, with that of the off-line phase fixed as 40 dBm. It is observed that the MSE in the general channel case of the proposed scheme is highest, while that of the benchmark scheme is lowest. The reason is that although the proposed scheme significantly reduces the training overhead (see Figs. 4(a) and 4(b)), the estimation error in 𝐇^bsuk{\bf\hat{H}}_{bsu_{k}} depends on both 𝐇^bs{\bf\hat{H}}_{bs} and 𝐡^suk{\bf\hat{h}}_{su_{k}}. Specifically, in the general channel case, the error in 𝐇^bs𝐇^bs{\bf\hat{H}}_{bs}\odot{\bf\hat{H}}_{bs} comes from 𝐇^bsa1{\bf\hat{H}}_{bsa_{1}}, 𝐇^bsa2{\bf\hat{H}}_{bsa_{2}} and 𝐡^a1sa2{\bf\hat{h}}_{a_{1}sa_{2}}, while in the special channel case the error in 𝐇^bs{\bf\hat{H}}_{bs} comes from 𝐇^bsa{\bf\hat{H}}_{bsa}. In contrast, for the benchmark scheme, 𝐇bsuk{\bf H}_{bsu_{k}} is estimated directly, which is thus less susceptible to noise/error. However, one can observe that the MSE of the proposed scheme is substantially reduced by increasing MM, and becomes even comparable when M=100M=100.

VI Conclusion

In this paper, we propose a new anchor-assisted channel estimation scheme for IRS-aided multiuser communications. By exploiting the fact that all BS-IRS-user cascaded channels share the same BS-IRS channel, the proposed scheme first estimates this common channel with only sign ambiguity via the anchor-assisted training. Then we show that the estimation of each cascaded BS-IRS-user channel is simplified to estimating each IRS-user channel with the number of unknowns significantly reduced from MNMN to NN, and the sign ambiguity in the estimated common channel does not affect the uniqueness of the recovered cascaded channels. Moreover, by exploring multi-antennas at the BS, the training overhead in estimating all users’ cascaded channels is reduced from NKNK to max(K,NKM)\text{max}(K,\left\lceil\frac{NK}{M}\right\rceil). Numerical results validate the effectiveness of the proposed scheme, especially when MM and/or KK is large. Considering the trend towards massive antenna arrays at the BS and massive connectivity with machine-type communications, our proposed scheme has the great potential of significantly improving the channel estimation efficiency in future IRS-aided wireless systems.

Appendix A

First, we show the following lemma.

Lemma 1: Given 𝐇^bsa1{\bf{\hat{H}}}_{bsa_{1}}, 𝐇^bs(m,n){\bf{\hat{H}}}_{bs}(m,n) can be estimated as 𝐇^bs(m,n)=αmn𝐇^bs(1,n){\bf{\hat{H}}}_{bs}(m,n)=\alpha_{mn}{\bf{\hat{H}}}_{bs}(1,n), where αmn=𝐇^bsa1(m,n)𝐇^bsa1(1,n)\alpha_{mn}=\frac{{\bf{\hat{H}}}_{bsa_{1}}(m,n)}{{\bf{\hat{H}}}_{bsa_{1}}(1,n)}, m\forall m.

Proof: Because 𝐇bsa1(1,n)=𝐇bs(1,n)𝐡sa1(n){\bf{H}}_{bsa_{1}}(1,n)={{\bf{H}}_{bs}(1,n)}{{\bf{h}}_{sa_{1}}(n)} and 𝐇bsa1(m,n)=𝐇bs(m,n)𝐡sa1(n){\bf{H}}_{bsa_{1}}(m,n)={{\bf{H}}_{bs}(m,n)}{{\bf{h}}_{sa_{1}}(n)}, we have 𝐇bs(m,n)=𝐇bsa1(m,n)𝐇bsa1(1,n)𝐇bs(1,n){\bf{H}}_{bs}(m,n)=\frac{{\bf{H}}_{bsa_{1}}(m,n)}{{\bf{H}}_{bsa_{1}}(1,n)}{\bf{H}}_{bs}(1,n), which thus completes the proof. \square

Lemma 1 reveals that for given 𝐇^bsa1{\bf{\hat{H}}}_{bsa_{1}}, there are NN rather than MNMN unknowns in 𝐇^bs{\bf{\hat{H}}}_{bs} and it can be rewritten as

𝐇^bs=[𝐇bs(1,1)𝐇bs(1,2)𝐇bs(1,N)α21𝐇bs(1,1)α22𝐇bs(1,2)α2N𝐇bs(1,N)......αM1𝐇bs(1,1)αM2𝐇bs(1,2)αMN𝐇bs(1,N)].\small{\bf{\hat{H}}}_{bs}\!=\!\left[\begin{aligned} &~{}~{}{{\bf{H}}_{bs}(1,1)}&{{\bf{H}}_{bs}(1,2)}~{}~{}&~{}~{}...&{{\bf{H}}_{bs}(1,N)}~{}~{}\\ &\alpha_{21}{{\bf{H}}_{bs}(1,1)}&\alpha_{22}{{\bf{H}}_{bs}(1,2)}&~{}~{}...&\alpha_{2N}{{\bf{H}}_{bs}(1,N)}\\ &~{}~{}~{}~{}~{}~{}~{}.&.~{}~{}~{}~{}~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}~{}~{}~{}\\ &~{}~{}~{}~{}~{}~{}~{}.&.~{}~{}~{}~{}~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}~{}~{}~{}\\ &\alpha_{M1}{{\bf{H}}_{bs}(1,1)}&\alpha_{M2}{{\bf{H}}_{bs}(1,2)}&~{}~{}...&\alpha_{MN}{{\bf{H}}_{bs}(1,N)}\end{aligned}\right].

Meanwhile, we can construct a candidate of 𝐇bs{\bf H}_{bs} as 𝐖{\bf W} in Proposition 1, given by

𝐖=[g11g12g1Nα21g11α22g12α2Ng1N......αM1g11αM2g12αMNg1N].\small{\bf{W}}=\left[\begin{aligned} &~{}~{}g_{11}&g_{12}~{}~{}&~{}~{}...&g_{1N}~{}~{}\\ &\alpha_{21}g_{11}&\alpha_{22}g_{12}&~{}~{}...&\alpha_{2N}g_{1N}\\ &~{}~{}~{}~{}.&.~{}~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}\\ &~{}~{}~{}~{}.&.~{}~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}\\ &\alpha_{M1}g_{11}&\alpha_{M2}g_{12}&~{}~{}...&\alpha_{MN}g_{1N}\end{aligned}\right]. (24)

Next, we consider the following two cases.

Case 1: 𝐇bs=𝐖{\bf H}_{bs}={\bf W}. In this case, we have 𝐇bs(1,n)=g1n,n{{\bf{H}}_{bs}(1,n)}=g_{1n},\forall n. Omitting the noise, we can express (18) as

𝐲¯b(k)=p𝐖𝐡suk.{\bf{\bar{y}}}_{b}^{(k)}=\sqrt{p}{\bf{W}}{\bf{h}}_{su_{k}}.\vspace{-1mm} (25)

Assuming that 𝐖{\bf W} is full-rank, 𝐡suk{\bf{h}}_{su_{k}} can be estimated as

𝐡^suk=1p(𝐖H𝐖)1𝐖H𝐲¯b(k),k=1,,K.{\bf{\hat{h}}}_{su_{k}}=\frac{1}{\sqrt{p}}({\bf{W}}^{H}{\bf{W}})^{-1}{\bf{W}}^{H}{\bf{\bar{y}}}_{b}^{(k)},k=1,...,K. (26)

Then, the cascaded BS-IRS-user channel is estimated as

𝐇^bsuk=𝐖diag(𝐡^suk)=𝐇^bsdiag(𝐡suk).{\bf{\hat{H}}}_{bsu_{k}}={\bf{W}}\text{diag}({\bf{\hat{h}}}_{su_{k}})={\bf{\hat{H}}}_{bs}\text{diag}({\bf{h}}_{su_{k}}). (27)

Case 2: 𝐇bs𝐖{\bf H}_{bs}\neq{\bf W}. Referring to (15), there must exist at least an nn such that 𝐇bs(1,n)=g1n{{\bf{H}}_{bs}(1,n)}=-g_{1n}. For illustration purpose, we assume that 𝐇bs(1,1)=g11{{\bf{H}}_{bs}(1,1)}=-g_{11} and 𝐇bs(1,n)=g1n{{\bf{H}}_{bs}(1,n)}=g_{1n} for n1\forall n\neq 1, while for all other possible 𝐇bs{{\bf{H}}_{bs}}’s, the result can be similarly proved. Then we have

𝐇^bs=[g11g12g1Nα21g11α22g12α2Ng1N......αM1g11αM2g12αMNg1N].\small{\bf{\hat{H}}}_{bs}=\left[\begin{aligned} &~{}-g_{11}&g_{12}~{}~{}&~{}~{}...&g_{1N}~{}~{}\\ &-\alpha_{21}g_{11}&\alpha_{22}g_{12}&~{}~{}...&\alpha_{2N}g_{1N}\\ &~{}~{}~{}~{}~{}~{}.&.~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}\\ &~{}~{}~{}~{}~{}~{}.&.~{}~{}~{}~{}&~{}~{}...&.~{}~{}~{}~{}~{}\\ &-\alpha_{M1}g_{11}&\alpha_{M2}g_{12}&~{}~{}...&\alpha_{MN}g_{1N}\end{aligned}\right]. (28)

By omitting the noise, (18) can be written as

𝐲¯b(k)=p𝐇^bs𝐡suk.{\bf{\bar{y}}}_{b}^{(k)}=\sqrt{p}{\bf{\hat{H}}}_{bs}{\bf{h}}_{su_{k}}. (29)

Accordingly, by using 𝐖{\bf W}, 𝐡suk{\bf h}_{su_{k}} is estimated as

𝐡^suk(2)=(𝐖H𝐖)1𝐖H𝐇^bs𝐡suk,k=1,,K.{\bf{\hat{h}}}_{su_{k}}^{(2)}=({\bf{W}}^{H}{\bf{W}})^{-1}{\bf{W}}^{H}{\bf{\hat{H}}}_{bs}{\bf{h}}_{su_{k}},k=1,...,K. (30)

Comparing 𝐇^bs{\bf{\hat{H}}}_{bs} in (28) with 𝐖{\bf{W}}, we observe that the only difference lies in the sign of the elements in the first column and thus the following equality holds

𝐇^bs=𝐖diag([1,1,1,,1]).\vspace{-1mm}{\bf{\hat{H}}}_{bs}={\bf{W}}\text{diag}([-1,1,1,...,1]). (31)

Since we have

[(𝐖^H𝐖)1𝐖H]𝐖=𝐈N,[({\bf{\hat{W}}}^{H}{\bf{W}})^{-1}{\bf{W}}^{H}]{\bf{W}}={\bf{I}}_{N},\vspace{-0mm} (32)

substituting (31) into (32) yields

[(𝐖H𝐖)1𝐖2H]𝐇^bs=diag([1,1,1,,1]).[({\bf{W}}^{H}{\bf{W}})^{-1}{\bf{W}}_{2}^{H}]{\bf{\hat{H}}}_{bs}=\text{diag}([-1,1,1,...,1]).\vspace{-0mm} (33)

Based on (30) and (33), the estimation of 𝐡suk{\bf{h}}_{su_{k}} obtained by using 𝐖{\bf W} can be written as

𝐡^suk(2)=diag([1,1,1,,1])𝐡suk,k=1,,K.{\bf{\hat{h}}}_{su_{k}}^{(2)}=\text{diag}([-1,1,1,...,1]){\bf{h}}_{su_{k}},k=1,...,K.\vspace{-0mm} (34)

Though 𝐡^suk(2){\bf{\hat{h}}}_{su_{k}}^{(2)} is not an exact estimation, the error only occurs in the sign of the first element of 𝐡suk{\bf{h}}_{su_{k}}. Finally, the cascaded BS-IRS-user channel is recovered by

𝐇^bsuk=𝐖diag(𝐡^suk(2))=𝐇^bsdiag(𝐡suk),k=1,,K,{\bf{\hat{H}}}_{bsu_{k}}={\bf{W}}\text{diag}({\bf{\hat{h}}}_{su_{k}}^{(2)})={\bf{\hat{H}}}_{bs}\text{diag}({\bf{h}}_{su_{k}}),k=1,...,K,\vspace{-1mm}

which is the same as (27).

Based on Cases 1 and 2, it is concluded that using 𝐖{\bf W} constructed in Proposition 1 to estimate 𝐇bsuk{\bf{H}}_{bsu_{k}} is always sufficient, regardless of whether 𝐖{\bf W} is exactly the same as 𝐇^bs{\bf\hat{H}}_{bs}, which thus completes the proof of Proposition 1.

References

  • [1] Q. Wu and R. Zhang, “Towards smart and reconfigurable environment: Intelligent reflecting surface aided wireless network,” IEEE Commun. Mag., vol. 58, no. 1, pp. 106-112, Jan. 2020.
  • [2] S. Gong et al., “Towards smart radio environment for wireless communications via intelligent reflecting surfaces: A comprehensive survey,” 2019. [Online]. Available: https://arxiv.org/abs/1912.07794.
  • [3] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network via joint active and passive beamforming,” IEEE Trans. Wireless Commun., vol. 18, no. 11, pp. 5394-5409, Nov. 2019.
  • [4] B. Zheng, Q. Wu, and R. Zhang, “Intelligent reflecting surface-assisted multiple access with user pairing: NOMA or OMA?” IEEE Commun. Lett., vol. 24, no. 4, pp. 753-757, Apr. 2020.
  • [5] Z. Ding and H. V. Poor, “A simple design of IRS-NOMA transmission,” IEEE Commun. Lett., vol. 24, no. 5, pp. 1119-1123, May 2020.
  • [6] Q. Wu and R. Zhang, “Joint active and passive beamforming optimization for intelligent reflecting surface assisted SWIPT under QoS constraints,” IEEE J. Sel. Areas Commun., 2019, to appear. [Online]. Available: https://arxiv.org/abs/1910.06220.
  • [7] C. Pan et al., “Intelligent reflecting surface enhanced MIMO broadcasting for simultaneous wireless information and power transfer,” 2019. [Online]. Available: http://arxiv.org/abs/1909.03272.
  • [8] X. Guan, Q. Wu, and R. Zhang, “Intelligent reflecting surface assisted secrecy communication: Is artificial noise helpful or not?” IEEE Wireless Commun. Lett., Jan. 2020, DOI: 10.1109/LWC.2020.2969629.
  • [9] H. Shen, W. Xu, S. Gong, Z. He, and C. Zhao, “Secrecy rate maximization for intelligent reflecting surface assisted multi-antenna communications,” IEEE Commun. Lett., vol. 23, no. 9, pp. 1488-1492, Sep. 2019.
  • [10] Y. Yang, B. Zheng, S. Zhang, and R. Zhang, “Intelligent reflecting surface meets OFDM: Protocol design and rate maximization,” 2019. [Online]. Available: https://arxiv.org/abs/1906.09956.
  • [11] D. Mishra and H. Johansson, “Channel estimation and low-complexity beamforming design for passive intelligent surface assisted MISO wireless energy transfer,” in Proc. IEEE ICASSP, May 2019.
  • [12] C. You, B. Zheng, and R. Zhang, “Channel estimation and passive beamforming for intelligent reflecting surface: discrete phase shift and progressive refinement,” IEEE J. Sel. Areas Commun., 2019, to appear. [Online] Available: https://arxiv.org/abs/1912.10646.
  • [13] B. Zheng and R. Zhang, “Intelligent reflecting surface enhanced OFDM: channel estimation and reflection optimization,” IEEE Wireless Commun. Lett., vol. 9, no. 4, pp. 518-522, Apr. 2020.