CN103117790A - 分层代数空时码在下行多用户多输入多输出系统中的仿射投影解码方法 - Google Patents

Abstract

本发明公开了分层代数空时码在下行多用户多输入多输出系统中的仿射投影解码方法：一：分别将接收到的信号Y，信道矩阵H，噪声矩阵W，写成列向量的形式；二：把信道矩阵作不同形式的分解；三：对H₁₁，H₁₂作极分解；四：将y₁沿着V₁投影到U₁上的仿射投影；五：将P₁ ^H乘以

Description

分层代数空时码在下行多用户多输入多输出系统中的仿射投影解码方法

技术领域

本发明属于无线通信技术领域，特别是多天线技术领域，具体是借助于仿射投影法把适用于点对点通信的分层代数空时码应用于下行多用户多输入多输出系统中。在这一应用中，基站并不需要任何有关信道的信息。

背景技术

在窄带无线通信系统中多径衰落是非常严重的问题，信号幅度的衰落会导致接收方误码率（BER）性能的降低。在无线移动通信中广泛使用了分集技术来减少多径衰落的影响，并且在不增加发射功率或牺牲通信带宽的情况下提高传输的可靠性。空间复用就是在发送端和接收端都配置多根天线，充分利用空间传播中的多径分量，在同一频带上使用多个数据通道来发射信号，从而使得容量随天线数量的增加而线性增加。

对角空时码是一种特殊的空时码，它的每一个码字是一个对角矩阵。对角空时码用于多天线系统时，每次发射信号时只有一根天线在工作，其他的天线都被闲置。分层代数空时码是一种每一层上都用一个对角空时码来填补的空时码。而利用分层空时码可以提高信息的传送率，充分利用多天线的优势。Mohamed Oussama Damen和Norman C.Beaulieu在IEEE Trans.Inform.Theory，vol.49,no.4,pp.1059-1063,April.2003发表的“On Two High-Rate AlgebraicSpace-Time Codes”中给出了如下空时码： $\{X=\frac{1}{\sqrt{2}}\left[\begin{array}{cc}{x}_1+x_2\mathrm{\θ}& x_3+x_4\mathrm{\θ}\\ \mathrm{\γ}(x_3-x_4\mathrm{\θ})& x_1-x_2\mathrm{\θ}\end{array}\right];x_1,x_2,x_3,x_4\∈Z[j\left]\right\},$ 其中，γ与θ都取e^jφ。这种码的构造是建立在代数数论中数域扩张理论的基础之上，他们的发现使得空时码的构造与代数数论联系在一起。当x₁,x₂,x₃,x₄可以取Ζ[j]中任意值时，这种空时码的非零最小行列式的值可能要趋向于零。那么是否存在一个空时码，使得非零的最小行列式不趋于零呢？回答是肯定的。而后，Belfiore J C和Rekaya G,Viterbo E在IEEE Trans.Inform.Theory，vol.51,no.4,pp.1432-1436,April.2005发表的“The Golden Code:A2×2Full-Rate Space-TimeCode WithNonvanishing Determinants”提出了黄金空时码，它是上述空时码规范化后得到的。黄金码使全速率，满秩且行列式下有非零界的完备空时码成为可能。Huiyong Liao,Haiquan Wang,和Xiang-Gen Xia在IEEE Trans.Inform.Theory，vol.55,no.2,pp.569-583,February.2009发表“Some Designs and Normalized Diversity Product Upper Bounds for Lattice-Based Diagonal andFull-Rate Space-Time Block Codes”提出了几乎是最优的2×2的完全空时码。

发明内容

以上所述的空时码都应用于点对点的无线通信，而本发明就是要把以上空时码应用于多用户多天线的下行信道中，提供一种分层代数空时码在下行多用户多输入多输出系统中的仿射投影解码方法。

