发明名称 用于在多小区多用户多天线系统中基于ICA的解码方法
摘要 本发明涉及一种用于在多小区多用户多天线系统中基于ICA的解码方法。本发明在接收方不知道信道状态信息的情况下,该解码方法基于ICA,采用少量的非正交导频序列,实现多小区多用户多天线系统中的解码。该发明可以提高数据传输速率且减少由导频污染带来的失真。
申请公布号 CN105207701A 申请公布日期 2015.12.30
申请号 CN201510471184.4 申请日期 2015.08.04
申请人 杭州电子科技大学 发明人 王海泉;刘杨懿;张婉
分类号 H04B7/04(2006.01)I;H04B7/08(2006.01)I;H04L5/00(2006.01)I 主分类号 H04B7/04(2006.01)I
代理机构 浙江杭州金通专利事务所有限公司 33100 代理人 王佳健
主权项 用于在多小区多用户多天线系统中基于ICA的解码方法,设定系统模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><msqrt><mi>&rho;</mi></msqrt><mi>H</mi><mi>B</mi><mi>X</mi><mo>+</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774356450000011.GIF" wi="982" he="88" /></maths>其中,该系统共有L个用户,基站侧天线数为M;再设定每个用户发送长为T的信息序列,其中每个元素均匀地取自某一QAM星座,且序列前端均包含长为τ的导频序列;L个用户的发送信号记为L×T的矩阵X=[Φ S],其中X的前τ列是导频序列矩阵Φ,它的第i行表示第i个用户所发送的导频序列,后(T‑τ)列为用户的信息矩阵S,它的第i行表示第i个用户所发送的信息;基站的接收信号可表示成M×T的接收信号矩阵Y;H是一个维度为M×L的信道矩阵;ρ为信噪比,<img file="FDA0000774356450000012.GIF" wi="574" he="97" />为对角矩阵,β<sub>l</sub>表示第l个用户到基站的路径损耗和阴影衰落增益系数,N表示随机噪声矩阵,1≤l≤L,令G=HB;其特征在于该方法包括以下步骤:步骤1.设计导频序列,具体是:根据该系统的用户数L,得到一个整数τ,要求τ满足条件4<sup>τ‑2</sup>&lt;L≤4<sup>τ‑1</sup>;构造一个L×τ的导频序列矩阵Φ,要求Φ的第一列的每个元素都是1,其它元素取自于4‑QAM,且保证Φ的任意两个行向量都不相同;步骤2.由Y估计出X的,具体如下:2‑1.基站接收到信号矩阵为Y=[y<sub>1</sub> y<sub>2</sub> … y<sub>T</sub>],y<sub>i</sub>是第i个列向量,1≤i≤T;计算Y的样本均值为<img file="FDA0000774356450000013.GIF" wi="286" he="148" />可得协方差矩阵E(yy<sup>H</sup>|G)的样本估计:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Sigma;</mi><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><mrow><mo>(</mo><msup><mi>YY</mi><mi>H</mi></msup><mo>-</mo><mi>T</mi><msup><mover><mrow><mi>y</mi><mi>y</mi></mrow><mo>&OverBar;</mo></mover><mi>H</mi></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774356450000014.GIF" wi="448" he="125" /></maths>其中,上标H表示矩阵的共轭转置;2‑2.对Σ进行特征值分解使得<img file="FDA0000774356450000015.GIF" wi="260" he="72" />其中<img file="FDA0000774356450000016.GIF" wi="58" he="70" />是一个M×M的对角矩阵,且对角线上的元素按降序排列,有<img file="FDA0000774356450000017.GIF" wi="518" he="105" /><img file="FDA0000774356450000018.GIF" wi="61" he="76" />是一个M×M的酉矩阵;令<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>D</mi><mo>^</mo></mover><mi>S</mi></msub><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><mrow><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mi>L</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000774356450000021.GIF" wi="538" he="109" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>D</mi><mo>^</mo></mover><mi>N</mi></msub><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><mrow><msub><mover><mi>d</mi><mo>^</mo></mover><mrow><mi>L</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mrow><mi>L</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mi>M</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000774356450000022.GIF" wi="634" he="108" /></maths>把<img file="FDA0000774356450000023.GIF" wi="49" he="74" />分解成<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mover><mi>U</mi><mo>^</mo></mover><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>U</mi><mo>^</mo></mover><mi>S</mi></msub></mtd><mtd><msub><mover><mi>U</mi><mo>^</mo></mover><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000774356450000024.GIF" wi="354" he="107" /></maths>其中<img file="FDA0000774356450000025.GIF" wi="78" he="84" />是<img file="FDA0000774356450000026.GIF" wi="59" he="75" />的前L列,<img file="FDA0000774356450000027.GIF" wi="84" he="87" />是剩下的M‑L列;由这些分解,有<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&Sigma;</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>U</mi><mo>^</mo></mover><mi>S</mi></msub></mtd><mtd><msub><mover><mi>U</mi><mo>^</mo></mover><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mover><mi>D</mi><mo>^</mo></mover><mi>S</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mover><mi>D</mi><mo>^</mo></mover><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msubsup><mover><mi>U</mi><mo>^</mo></mover><mi>S</mi><mi>H</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mover><mi>U</mi><mo>^</mo></mover><mi>N</mi><mi>H</mi></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000774356450000028.GIF" wi="697" he="182" /></maths>0表示零矩阵;利用PCA得到:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>V</mi><mi>H</mi></msup><mi>X</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msqrt><mrow><msub><mover><mi>D</mi><mo>^</mo></mover><mi>S</mi></msub><mo>-</mo><msub><mi>I</mi><mi>L</mi></msub></mrow></msqrt><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mover><mi>U</mi><mo>^</mo></mover><mi>S</mi><mi>H</mi></msubsup><mi>Y</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774356450000029.GIF" wi="1237" he="139" /></maths>其中,I<sub>L</sub>为L阶的单位矩阵,V<sup>H</sup>为L×L的酉矩阵;2‑3.在(2)式中V<sup>H</sup>和X未知,对(2)式运用复数的ICA算法,经迭代得到一个权重矩阵,记为L×L的酉矩阵W<sup>H</sup>;令<img file="FDA00007743564500000210.GIF" wi="539" he="159" />由此得到<img file="FDA00007743564500000211.GIF" wi="85" he="71" /><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mo>=</mo><msqrt><mi>&rho;</mi></msqrt><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mi>&rho;</mi></mfrac><msub><mi>I</mi><mi>L</mi></msub><mo>+</mo><msup><mover><mi>G</mi><mo>^</mo></mover><mi>H</mi></msup><mover><mi>G</mi><mo>^</mo></mover><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mover><mi>G</mi><mo>^</mo></mover><mi>H</mi></msup><mi>Y</mi></mrow>]]></math><img file="FDA00007743564500000212.GIF" wi="615" he="164" /></maths>令<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mover><mi>&Phi;</mi><mo>^</mo></mover></mtd><mtd><mover><mi>S</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00007743564500000213.GIF" wi="285" he="109" /></maths>其中<img file="FDA00007743564500000214.GIF" wi="58" he="75" />和<img file="FDA00007743564500000215.GIF" wi="42" he="73" />分别是L×τ和L×(T‑τ)的矩阵;将<maths num="0010" id="cmaths0010"><math><![CDATA[<mfenced open = '[' close = ']'><mtable><mtr><mtd><mover><mi>&Phi;</mi><mo>^</mo></mover></mtd><mtd><mover><mi>S</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00007743564500000216.GIF" wi="170" he="108" /></maths>每一行的第1个元素的幅度量化为1,然后记第l行第1个元素为<img file="FDA00007743564500000217.GIF" wi="109" he="71" />再将<maths num="0011" id="cmaths0011"><math><![CDATA[<mfenced open = '[' close = ']'><mtable><mtr><mtd><mover><mi>&Phi;</mi><mo>^</mo></mover></mtd><mtd><mover><mi>S</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00007743564500000218.GIF" wi="187" he="103" /></maths>第l行的每一个元素都乘上<img file="FDA00007743564500000219.GIF" wi="123" he="71" />每一行都执行此操作,操作完成后得到的矩阵记为<img file="FDA00007743564500000220.GIF" wi="210" he="100" />然后对<img file="FDA00007743564500000221.GIF" wi="172" he="101" />除第1列外其它元素,根据与4‑QAM和所采用的星座QAM距离最短的原则,分别进行判决得到的矩阵记为<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><mover><mi>&Phi;</mi><mo>&OverBar;</mo></mover></mtd><mtd><mover><mi>S</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA00007743564500000222.GIF" wi="202" he="98" /></maths>取Φ的第1行,与<img file="FDA00007743564500000223.GIF" wi="50" he="80" />的每一行做内积运算,内积最大那一行所在的行作为X第1行的解;再取Φ的第2行,与<img file="FDA00007743564500000224.GIF" wi="55" he="67" />的每一行做内积运算,内积最大那一行所在的行作为X第2行的解,由此类推,最终得到X每一行的解。
地址 310018 浙江省杭州市下沙高教园区2号大街