发明名称 用于在多小区多用户多天线系统中的解码方法
摘要 本发明涉及一种用于在多小区多用户多天线系统中的解码方法。本发明包括两大步骤,第一步是设计训练序列头,设<img file="2014106135712100004dest_path_image002.GIF" wi="24" he="24" />的前<img file="dest_path_image004.GIF" wi="13" he="16" />列是训练序列头,<img file="53509dest_path_image002.GIF" wi="24" he="24" />的第一列中的元素均为<img file="dest_path_image006.GIF" wi="8" he="16" />,<img file="908332dest_path_image002.GIF" wi="24" he="24" />的第2列至第<img file="688070dest_path_image004.GIF" wi="13" he="16" />列记为<img file="dest_path_image008.GIF" wi="22" he="24" />,其中<img file="270230dest_path_image008.GIF" wi="22" he="24" />中的元素可为<img file="22285dest_path_image006.GIF" wi="8" he="16" />或<img file="dest_path_image010.GIF" wi="21" he="16" />,且<img file="844747dest_path_image008.GIF" wi="22" he="24" />中各行向量均不相同。第二步是基于系统模型,用接收信号矩阵<img file="dest_path_image012.GIF" wi="18" he="18" />估计出发送信号矩阵<img file="dest_path_image014.GIF" wi="18" he="18" />。本发明在接收方不知道信道状态信息的情况下,该设计方法可以利用少量的导频信号实现多小区多用户多天线系统中对发送信号的估计。
申请公布号 CN104333406A 申请公布日期 2015.02.04
申请号 CN201410613571.2 申请日期 2014.11.04
申请人 杭州电子科技大学 发明人 王海泉;刘杨懿;张婉;许倩
分类号 H04B7/08(2006.01)I;H04L1/00(2006.01)I 主分类号 H04B7/08(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 用于在多小区多用户多天线系统中的解码方法,假定系统模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><msqrt><mi>&rho;</mi></msqrt><mi>HBX</mi><mo>+</mo><mi>W</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600877350000011.GIF" wi="1193" he="93" /></maths>其中,有L个来自不同小区的用户,接收端的基站天线数为M;再假定每个用户发送长为T的信息序列,且序列前端均包含有长为τ的训练序列头,则发送信号可表示成L×T的发送信号矩阵X,矩阵中每个元素均匀地取自于标准4‑QAM,某小区基站接收信号可表示成M×T的接收信号矩阵Y;H为M×L的信道矩阵,ρ为信噪比,B=diag{β<sub>1</sub>,β<sub>2</sub>,…β<sub>L</sub>}为对角矩阵,β<sub>i</sub>表示第i个小区到基站的大尺度衰落因子,W表示随机噪声;其中发送信号矩阵X规定为以下形式:X=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>T</sub>],x<sub>i</sub>是第i个列向量,x<sub>iR</sub>和x<sub>iI</sub>表示x<sub>i</sub>的实部和虚部,把x<sub>iR</sub>和x<sub>iI</sub>排列成2L×1的列向量<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>iR</mi></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>iI</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000600877350000012.GIF" wi="151" he="160" /></maths>则X可分解成<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>X</mi><mi>re</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>R</mi></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>R</mi></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>x</mi><mi>TR</mi></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>I</mi></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>I</mi></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>x</mi><mi>TI</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000600877350000013.GIF" wi="582" he="158" /></maths>其中前面的τ列是训练序列头,后(T‑τ)列为用户的有用信息序列,且τ满足2<sup>τ‑3</sup><L≤2<sup>τ‑2</sup>其特征在于该方法包括以下步骤:第一步是训练序列头的设计;具体如下:X<sub>re</sub>的前τ列是训练序列头,X<sub>re</sub>的第一列中的元素均为1,X<sub>re</sub>的第2列至第τ列记为X<sub>τ</sub>,其中X<sub>τ</sub>中的元素可为1或‑1,则可产生2<sup>τ‑1</sup>个不同的行向量,任取其中的2L个行向量构成X<sub>τ</sub>;第二步是利用系统模型,用接收信号矩阵Y估计出发送信号矩阵X具体如下:2.1基站接收到信号矩阵为Y=[y<sub>1</sub>,y<sub>2</sub>…y<sub>T</sub>],y<sub>i</sub>是第i个列向量,y<sub>iR</sub>和y<sub>iI</sub>表示y<sub>i</sub>的实部和虚部,把y<sub>iR</sub>和y<sub>iI</sub>排列成2M×1的列向量<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>y</mi><mi>iR</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>iI</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000600877350000014.GIF" wi="158" he="162" /></maths>则Y可分解成<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>Y</mi><mi>re</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>y</mi><mrow><mn>1</mn><mi>R</mi></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mn>2</mn><mi>R</mi></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>y</mi><mi>TR</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mn>1</mn><mi>I</mi></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mn>2</mn><mi>I</mi></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>Y</mi><mi>TI</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000600877350000021.GIF" wi="595" he="158" /></maths>计算Y<sub>re</sub>的样本均值为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>re</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>y</mi><mi>jR</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>jI</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000600877350000022.GIF" wi="380" he="170" /></maths>可得协方差矩阵E[y<sub>re</sub>y<sub>re</sub><sup>t</sup>]的样本估计:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>&Sigma;</mi><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><mrow><mo>(</mo><msub><mi>Y</mi><mi>re</mi></msub><msup><msub><mi>Y</mi><mi>re</mi></msub><mi>t</mi></msup><mo>-</mo><mi>T</mi><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>re</mi></msub><msup><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>re</mi></msub><mi>t</mi></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600877350000023.GIF" wi="510" he="127" /></maths>2.2将Σ进行特征值分解,Σ较大的2L个特征值组成的对角阵<img file="FDA0000600877350000024.GIF" wi="919" he="343" />U是2M×2L的矩阵,U中的第i个列向量是与D<sub>1</sub>对角线上第i个特征值相对应的特征矢量,可以得到:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msqrt><mn>2</mn><msub><mi>D</mi><mn>1</mn></msub><mo>-</mo><msub><mi>I</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub></msqrt><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>U</mi><mi>t</mi></msup><msub><mi>Y</mi><mi>re</mi></msub><mo>=</mo><msup><mi>V</mi><mi>t</mi></msup><msub><mi>X</mi><mi>re</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600877350000025.GIF" wi="1217" he="133" /></maths>2.3在(2)式中V<sup>t</sup>和X<sub>re</sub>未知,对(2)式运用快速ICA算法估计X<sub>re</sub>,其中快速ICA算法中负熵的近似,使用的非线性函数为g(y)=tanh(y)或g(y)=yexp(‑y<sup>2</sup>/2),最终得到估计结果<img file="FDA0000600877350000026.GIF" wi="110" he="84" />对<img file="FDA0000600877350000027.GIF" wi="78" he="86" />的每个元素进行判决,如果大于0,则该元素赋值为1,若小于0,则该元素赋值为‑1;2.4若解出<img file="FDA0000600877350000028.GIF" wi="78" he="83" />第1列中某一个元素为负,则该元素所在行符号发生了改变,<img file="FDA0000600877350000029.GIF" wi="84" he="86" />中该行乘上‑1,符号得以确定;记<img file="FDA00006008773500000210.GIF" wi="67" he="85" />是<img file="FDA00006008773500000211.GIF" wi="80" he="84" />的第2列至第τ列,取X<sub>τ</sub>的第1行,与<img file="FDA00006008773500000212.GIF" wi="68" he="87" />每一行做内积运算,内积最大那一行的所在行即为<img file="FDA00006008773500000213.GIF" wi="70" he="86" />正确排序的第1行,再取X<sub>τ</sub>的第2行,与<img file="FDA00006008773500000214.GIF" wi="74" he="81" />每一行做内积运算,内积最大那一行的所在行即为<img file="FDA00006008773500000215.GIF" wi="80" he="83" />正确排序的第2行,由此类推,最终确定每一行的顺序;在<img file="FDA00006008773500000216.GIF" wi="82" he="77" />的行序确定后,<img file="FDA00006008773500000217.GIF" wi="77" he="85" />的第1至L行为X实部的估计<img file="FDA00006008773500000218.GIF" wi="104" he="84" />第L+1至2L行为X虚部的估计<img file="FDA00006008773500000219.GIF" wi="104" he="86" />最后得到X的估计<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mi>X</mi><mo>~</mo></mover><mo>=</mo><msub><mover><mi>X</mi><mo>~</mo></mover><mi>R</mi></msub><mo>+</mo><mi>j</mi><msub><mover><mi>X</mi><mo>~</mo></mover><mi>I</mi></msub><mo>.</mo></mrow>]]></math><img file="FDA00006008773500000220.GIF" wi="320" he="92" /></maths>
地址 310018 浙江省杭州市下沙高教园区2号大街