假定基站配有M根发射天线、K个用户，且每个用户有N根接收天线的MIMO系统。假设发送时间间隔为T。要求N≥M和T=M并且K≤M。发送端发送给第i个用户的信息序列S_i=(s_i1，s_i2，...S_iM)^t，1≤i≤K其中的每个元素都是取自某一特定星座。假设K＜M，当i＞K时S_i=0。将S₁,S₂...S_M分别经过变换矩阵G，得到(x₁₁,x₁₂...x_1M)^t=G(s₁₁,s₁₂,...s_1M)^t，(x₂₁,x₂₂...x_2M)^t=G(s₂₁,s₂₂,...s_2M)^t，...(x_M1,x_M2...x_MM)^t=G(s_M1,s_M2,...s_MM)^t。把x₁₁,x₁₂...x_1M按顺序排列到主对角线上,x₂₁,x₂₂...x_2M按顺序排列到次对角线上，以此类推，x_M1,x_M2...x_MM按顺序排列到第M次对角线上。这样就得到了一个M×M维的发射信号矩阵第i层上的符号包含了发送给第i个用户的信息。第i个用户接收到的信息模型为 $Y=\sqrt{\frac{p}{M}}\mathrm{HX}+W.$

下面是介绍采用本发明技术方案对第1个用户的解码方法，其它用户的解码方法可以用矩阵的列置换及以下相同的方法得到：

第一步：分别将接收到的信号Y，信道矩阵H,噪声矩阵W，写成列向量的形式，

Y=[y₁ y₂ ... y_T]，H=[h₁ h₂ ... h_M]，W=[w₁ w₂ ... w_T].

则

第二步:把信道矩阵作一下不同形式的分解，H₁₁=[h₁],H₁₂=[h₂ h₃ ... h_M]；

H₂₁=[h₂]，H₂₂=[h₁ h₃ ... h_M]；...；H_M1=[h_M]，H_M2=[h₁ h₂ ... h_M-1]。

第三步：对H₁₁，H₁₂作极分解得H₁₁=P₁C₁₁,H₁₂=Q₁C₁₂，其中P₁（Q₁）是N×1（N×(M-1)）的酉矩阵，C₁₁（C₁₁）是1×1（(M-1)×(M-1)）的正定矩阵。令U₁是由P₁的列向量组成的子空间，V₁是由Q₁的列向量组成的子空间。

第四步：将y₁沿着V₁投影到U₁上的仿射投影为：

化简得

第五步：将P₁ ^H乘以

得到 ${\tilde{y}}_1=\sqrt{\frac{P}{M}}{C}_{11}{x}_{11}+{\tilde{w}}_1,$ ${\tilde{w}}_1={(1-{P_1}^H{Q}_1{Q_1}^H{P}_1)}^{-1}{P_1}^H(I_N-Q_1{Q_1}^H)w_1$

重复第三步、第四步、第五步，可以得到 ${\tilde{y}}_2=\sqrt{\frac{P}{M}}{C}_{21}{x}_{12}+{\tilde{w}}_2,$ ${\tilde{y}}_M=\sqrt{\frac{P}{M}}{C}_{M1}{x}_{1M}+{\tilde{w}}_M$

其中， ${\tilde{w}}_2={(1-{P_2}^H{Q_2{Q}_2}^H{P}_2)}^{-1}{P_2}^H(I_N-Q_2{Q_2}^H)w_2,$

${\tilde{w}}_M={(1-{P_M}^H{Q}_M{Q_M}^H{P}_M)}^{-1}{P_M}^H(I_N-Q_M{Q_M}^H)w_M$

第六步：记

x=[x₁₁ x₁₂ ... x_1M]^t,

通过最大似然解码，解出第1个用户接收到的信号信息为

其中

下面介绍一下有关仿射投影法：

已知U和V是复矢量空间或实矢量空间内积<，>的两个子空间，并且假设

令{u₁,u₂,…u_n}和{v₁,v₂,…v_m}分别是U和V的基向量。定义

$G_{11}\stackrel{\mathrm{\&Delta;}}{=}{(<u_i,u_j>)}_{n\&times;n},$ $G_{12}\stackrel{\mathrm{\&Delta;}}{=}{(<u_i,v_j>)}_{n\&times;m},$

