发明名称 宽带无线全双工MIMO通信系统回波自干扰自适应抑制方法
摘要 本发明公开了一种宽带无线全双工多天线MIMO通信系统回波自干扰自适应抑制方法,将自适应滤波法的迭代思想运用到回波自干扰抑制方法中,通过自适应迭代算法,实现对自干扰信号的准确估计及系统误码率性能的提升,并且充分考虑非理想衰落信道传输条件以及多天线MIMO系统配置,提出了基于归一化最小均方误差NLMS算法的宽带无线全双工MIMO通信系统回波自干扰抑制方法。本发明能在显著抑制高强度自干扰信号的同时,有效降低了系统误码率损失,提升了系统信道传输容量,以及降低了自干扰抑制的实现复杂度,综合提高了无线全双工系统的整体综合性能。
申请公布号 CN104468055A 申请公布日期 2015.03.25
申请号 CN201410840597.0 申请日期 2014.12.29
申请人 西安电子科技大学 发明人 王勇;袁璐;宫丰奎;葛建华
分类号 H04L1/06(2006.01)I 主分类号 H04L1/06(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 一种宽带无线全双工MIMO通信系统回波自干扰自适应抑制方法,其特征在于,该宽带无线全双工MIMO通信系统回波自干扰自适应抑制方法将迭代滤波法的迭代思想运用到回波自干扰抑制方法中,通过自适应迭代算法实现对自干扰信号的准确估计及系统误码率BER性能的提升;具体步骤如下:步骤一:近端通信节点的接收信号<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>t</mi><mi>R</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msubsup><mi>t</mi><msub><mi>N</mi><mn>1</mn></msub><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><msub><mi>N</mi><mi>R</mi></msub><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mo>]</mo><mi>T</mi></msup><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>N</mi><mi>R</mi></msub><mo>&times;</mo><msub><mi>N</mi><mi>f</mi></msub></mrow></msup></mrow>]]></math><img file="FDA0000646600170000011.GIF" wi="834" he="93" /></maths>为:t<sup>R</sup>(n)=H<sup>FE</sup>(n)t<sup>FE</sup>(n)+H<sup>NE</sup>(n)t<sup>NE</sup>(n)+w(n);其中,<img file="FDA0000646600170000012.GIF" wi="812" he="93" />为来自远端节点的有用目标接收信号;而<img file="FDA0000646600170000013.GIF" wi="786" he="92" />为近端节点自身发射信号,即回波自干扰信号;<img file="FDA0000646600170000014.GIF" wi="324" he="84" />分别表示近端和远端第j(j=N<sub>1</sub>,…,N<sub>T</sub>)条天线上的发送信号;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>H</mi><mi>FE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msubsup><mi>h</mi><msub><mi>N</mi><mn>1</mn></msub><mi>FE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>h</mi><mi>j</mi><mi>FE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>h</mi><msub><mi>N</mi><mi>T</mi></msub><mi>FE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>N</mi><mi>R</mi></msub><mo>&times;</mo><msub><mi>N</mi><mi>T</mi></msub></mrow></msup></mrow>]]></math><img file="FDA0000646600170000015.GIF" wi="954" he="87" /></maths>与<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>H</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msubsup><mi>h</mi><msub><mi>N</mi><mn>1</mn></msub><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>h</mi><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>h</mi><msub><mi>N</mi><mi>T</mi></msub><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>N</mi><mi>R</mi></msub><mo>&times;</mo><msub><mi>N</mi><mi>T</mi></msub></mrow></msup></mrow>]]></math><img file="FDA0000646600170000016.GIF" wi="974" he="84" /></maths>分别为远端和近端发射信号的信道转移函数;w(n)为信道加性高斯白噪声;其中,N<sub>T</sub>表示通信节点发射天线数目,N<sub>R</sub>是接收天线数目,N<sub>f</sub>是信号每帧长度,(·)<sup>T</sup>表示对矩阵或矢量的转置运算符号;步骤二:在接收端通过利用归一化最小均方误差NLMS算法对混有自干扰、信道噪声的接收信号进行自干扰抑制,定义算法的代价函数为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>Min</mi><mo>{</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>(</mo><msup><mi>e</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>}</mo><mo>=</mo><mi>Min</mi><mo>{</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>[</mo><msup><mi>t</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>t</mi><mo>^</mo></mover><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000017.GIF" wi="1109" he="139" /></maths>其中,Min表示取最小值,n表示第n时刻,E[(e<sup>NE</sup>(n))<sup>2</sup>]表示近端误差信号<img file="FDA0000646600170000018.GIF" wi="518" he="87" />的平均功率,E[·]表示期望运算符,t<sup>NE</sup>(n)表示近端发送天线的实际发送信号,<img file="FDA0000646600170000019.GIF" wi="162" he="90" />表示对近端总接收信号滤波后,获得的对近端发送信号t<sup>NE</sup>(n)的估计值;步骤三:设置采用归一化最小均方误差NLMS算法进行自干扰抑制的相关初始值:令初始迭代次数k=1,并设置最大迭代次数K及根据近端输入信号的自相关矩阵设置收敛步长因子μ<sup>NE</sup>,自适应滤波器的初始化权值矢量α<sup>NE</sup>(0)以及滤波器的长度M,开始迭代过程,分别设置K、M、<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>&alpha;</mi><mi>NE</mi></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>&alpha;</mi><mn>1</mn><mi>NE</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>M</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mrow><mo>[</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>0</mn><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000021.GIF" wi="1294" he="102" /></maths>μ<sup>NE</sup>;步骤四:根据公式<img file="FDA0000646600170000022.GIF" wi="555" he="93" />按照以下公式求出近端的估计信号<img file="FDA0000646600170000023.GIF" wi="182" he="92" />具体过程如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msubsup><mover><mi>t</mi><mo>^</mo></mover><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>&alpha;</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><msup><mrow><mo>(</mo><msup><mrow><mo>[</mo><msubsup><mi>&alpha;</mi><mn>1</mn><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>M</mi><mi>NE</mi></msubsup><mmultiscripts><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mmultiscripts><mo>]</mo></mrow><mi>T</mi></msup><mo>)</mo></mrow><mi>T</mi></msup><mo>&CenterDot;</mo><msup><mrow><mo>[</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup></mtd></mtr><mtr><mtd><mo>=</mo><mo>[</mo><msubsup><mi>&alpha;</mi><mn>1</mn><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&alpha;</mi><mi>M</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><msup><mrow><mo>[</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000024.GIF" wi="1623" he="493" /></maths>其中j=N<sub>1</sub>,…,N<sub>T</sub>,N<sub>T</sub>表示发送天线总数目,M为自适应滤波器的长度,α<sup>NE</sup>(n)在表示n时刻的权值矢量,<img file="FDA0000646600170000025.GIF" wi="161" he="86" />为n时刻第j条接收天线经自适应滤波后获得的近端误差信号,<img file="FDA0000646600170000026.GIF" wi="134" he="90" />为第j条近端接收天线上的接收信号;j<N<sub>T</sub>,则令j=j+1,估计下一接收天线上的估计信号<img file="FDA0000646600170000027.GIF" wi="222" he="91" />j=N<sub>T</sub>,则前进至步骤五;步骤五:根据下式,更新n时刻的权值矢量<img file="FDA0000646600170000028.GIF" wi="727" he="109" />并根据迭代结果输出近端发送信号t<sup>NE</sup>(n)的估计信号<img file="FDA0000646600170000029.GIF" wi="185" he="90" />具体过程如下:如下式更新下一时刻的权值矢量:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>&alpha;</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mi>&alpha;</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msup><mi>&mu;</mi><mi>NE</mi></msup><mrow><msup><mrow><mo>(</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><msubsup><mi>e</mi><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>t</mi><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00006466001700000210.GIF" wi="1228" he="174" /></maths>其中,j=1,…,N<sub>T</sub>,N<sub>T</sub>表示发送天线总数目,ε表示的是权值矢量α<sup>NE</sup>(n)在迭代过程中的调整因数,<img file="FDA00006466001700000211.GIF" wi="135" he="87" />为第j条近端接收天线上的接收信号,<img file="FDA00006466001700000212.GIF" wi="161" he="86" />为n时刻第j条接收天线经NLMS自适应滤波后获得的近端误差信号,μ<sup>NE</sup>表示收敛步长因子,(·)<sup>T</sup>表示对矩阵或矢量的转置运算符;步骤六:根据最佳权值矢量α<sup>NE</sup>(n)以及公式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mover><mi>t</mi><mo>^</mo></mover><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>&alpha;</mi><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>t</mi><mi>j</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000031.GIF" wi="1079" he="140" /></maths>由下式得近端估计信号<img file="FDA0000646600170000032.GIF" wi="152" he="90" />的最终表达式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msup><mover><mi>t</mi><mo>^</mo></mover><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mover><mi>t</mi><mo>^</mo></mover><msub><mi>N</mi><mn>1</mn></msub><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mover><mi>t</mi><mo>^</mo></mover><mi>j</mi><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mover><mi>t</mi><mo>^</mo></mover><msub><mi>N</mi><mi>T</mi></msub><mi>NE</mi></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000033.GIF" wi="876" he="90" /></maths>其中,j=1,…,N<sub>T</sub>,N<sub>T</sub>表示发送天线总数目,α<sup>NE</sup>(n)表示n时刻权值矢量,<img file="FDA0000646600170000034.GIF" wi="172" he="84" />表示n时刻的权值,其中i=1,…,M,M表示滤波器的长度,<img file="FDA0000646600170000035.GIF" wi="120" he="84" />表示第j条接收天线上的接收信号;步骤七:从总体接收信号t<sup>R</sup>(n)中滤除估计出的回波自干扰信号,以获得来自远端节点的有用传输信号,将该信号送入后续的MIMO译码检测单元,以获得对远端发送信号的准确估计,具体包括:第一步,从接收信号t<sup>R</sup>(n)中减去回波自干扰估计信号<img file="FDA0000646600170000036.GIF" wi="176" he="93" />得到来自远端节点的有用传输信号t<sup>ES</sup>(n),即:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msup><mi>t</mi><mi>ES</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>H</mi><mi>FE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>t</mi><mi>FE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>t</mi><mi>R</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>t</mi><mo>^</mo></mover><mi>NE</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000646600170000037.GIF" wi="880" he="90" /></maths>第二步,将信号t<sup>ES</sup>(n)送入后续的MIMO译码检测单元,以获得对远端发送信号t<sup>FE</sup>(n)的估计。
地址 710071 陕西省西安市太白南路2号西安电子科技大学