发明名称 一种OFDM信号的多级时域盲均衡方法
摘要 本发明方法涉及一种OFDM信号的多级时域盲均衡方法,包括以下步骤:首先对接收序列利用相关峰搜索方法确定接收信号的帧格式;然后利用OFDM信号的帧同步符号信息完成对接收信号的第一级时域均衡;接下来利用OFDM信号的循环前缀序列完成对接收信号的第二级时域均衡,经该级时域均衡后的信号即是输出信号,可用于信号的后续处理。该方法在确定均衡器系数的过程中充分利用了OFDM信号帧格式中固有的重复性,在没有接收信号的先验信息的条件下,实现了信号的盲均衡,均衡效果良好;并且该方法对于有很多子载波的OFDM信号,具有计算简单、处理迅速的特点。
申请公布号 CN103957177B 申请公布日期 2017.04.05
申请号 CN201410214566.4 申请日期 2014.05.21
申请人 北京理工大学 发明人 职如昕;杨杰;袁莹莹;张序绮
分类号 H04L25/03(2006.01)I;H04L27/26(2006.01)I 主分类号 H04L25/03(2006.01)I
代理机构 代理人
主权项 一种OFDM信号的多级时域盲均衡方法,其特征在于,包括以下步骤:步骤1:对接收序列利用相关峰搜索方法确定接收信号的帧格式,过程如下:在接收端对接收信号利用相关峰搜索方法根据重复信号能量差最小的准则,统计获得信号的帧长L<sub>frm</sub>,OFDM符号的长度L<sub>S</sub>,每个符号内循环前缀长度L<sub>CP</sub>,并确定FFT的长度L<sub>fft</sub>=L<sub>S</sub>‑L<sub>CP</sub>;步骤2:在步骤1中确定了OFDM信号的帧格式后,利用OFDM信号的帧同步符号信息完成对接收信号的第一级时域均衡,均衡过程包括以下步骤:(1)第一级时域均衡器采用的级数为L<sub>1</sub>,初始化该级时域均衡器系数:第一阶系数为1,其余系数为0,即w<sub>1</sub>=[w<sub>1</sub>(1),w<sub>1</sub>(2),w<sub>1</sub>(3),…,w<sub>1</sub>(L<sub>1</sub>)]<sup>T</sup>=[1,0,0,…,0]<sup>T</sup>;(2)计算迭代误差,此误差由相邻两个帧同步符号经过矩阵重组之后,与均衡器权值向量进行相乘运算得到的结果的差值来表示,形式化描述如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>e</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>w</mi><mn>1</mn><mi>H</mi></msubsup><msub><mi>v</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>w</mi><mn>1</mn><mi>H</mi></msubsup><msub><mi>v</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000011.GIF" wi="750" he="79" /></maths>其中,e<sub>1</sub>(i)为当前的误差信号,v<sub>n</sub>(j)=[u<sub>n</sub>(j+L<sub>1</sub>‑1),u<sub>n</sub>(j+L<sub>1</sub>‑2),…,u<sub>n</sub>(j)]<sup>T</sup>,u<sub>n</sub>(j)表示第n帧信号数据的帧同步符号的第j个采样数据,u<sub>n</sub>为接收信号的第n帧信号数据的帧同步符号,v<sub>n+1</sub>(j)=[u<sub>n+1</sub>(j+L<sub>1</sub>‑1),u<sub>n+1</sub>(j+L<sub>1</sub>‑2),…,u<sub>n+1</sub>(j)]<sup>T</sup>,u<sub>n+1</sub>(j)表示第n+1帧信号数据的帧同步符号的第j个采样数据,u<sub>n+1</sub>为第n+1帧信号数据的帧同步符号,j=1,2,…,L<sub>S</sub>‑L<sub>1</sub>+1;(3)根据迭代误差和均衡器的当前权值向量按照相应的关系计算新的均衡器权值向量,作为下一次均衡过程所用的权值向量,新的均衡器权值向量w′<sub>1</sub>的计算公式如下:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>w</mi><mn>1</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><mo>-</mo><mn>2</mn><msub><mi>&mu;</mi><mn>1</mn></msub><msubsup><mi>e</mi><mn>1</mn><mo>*</mo></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>v</mi><mi>n</mi></msub><mo>(</mo><mi>j</mi><mo>)</mo><mo>-</mo><msub><mi>v</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>(</mo><mi>j</mi><mo>)</mo><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000012.GIF" wi="934" he="71" /></maths>其中w<sub>1</sub>为当前均衡器权值向量,e<sub>1</sub>(i)为步骤(2)中的迭代误差,*表示共轭操作,μ<sub>1</sub>为均衡器调整步长;(4)重复执行上述步骤(2)和(3)共计L<sub>S</sub>‑L<sub>1</sub>+1次,得到第一级时域均衡器权值向量w<sub>1</sub>;(5)根据均衡器权值及均衡器的输入信号确定此级均衡器的输出信号,输出信号的计算过程以及后续的处理如下:将均衡器权值向量与对应的时域信号相乘,得到经过第一级时域均衡的信号,即第一级时域均衡器输出信号b<sub>n</sub>,b<sub>n</sub>的计算公式如下:<maths num="0003"><math><![CDATA[<mrow><msub><mi>b</mi><mi>n</mi></msub><mo>=</mo><msubsup><mi>w</mi><mn>1</mn><mi>H</mi></msubsup><msub><mi>a</mi><mi>n</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000021.GIF" wi="308" he="78" /></maths>其中,w<sub>1</sub>为经过步骤(4)得到的第一级时域均衡器权值向量,上标H表示共轭转置,a<sub>n</sub>为当前处理的帧同步信号所在的OFDM帧信号序列经过矩阵重组的序列矩阵,矩阵中每一列元素组成如下a<sub>n</sub>(i)=[x(i),x(i‑1),x(i‑2),…,x(i‑L<sub>1</sub>+1)]<sup>T</sup>,x(i)为接收信号的第i个采样点;当对下一帧信号进行均衡处理时,若检测到后续还有其他帧数据,则均衡过程从步骤(1)重新开始;若是最后一帧数据,则均衡器权值向量不进行更新,使用上一帧数据得到的均衡器权值向量完成对最后一帧数据的均衡;步骤3:对经过步骤2均衡后的信号,利用OFDM信号的循环前缀序列完成对接收信号的第二级时域均衡,经该级时域均衡后的信号即是输出信号,可用于信号的后续处理,均衡过程包括以下步骤:(1)第二级时域均衡器采用的级数为L<sub>2</sub>,初始化该级时域均衡器系数:第一阶系数为1,其余系数为0,即w<sub>2</sub>=[w<sub>2</sub>(1),w<sub>2</sub>(2),w<sub>2</sub>(3),…,w<sub>2</sub>(L<sub>2</sub>)]<sup>T</sup>=[1,0,0,…,0]<sup>T</sup>;(2)计算迭代误差,此误差由一个OFDM符号内循环前缀和用来复制得到循环前缀的序列经过矩阵重组之后,与均衡器权值向量进行相乘运算得到的结果的差值来表示,形式化描述如下:<maths num="0004"><math><![CDATA[<mrow><msub><mi>e</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>w</mi><mn>2</mn><mi>H</mi></msubsup><msub><mi>r</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>w</mi><mn>2</mn><mi>H</mi></msubsup><msub><mi>r</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>f</mi><mi>f</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000031.GIF" wi="837" he="79" /></maths>其中,e<sub>2</sub>(i)为当前的误差信号,r<sub>n</sub>(j)=[y<sub>n</sub>(j+L<sub>2</sub>‑1),y<sub>n</sub>(j+L<sub>2</sub>‑2),…,y<sub>n</sub>(j)]<sup>T</sup>,y<sub>n</sub>(j)表示第n个OFDM符号的循环前缀的第j个采样数据,j=1,2,…,L<sub>CP</sub>‑L<sub>2</sub>+1;(3)根据迭代误差和均衡器的当前权值向量按照相应的关系计算新的均衡器权值向量,作为下一次均衡过程所用的权值向量,新的均衡器权值向量w′<sub>2</sub>的计算公式如下:<maths num="0005"><math><![CDATA[<mrow><msubsup><mi>w</mi><mn>2</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>w</mi><mn>2</mn></msub><mo>-</mo><mn>2</mn><msub><mi>&mu;</mi><mn>2</mn></msub><msubsup><mi>e</mi><mn>2</mn><mo>*</mo></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>r</mi><mi>n</mi></msub><mo>(</mo><mi>j</mi><mo>)</mo><mo>-</mo><msub><mi>r</mi><mi>n</mi></msub><mo>(</mo><mrow><mi>j</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>f</mi><mi>f</mi><mi>t</mi></mrow></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000032.GIF" wi="1022" he="71" /></maths>其中w<sub>2</sub>为当前均衡器权值向量,e<sub>2</sub>(i)为步骤(2)中的迭代误差,μ<sub>2</sub>为均衡器调整步长;(4)重复执行上述步骤(2)和(3)共计L<sub>CP</sub>‑L<sub>2</sub>+1次,得到第二级时域均衡器权值向量w<sub>2</sub>;(5)根据均衡器权值及均衡器的输入信号确定此级均衡器的输出信号,完成对接收信号的时域均衡;输出信号的计算过程如下:将均衡器权值向量与对应的输入时域信号相乘,得到经过第二级时域均衡的信号,即第二级时域均衡器输出信号d<sub>n</sub>,d<sub>n</sub>的计算公式如下:<maths num="0006"><math><![CDATA[<mrow><msub><mi>d</mi><mi>n</mi></msub><mo>=</mo><msubsup><mi>w</mi><mn>2</mn><mi>H</mi></msubsup><msub><mi>c</mi><mi>n</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001113471830000033.GIF" wi="302" he="78" /></maths>其中,w<sub>2</sub>为第二级时域均衡器权值向量,c<sub>n</sub>为当前处理的OFDM符号经过矩阵重组的序列矩阵,矩阵中每一列元素组成如下c<sub>n</sub>(i)=[b(i),b(i‑1),b(i‑2),…,b(i‑L<sub>2</sub>+1)]<sup>T</sup>,b(i)为从第一级时域均衡器输出到第二级时域均衡器信号的第i个采样点;当对下一个OFDM符号进行均衡处理时,从步骤(1)重新开始均衡过程。
地址 100081 北京市海淀区中关村南大街5号