$G_{21}\stackrel{\mathrm{\&Delta;}}{=}{(<v_i,u_j>)}_{m\&times;n},$ $G_{22}\stackrel{\mathrm{\&Delta;}}{=}{(<v_i,v_j>)}_{m\&times;m}$

并且

$G=\left(\begin{array}{cc}{G}_{11}& G_{12}\\ G_{21}& G_{22}\end{array}\right)$

因为u_j，j=1,2,…,n，v_j，j=1,2,…,m分别是相互独立的，因此G₁₁和G₂₂是可逆的。矩阵G也是可逆的。因为W是U和V的直和，那么{u₁,u₂,…u_n,v₁,v₂,…v_m}就是W的基向量。

设是x空间W中的任意向量，存在特殊的分解x=x₁+x₂，其中x₁∈U，x₂∈V。x₁可以用{u₁,u₂,…u_n}的线性组合表示，x₂可以用{v₁,v₂,…v_m}的线性组合表示。接下来，就先计算线性组合的系数。

假设 $x_1=\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}{x}_i{u}_i,$ $x_2=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{x}_{n+i}{v}_i,$ 其中x_j，j=1,2,…,n+m，是系数。令

a=(x₁,x₂,…,x_n)^t，b=(x_n+1,x_n+2,…,x_n+m)^t

那么就容易得到

$\left(\begin{array}{c}a\\ b\end{array}\right)=G^{-1}\left(\begin{array}{c}<u_1,x>\\ ...\\ <u_n,x>\\ <v_1,x>\\ ...\\ <v_m,x>\end{array}\right)$

为了简化，定义

<u,x>=(<u₁,x>,…,<u_n,x>)^t，<v,x>=(<v₁,x>,…,<v_m,x>)^t。那么根据上面一个等式，a和b可计算如下：

$a=G_{11}^{-1}(I_n+G_{12}{G}_{\mathrm{22,1}}^{-1}{G}_{21}{G}_{11}^{-1})<u,x>$

$-G_{11}^{-1}{G}_{12}{G}_{\mathrm{22,1}}^{-1}<v,x>$

$b=-G_{22}^{-1}{G}_{21}{G}_{\mathrm{11,2}}^{-1}<u,x>+G_{\mathrm{22,1}}^{-1}<v,x>$

在上式中 $G_{\mathrm{22,1}}=G_{22}-G_{21}{G}_{11}^{-1}{G}_{12},$ $G_{\mathrm{11,2}}=G_{11}-G_{12}{G}_{22}^{-1}{G}_{21},$ 因此，x₁和x₂可表示为x₁=(u₁,…,u_n)a，x₂=(v₁,…,v_n)b。根据投影的定义，把向量x₁(或x₂)叫做x沿着空间V(或U)在U(或V)空间上的仿射投影。他们的表达式是这样的：

$P_{\mathrm{uv}}\left(x\right)\stackrel{\mathrm{\&Delta;}}{=}{x}_1=(u_1,\&CenterDot;\&CenterDot;\&CenterDot;,u_n)a,$ $P_{\mathrm{vu}}\left(x\right)\stackrel{\mathrm{\&Delta;}}{=}{x}_2=(v_1,\&CenterDot;\&CenterDot;\&CenterDot;,v_n)b$

具体如图4所示。

本发明的优点是：在发送端不知道信道状态信息的情况下，用仿射投影法可以消除下行多用户MIMO系统中多用户间的干扰信息，从而恢复出发送端发送信号的信息。而所有采用预编码技术来消除信道干扰都要求知道信道的状态信息。

附图说明

图1是基于仿射投影法的分层代数空时码的误比特率分析。

图2是基于仿射投影法的黄金码的误比特率分析。

图3是基于仿射投影法的Alamouti码的误比特率分析。

图4是仿射投影图。

具体实施方式

下面对本发明实施例作详细说明。

实施例1：

设定系统配有4根发射天线，4个用户每个用户有4根接收天线。星座用4-QAM,总能量为4。传送率为8比特pcu。发射天线系统采用四阶的分层空时码。

