Introduction
在做 system model 和 problem formulation 的时候,对于信道的构成以及标准表达格式理解的不太清楚,是以想要完成这篇总结,后续随着理解加深会不断补充其他情况。每个人对于信道模型的理解有所不同,但是对于处于研一阶段的我来说,由于线性代数以及矩阵论的学习还不深入,在对于所看论文里的系统模型以及传输协议进行理解的过程中,我认为最重要的就是保证维度的匹配,这也是在我完成第一个 system model 时耗时甚久的重要原因。
Notation and Identities
首先,需要先对不同的符号以及是否加粗进行定义:(在初学阶段,格式的统一很重要)
- 标量(scalar)不需要加粗,用小写字母表示;
- column vectors 需要加粗表示,用小写字母表示;
- matrix 用大写字母加粗表示;
- ⊗ \otimes ⊗ 代表Kronecker product;
- A T A^T AT代表矩阵的转置,在信道的分解中主要目的是匹配维度,不过同时还有 A H A^H AH代表Hermitian vector, 目前我没搞懂这俩对于信号的实际意义有什么区别。
Channel model without IRSs
考虑一个MIMO系统,输入和输出的信号遵循以下关系:
y
=
H
x
+
v
∈
C
M
r
×
1
\begin{equation} \mathbf{y}=\mathbf{H x}+\mathbf{v} \in \mathbb{C}^{M_r \times 1} \end{equation}
y=Hx+v∈CMr×1
其中,
x
∈
C
M
t
×
1
\mathbf{x} \in \mathbb{C}^{M_t \times 1}
x∈CMt×1 是传输的信号,
v
v
v 是RX端的Additive White Gaussian Noise(AWGN),
H
∈
C
M
r
×
M
t
\mathbf{H} \in \mathbb{C}^{M_r \times M_t}
H∈CMr×Mt 可以被modeled 成如下形式:
H
=
∑
p
=
1
P
γ
p
a
p
(
r
x
)
b
p
(
t
x
)
T
∈
C
M
r
×
M
t
\begin{equation} \mathbf{H}=\sum_{p=1}^P \gamma_p \mathbf{a}_p^{(\mathrm{rx})} \mathbf{b}_p^{(\mathrm{tx}) \mathrm{T}} \in \mathbb{C}^{M_r \times M_t} \end{equation}
H=p=1∑Pγpap(rx)bp(tx)T∈CMr×Mt
其中,
γ
p
\gamma_p
γp 是path loss 和fast-fading components,
a
p
(
r
x
)
\mathbf{a}_p^{(\mathrm{rx})}
ap(rx) 和
b
p
(
t
x
)
\mathbf{b}_p^{(\mathrm{tx})}
bp(tx) 分别是
p
p
pth steering vectors,它们的定义如下:
a
p
(
r
x
)
=
[
1
,
…
,
e
−
j
π
(
M
r
−
1
)
cos
(
α
p
)
]
T
∈
C
M
r
×
1
\begin{equation} \mathbf{a}_p^{(\mathrm{rx})}=\left[1, \ldots, e^{-j \pi\left(M_r-1\right) \cos \left(\alpha_p\right)}\right]^{\mathrm{T}} \in \mathbb{C}^{M_r \times 1} \end{equation}
ap(rx)=[1,…,e−jπ(Mr−1)cos(αp)]T∈CMr×1
b p ( t x ) = [ 1 , … , e − j π ( M t − 1 ) cos ( β p ) ] T ∈ C M t × 1 \begin{equation} \mathbf{b}_p^{(\mathrm{tx})}=\left[1, \ldots, e^{-j \pi\left(M_t-1\right) \cos \left(\beta_p\right)}\right]^{\mathrm{T}} \in \mathbb{C}^{M_t \times 1} \end{equation} bp(tx)=[1,…,e−jπ(Mt−1)cos(βp)]T∈CMt×1
其中,
α
p
\alpha_p
αp and
β
p
\beta_p
βp 分别表示的是
p
p
p 条路径的Angle of Arrival(AoA)和Angle of Departure(AoD)。以上的信道模型还未考虑IRS,接下来让我们研究加入多个IRS后的信道模型。
Channel model with IRSs
y k ( l ) = ∑ p = 1 P γ p a p ( r x ) a p ( i r s ) T diag k ( S ( p , l ) ) b p ( i r s ) b p ( t x ) T x k ( l ) + v k ( l ) ∈ C M r × 1 \begin{equation} \begin{aligned} \mathbf{y}_k^{(l)} &=\sum_{p=1}^P \gamma_p \mathbf{a}_p^{(\mathrm{rx})} \mathbf{a}_p^{(\mathrm{irs}) \mathrm{T}} \operatorname{diag}_k\left(\mathbf{S}^{(p, l)}\right)\mathbf{b}_p^{(\mathrm{irs})} \mathbf{b}_p^{(\mathrm{tx}) \mathrm{T}} \mathbf{x}_k^{(l)}\\ &+\mathbf{v}_k^{(l)} \in \mathbb{C}^{M_r \times 1} \end{aligned} \end{equation} yk(l)=p=1∑Pγpap(rx)ap(irs)Tdiagk(S(p,l))bp(irs)bp(tx)Txk(l)+vk(l)∈CMr×1
其中
x
k
(
l
)
∈
C
M
r
×
1
\mathbf{x}_k^{(l)} \in \mathbb{C}^{M_r \times 1}
xk(l)∈CMr×1 是所要传输的pilot symbol,
v
k
(
l
)
∈
C
M
r
×
1
\mathbf{v}_k^{(l)} \in \mathbb{C}^{M_r \times 1}
vk(l)∈CMr×1 是AWGN的部分,除此之外,
a
p
(
i
r
s
)
∈
C
N
×
1
\mathbf{a}_p^{(\mathrm{irs})} \in \mathbb{C}^{N \times 1}
ap(irs)∈CN×1 和
b
p
(
i
r
s
)
∈
C
N
×
1
\mathbf{b}_p^{(\mathrm{irs})} \in \mathbb{C}^{N \times 1}
bp(irs)∈CN×1 是IRS array steering vectors。他们的定义如下:
b
p
(
irs
)
=
b
p
v
(
irs
)
⊗
b
p
h
(
irs
)
∈
C
N
h
N
v
×
1
\mathbf{b}_p^{(\text {irs })}=\mathbf{b}_p^{v(\text { irs })} \otimes \mathbf{b}_p^{h(\text { irs })} \in \mathbb{C}^{N_h N_v \times 1}
bp(irs )=bpv( irs )⊗bph( irs )∈CNhNv×1
这是AoA response 的定义,with
N
h
N
v
=
N
N_h N_v=N
NhNv=N,
b
p
v
(
irs
)
\mathbf{b}_p^{v(\text { irs })}
bpv( irs )和
b
p
h
(
irs
)
\mathbf{b}_p^{h(\text { irs })}
bph( irs ) 分别是它的elevation and azimuth steering vectors。
n
n
nth element of
b
p
(
irs
)
\mathbf{b}_p^{(\text {irs })}
bp(irs ) 计算如下:
[
b
p
(
i
r
s
)
]
n
h
+
(
n
v
−
1
)
N
h
=
e
j
π
(
(
n
h
−
1
)
cos
ϕ
p
b
sin
ψ
p
b
+
(
n
v
−
1
)
cos
ϕ
p
b
)
\left[\mathbf{b}_p^{(\mathrm{irs})}\right]_{n_h+\left(n_v-1\right) N_h}=e^{j \pi\left(\left(n_h-1\right) \cos \phi_p^b \sin \psi_p^b+\left(n_v-1\right) \cos \phi_p^b\right)}
[bp(irs)]nh+(nv−1)Nh=ejπ((nh−1)cosϕpbsinψpb+(nv−1)cosϕpb)
with
n
h
=
{
1
,
…
,
N
h
}
n_h=\left\{1, \ldots, N_h\right\}
nh={1,…,Nh},
n
v
=
{
1
,
…
,
N
v
}
n_v=\left\{1, \ldots, N_v\right\}
nv={1,…,Nv}, and
n
=
n
h
+
(
n
v
−
1
)
N
h
n=n_h+(n_v-1)N_h
n=nh+(nv−1)Nh, 其中
ϕ
p
b
\phi_p^b
ϕpb和
ψ
p
b
\psi_p^b
ψpb 分别是AoAs的 azimuth 和 elevation。
同理,AoD response的定义如下:
a
p
(
i
r
s
)
=
a
p
v
(
i
r
s
)
⊗
a
p
h
(
i
r
s
)
∈
C
N
h
N
v
×
1
\mathbf{a}_p^{(\mathrm{irs})}=\mathbf{a}_p^{v(\mathrm{irs})} \otimes \mathbf{a}_p^{h(\mathrm{irs})} \in \mathbb{C}^{N_h N_v \times 1}
ap(irs)=apv(irs)⊗aph(irs)∈CNhNv×1
考虑到
n
n
nth element,进而可被表示为:
[
a
p
(
i
r
s
)
]
n
h
+
(
n
v
−
1
)
N
h
=
e
j
π
(
(
n
h
−
1
)
cos
ϕ
p
a
sin
ψ
p
a
+
(
n
v
−
1
)
cos
ϕ
p
a
)
\left[\mathbf{a}_p^{(\mathrm{irs})}\right]_{n_h+\left(n_v-1\right) N_h}=e^{j \pi\left(\left(n_h-1\right) \cos \phi_p^a \sin \psi_p^a+\left(n_v-1\right) \cos \phi_p^a\right)}
[ap(irs)]nh+(nv−1)Nh=ejπ((nh−1)cosϕpasinψpa+(nv−1)cosϕpa)
其中所涉及的角度为AoD对应的azimuth和elevation。在(5)式中,
S
(
p
,
l
)
∈
C
N
×
K
\mathbf{S}^{(p, l)} \in \mathbb{C}^{N \times K}
S(p,l)∈CN×K 是一个包含了IRS的phase shift的矩阵。可被表示为:
s
k
(
p
,
l
)
=
[
β
1
,
k
(
p
,
l
)
e
j
θ
1
,
k
(
p
)
,
…
,
β
N
,
k
(
p
,
l
)
e
j
θ
N
,
k
(
p
)
]
T
∈
C
N
×
1
\begin{equation} \mathbf{s}_k^{(p, l)}=\left[\beta_{1, k}^{(p, l)} e^{j \theta_{1, k}^{(p)}}, \ldots, \beta_{N, k}^{(p, l)} e^{j \theta_{N, k}^{(p)}}\right]^{\mathrm{T}} \in \mathbb{C}^{N \times 1} \end{equation}
sk(p,l)=[β1,k(p,l)ejθ1,k(p),…,βN,k(p,l)ejθN,k(p)]T∈CN×1
其中,
β
n
,
k
(
p
,
l
)
=
{
0
,
1
}
\beta_{n, k}^{(p, l)}=\{0,1\}
βn,k(p,l)={0,1} 和
θ
n
,
k
\theta_{n, k}
θn,k 分别是reflection coefficient 和designed phase-shift for the
n
n
nth IRS element, for the
k
k
kth time slot。
Conclusion
我按照自己的理解分别介绍了有无IRS参与的两种信道模型,在这个过程中,我仍然想要提醒像我一样的初学者:维度的匹配十分重要!!!这篇总结主要想加深自己的理解,如果能同时帮助到大家理解信道模型的分解,我会更加高兴。
K. B. d. A. Benício, B. Sokal and A. L. F. de Almeida, “Channel Estimation and Performance Evaluation of Multi-IRS Aided MIMO Communication System,” 2021 Workshop on Communication Networks and Power Systems (WCNPS), 2021, pp. 1-6, doi: 10.1109/WCNPS53648.2021.9626243.