发明名称 多天线系统中基于迫零判决反馈检测的功率控制方法
摘要 本发明多天线系统中基于迫零判决反馈检测的功率控制方法,特征是在最小化误块率的准则下,计算功率分配矩阵,通过部分反馈信道信息即只反馈功率信息,把信道矩阵与功率矩阵的乘积分解成一个酉阵和一个上三角阵,通过用功率矩阵加权发送信号,使分解得到的上三角矩阵具有等对角线的特点,以优化迫零判决反馈检测的检测顺序和检测性能,从而降低了平均误比特率和误块率。这种检测算法同时具有较低的复杂度和较低的系统平均误比特率和误块率,在实际系统中较易于实现。
申请公布号 CN1909398B 申请公布日期 2011.04.20
申请号 CN200510041397.X 申请日期 2005.08.03
申请人 中国科学技术大学 发明人 邱玲;张璇;赵昆;朱近康
分类号 H04B7/005(2006.01)I;H04L1/06(2006.01)I 主分类号 H04B7/005(2006.01)I
代理机构 安徽合肥华信知识产权代理有限公司 34112 代理人 余成俊
主权项 1.一种多天线系统中基于迫零判决反馈检测的功率控制方法,首先将M个发射天线和N个接收天线的迫零判决反馈检测系统中的接收信号表示为r=Hx+n的形式,其中,H是N×M信道增益矩阵,发送信号向量x=[x<sub>1</sub> x<sub>2</sub> …x<sub>M</sub>]<sup>T</sup>,噪声向量n=[n<sub>1</sub> n<sub>2</sub> …n<sub>N</sub>]<sup>T</sup>;在接收端把矩阵H分解成一个酉阵Q和一个上三角阵R,即:做QR分解H=[Q Q′][R 0]<sup>T</sup>,左乘Q<sup>H</sup>,变换后的接收信号<maths num="0001"><![CDATA[<math><mrow><mover><mi>r</mi><mo>~</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mn>1</mn></msub></mtd><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mn>2</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mi>M</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><mi>Q</mi><mi>H</mi></msup><mi>r</mi><mo>,</mo></mrow></math>]]></maths>变换后的噪声向量<maths num="0002"><![CDATA[<math><mrow><mover><mi>n</mi><mo>~</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mn>1</mn></msub></mtd><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mn>2</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mi>M</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><mi>Q</mi><mi>H</mi></msup><mi>n</mi><mo>,</mo></mrow></math>]]></maths>得到<img file="FA20181857200510041397X01C00013.GIF" wi="19" he="32" />中的第k个元素<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mover><mi>n</mi><mo>~</mo></mover><mi>k</mi></msub></mrow></math>]]></maths>k=1,2,...,M,其中,R的对角线元素R<sub>k,k</sub> k=1,2,...,M,R的第k行第i列元素R<sub>k,i</sub> k=1,2,...,M,i=k+1,k+2,...,M;然后进行初始化检测,检测发送信号向量的第M个元素x<sub>M</sub>,定义<maths num="0004"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mfrac><msub><mover><mi>r</mi><mo>~</mo></mover><mi>M</mi></msub><msub><mi>R</mi><mrow><mi>M</mi><mo>,</mo><mi>M</mi></mrow></msub></mfrac><mo>,</mo></mrow></math>]]></maths>通过硬判决<maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>M</mi></msub><mo>=</mo><mi>q</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></math>]]></maths>得到x<sub>M</sub>的估计值<img file="FA20181857200510041397X01C00017.GIF" wi="72" he="47" />再进行干扰抵消检测,当k=M-1时,将<img file="FA20181857200510041397X01C00018.GIF" wi="46" he="47" />代入公式<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mover><mi>n</mi><mo>~</mo></mover><mi>k</mi></msub></mrow></math>]]></maths>作干扰抵消,得到<maths num="0007"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mfrac><mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>M</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>M</mi></mrow></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>M</mi></msub></mrow><msub><mi>R</mi><mrow><mi>M</mi><mo>,</mo><mi>M</mi></mrow></msub></mfrac><mo>,</mo></mrow></math>]]></maths>并硬判决得到发送信号向量的第M-1个元素x<sub>M-1</sub>的估计值<img file="FA20181857200510041397X01C000111.GIF" wi="131" he="61" />循环该操作直到检测出发送信号向量的第1个元素x<sub>1</sub>的估计值<img file="FA20181857200510041397X01C000112.GIF" wi="65" he="54" />恢复出原始发送数据;其特征在于:将所述接收信号表示为考虑了功率控制的r=HPx+n的形式,其中,P是发送功率矩阵,在接收端对矩阵HP作QR分解HP=[Q Q′][R 0]<sup>T</sup>,进行左乘Q<sup>H</sup>的操作,其中Q和Q′分别是N×M和N×(N-M)酉阵,0是(N-M)×M全零矩阵,R是M×M上三角方阵,变换后的接收信号<maths num="0008"><![CDATA[<math><mrow><mover><mi>r</mi><mo>~</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mn>1</mn></msub></mtd><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mn>2</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mover><mi>r</mi><mo>~</mo></mover><mi>M</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><mi>Q</mi><mi>H</mi></msup><mi>r</mi><mo>,</mo></mrow></math>]]></maths>变换后的噪声向量<maths num="0009"><![CDATA[<math><mrow><mover><mi>n</mi><mo>~</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mn>1</mn></msub></mtd><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mn>2</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mover><mi>n</mi><mo>~</mo></mover><mi>M</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><mi>Q</mi><mi>H</mi></msup><mi>n</mi><mo>,</mo></mrow></math>]]></maths>得到<img file="FA20181857200510041397X01C000115.GIF" wi="19" he="32" />中的第k个元素<maths num="0010"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mover><mi>n</mi><mo>~</mo></mover><mi>k</mi></msub></mrow></math>]]></maths>k=1,2,...,M,其中,R的对角线元素R<sub>k,k</sub> k=1,2,...,M,R的第k行第i列元素R<sub>k,i</sub> k=1,2,...,M,i=k+1,k+2,...,M;采用无编码的q个星座点的正交幅度调制,将第k个被检测的符号的误比特率表示为P<sub>e,k</sub>,信噪比表示为ρ<sub>k</sub>,在高信噪比范围下有<maths num="0011"><![CDATA[<math><mrow><msub><mi>P</mi><mrow><mi>e</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>&ap;</mo><mn>0.2</mn><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msub><mrow><mn>1.6</mn><mi>&rho;</mi></mrow><mi>k</mi></msub><mrow><mi>q</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><mo>=</mo><mn>0.2</mn><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1.6</mn><mrow><mi>q</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mfrac><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi></mrow><mn>2</mn></msubsup><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>k=1,2,...,M,对某一确定的总发射功率,在接收端计算最小化误块率的优化的发送功率矩阵<maths num="0012"><![CDATA[<math><mrow><mi>P</mi><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mrow><mi>tr</mi><mrow><mo>(</mo><msup><mi>P</mi><mi>H</mi></msup><mi>P</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>t</mi></msub></mrow></munder><mo>{</mo><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>e</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></math>]]></maths>使得R矩阵是等对角矩阵;令R=RP,其中R是M×M的上三角方阵,且满足H=QR,求出各个天线优化的发送功率的解:<maths num="0013"><![CDATA[<math><mrow><msub><mi>P</mi><mi>k</mi></msub><mo>=</mo><msub><mi>P</mi><mi>t</mi></msub><mfrac><msubsup><mover><mi>R</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi></mrow><mrow><mo>-</mo><mn>2</mn></mrow></msubsup><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msubsup><mover><mi>R</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>-</mo><mn>2</mn></mrow></msubsup></mrow></mfrac></mrow></math>]]></maths>k=1,2,...,M,<maths num="0014"><![CDATA[<math><mrow><mi>P</mi><mo>=</mo><mi>diag</mi><mrow><mo>(</mo><msqrt><msub><mi>P</mi><mn>1</mn></msub></msqrt><mo>,</mo><msqrt><msub><mi>P</mi><mn>2</mn></msub></msqrt><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msqrt><msub><mi>P</mi><mi>M</mi></msub></msqrt><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>将该优化的发射功率矩阵P通过反馈信道从接收端反馈给发送端,各天线在发送端对信号乘以优化的发送功率矩阵P,通过各天线发送信号;接收端对优化的功率控制的发送符号进行迫零判决反馈检测,检测出原始发送数据。
地址 230026 安徽省合肥市金寨路96号