步骤1：根据以上的解码方案。首先将发送接收方程写成

$\left[\begin{array}{cccc}{y}_1& y_2& y_3& y_4\end{array}\right]=\sqrt{\frac{p}{M}}\left[\begin{array}{cccc}{h}_1& h_2& h_3& h_4\end{array}\right]\left[\begin{array}{cccc}{x}_{11}& x_{21}& x_{31}& x_{41}\\ x_{44}& x_{12}& x_{22}& x_{32}\\ x_{34}& x_{43}& x_{13}& x_{23}\\ x_{24}& x_{33}& x_{42}& x_{14}\end{array}\right]+\left[\begin{array}{cccc}{w}_1& w_2& w_3& w_4\end{array}\right]$

则， $y_1=\sqrt{\frac{p}{M}}(h_1{x}_{11}+h_2{x}_{44}+h_3{x}_{34}+h_4{x}_{24})+w_1,$ $y_2=\sqrt{\frac{p}{M}}(h_1{x}_{21}+h_2{x}_{12}+h_3{x}_{43}+h_4{x}_{33})+w_2,$

$y_3=\sqrt{\frac{p}{M}}(h_1{x}_{31}+h_2{x}_{22}+h_3{x}_{13}+h_4{x}_{42})+w_3,$ $y_4=\sqrt{\frac{p}{M}}(h_1{x}_{41}+h_2{x}_{32}+h_3{x}_{23}+h_4{x}_{14})+w_4$

步骤2：把信道矩阵分为两部分H₁₁=[h₁]，H₁₂=[h₂,h₃,h₄]，同理H₂₁=[h₂]，H₂₂=[h₁,h₃,h₄]；

H₃₁=[h₃]，H₃₂=[h₁,h₂,h₄]；H₄₁=[h₄]，H₃₂=[h₁,h₂,h₃]。

步骤3：作H₁₁，H₁₂的极分解得H₁₁=P₁C₁₁,H₁₂=Q₁C₁₂。其中P₁（Q₁)是4×1（4×3）的酉矩阵，C₁₁（C₁₂）是1×1（3×3）的正定矩阵。设U₁是由P₁的列向量组成的子空间,V₁是由Q₁的列向量组成的子空间。

步骤4：将y₁沿着V₁投影到U₁上的仿射投影为：

步骤5：令

${\tilde{w}}_1={(1-{P_1}^H{Q}_1{Q_1}^H{P}_1)}^{-1}{P_1}^H(I_4-Q_1{Q_1}^H)w_1$

重复步骤3、4、5得

${\tilde{y}}_2=\sqrt{\frac{P}{M}}{C}_{21}{x}_{12}+{\tilde{w}}_2,$ ${\tilde{y}}_3=\sqrt{\frac{P}{M}}{C}_{31}{x}_{13}+{\tilde{w}}_3,$ ${\tilde{y}}_4=\sqrt{\frac{P}{M}}{C}_{41}{x}_{14}+{\tilde{w}}_4,$ 其中，

${\tilde{w}}_2={(1-{P_2}^H{Q_2{Q}_2}^H{P}_2)}^{-1}{P_2}^H(I_4-Q_2{Q_2}^H)w_2,$ ${\tilde{w}}_3={(1-{P_3}^H{Q_3{Q}_3}^H{P}_3)}^{-1}{P_3}^H(I_4-Q_3{Q_3}^H)w_3,$

${\tilde{w}}_4={(1-{P_4}^H{Q}_4{Q_4}^H{P}_4)}^{-1}{P_4}^H(I_4-Q_4{Q_4}^H)w_4$

步骤6：记

