发明名称 在移动环境下高阶调制MIMO系统的检测方法
摘要 本发明提出了一种在移动环境下高阶调制MIMO系统的检测方法,主要解决现有高阶4q-QAM调制多输入多输出MIMO复用系统检测复杂度高的问题。其检测步骤是:利用4q-QAM的分层特性,对4q-QAM的星座点进行区域划分;根据权重系数较大的调制子层数据确定侯选区域;在侯选区域中按星座点与对应的接收数据间的最小距离准则进行星座点优选;采用球形检测SD算法对各层优选后的星座点进行搜索,确定发送数据。本发明与多状态球形检测MSD方法相比,其复杂度有明显下降,且误码率性能有改善;与球形检测SD相比,本发明的复杂度在高信噪比SNR或低信噪比SNR下均有大幅下降,可用于移动环境下高阶调制MIMO系统中接收机的设计与实现。
申请公布号 CN101958875A 申请公布日期 2011.01.26
申请号 CN201010289232.5 申请日期 2010.09.24
申请人 西安电子科技大学 发明人 任光亮;段昕利;杨丽花;张会宁
分类号 H04L27/36(2006.01)I;H04B7/06(2006.01)I;H04B17/00(2006.01)I 主分类号 H04L27/36(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种在移动环境下高阶调制MIMO系统的检测方法,包括如下步骤:(1)根据接收天线接收的数据及信道值,构建检测的发射天线数据公式:<maths num="0001"><![CDATA[<math><mrow><mover><mi>s</mi><mo>^</mo></mover><mo>=</mo><mo>[</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>3</mn></mrow></msub><mo>,</mo><mi>L</mi><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><msub><mi>N</mi><mi>T</mi></msub></mrow></msub><mo>]</mo><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><msub><mi>s</mi><mi>ti</mi></msub><mo>&Element;</mo><msub><mi>q</mi><mrow><msup><mn>4</mn><mi>q</mi></msup><mo>-</mo><mi>QAM</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><msub><mi>N</mi><mi>T</mi></msub><mo>)</mo></mrow></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><msub><mi>N</mi><mi>T</mi></msub></mrow><mn>1</mn></munderover><msup><mrow><mo>|</mo><mo>|</mo><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>s</mi><mi>ti</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>s</mi><mi>tj</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>&le;</mo><mi>C</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000026914560000012.GIF" wi="21" he="40" />代表最终确定的发射向量,<img file="FDA0000026914560000013.GIF" wi="36" he="52" />为最终检测确定的第i层对应的发射数据,1≤i≤N<sub>T</sub>,N<sub>T</sub>代表发射天线数,  y′=Q<sup>H</sup>y,H代表共轭转置,y为接收向量,y′<sub>i</sub>代表y′向量的第i个数据,H<sub>Ch</sub>=QR,H<sub>Ch</sub>代表信道矩阵,Q为列正交矩阵,R为上三角矩阵,r<sub>i,j</sub>是R的第i行和第j列的元素,s<sub>ti</sub>代表第i根发射天线发射的数据,1≤i≤N<sub>T</sub>,s<sub>ti</sub>∈4<sup>q</sup>-QAM,4<sup>q</sup>-QAM代表调制方式,q代表该4<sup>q</sup>-QAM的调制子层数,s<sub>tj</sub>代表第j根发射天线发射的数据,(i+1)≤j≤N<sub>T</sub>,C为球半径;(2)将4<sup>q</sup>-QAM调制数据用各个子层调制数据的线性组合表示:<maths num="0002"><![CDATA[<math><mrow><msub><mi>s</mi><mi>ti</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>q</mi></munderover><msup><mn>2</mn><mrow><mi>q</mi><mo>-</mo><mi>k</mi></mrow></msup><mrow><mo>(</mo><msqrt><mn>2</mn></msqrt><mo>/</mo><mn>2</mn><mo>)</mo></mrow><msub><mi>s</mi><mi>ik</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,s<sub>ik</sub>∈Q<sub>QPSK</sub>代表第i层第k个调制子层数据,Q<sub>QPSK</sub>代表QPSK调制方式所包含的星座点,<img file="FDA0000026914560000015.GIF" wi="246" he="84" />代表各调制子层权重系数,q代表该4<sup>q</sup>-QAM的调制子层数;(3)将1)式中的s<sub>ti</sub>用2)式代替,得到3)式,构建检测各层各调制子层的数据公式:<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mi>ik</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><msub><mi>s</mi><mi>ik</mi></msub><mo>&Element;</mo><msub><mi>Q</mi><mi>QPSK</mi></msub><mo>,</mo><mrow><mo>(</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><msub><mi>N</mi><mi>T</mi></msub><mo>)</mo></mrow></mrow></munder><mo>|</mo><mo>|</mo><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi></mrow></msub><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mn>2</mn><mrow><mi>q</mi><mo>-</mo><mi>p</mi></mrow></msup><mrow><mo>(</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><mo>)</mo></mrow><msub><mi>s</mi><mi>ip</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>s</mi><mi>tj</mi></msub><mo>|</mo><mo>|</mo><mo>&le;</mo><mi>C</mi><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mi>q</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000026914560000017.GIF" wi="43" he="52" />代表第i层第k个调制子层检测确定的数据;(4)利用3)式逐子层确定第i层中各调制子层数据<img file="FDA0000026914560000018.GIF" wi="66" he="52" />直到k=q-1,(5)将第i层中确定的前q-1个调制子层数据<img file="FDA0000026914560000019.GIF" wi="58" he="52" />(1≤k≤q)按式2)进行线性组合,得到用于第i层最终区域范围确定的值s<sub>ti(q-1)</sub>:<maths num="0004"><![CDATA[<math><mrow><msub><mi>s</mi><mrow><mi>ti</mi><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>q</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mn>2</mn><mrow><mi>q</mi><mo>-</mo><mi>k</mi></mrow></msup><mrow><mo>(</mo><msqrt><mn>2</mn></msqrt><mo>/</mo><mn>2</mn><mo>)</mo></mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mi>ik</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>(6)利用4)式进行第i层检测星座点范围的确定:当|r<sub>i,i</sub>|<η<sub>1</sub>时,在4<sup>q</sup>-QAM调制星座图中选取距离s<sub>ti(q-1)</sub>最近的十六个星座点作为第i层检测星座点范围,并称此范围为扩展区域二;当η<sub>1</sub>≤|r<sub>i,i</sub>|<η<sub>2</sub>时,在4<sup>q</sup>-QAM调制星座图中选取距离s<sub>ti(q-1)</sub>最近的九个星座点作为第i层检测星座点范围,并称此范围为扩展区域一;当η<sub>2</sub>≤|r<sub>i,i</sub>|时,在4<sup>q</sup>-QAM调制星座图中选取距离s<sub>ti(q-1)</sub>最近的四个星座点作为第i层检测星座点范围,并称此范围为基本区域,其中,η<sub>1</sub>和η<sub>2</sub>是用来确定不同区域划分的两个不同的门限值,  0≤η<sub>1</sub>≤1,η<sub>2</sub>>1,其具体值根据系统性能要求通过信道统计特性确定,|·|代表取绝对值;(7)在步骤(6)中确定的第i层检测星座点范围基础上,按照区域中星座点与第i层接收信号距离大小进行第i层星座点的自适应优选,即删除与第i层接收信号距离较大的区域中的星座点:对于基本区域,将φ<sub>1</sub>与|r<sub>i,i</sub>|和η<sub>2</sub>进行比较:当|r<sub>i,i</sub>|≥φ<sub>1</sub>,优选星座点的数目为Z<sub>1</sub>;当η<sub>2</sub>≤|r<sub>i,i</sub>|<φ<sub>1</sub>,优选星座点的数目为Z<sub>2</sub>;对于扩展区域一,将φ<sub>2</sub>与|r<sub>i,i</sub>|,η<sub>1</sub>和η<sub>2</sub>进行比较:当φ<sub>2</sub>≤|r<sub>i,i</sub>|<η<sub>2</sub>,优选星座点的数目为Z<sub>3</sub>;当η<sub>1</sub>≤|r<sub>i,i</sub>|<φ<sub>2</sub>,优选星座点的数目为Z<sub>4</sub>;对于扩展区域二,将φ<sub>3</sub>与|r<sub>i,i</sub>|和η<sub>1</sub>进行比较:当φ<sub>3</sub>≤|r<sub>i,i</sub>|<η<sub>1</sub>,优选星座点的数目为Z<sub>5</sub>;当|r<sub>i,i</sub>|<φ<sub>3</sub>,优选星座点的数目为Z<sub>6</sub>;其中φ<sub>1</sub>,φ<sub>2</sub>和φ<sub>3</sub>是针对不同区域划分确定优选星座点数目的三个不同的门限值,φ<sub>1</sub>≥η<sub>2</sub>,η<sub>1</sub>≤φ<sub>2</sub><η<sub>2</sub>,0≤φ<sub>3</sub><η<sub>1</sub>,其具体值和优选星座点的数目Z<sub>1</sub>,Z<sub>2</sub>,Z<sub>3</sub>,Z<sub>4</sub>,Z<sub>5</sub>和Z<sub>6</sub>基于不同的η<sub>1</sub>,η<sub>2</sub>及性能需求直接给出;(8)将第i层优选后星座点构成的集合记为L<sub>i</sub>,将3)式中的各子层数据的线性组合项<img file="FDA0000026914560000022.GIF" wi="359" he="126" />用各层数据s<sub>ti</sub>代替,且用优选后的星座点集合L<sub>i</sub>,1≤i≤N<sub>T</sub>给出s<sub>ti</sub>的检测范围,按5)式得第i层检测的发射数据取值集合<img file="FDA0000026914560000023.GIF" wi="62" he="47" /><maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>ti</mi></msub><mo>=</mo><mo>{</mo><msub><mi>s</mi><mi>ti</mi></msub><mo>&Element;</mo><msub><mi>L</mi><mi>i</mi></msub><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><msub><mi>N</mi><mi>T</mi></msub></mrow><mi>i</mi></munderover><msup><mrow><mo>|</mo><mo>|</mo><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>s</mi><mi>ti</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>tj</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>&le;</mo><mi>C</mi><mo>,</mo><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>tj</mi></msub><mo>&Element;</mo><msub><mi>L</mi><mi>j</mi></msub><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000026914560000032.GIF" wi="37" he="47" />为确定的第i层检测的发射数据取值集合,L<sub>i</sub>为第i层优选后的星座点集合,<img file="FDA0000026914560000033.GIF" wi="38" he="53" />为确定的第j层检测的发射数据范围,L<sub>j</sub>为第j层优选后的星座点集合,令i=i-1,重复步骤(4)-(8),直到确定出所有层检测的发射数据取值集合。(9)利用步骤(8)中确定的各层检测的发射数据取值集合<img file="FDA0000026914560000034.GIF" wi="63" he="47" />1≤i≤N<sub>T</sub>给出s<sub>ti</sub>的取值范围,找出与y′距离最近的各层候选星座点组合作为最终检测的各层发射数据:<maths num="0006"><![CDATA[<math><mrow><mover><mi>s</mi><mo>^</mo></mover><mo>=</mo><mo>[</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><mn>3</mn></mrow></msub><mi>L</mi><mo>,</mo><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>t</mi><msub><mi>N</mi><mi>T</mi></msub></mrow></msub><mo>]</mo><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><msub><mi>s</mi><mi>ti</mi></msub><mo>&Element;</mo><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mi>ti</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><msub><mi>N</mi><mi>T</mi></msub><mo>)</mo></mrow></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><msub><mi>N</mi><mi>T</mi></msub></mrow><mn>1</mn></munderover><msup><mrow><mo>|</mo><mo>|</mo><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>s</mi><mi>ti</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><msub><mi>r</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>s</mi><mi>tj</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>.</mo></mrow></math>]]></maths>其中,<img file="FDA0000026914560000036.GIF" wi="21" he="40" />代表最终确定的发射向量,<img file="FDA0000026914560000037.GIF" wi="36" he="52" />为最终检测确定的第i层对应的发射数据,1≤i≤N<sub>T</sub>,N<sub>T</sub>代表发射天线数,<img file="FDA0000026914560000038.GIF" wi="37" he="47" />代表确定的第i层检测的发射数据范围,1≤i≤N<sub>T</sub> y′=Q<sup>H</sup>y,H代表共轭转置,y为接收向量,y′<sub>i</sub>代表y′向量的第i个数据,H<sub>Ch</sub>=QR,H<sub>Ch</sub>代表信道矩阵,Q为列正交矩阵,R为上三角矩阵,r<sub>i,j</sub>是R的第i行和第j列的元素,s<sub>ti</sub>代表第i根发射天线发射的数据,1≤i≤N<sub>T</sub><img file="FDA0000026914560000039.GIF" wi="150" he="48" /><img file="FDA00000269145600000310.GIF" wi="38" he="48" />为第i层检测确定的发射数据取值集合,s<sub>tj</sub>代表第j根发射天线发射的数据,(i+1)≤j≤N<sub>T</sub>,<img file="FDA00000269145600000311.GIF" wi="151" he="53" /><img file="FDA00000269145600000312.GIF" wi="38" he="53" />为第j层检测确定的发射数据取值集合,C为球半径。
地址 710071 陕西省西安市太白南路2号