发明名称 空间相关MISO系统中基于不完全信道信息预编码方法
摘要 本发明公开了一种空间相关MISO系统中基于步完全信道信息的预编码方法,涉及无线通信中自适应资源分配技术领域,本发明包括建立MISO系统和信道模型、确定系统有效信噪比、确定条件相关矩阵及预编码波束成型矩阵、确定预编码中功率分配、最小化系统平均比特率、确定系统平均误比特率等步骤;本发明的预编码方法可获得闭式功率分配,避免了已有最优方法的迭代计算;所提供的功率分配表达式较为简单,可直接计算;考虑了实际中反馈错误,有着很好的实用性;整个方法的复杂度低,代价小,时延少,易于实现。
申请公布号 CN103138815B 申请公布日期 2016.04.20
申请号 CN201110392698.2 申请日期 2011.12.01
申请人 南京航空航天大学 发明人 虞湘宾;陈小敏;刘晓帅;殷馨;辛瑜玉
分类号 H04B7/06(2006.01)I;H04B7/08(2006.01)I 主分类号 H04B7/06(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 一种空间相关MISO系统中基于不完全信道信息的预编码方法,其特征在于,包括以下步骤:步骤一:建立系统模型和相应的信道模型;建立空间相关无线MISO通信系统模型,该系统有着M个发射天线和单个接收天线,其对应信道模型可表示为信道矩阵h={h<sub>m</sub>},建模为:<img file="FDA0000856765870000011.GIF" wi="271" he="80" />其中h<sub>w</sub>为信道的随机变化部分,R<sub>t</sub>为发送天线相关矩阵,表示信道协方差矩阵;步骤二:确定系统的有效信噪比;由接收信号矩阵,求得接收端的有效信噪比<img file="FDA0000856765870000012.GIF" wi="382" he="104" />其中:<img file="FDA0000856765870000013.GIF" wi="352" he="108" />为对第m个发送天线进行波束成型后信道功率值;<img file="FDA0000856765870000014.GIF" wi="246" he="82" />为平均SNR,u<sub>im</sub>为波束成形矩阵U的元素,r是空时编码的码率;{P<sub>m</sub>,m=1,…,M}为对M个特征波束的功率控制,满足下列功率约束:<math><![CDATA[<mrow><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000856765870000015.GIF" wi="283" he="111" /></maths>P<sub>m</sub>≥0;步骤三:确定条件相关矩阵和波束成型矩阵;利用预白化<img file="FDA0000856765870000016.GIF" wi="236" he="85" />和变量变换,<img file="FDA0000856765870000017.GIF" wi="64" he="78" />为信道白化矢量,<img file="FDA0000856765870000018.GIF" wi="62" he="78" />为条件相关矩阵;来自步骤二的接收信噪比可进一步表示为:<math><![CDATA[<mrow><mi>&gamma;</mi><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>P</mi><mi>m</mi></msub><msub><mi>&zeta;</mi><mi>m</mi></msub><mo>|</mo><msub><mover><mi>h</mi><mo>&OverBar;</mo></mover><mi>m</mi></msub><msup><mo>|</mo><mn>2</mn></msup><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mo>/</mo><mi>r</mi><mo>,</mo></mrow>]]></math><img file="FDA0000856765870000019.GIF" wi="526" he="112" /></maths>其中:<math><![CDATA[<mrow><mover><mi>h</mi><mo>&OverBar;</mo></mover><mo>=</mo><msub><mover><mi>h</mi><mo>~</mo></mover><mi>w</mi></msub><msub><mi>V</mi><mi>h</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA00008567658700000110.GIF" wi="237" he="83" /></maths>与<img file="FDA00008567658700000111.GIF" wi="63" he="85" />有着相同的分布,即其元素为服从均值为0和方差为1的复高斯分布,<img file="FDA00008567658700000112.GIF" wi="60" he="79" />为矢量<img file="FDA00008567658700000113.GIF" wi="45" he="60" />中第m个元素,<img file="FDA00008567658700000114.GIF" wi="51" he="60" />为平均信噪比;<img file="FDA00008567658700000115.GIF" wi="81" he="84" />为条件相关矩阵,<math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>~</mo></mover><mi>h</mi></msub><mo>=</mo><mi>E</mi><mo>{</mo><msup><mi>h</mi><mi>H</mi></msup><mi>h</mi><mo>}</mo><mo>=</mo><msup><mi>c</mi><mn>2</mn></msup><msup><mover><mi>h</mi><mo>^</mo></mover><mi>H</mi></msup><mover><mi>h</mi><mo>^</mo></mover><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>c</mi><mn>2</mn></msup><mo>)</mo></mrow><msub><mi>R</mi><mi>t</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA00008567658700000116.GIF" wi="685" he="87" /></maths>其中:<img file="FDA00008567658700000117.GIF" wi="50" he="71" />为信道均值反馈,c为用于测量MMSE估计精度的归一化相关系数;进行特征值分解为:<img file="FDA00008567658700000118.GIF" wi="364" he="85" />其中波束成形矩阵U=V<sub>h</sub>,V<sub>h</sub>为对<img file="FDA00008567658700000119.GIF" wi="61" he="76" />特征值分解后由特征矢量构成的特征矩阵,Λ<sub>h</sub>为特征值分解后由特征值构成的对角矩阵;步骤四:确定系统平均误比特率;根据获得的均值反馈,利用随机变量变换,通过平均化卡方分布,确定出系统平均误比特率(BER):<img file="FDA00008567658700000120.GIF" wi="571" he="145" />其中λ=gE<sub>s</sub>/(rN<sub>0</sub>)为等效符号信噪比,g为调制参数;步骤五:最小化平均误比特率;在固定功率约束下,以最小化步骤四中的系统BER为优化目标,建立相应的目标优化函数<math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>P</mi><mn>1</mn></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mi>P</mi><mi>M</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>&lsqb;</mo><mi>log</mi><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msub><mi>&lambda;P</mi><mi>m</mi></msub><msub><mi>&zeta;</mi><mi>m</mi></msub></mrow><mo>)</mo><mo>&rsqb;</mo><mo>+</mo><msub><mi>&eta;</mi><mn>1</mn></msub><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>P</mi><mi>m</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000856765870000021.GIF" wi="1076" he="154" /></maths>其中:η<sub>1</sub>为拉格朗日乘子;步骤六:确定功率分配;对步骤五中目标函数关于{P<sub>m</sub>}求偏导,并令偏导值等于零,利用功率约束条件,得到闭式功率分配<math><![CDATA[<mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>&lambda;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msubsup><mi>&zeta;</mi><mi>l</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>&zeta;</mi><mi>m</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>+</mo></msup><mo>;</mo></mrow>]]></math><img file="FDA0000856765870000022.GIF" wi="837" he="112" /></maths>其中:<math><![CDATA[<mrow><msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo></msup><munder><mi>&Delta;</mi><mo>=</mo></munder><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000856765870000023.GIF" wi="358" he="85" /></maths>
地址 210016 江苏省南京市白下区御道街29号