$\tilde{Y}=\left[\begin{array}{c}{\tilde{y}}_1\\ {\tilde{y}}_2\\ {\tilde{y}}_3\\ {\tilde{y}}_4\end{array}\right]=\sqrt{\frac{P}{M}}\left[\begin{array}{c}{C}_{11}{x}_{11}\\ C_{21}{x}_{12}\\ C_{31}{x}_{13}\\ C_{41}{x}_{14}\end{array}\right]+\left[\begin{array}{c}{\tilde{w}}_1\\ {\tilde{w}}_2\\ .\\ .\\ .\\ {\tilde{w}}_T\end{array}\right]=\sqrt{\frac{P}{M}}\left[\begin{array}{cccc}{C}_{11}& & & \\ & C_{21}& & \\ & & C_{31}& \\ & & & C_{41}\end{array}\right]\left[\begin{array}{c}{x}_{11}\\ x_{12}\\ x_{13}\\ x_{14}\end{array}\right]+\left[\begin{array}{c}{\tilde{w}}_1\\ {\tilde{w}}_2\\ {\tilde{w}}_3\\ {\tilde{w}}_4\end{array}\right],$

x=[x₁₁ x₁₂ x₁₃ x₁₄]^t

通过最大似然解码，解出第1个用户接收到的信号信息。

$\hat{x}=\arg \min {\left((\tilde{Y}-\sqrt{\frac{P}{M}}\left[\begin{array}{cccc}{C}_{11}& & & \\ & C_{21}& & \\ & & C_{31}& \\ & & & C_{41}\end{array}\right]x\right)}^H{\mathrm{\&Sigma;}}^{-1}(\tilde{Y}-\sqrt{\frac{P}{M}}\left[\begin{array}{cccc}{C}_{11}& & & \\ & C_{21}& & \\ & & C_{31}& \\ & & & C_{41}\end{array}\right]x))$

其中， $\mathrm{\&Sigma;}=\left[\begin{array}{cccc}{(1-{P_1}^H{Q}_1{Q_1}^H{P}_1)}^{-1}& & & \\ & {(1-{P_2}^H{Q}_2{Q_2}^H{P}_2)}^{-1}& & \\ & & {(1-{P_3}^H{Q}_3{Q_3}^H{P}_3)}^{-1}& \\ & & & {(1-{P_4}^H{Q}_4{Q_4}^H{P}_4)}^{-1}\end{array}\right].$

得到图1的仿真结果。

实施例2：

设定系统有2根发射天线，2个用户每个用户有4根接收天线，发送时间间隔T=2。星座用4-QAM,总能量是4则传送率为8比特pcu。系统模型为

采用黄金码作为发送信号矩阵X，其中x₁，x₂是发送给第1个用户的信息，x₃，x₄是发送给第2个用户的信息，则发送接收方程写成

$\left[\begin{array}{cc}{y}_1& y_2\end{array}\right]=\sqrt{\frac{p}{M}}\left[\begin{array}{cc}{h}_1& h_2\end{array}\right]\frac{1}{\sqrt{5}}\left[\begin{array}{cc}\mathrm{\&alpha;}(x_1+{\mathrm{\&theta;x}}_2)& \mathrm{\&alpha;}(x_3+{\mathrm{\&theta;x}}_4)\\ \mathrm{i\&delta;}\left(\mathrm{\&alpha;}\right)(x_3+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_4)& \mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2)\end{array}\right]$ 则，

$+\left[\begin{array}{cc}{w}_1& w_2\end{array}\right]$

$y_1=\sqrt{\frac{p}{5M}}(h_1\mathrm{\&alpha;}(x_1+{\mathrm{\&theta;x}}_2)+h_2\mathrm{i\&delta;}\left(\mathrm{\&alpha;}\right)(x_3+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_4))+w_1$

$y_2=\sqrt{\frac{p}{5M}}(h_1\mathrm{\&alpha;}(x_3+{\mathrm{\&theta;x}}_4)+h_2\mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2))+w_2.$

根据以上的解码方案：

步骤1：将h₁h₂进行极分解得h₁=PC₁，h₂=QC₂

步骤2：将y₁沿着h₂在h₁（y₂沿着h₁在h₂）上做仿射投影，则

将y₁，y₂分别带入上式，化简得

步骤3：将P^H，Q^H分别乘以

得

