发明名称 MIMO认知无线电干扰网络中基于信道学的频谱共享方法
摘要 本发明涉及MIMO认知无线电干扰网络中基于信道学的频谱共享方法,通过分析接收的数据的二阶统计量获取需要的主用户干扰信道的空间特征,融合MIMO多子流传输和认知多信道切换技术,设计发送策略,避免对主用户的干扰,提高频谱共享效率;针对干扰信道状态信息缺失下次用户干扰网络与主用户系统共存问题,兼顾MIMO认知无线电网络的空间传输能力,融合主用户的信道切换能力,联合利用空间和频率资源的高效频谱共享方法,比以往单一资源利用方法的效率更高;分析了次用户通过设置信道学时间对网络性能的影响,说明通过控制时间可以有效控制对主用户的干扰功率;从自由度检测性能、主用户误码率、算法收敛性和网络容量等多方面分析。
申请公布号 CN104918261A 申请公布日期 2015.09.16
申请号 CN201510307075.9 申请日期 2015.06.08
申请人 中国人民解放军信息工程大学 发明人 任修坤;朱世磊;胡捍英;郑娜娥;赵远;陈松;王盛;李玉翔;范立岩
分类号 H04W16/14(2009.01)I 主分类号 H04W16/14(2009.01)I
代理机构 郑州大通专利商标代理有限公司 41111 代理人 陈大通
主权项 一种MIMO认知无线电干扰网络中基于信道学习的频谱共享方法,其特征在于:包含如下步骤:步骤1.根据次用户通信需求和天线配置,确定各个用户的发送自由度d<sub>k</sub>,设置次用户内部干扰门限Γ<sub>k</sub>,选择检测时间长度;步骤2.次用户选择通信频段和信道学习时间,实时监听信道状态,利用接收到的信号y<sub>1k</sub>(n),计算采样协方差矩阵<img file="FDA0000733538390000011.GIF" wi="539" he="149" />对采样协方差矩阵<img file="FDA0000733538390000012.GIF" wi="76" he="85" />进行特征值分解<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>^</mo></mover><mrow><mn>1</mn><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub><msub><mi>&Lambda;</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub><msubsup><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow><mi>H</mi></msubsup><mo>=</mo><mo>[</mo><msubsup><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></msubsup><mo>]</mo><msub><mi>&Lambda;</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub><msup><mrow><mo>[</mo><msubsup><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>V</mi><mrow><mn>1</mn><mi>k</mi></mrow><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></msubsup><mo>]</mo></mrow><mi>H</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000733538390000013.GIF" wi="922" he="86" /></maths>得到特征值<img file="FDA0000733538390000014.GIF" wi="273" he="77" />和干扰信道的正交矩阵<img file="FDA0000733538390000015.GIF" wi="242" he="81" />并根据计算结果对所选频段做调整,使得总子流数能够满足<img file="FDA0000733538390000016.GIF" wi="619" he="77" />其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mi>k</mi></munder><mrow><mo>(</mo><msub><mi>M</mi><mn>0</mn></msub><mo>-</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>N</mi><mi>L</mi></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>0</mn></msub></munderover><msup><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mrow><mo>(</mo><msub><mi>M</mi><mn>0</mn></msub><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mrow></msup></mrow><mrow><mfrac><mn>1</mn><mrow><msub><mi>M</mi><mn>0</mn></msub><mo>-</mo><mi>k</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>0</mn></msub></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>k</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>M</mi><mn>0</mn></msub><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mi>ln</mi><msub><mi>N</mi><mi>L</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000733538390000017.GIF" wi="1267" he="313" /></maths>M<sub>0</sub>为主用户的天线个数;步骤3.通过步骤2得到的干扰信道的正交矩阵<img file="FDA0000733538390000018.GIF" wi="273" he="94" />计算等效信道<img file="FDA0000733538390000019.GIF" wi="758" he="86" />初始化等效预编码矩阵A<sub>k</sub>(k=1,2,…,K),设置<img file="FDA00007335383900000110.GIF" wi="440" he="112" />其中e<sub>i</sub>(i=1,2,…,d<sub>k</sub>)为基本向量,H<sub>lk</sub>为次用户干扰信道,计算等效发送<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>M</mi><mi>eq</mi></msup><mo>=</mo><mi>M</mi><mo>-</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>0</mn></msub></mrow>]]></math><img file="FDA00007335383900000111.GIF" wi="288" he="83" /></maths>和接收天线数<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>N</mi><mi>eq</mi></msup><mo>=</mo><mi>N</mi><mo>-</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>0</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA00007335383900000112.GIF" wi="289" he="98" /></maths>步骤4.通过等效信道<img file="FDA00007335383900000113.GIF" wi="88" he="76" />及等效预编码矩阵A<sub>k</sub>,计算<img file="FDA00007335383900000114.GIF" wi="615" he="149" />对<img file="FDA00007335383900000115.GIF" wi="75" he="80" />进行特征值分解,取N<sup>eq</sup>‑d<sub>l</sub>个最大特征值对应的特征向量构造等效干扰子空间矩阵C<sub>l</sub>;步骤5.利用等效干扰子空间矩阵C<sub>l</sub>下标k和l的对称性,计算<img file="FDA00007335383900000116.GIF" wi="779" he="154" />对<img file="FDA00007335383900000117.GIF" wi="80" he="83" />进行特征值分解,取d<sub>k</sub>个最小特征值对应的特征向量构造等效预编码矩阵A<sub>k</sub>;步骤6.通过步骤5得到的等效预编码矩阵A<sub>k</sub>,计算用户的干扰协方差矩阵<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>l</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>&NotEqual;</mo><mi>l</mi></mrow><mi>K</mi></munderover><msub><mi>H</mi><mi>lk</mi></msub><msub><mi>A</mi><mi>k</mi></msub><msubsup><mi>A</mi><mi>k</mi><mi>H</mi></msubsup><msubsup><mi>H</mi><mi>lk</mi><mi>H</mi></msubsup><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>K</mi><mo>,</mo></mrow>]]></math><img file="FDA0000733538390000021.GIF" wi="800" he="146" /></maths>若trace(Q<sub>k</sub>)<Γ<sub>k</sub>对<img file="FDA0000733538390000022.GIF" wi="67" he="62" />成立,则<img file="FDA0000733538390000023.GIF" wi="193" he="81" />计算预编码矩阵<img file="FDA0000733538390000024.GIF" wi="238" he="82" />和接收干扰抑制矩阵<img file="FDA0000733538390000025.GIF" wi="275" he="83" />结束;否则,返回步骤4,重新迭代计算。
地址 450052 河南省郑州市高新区科学大道62号