发明名称 基于基扩展模型的OFDM系统信道估计与信号检测方法
摘要 基于基扩展模型的OFDM系统信道估计与信号检测方法,属于无线与移动通信技术领域。具体包括以下步骤:发送端发送OFDM信号、接收端信道建模、初始化、估计基扩展模型系数、信道均衡与信号检测、迭代、输出。接收端信道建模对接收端频域信号进行干扰和有用信息的分离,初始化假设无循环前缀限制,估计出基扩展模型系数并且检测得到当前发送数据符号,通过迭代求出符号间干扰频响以及循环前缀重构部分的频响,进而消除符号间干扰带来的影响。本发明提出的估计方法,一方面,可以估计时变性明显的信道,又能有效地消除多径时延引入的符号间干扰,完成对双选信道下OFDM系统的信道估计以及信号检测,提高系统性能。
申请公布号 CN103441967B 申请公布日期 2017.03.15
申请号 CN201310389937.8 申请日期 2013.08.31
申请人 电子科技大学 发明人 雷霞;曹海波;宋阳;罗阳;陈晓
分类号 H04L25/02(2006.01)I;H04L25/03(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李顺德;王睿
主权项 基于基扩展模型的OFDM系统信道估计与信号检测方法,包括如下步骤:步骤一:设当前OFDM系统发送端发送的导频和数据混合的频域OFDM符号为X<sup>i</sup>,记为<img file="FDA0000941488670000011.GIF" wi="422" he="83" />其中N表示OFDM符号的长度,频域OFDM符号X<sup>i</sup>中一些确定位置的子载波分配导频,其它为数据子载波;则发送端当前发送的相应时域OFDM符号x<sup>i</sup>可记为<img file="FDA0000941488670000012.GIF" wi="388" he="83" />且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符号为<img file="FDA0000941488670000013.GIF" wi="446" he="83" />利用复指数基扩展模型描述快时变和多径效应叠加的双选信道<img file="FDA00009414886700000112.GIF" wi="418" he="64" />其中<img file="FDA00009414886700000113.GIF" wi="477" he="54" />表示第<img file="FDA00009414886700000114.GIF" wi="19" he="41" />条径冲激响应,如果用<img file="FDA00009414886700000115.GIF" wi="97" he="59" />表示信道第<img file="FDA00009414886700000116.GIF" wi="23" he="37" />条径对应的BEM系数,b<sub>q</sub>表示基函数,则有:<maths num="0001"><math><![CDATA[<mrow><msub><mi>h</mi><mi>l</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><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><img file="FDA0000941488670000014.GIF" wi="998" he="123" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>b</mi><mi>q</mi></msub><mo>=</mo><mo>&lsqb;</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><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>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000941488670000015.GIF" wi="1170" he="143" /></maths><img file="FDA0000941488670000016.GIF" wi="995" he="66" />其中,<img file="FDA0000941488670000017.GIF" wi="57" he="65" />表示向上取整,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="0003"><math><![CDATA[<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><mrow><mi>I</mi><mi>S</mi><mi>I</mi></mrow><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>A</mi><mi>d</mi><mi>d</mi></mrow><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000941488670000018.GIF" wi="1035" he="62" /></maths>其中,<img file="FDA0000941488670000019.GIF" wi="47" he="57" />表示OFDM系统中循环前缀的长度大于信道的最大多径时延时接收端接收到第i个OFDM符号,<img file="FDA00009414886700000110.GIF" wi="68" he="59" />代表前一个OFDM发送符号x<sup>i‑1</sup>对当前OFDM接收符号y<sup>i</sup>的干扰,<img file="FDA00009414886700000111.GIF" wi="81" he="59" />代表循环前缀的重构部分;将式子(4)变换到频域,则:<maths num="0004"><math><![CDATA[<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><mrow><mi>I</mi><mi>S</mi><mi>I</mi></mrow><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>Y</mi><mrow><mi>A</mi><mi>d</mi><mi>d</mi></mrow><mi>i</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000941488670000021.GIF" wi="1061" he="73" /></maths><maths num="0005"><math><![CDATA[<mrow><msubsup><mi>Y</mi><mrow><mi>I</mi><mi>S</mi><mi>I</mi></mrow><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><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><img file="FDA0000941488670000022.GIF" wi="1081" he="85" /></maths><maths num="0006"><math><![CDATA[<mrow><msubsup><mi>Y</mi><mrow><mi>A</mi><mi>d</mi><mi>d</mi></mrow><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><img file="FDA0000941488670000023.GIF" wi="1074" he="91" /></maths>步骤三:初始化;设置迭代次数I=0,令<img file="FDA0000941488670000024.GIF" wi="300" he="69" />这样<img file="FDA0000941488670000025.GIF" wi="170" he="64" />步骤四:估计相应的信道基扩展模型系数;接收端信号Y中导频位置观测量Y<sup>p</sup>可以表示为:<maths num="0007"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msup><mi>Y</mi><mi>p</mi></msup><mo>=</mo><msup><mover><mi>Y</mi><mo>~</mo></mover><mi>p</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msup><mi>D</mi><mi>p</mi></msup><msup><mi>S</mi><mi>p</mi></msup><msub><mi>h</mi><mrow><mi>e</mi><mi>q</mi><mi>v</mi></mrow></msub><mo>+</mo><msup><mi>D</mi><mi>d</mi></msup><msup><mi>S</mi><mi>d</mi></msup><msub><mi>h</mi><mrow><mi>e</mi><mi>q</mi><mi>v</mi></mrow></msub><mo>+</mo><msup><mi>W</mi><mi>p</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>Ph</mi><mrow><mi>e</mi><mi>q</mi><mi>v</mi></mrow></msub><mo>+</mo><msup><mi>D</mi><mi>d</mi></msup><msup><mi>S</mi><mi>d</mi></msup><msub><mi>h</mi><mrow><mi>e</mi><mi>q</mi><mi>v</mi></mrow></msub><mo>+</mo><msup><mi>W</mi><mi>p</mi></msup></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000941488670000026.GIF" wi="1135" he="209" /></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="0008"><math><![CDATA[<mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>e</mi><mi>q</mi><mi>v</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>P</mi><mi>H</mi></msup><mi>P</mi><mo>+</mo><mi>&alpha;</mi><mi>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><img file="FDA0000941488670000027.GIF" wi="1079" he="81" /></maths>考虑对矩阵P求逆运算的准确性,加上一个小的扰动,其中α为扰动因子;将估计得到的BEM系数<img file="FDA0000941488670000028.GIF" wi="68" he="74" />代入等式(1),即可求出信道时域冲激响应<img file="FDA0000941488670000029.GIF" wi="52" he="73" />步骤五:信道均衡与信号检测;利用估计得到的信道时域冲激响应<img file="FDA00009414886700000210.GIF" wi="34" he="57" />构造矩阵<img file="FDA00009414886700000211.GIF" wi="82" he="60" />求出信道频率响应矩阵<img file="FDA00009414886700000212.GIF" wi="81" he="68" />利用MMSE准则构造第i帧OFDM数据单元的频域均衡器有<img file="FDA00009414886700000213.GIF" wi="516" he="112" />其中,σ<sup>2</sup>是噪声方差,I<sub>N</sub>是N×N的单位矩阵,然后检测得到发端当前频域符号矢量<img file="FDA00009414886700000214.GIF" wi="75" he="64" />步骤六:判断;令:I=I+1                 (10)如果I≤M,M为取值范围为自然数的迭代次数,执行步骤七,否则,转到步骤八;步骤七:迭代;利用估计得到的信道时域冲激响应<img file="FDA0000941488670000031.GIF" wi="24" he="46" />构造矩阵<img file="FDA0000941488670000032.GIF" wi="80" he="59" />再加上前一时刻估计得到的频域符号矢量<img file="FDA0000941488670000033.GIF" wi="89" he="62" />和当前频域符号矢量<img file="FDA0000941488670000034.GIF" wi="58" he="62" />时域代入式子(6)(7),求出<img file="FDA0000941488670000035.GIF" wi="72" he="75" />和<img file="FDA0000941488670000036.GIF" wi="81" he="75" />后更新等式(5),然后转到步骤四;步骤八:检测得到发端当前频域符号矢量<img file="FDA0000941488670000037.GIF" wi="86" he="64" />计算最后的估计量<img file="FDA0000941488670000038.GIF" wi="145" he="75" />
地址 611731 四川省成都市高新区(西区)西源大道2006号