${\tilde{y}}_1=\sqrt{\frac{p}{M}}\mathrm{\&alpha;}{C}_1(x_1+{\mathrm{\&theta;x}}_2)+{\tilde{w}}_1,$ ${\tilde{y}}_2=\sqrt{\frac{p}{M}}\mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)C_2(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2)+{\tilde{w}}_2$

其中 ${\tilde{w}}_1={(1-P^{H}Q{Q}^{H}P)}^{-1}{P}^H(I_2-Q{Q}^H)w_1$

${\tilde{w}}_2={(1-Q^{H}P{P}^{H}Q)}^{-1}{Q}^H(I_2-P{P}^H)w_2$

记 $\tilde{Y}=\left[\begin{array}{c}{\tilde{y}}_1\\ {\tilde{y}}_2\end{array}\right]=\sqrt{\frac{p}{5M}}\left[\begin{array}{cc}{C}_1& \\ & C_2\end{array}\right]\left[\begin{array}{c}\mathrm{\&alpha;}(x_1+{\mathrm{\&theta;x}}_2)\\ \mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2)\end{array}\right]+\left[\begin{array}{c}{\tilde{w}}_1\\ {\tilde{w}}_2\end{array}\right]$

步骤4：通过最大似然解码，得到第1个用户的信息

$\hat{x}=\arg \min {\left((\tilde{Y}-\sqrt{\frac{p}{5M}}\left[\begin{array}{cc}{C}_1& \\ & C_2\end{array}\right]\left[\begin{array}{c}\mathrm{\&alpha;}(x_1+\mathrm{\&theta;}{x}_2)\\ \mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2)\end{array}\right]\right)}^H$

$\left[\begin{array}{cc}(1-P^{H}Q{Q}^{H}P)& \\ & (1-Q^{H}P{P}^{H}Q)\end{array}\right](\tilde{Y}-\sqrt{\frac{p}{5M}}\left[\begin{array}{cc}{C}_1& \\ & C_2\end{array}\right]\left[\begin{array}{c}\mathrm{\&alpha;}(x_1+{\mathrm{\&theta;x}}_2)\\ \mathrm{\&delta;}\left(\mathrm{\&alpha;}\right)(x_1+\mathrm{\&delta;}\left(\mathrm{\&theta;}\right)x_2)\end{array}\right]))$

得到图2的仿真结果。

实施例3：

设定系统有2根发射天线，2个用户每个用户有2根接收天线，发送时间间隔T=2。星座用4-QAM,总能量是4，则传送率为4比特pcu。系统模型为

采用Alamouti码作为发送信号矩阵X，其中x₁，是发送给第1个用户的信息，x₂是发送给第2个用户的信息，则发送接收方程写成 $\left[\begin{array}{cc}{y}_1& y_2\end{array}\right]=\sqrt{\frac{p}{M}}\left[\begin{array}{cc}{h}_1& h_2\end{array}\right]\left[\begin{array}{cc}{x}_1& x_2\\ {-x}_2^{*}& x_1^{*}\end{array}\right]+\left[\begin{array}{cc}{w}_1& w_2\end{array}\right]$

则， $y_1=\sqrt{\frac{p}{M}}(h_1{x}_1+h_2(-x_2^{*}))+w_1,$ $y_2=\sqrt{\frac{p}{M}}(h_1{x}_2+h_2{x}_1^{*})+w_2.$

根据以上的解码方案：

步骤1：将h₁h₂进行极分解得h₁=PC₁，h₂=QC₂

步骤2：将y₁沿着h₂在h₁（y₂沿着h₁在h₂）上做仿射投影，

将y₁，y₂分别带入上式，化简得

步骤3：将PH，QH分别乘以

得

${\tilde{y}}_1=\sqrt{\frac{p}{M}}P{C}_1{x}_1+{\tilde{w}}_1,$ 其中 ${\tilde{w}}_1={(1-P^{H}Q{Q}^{H}P)}^{-1}{P}^H\left(I_{2}Q{Q}^H\right)w_1$

