参考书目:《MIMO-OFDM无线通信技术及MATLAB实现》第13章
1. 数学模型
N R x × N T x N_{\mathrm{Rx}} \times N_{\mathrm{Tx}} NRx×NTx 单用户 M I M O \mathrm{MIMO} MIMO 系统的信道容量与 N min = min ( N T x , N R x ) N_{\min }=\min \left(N_{\mathrm{Tx}}, N_{\mathrm{Rx}}\right) Nmin=min(NTx,NRx)成正比。事实上, 在高信噪比情况下, MIMO 技术是一种增加容量的必要方法, 它能够提供最多 N min N_{\min } Nmin 个空间自由度。在单用户 MIMO 系统中, 空间复用技术能够支持高数据速率的点对点传输, 同时还能提供空间分集增益(使用多个时频空资源等发送信号副本)。
一个典型的多用户通信环境: 在蜂窝系统中, 一个基站 (BS) 同时服务多个移动台 (MS), 并为他们分配通信资源, 如时间、频率和空间流。假设 BS 和每个 MS 分别有 N B N_{\mathrm{B}} NB 和 N M N_{\mathrm{M}} NM 根天线。因为 K K K 个独立的用户可以形成一个 K ⋅ N M K \cdot N_{\mathrm{M}} K⋅NM的虚天线集, 当与一个具有 N B N_{\mathrm{B}} NB 根天线的 BS 通信时, 这种端到端配置可以被认为是一个 ( K ⋅ N M ) × N B \left(K \cdot N_{\mathrm{M}}\right) \times N_{\mathrm{B}} (K⋅NM)×NB 的 MIMO 下行系统, 或者被认为是一个 N B × ( K ⋅ N M ) N_{\mathrm{B}} \times\left(K \cdot N_{\mathrm{M}}\right) NB×(K⋅NM) 的 MIMO 上行系统。在这个多用户通信系统中, 多天线允许独立的用户在上行链路 (多对一) 同时发送各自的数据流, 或者 BS 在下行链路 (一对多) 向多个用户发送数据流, 其中每个用户独立解码。这归因于多天线增加了自由度, 如同在单用户 MIMO 系统中一样。
在多用户MIMO系统,下行和上行信道分别称为广播信道(Broadcast Channel, BC)和多址接入信道(Multiple Access Channel, MAC)。在MAC中,基站的接收机可以获得K个独立用户的所有数据流,此时多用户MIMO系统相当于一个 ( K ⋅ N M ) × N B \left(K \cdot N_{\mathrm{M}}\right) \times N_{\mathrm{B}} (K⋅NM)×NB 的上行单用户MIMO系统。多用户MIMO系统上行容量与 min ( N B , K ⋅ N M ) \min \left(N_{\mathrm{B}}, K \cdot N_{\mathrm{M}}\right) min(NB,K⋅NM)成正比。
考虑有
K
K
K 个独立用户的多用户 MIMO 系统。假设 BS 和每个 MS 分别配置
N
B
N_{\mathrm{B}}
NB 和
N
M
N_{\mathrm{M}}
NM 根天线。图 13.2 显示了具有
K
K
K 个独立用户的上行信道。令
x
u
∈
C
N
M
×
1
,
u
=
1
,
⋯
,
K
x_u \in \mathbb{C}^{N_{\mathrm{M}} \times 1}, u=1, \cdots, K
xu∈CNM×1,u=1,⋯,K, 表示来自第
u
u
u 个用户的发送信号;
y
M
A
C
∈
C
N
B
×
1
\boldsymbol{y}_{\mathrm{MAC}} \in \mathbb{C}^{N_{\mathrm{B}} \times 1}
yMAC∈CNB×1 表示
B
S
\mathrm{BS}
BS 的接收信号;
H
u
U
L
∈
C
N
B
×
N
M
,
u
=
1
,
⋯
,
K
\boldsymbol{H}_u^{\mathrm{UL}} \in \mathbb{C}^{N_{\mathrm{B}} \times N_{\mathrm{M}}}, u=1, \cdots, K
HuUL∈CNB×NM,u=1,⋯,K, 表示第
u
u
u 个用户和基站之间的信道增益。那么, 接收信号可以表示为
y
M
A
C
=
H
1
U
L
x
1
+
H
2
U
L
x
2
+
⋯
+
H
K
U
L
x
K
+
z
=
[
H
1
U
L
H
2
U
L
⋯
H
K
U
L
]
⏟
H
U
L
[
x
1
⋮
x
K
]
+
z
=
H
U
L
[
x
1
⋮
x
K
]
+
z
\begin{aligned} \boldsymbol{y}_{\mathrm{MAC}} & =\boldsymbol{H}_1^{\mathrm{UL}} \boldsymbol{x}_1+\boldsymbol{H}_2^{\mathrm{UL}} \boldsymbol{x}_2+\cdots+\boldsymbol{H}_K^{\mathrm{UL}} \boldsymbol{x}_K+\boldsymbol{z} \\ & =\underbrace{\left[\boldsymbol{H}_1^{\mathrm{UL}} \boldsymbol{H}_2^{\mathrm{UL}} \cdots \boldsymbol{H}_K^{\mathrm{UL}}\right]}_{\boldsymbol{H}^{\mathrm{UL}}}\left[\begin{array}{c} \boldsymbol{x}_1 \\ \vdots \\ \boldsymbol{x}_K \end{array}\right]+\boldsymbol{z}=\boldsymbol{H}^{\mathrm{UL}}\left[\begin{array}{c} \boldsymbol{x}_1 \\ \vdots \\ \boldsymbol{x}_K \end{array}\right]+\boldsymbol{z} \end{aligned}
yMAC=H1ULx1+H2ULx2+⋯+HKULxK+z=HUL
[H1ULH2UL⋯HKUL]
x1⋮xK
+z=HUL
x1⋮xK
+z
其中,
z
∈
C
N
B
×
1
z \in \mathbb{C}^{N_{\mathrm{B}} \times 1}
z∈CNB×1 为接收机的加性噪声, 将其建模为一个 ZMCSCG 随机向量。
另一方面, 图 13.3 显示了下行信道, 其中
x
∈
C
N
B
×
1
\boldsymbol{x} \in \mathbb{C}^{N_{\mathrm{B}} \times 1}
x∈CNB×1 为
B
S
\mathrm{BS}
BS 的发射信号,
y
u
∈
C
N
M
×
1
\boldsymbol{y}_u \in \mathbb{C}^{N_{\mathrm{M}} \times 1}
yu∈CNM×1 为第
u
u
u 个用户的接收信号,
u
=
1
,
2
,
⋯
,
K
u=1,2, \cdots, K
u=1,2,⋯,K 。令
H
u
D
L
∈
C
N
M
×
N
B
\boldsymbol{H}_u^{\mathrm{DL}} \in \mathbb{C}^{N_{\mathrm{M}} \times N_{\mathrm{B}}}
HuDL∈CNM×NB 表示
B
S
\mathrm{BS}
BS 和第
u
u
u 个用户之间的信道增益。在BC中, 第
u
u
u 个用户的接收信号可以表示为
y
y
=
H
u
D
L
x
+
z
u
,
u
=
1
,
2
,
⋯
,
K
\boldsymbol{y}_{\boldsymbol{y}}=\boldsymbol{H}_u^{\mathrm{DL}} \boldsymbol{x}+\boldsymbol{z}_u, u=1,2, \cdots, K
yy=HuDLx+zu,u=1,2,⋯,K
其中,
z
u
∈
C
N
M
×
1
z_u \in \mathbb{C}^{N_{\mathrm{M}} \times 1}
zu∈CNM×1 为第
u
u
u 个用户的加性 ZMCSCG 噪声。用单个向量表示所有用户的信号, 则整个系统可以表示为
[
y
1
y
2
⋮
y
K
]
⏟
y
B
C
=
[
H
1
D
L
H
2
D
L
⋮
H
K
D
L
]
⏟
H
D
L
x
+
[
z
1
z
2
⋮
z
K
]
⏟
z
\begin{equation} \underbrace{\left[\begin{array}{c} \boldsymbol{y}_1 \\ \boldsymbol{y}_2 \\ \vdots \\ \boldsymbol{y}_K \end{array}\right]}_{\boldsymbol{y}_{\mathrm{BC}}}=\underbrace{\left[\begin{array}{c} \boldsymbol{H}_1^{\mathrm{DL}} \\ \boldsymbol{H}_2^{\mathrm{DL}} \\ \vdots \\ \boldsymbol{H}_K^{\mathrm{DL}} \end{array}\right]}_{\boldsymbol{H}_{\mathrm{DL}}} \boldsymbol{x}+\underbrace{\left[\begin{array}{c} \boldsymbol{z}_1 \\ \boldsymbol{z}_2 \\ \vdots \\ z_K \end{array}\right]}_{\boldsymbol{z}} \end{equation}
yBC
y1y2⋮yK
=HDL
H1DLH2DL⋮HKDL
x+z
z1z2⋮zK
2. BC信道容量
高斯广播信道的容量区域仍然是一个没有解决的问题。本节将讨论在
N
B
=
2
、
N
M
=
1
、
K
=
2
N_{\mathrm{B}}=2 、 N_{\mathrm{M}}=1 、 K=2
NB=2、NM=1、K=2的情况下, 使用脏纸编码 (Dirty Paper Coding, DPC) 可以获得的下行信道容量。在这种情况下, 式 (1) 中的接收信号可以表示为
[
y
1
y
2
]
⏟
D
B
C
=
[
H
1
D
L
H
2
D
L
]
⏟
H
D
L
[
x
1
x
2
]
+
[
z
1
z
2
]
\underbrace{\left[\begin{array}{l} \boldsymbol{y}_1 \\ \boldsymbol{y}_2 \end{array}\right]}_{\boldsymbol{D}_{\mathrm{BC}}}=\underbrace{\left[\begin{array}{l} \boldsymbol{H}_1^{\mathrm{DL}} \\ \boldsymbol{H}_2^{\mathrm{DL}} \end{array}\right]}_{\boldsymbol{H}^{\mathrm{DL}}}\left[\begin{array}{l} x_1 \\ x_2 \end{array}\right]+\left[\begin{array}{l} \boldsymbol{z}_1 \\ z_2 \end{array}\right]
DBC
[y1y2]=HDL
[H1DLH2DL][x1x2]+[z1z2]
其中, H u D L ∈ C 1 × 2 \boldsymbol{H}_u^{\mathrm{DL}} \in \mathbb{C}^{1 \times 2} HuDL∈C1×2 表示基站和第 u u u 个用户之间的信道矩阵, u = 1 , 2 ; x i u=1,2 ; \quad x_i u=1,2;xi 是第 i i i 根发射天线的发送信号, i = 1 , 2 i=1,2 i=1,2 。
如果基站能够获得完全的信道信息, 那么可以对整个信道进行 LQ (下三角
×
\times
×正交矩阵)分解:
H
D
L
=
[
l
11
0
l
21
l
22
]
⏟
L
(
下三角矩阵
)
[
q
1
q
2
]
⏟
Q
(
正交矩阵
)
\boldsymbol{H}^{\mathrm{DL}}=\underbrace{\left[\begin{array}{ll} l_{11} & 0 \\ l_{21} & l_{22} \end{array}\right]}_{\boldsymbol{L(下三角矩阵)}} \underbrace{\left[\begin{array}{l} \boldsymbol{q}_1 \\ \boldsymbol{q}_2 \end{array}\right]}_{\boldsymbol{Q(正交矩阵)}}
HDL=L(下三角矩阵)
[l11l210l22]Q(正交矩阵)
[q1q2]
其中
l
11
=
∥
H
1
D
L
∥
q
1
=
1
l
11
H
1
D
L
l
21
=
q
1
⋅
(
H
2
D
L
)
H
l
22
=
∥
H
2
D
L
−
l
21
q
1
∥
q
2
=
1
l
22
(
H
2
D
L
−
l
21
q
1
)
\begin{gathered} l_{11}=\left\|\boldsymbol{H}_1^{\mathrm{DL}}\right\| \\ \boldsymbol{q}_1=\frac{1}{l_{11}} \boldsymbol{H}_1^{\mathrm{DL}} \\ l_{21}=\boldsymbol{q}_1 \cdot\left(\boldsymbol{H}_2^{\mathrm{DL}}\right)^{\mathrm{H}} \\ l_{22}=\left\|\boldsymbol{H}_2^{\mathrm{DL}}-l_{21} \boldsymbol{q}_1\right\| \\ \boldsymbol{q}_2=\frac{1}{l_{22}}\left(\boldsymbol{H}_2^{\mathrm{DL}}-l_{21} \boldsymbol{q}_1\right) \end{gathered}
l11=
H1DL
q1=l111H1DLl21=q1⋅(H2DL)Hl22=
H2DL−l21q1
q2=l221(H2DL−l21q1)
其中,
{
q
i
}
i
=
1
2
\left\{\boldsymbol{q}_i\right\}_{i=1}^2
{qi}i=12 是正交行向量。给定
{
q
i
}
i
=
1
2
\left\{\boldsymbol{q}_i\right\}_{i=1}^2
{qi}i=12, 可以对发送信号进行预编码:
[
x
1
x
2
]
⏟
x
=
Q
H
[
x
~
1
x
~
2
−
1
l
22
l
21
x
~
1
]
\underbrace{\left[\begin{array}{l} x_1 \\ x_2 \end{array}\right]}_{\boldsymbol{x}}=\boldsymbol{Q}^{\mathrm{H}}\left[\begin{array}{c} \tilde{x}_1 \\ \tilde{x}_2-\frac{1}{l_{22}} l_{21} \tilde{x}_1 \end{array}\right]
x
[x1x2]=QH[x~1x~2−l221l21x~1]
对于上述的预掮码信号, 接收信号可以表示为
y
B
C
=
H
D
L
x
+
z
=
[
l
11
0
l
21
l
22
]
[
q
1
q
2
]
[
q
1
H
q
2
H
]
[
x
~
1
x
~
2
−
1
l
22
l
21
x
~
1
]
+
z
=
[
l
11
0
0
l
22
]
[
x
~
1
x
~
2
]
+
z
(
effective channel
)
=
[
∥
H
1
D
L
∥
0
0
∥
H
1
D
L
−
l
21
q
1
∥
]
[
x
~
1
x
~
2
]
+
z
\begin{aligned} \boldsymbol{y}_{\mathrm{BC}} & =\boldsymbol{H}^{\mathrm{DL}} \boldsymbol{x}+\boldsymbol{z} \\ & =\left[\begin{array}{ll} l_{11} & 0 \\ l_{21} & l_{22} \end{array}\right]\left[\begin{array}{l} \boldsymbol{q}_1 \\ \boldsymbol{q}_2 \end{array}\right]\left[\begin{array}{ll} \boldsymbol{q}_1^{\mathrm{H}} & \boldsymbol{q}_2^{\mathrm{H}} \end{array}\right]\left[\begin{array}{c} \tilde{x}_1 \\ \tilde{x}_2-\frac{1}{l_{22}} l_{21} \tilde{x}_1 \end{array}\right]+\boldsymbol{z} \\ & =\left[\begin{array}{cc} l_{11} & 0 \\ 0 & l_{22} \end{array}\right]\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \end{array}\right]+\boldsymbol{z}\qquad (\text{effective channel})\\ & =\left[\begin{array}{cc} \left\|\boldsymbol{H}_1^{\mathrm{DL}}\right\| & 0 \\ 0 & \left\|\boldsymbol{H}_1^{\mathrm{DL}}-l_{21} \boldsymbol{q}_1\right\| \end{array}\right]\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \end{array}\right]+\boldsymbol{z} \end{aligned}
yBC=HDLx+z=[l11l210l22][q1q2][q1Hq2H][x~1x~2−l221l21x~1]+z=[l1100l22][x~1x~2]+z(effective channel)=
H1DL
00
H1DL−l21q1
[x~1x~2]+z
从上式可以看到,产生了两个无干扰的虚拟信道。
3. BC信道的四种预编码方案
在BC中,数据传输的主要困难是不能在接收机间直接进行协同信号检测,因此需要在BS消除数据间干扰。
3.1 信道反转
单天线用户,令
x
~
u
\tilde{x}_u
x~u表示第
u
u
u个用户的信号,
H
u
DL
∈
C
1
×
K
\boldsymbol{H}_u^{\text{DL}}\in \mathbb{C}^{1\times K}
HuDL∈C1×K表示BS和第
u
u
u 个用户之间的信道矩阵。所有用户的接收信号可以表示为:
[
y
1
y
2
⋮
y
K
]
⏟
y
B
C
=
[
H
1
D
L
H
2
D
L
⋮
H
K
D
L
]
⏟
H
D
L
[
x
~
1
x
~
2
⋮
x
~
K
]
⏟
x
+
[
z
1
z
2
⋮
z
K
]
⏟
z
\underbrace{\left[\begin{array}{c} y_1 \\ y_2 \\ \vdots \\ y_K \end{array}\right]}_{\boldsymbol{y}_{\mathrm{BC}}}=\underbrace{\left[\begin{array}{c} \boldsymbol{H}_1^{\mathrm{DL}} \\ \boldsymbol{H}_2^{\mathrm{DL}} \\ \vdots \\ \boldsymbol{H}_K^{\mathrm{DL}} \end{array}\right]}_{\boldsymbol{H}^{\mathrm{DL}}} \underbrace{\left[\begin{array}{c} \tilde{x}_1 \\ \tilde{x}_2 \\ \vdots \\ \tilde{x}_K \end{array}\right]}_{\boldsymbol{x}}+\underbrace{\left[\begin{array}{c} z_1 \\ z_2 \\ \vdots \\ z_K \end{array}\right]}_{\boldsymbol{z}}
yBC
y1y2⋮yK
=HDL
H1DLH2DL⋮HKDL
x
x~1x~2⋮x~K
+z
z1z2⋮zK
在上式中, 每个用户的接收信号是一个标量。而在用户多天线系统中, 每个用户的接收信号是一个向量。
由于每个用户配置一根天线(无空间复用和天线分集的能力), 不能消除由其他信号引起的干扰。因而, 可以考虑信道反转和规则化信道预编码技术。
在多用户MIMO方案中,信道反转过程与ZF预均衡类似,需要将这里的 H D L \boldsymbol{H}^{\mathrm{DL}} HDL替换掉 W ZF = β H − 1 \mathbf{W}_{\text{ZF}}=\beta \boldsymbol{H}^{-1} WZF=βH−1 中的 H \boldsymbol{H} H。
使用MMSE准则可以减轻噪声增强的影响,将将这里的 H D L \boldsymbol{H}^{\mathrm{DL}} HDL替换掉 W MMSE = β H H ( H H H + σ z 2 σ x 2 I ) − 1 \mathbf{W}_{\text{MMSE}}=\beta \boldsymbol{H}^{H}(\boldsymbol{H}\boldsymbol{H}^{H}+\frac{\sigma_z^2}{\sigma_x^2}\boldsymbol{I})^{-1} WMMSE=βHH(HHH+σx2σz2I)−1 中的 H \boldsymbol{H} H。 β = N Tx Tr ( H − 1 ( H − 1 ) H ) \beta=\sqrt {\frac{N_{\text{Tx}}}{\text{Tr}(\boldsymbol{H}^{-1}(\boldsymbol{H}^{-1})^{H})}} β=Tr(H−1(H−1)H)NTx
3.2 块对角化(Block Diagonalization, BD)
源自同一信号的天线间干扰和其他用户的干扰都能在信道求逆过程中被消除或者减轻,但从目标用户的角度看,噪声增强变得更严重。
BD方法在预编码的过程中只消除来自其他用户的信号干扰。对于同一用户的天线间干扰,可以用各种信号检测方法消除。
将一个信号以 V ~ u zero \tilde{\boldsymbol{V}}_u^{\text{zero}} V~uzero方向发送,第 u u u个用户以外的其他用户根本接收不到该信号。其中, V ~ u zero \tilde{\boldsymbol{V}}_u^{\text{zero}} V~uzero是 H ~ u DL \tilde{\boldsymbol{H}}_u^{\text{DL}} H~uDL的零空间。 H ~ u DL \tilde{\boldsymbol{H}}_u^{\text{DL}} H~uDL是构建的除第 u u u个用户外的其他所有用户的信道矩阵 H ~ u D L = [ ( H 1 D L ) H ⋯ ( H u − 1 D L ) H ( H u + 1 D L ) H ⋯ ( H K D L ) H ] H \tilde{\boldsymbol{H}}_u^{\mathrm{DL}}=\left[\left(\boldsymbol{H}_1^{\mathrm{DL}}\right)^{\mathrm{H}} \cdots\left(\boldsymbol{H}_{u-1}^{\mathrm{DL}}\right)^{\mathrm{H}}\left(\boldsymbol{H}_{u+1}^{\mathrm{DL}}\right)^{\mathrm{H}} \cdots\left(\boldsymbol{H}_K^{\mathrm{DL}}\right)^{\mathrm{H}}\right]^{\mathrm{H}} H~uDL=[(H1DL)H⋯(Hu−1DL)H(Hu+1DL)H⋯(HKDL)H]H。
因此,第
u
u
u个用户的接收信号为
y
u
=
H
u
D
L
V
~
u
zero
x
~
u
⏟
预编码信号
x
u
+
z
u
\boldsymbol{y}_u=\boldsymbol{H}_u^{\mathrm{DL}}\underbrace{\tilde{\boldsymbol{V}}_u^{\text{zero}}\tilde{\boldsymbol{x}}_u}_{预编码信号\boldsymbol{x}_u}+\boldsymbol{z}_u
yu=HuDL预编码信号xu
V~uzerox~u+zu
3.3 脏纸编码(DPC)
DPC 是一种预编码方法,能够将发射机已知的一些干扰删除。更具体地说,对第 k k k个用户的信号进行预编码时,能够删除第一个到第 k − 1 k-1 k−1 个用户信号所造成的干扰。
这里只考虑
N
B
=
3
,
K
=
3
,
N
M
,
u
=
1
,
u
=
1
,
2
,
3
N_{\mathrm{B}}=3, K=3, N_{\mathrm{M}, u}=1, u=1,2,3
NB=3,K=3,NM,u=1,u=1,2,3 的情况。如果给定第
u
u
u 个用户的信号
x
~
u
∈
C
\tilde{x}_u \in \mathbb{C}
x~u∈C, 那么接收信号可以表示为
[
y
1
y
2
y
3
]
=
[
H
1
D
L
H
2
D
L
H
3
D
L
]
⏟
H
D
L
[
x
~
1
x
~
2
x
~
3
]
+
[
z
1
z
2
z
3
]
\left[\begin{array}{l} y_1 \\ y_2 \\ y_3 \end{array}\right]=\underbrace{\left[\begin{array}{l} \boldsymbol{H}_1^{\mathrm{DL}} \\ \boldsymbol{H}_2^{\mathrm{DL}} \\ \boldsymbol{H}_3^{\mathrm{DL}} \end{array}\right]}_{\boldsymbol{H}^{\mathrm{DL}}}\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \\ \tilde{x}_3 \end{array}\right]+\left[\begin{array}{l} z_1 \\ z_2 \\ z_3 \end{array}\right]
y1y2y3
=HDL
H1DLH2DLH3DL
x~1x~2x~3
+
z1z2z3
其中,
H
u
D
L
∈
C
l
×
3
\boldsymbol{H}_u^{\mathrm{DL}} \in \mathbb{C}^{\mathrm{l} \times 3}
HuDL∈Cl×3 为
B
S
\mathrm{BS}
BS 和第
u
u
u 个用户之间的信道增益。对信道增益矩阵
H
D
L
\boldsymbol{H}^{\mathrm{DL}}
HDL 进行
L
Q
L Q
LQ 分解:
H
D
L
=
[
l
11
0
0
l
21
l
22
0
l
31
l
32
l
33
]
⏟
L
[
q
1
q
2
q
3
]
⏟
Q
\boldsymbol{H}^{\mathrm{DL}}=\underbrace{\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \end{array}\right]}_{\boldsymbol{L}} \underbrace{\left[\begin{array}{l} \boldsymbol{q}_1 \\ \boldsymbol{q}_2 \\ \boldsymbol{q}_3 \end{array}\right]}_{\boldsymbol{Q}}
HDL=L
l11l21l310l22l3200l33
Q
q1q2q3
其中,
{
q
i
}
i
=
1
3
∈
C
1
×
3
\left\{\boldsymbol{q}_i\right\}_{i=1}^3 \in \mathbb{C}^{1 \times 3}
{qi}i=13∈C1×3 为正交行向量。令
x
=
[
x
1
x
2
x
3
]
T
\boldsymbol{x}=\left[x_1~x_2~x_3\right]^{\mathrm{T}}
x=[x1 x2 x3]T 表示对已知信号
x
~
=
[
x
~
1
x
~
2
x
~
3
]
T
\tilde{\boldsymbol{x}}=\left[\tilde{x}_1 ~ \tilde{x}_2 ~ \tilde{x}_3\right]^{\mathrm{T}}
x~=[x~1 x~2 x~3]T 的预编码信号。通过发射
Q
H
x
\boldsymbol{Q}^{\mathrm{H}} \boldsymbol{x}
QHx, 在它通过信道后,能够删除
Q
\boldsymbol{Q}
Q 的作用,因此传输后只剩下三角矩阵,接收信号可以表示为
[
y
1
y
2
y
3
]
=
[
H
1
D
L
H
2
D
L
H
3
D
L
]
⏟
H
D
L
Q
H
x
+
[
z
1
z
2
z
3
]
=
[
l
11
0
0
l
21
l
22
0
l
31
l
32
l
33
]
[
x
1
x
2
x
3
]
+
[
z
1
z
2
z
3
]
\begin{equation}\scriptsize \left[\begin{array}{l} y_1 \\ y_2 \\ y_3 \end{array}\right]=\underbrace{\left[\begin{array}{l} \boldsymbol{H}_1^{\mathrm{DL}} \\ \boldsymbol{H}_2^{\mathrm{DL}} \\ \boldsymbol{H}_3^{\mathrm{DL}} \end{array}\right]}_{\boldsymbol{H}^{\mathrm{DL}}} \boldsymbol{Q}^{\mathrm{H}} \boldsymbol{x}+\left[\begin{array}{l} z_1 \\ z_2 \\ z_3 \end{array}\right]=\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \end{array}\right]\left[\begin{array}{l} x_1 \\ x_2 \\ x_3 \end{array}\right]+\left[\begin{array}{l} z_1 \\ z_2 \\ z_3 \end{array}\right] \end{equation}
[y1y2y3]=HDL
H1DLH2DLH3DL
QHx+[z1z2z3]=[l11l21l310l22l3200l33][x1x2x3]+[z1z2z3]
由式 (2) 可知第一个用户的接收信号为
y
1
=
l
11
x
1
+
z
1
y_1=l_{11} x_1+z_1
y1=l11x1+z1
因此, 从第一个用户角度来看, 无干扰的数据传输条件是
x
1
=
x
~
1
\begin{equation} x_1=\tilde{x}_1 \end{equation}
x1=x~1
从式(3)中可以看到预编码信号
x
1
x_1
x1 只由第一个用户的信号
x
~
1
\tilde{x}_1
x~1 构成。从式 (2) 和式 (3) 可知第二个用户的接收信号为
y
2
=
l
21
x
1
+
l
22
x
2
+
z
2
=
l
21
x
~
1
+
l
22
x
2
+
z
2
y_2=l_{21} x_1+l_{22} x_2+z_2=l_{21} \tilde{x}_1+l_{22} x_2+z_2
y2=l21x1+l22x2+z2=l21x~1+l22x2+z2
由上式可知(通过求解方程
l
21
x
~
1
+
l
22
x
2
=
x
~
2
l_{21} \tilde{x}_1+l_{22} x_2=\tilde{x}_2
l21x~1+l22x2=x~2)下式中的预编码在发射端消除了干扰项
l
21
x
1
l_{21} x_1
l21x1 或
l
21
x
~
1
l_{21} \tilde{x}_1
l21x~1 :
x
2
=
x
~
2
−
l
21
l
22
x
1
=
x
~
2
−
l
21
l
22
x
~
1
\begin{equation} x_2=\tilde{x}_2-\frac{l_{21}}{l_{22}} x_1=\tilde{x}_2-\frac{l_{21}}{l_{22}} \tilde{x}_1 \end{equation}
x2=x~2−l22l21x1=x~2−l22l21x~1
由上式可知预编码信号
x
2
x_2
x2 由用户的信号
x
~
1
\tilde{x}_1
x~1 和
x
~
2
\tilde{x}_2
x~2 组成。最后, 第三个用户的接收信号为
y
3
=
l
31
x
1
+
l
32
x
2
+
l
33
x
3
+
z
3
\begin{equation} y_3=l_{31} x_1+l_{32} x_2+l_{33} x_3+z_3 \end{equation}
y3=l31x1+l32x2+l33x3+z3
其中, 预编码信号
x
1
x_1
x1 和
x
2
x_2
x2 由已知的用户信号
x
~
1
\tilde{x}_1
x~1 和
x
~
2
\tilde{x}_2
x~2 组成, 它们分别由式 (3) 和式 (4)给出。从第三个用户的角度来看, 式 (5) 中的预编码信号
x
1
x_1
x1 和
x
2
x_2
x2 是干扰项, 能够通过下面的预编码在发射端将其消除:
x
3
=
x
~
3
−
l
31
l
33
x
1
−
l
32
l
33
x
2
\begin{equation} x_3=\tilde{x}_3-\frac{l_{31}}{l_{33}} x_1-\frac{l_{32}}{l_{33}} x_2 \end{equation}
x3=x~3−l33l31x1−l33l32x2
结合式(3)(4)(6),可以用下面的矩阵形式表示DPC:
[
x
1
x
2
x
3
]
=
[
1
0
0
0
1
0
−
l
31
l
33
−
l
32
l
33
1
]
[
1
0
0
−
l
21
l
22
1
0
0
0
1
]
[
1
0
0
0
1
0
0
0
1
]
[
x
~
1
x
~
2
x
~
3
]
=
[
1
0
0
−
l
21
l
22
1
0
−
l
31
l
33
+
l
32
l
33
l
21
l
22
−
l
32
l
33
1
]
[
x
~
1
x
~
2
x
~
3
]
\scriptsize \begin{aligned} {\left[\begin{array}{l} x_1 \\ x_2 \\ x_3 \end{array}\right] } & =\left[\begin{array}{ccc} 1 & 0 & 0 \\ 0 & 1 & 0 \\ -\frac{l_{31}}{l_{33}} & -\frac{l_{32}}{l_{33}} & 1 \end{array}\right]\left[\begin{array}{ccc} 1 & 0 & 0 \\ -\frac{l_{21}}{l_{22}} & 1 & 0 \\ 0 & 0 & 1 \end{array}\right]\left[\begin{array}{lll} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{array}\right]\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \\ \tilde{x}_3 \end{array}\right] \\ & =\left[\begin{array}{ccc} 1 & 0 & 0 \\ -\frac{l_{21}}{l_{22}} & 1 & 0 \\ -\frac{l_{31}}{l_{33}}+\frac{l_{32}}{l_{33}} \frac{l_{21}}{l_{22}} & -\frac{l_{32}}{l_{33}} & 1 \end{array}\right]\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \\ \tilde{x}_3 \end{array}\right] \end{aligned}
[x1x2x3]=[10−l33l3101−l33l32001][1−l22l210010001][100010001][x~1x~2x~3]=[1−l22l21−l33l31+l33l32l22l2101−l33l32001][x~1x~2x~3]
使用上面的预编码矩阵, 式 (2) 能够被重写为
[
y
1
y
2
y
3
]
=
[
l
11
0
0
l
21
l
22
0
l
31
l
32
l
33
]
[
x
1
x
2
x
3
]
+
[
z
1
z
2
z
3
]
=
[
l
11
0
0
l
21
l
22
0
l
31
l
32
l
33
]
[
1
0
0
−
l
21
l
22
1
0
]
[
x
~
1
x
~
2
x
~
3
]
+
[
z
1
z
33
z
2
z
3
]
=
[
l
11
0
0
0
l
22
0
0
0
l
33
]
[
l
~
12
x
~
2
x
~
2
x
~
3
]
+
[
z
1
z
2
z
3
]
\begin{aligned} {\left[\begin{array}{l} y_1 \\ y_2 \\ y_3 \end{array}\right] } & =\left[\begin{array}{lll} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \end{array}\right]\left[\begin{array}{l} x_1 \\ x_2 \\ x_3 \end{array}\right]+\left[\begin{array}{l} z_1 \\ z_2 \\ z_3 \end{array}\right] \\ & =\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \end{array}\right]\left[\begin{array}{ccc} 1 & 0 & 0 \\ -\frac{l_{21}}{l_{22}} & 1 & 0 \end{array}\right]\left[\begin{array}{l} \tilde{x}_1 \\ \tilde{x}_2 \\ \tilde{x}_3 \end{array}\right]+\left[\begin{array}{l} z_1 \\ z_{33} \\ z_2 \\ z_3 \end{array}\right] \\ & =\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ 0 & l_{22} & 0 \\ 0 & 0 & l_{33} \end{array}\right]\left[\begin{array}{l} \tilde{l}_{12} \\ \tilde{x}_2 \\ \tilde{x}_2 \\ \tilde{x}_3 \end{array}\right]+\left[\begin{array}{l} z_1 \\ z_2 \\ z_3 \end{array}\right] \end{aligned}
y1y2y3
=
l11l21l310l22l3200l33
x1x2x3
+
z1z2z3
=
l11l21l310l22l3200l33
[1−l22l210100]
x~1x~2x~3
+
z1z33z2z3
=
l11000l22000l33
l~12x~2x~2x~3
+
z1z2z3
由上式可知,可以对每个用户进行无干扰检测;而且 DPC 的预编码矩阵是一个被缩放的下三角矩阵的逆矩阵,其中下三角矩阵由信道增益矩阵得到。也就是说,
[
1
0
0
−
l
21
l
22
1
0
−
l
31
l
33
+
l
32
l
33
l
21
l
22
−
l
32
l
33
1
]
=
[
l
11
0
0
l
21
l
22
0
l
31
l
32
l
33
]
−
1
[
l
11
0
0
0
l
22
0
0
0
l
33
]
\small \left[\begin{array}{ccc} 1 & 0 & 0 \\ -\frac{l_{21}}{l_{22}} & 1 & 0 \\ -\frac{l_{31}}{l_{33}}+\frac{l_{32}}{l_{33}} \frac{l_{21}}{l_{22}} & -\frac{l_{32}}{l_{33}} & 1 \end{array}\right]=\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ l_{21} & l_{22} & 0 \\ l_{31} & l_{32} & l_{33} \end{array}\right]^{-1}\left[\begin{array}{ccc} l_{11} & 0 & 0 \\ 0 & l_{22} & 0 \\ 0 & 0 & l_{33} \end{array}\right]
1−l22l21−l33l31+l33l32l22l2101−l33l32001
=[l11l21l310l22l3200l33]−1[l11000l22000l33]
3.4 TH预编码
在发射端的DPC与在接收端的判决反馈均衡 (Decision Feedback Equalization,DFE)非常相似。实际上,已经证明了将对称模运算与 DPC结合同Tomlinson-Harashima(TH)预编码是等价的。DFE存在差错传播的问题,最初发明 TH预编码是为了减少DFE中的峰值功率或平均功率。在 DFE 中进行 TH 预编码的原始想法是消除发射机的后体(Post-cursor)ISI,其中无差错地已知之前发送的符号。实际上,TH 预编码需要完全的信道脉冲响应知识,对于时不变或时间慢变的信道,只有通过接收机的反馈才能获得。为了方便阐述这个想法,考虑一维情