发明名称 基于基扩展模型的OFDM系统信道估计与信号检测方法
摘要 基于基扩展模型的OFDM系统信道估计与信号检测方法,属于无线与移动通信技术领域。具体包括以下步骤:发送端发送OFDM信号、接收端信道建模、初始化、估计基扩展模型系数、信道均衡与信号检测、迭代、输出。接收端信道建模对接收端频域信号进行干扰和有用信息的分离,初始化假设无循环前缀限制,估计出基扩展模型系数并且检测得到当前发送数据符号,通过迭代求出符号间干扰频响以及循环前缀重构部分的频响,进而消除符号间干扰带来的影响。本发明提出的估计方法,一方面,可以估计时变性明显的信道,又能有效地消除多径时延引入的符号间干扰,完成对双选信道下OFDM系统的信道估计以及信号检测,提高系统性能。
申请公布号 CN103441967A 申请公布日期 2013.12.11
申请号 CN201310389937.8 申请日期 2013.08.31
申请人 电子科技大学 发明人 雷霞;曹海波;宋阳;罗阳;陈晓
分类号 H04L25/02(2006.01)I;H04L25/03(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李顺德;王睿
主权项 1.基于基扩展模型的OFDM系统信道估计与信号检测方法,包括如下步骤:步骤一:设当前OFDM系统发送端发送的导频和数据混合的频域OFDM符号为X<sup>i</sup>,记为<img file="FDA0000375531660000011.GIF" wi="434" he="80" />其中N表示OFDM符号的长度,频域OFDM符号X<sup>i</sup>中一些确定位置的子载波分配导频,其它为数据子载波;则发送端当前发送的相应时域OFDM符号x<sup>i</sup>可记为<img file="FDA0000375531660000012.GIF" wi="394" he="91" />且x<sup>i</sup>=F<sup>H</sup>X<sup>i</sup>,其中F<sup>H</sup>表示N点的IFFT变换矩阵,i表示第i个OFDM符号,即OFDM符号的索引号;步骤二:假设时域OFDM符号x<sup>i</sup>经快时变和多径效应叠加的双选信道后,接收端接收到的时域OFDM符号为<maths num="0001"><![CDATA[<math><mrow><msup><mi>y</mi><mi>i</mi></msup><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>y</mi><mn>0</mn><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mn>1</mn><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mn>2</mn><mi>i</mi></msubsup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msubsup><mi>y</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>利用复指数基扩展模型描述快时变和多径效应叠加的双选信道h=[h<sub>0</sub>,h<sub>1</sub>,...,h<sub>l</sub>,...,h<sub>L</sub>],其中h<sub>l</sub>=[h<sub>0,l</sub>,h<sub>1,l</sub>,…,h<sub>n,l</sub>,…,h<sub>N-1,l</sub>]表示第l条径冲激响应,如果用h<sub>q</sub>(l)表示信道第l条径对应的BEM系数,b<sub>q</sub>表示基函数,则有:<maths num="0002"><![CDATA[<math><mrow><msub><mi>h</mi><mi>l</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>b</mi><mi>q</mi></msub><msub><mi>h</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>b</mi><mi>q</mi></msub><mo>=</mo><mo>[</mo><mn>1</mn><mo>,</mo><msup><mi>e</mi><mrow><mfrac><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mi>Q</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow></msup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msup><mi>e</mi><mrow><mfrac><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mi>Q</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow></msup><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA00003755316600000111.GIF" wi="992" he="71" />其中,<img file="FDA0000375531660000016.GIF" wi="63" he="73" />表示向上取整,f<sub>d</sub>表示最大多普勒频移、且f<sub>d</sub>=f<sub>c</sub>v/c,而f<sub>c</sub>表示OFDM系统载波的中心频率,v表示发送端和接收端之间的相对运动速度,c表示光速,T<sub>s</sub>表示采样周期;时域OFDM符号x<sup>i</sup>经快时变和多径效应叠加的双选信道后,接收端接收到的时域OFDM符号y<sup>i</sup>可表示为:<maths num="0004"><![CDATA[<math><mrow><msup><mover><mi>y</mi><mo>~</mo></mover><mi>i</mi></msup><mo>=</mo><msup><mi>y</mi><mi>i</mi></msup><mo>-</mo><msubsup><mi>y</mi><mi>ISI</mi><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>Add</mi><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000375531660000018.GIF" wi="48" he="66" />表示OFDM系统中循环前缀的长度大于信道的最大多径时延时接收端接收到第i个OFDM符号,<img file="FDA0000375531660000019.GIF" wi="67" he="67" />代表前一个OFDM发送符号x<sup>i-1</sup>对当前OFDM接收符号y<sup>i</sup>的干扰,<img file="FDA00003755316600000110.GIF" wi="85" he="67" />代表循环前缀的重构部分;将式子(4)变换到频域,则:<maths num="0005"><![CDATA[<math><mrow><msup><mover><mi>Y</mi><mo>~</mo></mover><mi>i</mi></msup><mo>=</mo><msup><mi>Y</mi><mi>i</mi></msup><mo>-</mo><msubsup><mi>Y</mi><mi>ISI</mi><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>Y</mi><mi>Add</mi><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>Y</mi><mi>ISI</mi><mi>i</mi></msubsup><mo>=</mo><mi>F</mi><mrow><mo>(</mo><msubsup><mi>F</mi><mn>1</mn><mi>i</mi></msubsup><mo>)</mo></mrow><msup><mi>F</mi><mi>H</mi></msup><msup><mi>X</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>Y</mi><mi>Add</mi><mi>i</mi></msubsup><mo>=</mo><mi>F</mi><mrow><mo>(</mo><msubsup><mi>H</mi><mn>1</mn><mi>i</mi></msubsup><mo>)</mo></mrow><msup><mi>F</mi><mi>H</mi></msup><msup><mi>X</mi><mi>i</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三:初始化;设置迭代次数I=0,令<img file="FDA0000375531660000024.GIF" wi="299" he="78" />这样<img file="FDA0000375531660000025.GIF" wi="180" he="73" />步骤四:估计相应的信道基扩展模型系数;接收端信号Y中导频位置观测量Y<sup>p</sup>可以表示为:<maths num="0008"><![CDATA[<math><mrow><msup><mi>Y</mi><mi>p</mi></msup><mo>=</mo><msup><mover><mi>Y</mi><mo>~</mo></mover><mi>p</mi></msup></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><mo>=</mo><msup><mi>D</mi><mi>p</mi></msup><msup><mi>S</mi><mi>p</mi></msup><msub><mi>h</mi><mi>eqv</mi></msub><mo>+</mo><msup><mi>D</mi><mi>d</mi></msup><msup><mi>S</mi><mi>d</mi></msup><msub><mi>h</mi><mi>eqv</mi></msub><mo>+</mo><msup><mi>W</mi><mi>p</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><mo>=</mo><msub><mi>Ph</mi><mi>eqv</mi></msub><mo>+</mo><msup><mi>D</mi><mi>d</mi></msup><msup><mi>S</mi><mi>d</mi></msup><msub><mi>h</mi><mi>eqv</mi></msub><mo>+</mo><msup><mi>W</mi><mi>p</mi></msup></mrow></math>]]></maths>其中,Y<sup>p</sup>表示接收端导频位置处的频响,Ph<sub>eqv</sub>表示导频产生的频响,D<sup>d</sup>S<sup>d</sup>h<sub>eqv</sub>表示OFDM符号数据子载波对导频子载波的干扰部分,W<sup>p</sup>表示接收端导频位置处的频域噪声;利用LS估计准则对h<sub>eqv</sub>进行估计,则估计方程为:<maths num="0011"><![CDATA[<math><mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>eqv</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>P</mi><mi>H</mi></msup><mi>P</mi><mo>+</mo><mi>&alpha;I</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>P</mi><mi>H</mi></msup><msup><mi>Y</mi><mi>p</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>考虑对矩阵P求逆运算的准确性,常常加上一个小的扰动,其中α为扰动因子。将估计得到的BEM系数<img file="FDA00003755316600000210.GIF" wi="72" he="82" />代入等式(1),即可求出信道时域冲激响应<img file="FDA00003755316600000211.GIF" wi="57" he="81" />步骤五:信道均衡与信号检测;利用估计得到的信道时域冲激响应<img file="FDA00003755316600000212.GIF" wi="36" he="63" />构造矩阵<img file="FDA00003755316600000213.GIF" wi="91" he="63" />求出信道频率响应矩阵<img file="FDA00003755316600000214.GIF" wi="89" he="78" />利用MMSE准则构造第i帧OFDM数据单元的频域均衡器有<img file="FDA00003755316600000215.GIF" wi="517" he="117" />其中,σ<sup>2</sup>是噪声方差,I<sub>N</sub>是N×N的单位矩阵,然后检测得到发端当前频域符号矢量<img file="FDA00003755316600000216.GIF" wi="84" he="82" />步骤六:判断;令:I=I+1                (10)如果I≤M,执行步骤七,否则,转到步骤八步骤七:迭代;利用估计得到的信道时域冲激响应<img file="FDA0000375531660000031.GIF" wi="38" he="64" />构造矩阵<img file="FDA0000375531660000032.GIF" wi="84" he="68" />再加上前一时刻估计得到的频域符号矢量<img file="FDA0000375531660000033.GIF" wi="86" he="69" />和当前频域符号矢量<img file="FDA0000375531660000034.GIF" wi="64" he="69" />时域代入式子(6)(7),求出<img file="FDA0000375531660000035.GIF" wi="70" he="82" />和<img file="FDA0000375531660000036.GIF" wi="88" he="82" />后更新等式(5),然后转到步骤四;步骤八:检测得到发端当前频域符号矢量<img file="FDA0000375531660000037.GIF" wi="88" he="71" />计算最后的估计量<img file="FDA0000375531660000038.GIF" wi="146" he="83" />
地址 611731 四川省成都市高新区(西区)西源大道2006号