${\tilde{y}}_2=\sqrt{\frac{p}{M}}{C}_2{x}_1^{*}+{\tilde{w}}_2,$ 其中 ${\tilde{w}}_2={(1-Q^{H}P{P}^{H}Q)}^{-1}{Q}^H(I_2-P{P}^H)w_2$

记 $\tilde{Y}=\left[\begin{array}{c}{\tilde{y}}_1\\ {\tilde{y}}_2\end{array}\right]=\sqrt{\frac{p}{M}}\left[\begin{array}{c}{C}_1\\ C_2^{*}\end{array}\right]x_1+\left[\begin{array}{c}{\tilde{w}}_1\\ {\tilde{w}}_2^{*}\end{array}\right]$

步骤4：通过最大似然解码，得到第1个用户的信息

${\hat{x}}_1=\arg \min {\left((\tilde{Y}-\sqrt{\frac{p}{M}}\left[\begin{array}{c}{C}_1\\ C_2^{*}\end{array}\right]x_1\right)}^H\left[\begin{array}{cc}(1-P^{H}Q{Q}^{H}P)& \\ & (1-Q^{H}P{P}^{H}Q)\end{array}\right](\tilde{Y}-\sqrt{\frac{p}{M}})\left[\begin{array}{c}{C}_1\\ C_2^{*}\end{array}\right]x_1\left)\right)$

得到图3的仿真结果。

当然，本发明还可有其他多种实施例，在不背离发明精神及其实质的情况下，本领域的技术人员可根据本发明做出各种相应的改变和变形，但这些相应的改变和变形都落入本发明的保护范围。

Claims (1)

1.分层代数空时码在下行多用户多输入多输出系统中的仿射投影解码方法，其特征是按如下步骤：

第一步：分别将接收到的信号Y，信道矩阵H，噪声矩阵W，写成列向量的形式，Y=[y₁ y₂ ... y_T]，H=[h₁ h₂ ... h_M]，W=[w₁ w₂ ... w_T].

则

第二步:把信道矩阵作不同形式的分解，H₁₁=[h₁],H₁₂=[h₂ h₃ ... h_M]；H₂₁=[h₂]，

H₂₂=[h₁ h₃ ... h_M]；...；H_M1=[h_M]，H_M2=[h₁ h₂ ... h_M-1]

第三步：对H₁₁，H₁₂作极分解得H₁₁=P₁C₁₁,H₁₂=Q₁C₁₂，其中P₁（Q₁)是N×1（N×(M-1)）的酉矩阵，C₁₁（C₁₂）是1×1（(M-1)×(M-1)）的正定矩阵；令U₁是由P₁的列向量组成的子空间，V₁是由Q₁的列向量组成的子空间；

第四步：将y₁沿着V₁投影到U₁上的仿射投影为：

化简得

第五步：将

乘以

得到 ${\tilde{y}}_1=\sqrt{\frac{P}{M}}{C}_{11}{x}_{11}+{\tilde{w}}_1,$ ${\tilde{w}}_1={(1-{P_1}^H{Q}_1{Q_1}^H{P}_1)}^{-1}{P_1}^H(I_N-Q_1{Q_1}^H)w_1$

重复第三步、第四步、第五步，得到 ${\tilde{y}}_2=\sqrt{\frac{P}{M}}{C}_{21}{x}_{12}+{\tilde{w}}_2,$ ${\tilde{y}}_M=\sqrt{\frac{P}{M}}{C}_{M1}{x}_{1M}+{\tilde{w}}_M$

其中， ${\tilde{w}}_2{={(1-{P_2}^H{Q}_2{Q_2}^H{P}_2)}^{-1}{P}_2}^H(I_N-Q_2{Q_2}^H)w_2,$

${\tilde{w}}_M={(1-{P_M}^H{Q}_M{Q_M}^H{P}_M)}^{-1}{P_M}^H(I_N-Q_M{Q_M}^H)w_M$

第六步：记

x=[x₁₁ x₁₂ ... x_1M]^t,

通过最大似然解码，解出第1个用户接收到的信号信息为

其中

Images