一种基于粒子滤波器和神经网络的双选择性信道跟踪方法

专利名称：一种基于粒子滤波器和神经网络的双选择性信道跟踪方法
技术领域：
本发明涉及多输入多输出正交频分复用(MIMO-OFDM)系统中的信道跟踪方法，确切 的说是MIMO-OFDM系统中基于粒子滤波器和神经网络的时间-频率双选择性信道跟踪方 法，特别涉及一种基于粒子滤波器和神经网络的双选择性信道跟踪方法。
背景技术：
众所周知，MIMO和OFDM技术将是下一代移动通信的关键技术。两者的结合充分利 用了空间分集、频率分集和时间分集来改善系统的容量和抗噪声性能。然而MIMO-OFDM的 高效性能依赖于接收端的信道估计的准确与否。传统的信道估计的方法为基于导频的信道估 计算法，但是为了获得可信的信道估计，信道带宽的相当一部分被训练序列所占用。1996 年，Stuber在他的著作"Principles of mobile communication"中指出，在快衰落、高多普勒 条件下，带宽被占用的比例可高达50%。2002年，Liu等在文章"Space time coding and Kalman filtering for time-selective fading channel"中提出利用卡尔曼滤波来进行信道估计和跟踪，以 此来减少训练信号的长度。对于时变信道的跟踪，在线性高斯系统中，没有一种算法能优于 卡尔曼滤波。但在无线通信系统中，接收信号经常被非高斯噪声所污染。在涉及到非线性状 态转移与非高斯噪声的信道跟踪问题上，序贯蒙特卡罗滤波比扩展卡尔曼滤波更优越。
粒子滤波算法是序贯蒙特卡罗滤波的一种基础方法，其主要思想是利用一组相关权值的 随机变量样本计算完成估计，以逼近实际的后验概率密度，当样本数目很大时，这种概率估 计可近拟于真实后验概率。
设<formula>formula see original document page 3</formula>表示随机粒子，其含有后验概率<formula>formula see original document page 3</formula>特性，其中<formula>formula see original document page 3</formula>是与权值<formula>formula see original document page 3</formula> 相关的状态集，且权值满足<formula>formula see original document page 3</formula>，贝lh时刻的
后验概率密度可近似表示<formula>formula see original document page 3</formula>.
对于粒子滤波跟踪算法，跟踪的初值是影响跟踪性能的主要因素之一。在用粒子滤波算 法进行跟踪时，大多假定信道的初值已知，或者采用训练序列的方法得到信道的初值。本发 明中我们采用径向基神经网络对信道进行学习训练以得到粒子滤波跟踪的精确初值。
径向基神经网络是一种三层前向局部网络输入层由信号源节点构成，第二层为隐含层， 第三层为输出层，它对输入模式的作用做出响应。从输入空间到隐含空间的变换是非线性的，
而隐含空间到输出空间的变换是线性的。隐节点的激活函数是径向基函数(RBF)，它是一
种关于中心点对称的非负非线性函数，通常我们采用高斯函数作为径向基函数。

发明内容
为在提高频谱利用率的前提下，使系统得到良好的误码率性能，本发明提供一种
MMO-OFDM双选择性信道的跟踪方法。为提供该跟踪方法，首先给出MIMO-OFDM双选
择信道的信道模型。
设MIMO-OFDM系统中发送天线Z和接收天线y'间的双选择信道的冲激响应设为 h,/","，其中/ = 0，l，...,Mr， / = 0,l，...，Mfl， Mr和A^分别为发送天线和接收天线的数 目。则参考Proakis在"Digital Communication (Forth Edition)"所述的抽头延迟线模型 <formula>formula see original document page 4</formula> (1 ) 其中丄代表多径数，A,和r。分别为各径的路径增益和延迟扩展，w,J")为复噪声随机 序列。单径信道h。(")建模为一阶自回归AR过程，即 <formula>formula see original document page 4</formula> (2) 参考1996年Tsatasnis等发表的"Estimation and equalization of fading channels with random coefficients "将系数a确定为 <formula>formula see original document page 4</formula> (3) 其中， /。(.)是零阶Bessel函数，^为最大多普勒频移，/。是载波频偏，7;为符号时间 一种基于粒子滤波器和神经网络的MIMO-OFDM双选择性信道跟踪方法，步骤如下 1) 利用径向基神经网络对双选信道进行训练得到粒子滤波跟踪的初始值<formula>formula see original document page 4</formula>似，权值w[(0)均置l/M; 2) 得到一个观测值后，利用公式(1) (2)可由h^("-l)得到新的粒子集h^(n); 3) 对于每一个粒子，应用似然函数计算粒子权值， (<formula>formula see original document page 4</formula> (4) 4) 计算归一化权值W^("); 5) 重采样具有大权值的样本，重复采样；同时小权值样本，尽量抛弃，如此获得M 个随机样本，其权值设为w;:,(")-i/M，样本分布渐趋于后验分布， <formula>formula see original document page 4</formula>(5) 6) 由后验分布得到信道估计值为 <formula>formula see original document page 4</formula>( 6) 本发明方法的优点可以对无线通信中的时间-频率双选择性信道进行良好的信道跟踪， 跟踪过程不再需要插入训练序列或者导频符号， 可提高系统的频谱利用率， 得到高的误码率性能。


图l为本发明方法采用的基于空时分组编码的两发两收MIMO-OFDM系统方框其中，l.STBC编码器，2.STBC译码器，3.神经网络结合粒子滤波信道跟踪模块。
图2是两发两收MIMO-OFDM系统空时编码过程示意其中4.频率，5.时间，6.空时编码。
图3是径向基神经网络训练参数表；
图4是粒子滤波跟踪算法参数表；
图5是粒子滤波算法在不同跟踪初值时的系统BER性能比较。
具体实施方式
实施例
以发射天线数为2，接受天线数为2的MIMO-OFDM系统为例来介绍这种算法。系统 模型如图1所示信息序列经过调制器调制，然后把这些调制符号通过空时编码器1进行 空时编码(STBC)，经串并变换后经过快速傅立叶反变换(IFFT)，最后把码字通过发射天 线发射出去。在接收端，首先对接收信号进行快速傅立叶变换(FFT)，然后经过空时译码 器2译码和信道跟踪模块3，并最终译出信息序列。
图l所示系统中空时编码的过程如下设OFDM的子载波数为A^，输入的连续信息比
特流经过调制后，以A^为单位进行分组。在每两个连续OFDM符号间隔的开始，将两个连 续的分组经串/并变换后得到符号向量X(") = [X(0,"),X(1，")，...，X(tV, -1，")]和
X(" + l) = [X(0," + l),X(l," + lV..,X(iV, -l,"十l)]。将符号向量对应的分量进行Alamouti空
时分组编码，如图2所示，其中纵向4表示频率，横向5表示时间，经空时编码6编码后得 到4个向量X。(w)、 Xj(")、 X。(" + l)、 X,(" + l)，即4个OFDM符号，其中n代表第"个
OFDM符号间隔，而下标代表天线索引。在第"个OFDM符号间隔内，编码结果X。(")，
X,(w)分别经过IFFT变换后由第1和第2个天线同时发射出去，在第n + l个OFDM符号间
隔内，X。(n + 1)、 X,(n + 1)分别经过IFFT变换后由第1和第2个天线同时发射出去。
假设信道在一个OFDM符号间隔内不变，则接收端经过解调(FFT变换)后，每个子载 波上的接收信号是两个畸变发送信号的迭加，可表示为
<formula>formula see original document page 5</formula>
其中X,("A)， Y,(",yfc)分别是第"个OFDM符号的第A个子载波的在第/个发送天线 和第y个接收天线上的的发送信号和接收信号；W/n,A:)是第j'个接收天线在第"个OFDM 符号期间第A:个子载波上的随机噪声。Hy("，"是从第/个发射天线到第j'个接收天线在第"
个OFDM符号期间第A:个子载波上的信道频率响应系数，
<formula>formula see original document page 6</formula>
其中11,/("，0 / = 0，1; j'-O,l为两发两收MIMO-OFDM系统中双选择信道的冲激响应。
对于该信道的跟踪算法步骤如下
1) 利用径向基神经网络对双选信道进行训练得到粒子滤波跟踪的初始值<formula>formula see original document page 6</formula>权值(0)均置1 /M 。
2) 得到一个观测值后，利用公式(i) (2)可由^("-1)得到新的粒子集h^(")。
3) 对于每一个粒子，应用似然函数计算粒子权值，
<formula>formula see original document page 6</formula>
4) 计算归一化权值w^(");
5) 重采样具有大权值的样本，重复采样；同时小权值样本，尽量抛弃。如此获得 M个随机样本，其权值设为W:.("):1/M 。样本分布渐趋于后验分布，
<formula>formula see original document page 6</formula>
6) 由后验分布得到信道估计值为
<formula>formula see original document page 6</formula>
图3给出了径向基神经网络对双选信道进行学习训练时的参数列表G04丄为网络输出
值与目标值的均方误差，CT,为径向基函数的宽度，C为隐层神经元的个数。参数的取得根据
经验值得到。图4给出了粒子滤波算法对MIM0-0FDM双选择性信道进行跟踪时的参数列 表。在采用粒子滤波算法进行信道跟踪时，粒子数的多少会影响跟踪的性能，本发明中采用 粒子数^ = 50对信道进行跟踪。图5给出了粒子滤波算法在不同跟踪初值时的系统误码率性 能比较。从图中可以看出用径向基神经网络获得跟踪初值的系统误码率性能与跟踪初值已知 时的误码率性能非常接近。证明了径向基神经网络结合粒子滤波算法的高效性。
权利要求
1.一种基于粒子滤波器和神经网络的双选择性信道跟踪方法，步骤如下1)利用径向基神经网络对双选信道进行训练得到粒子滤波跟踪的初始值hi，jm(0)，i＝0，1，...，MT；j＝0，1，...，MR；m＝1，...，M，权值wi，jm(0)均置1/M；2)得到一个观测值后，利用公式<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>n</mi> <mo>,</mo> <mi>&tau;</mi> <mo>)</mo></mrow><mo>=</mo><msubsup> <mi>&Sigma;</mi> <mrow><mi>l</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>L</mi><mo>-</mo><mn>1</mn> </mrow></msubsup><msub> <mi>A</mi> <mi>l</mi></msub><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>l</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <msubsup><mi>&tau;</mi><mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi></mrow><mi>l</mi> </msubsup> <mo>)</mo></mrow><mo>+</mo><msub> <mi>w</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow> </mrow>]]></math> id="icf0001" file="A2008100167480002C1.tif" wi="67" he="6" top= "60" left = "31" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>i＝0，1，...，MT；j＝0，1，...，MR和公式<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>l</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>=</mo><mi>&alpha;</mi><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>l</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow> </mrow>]]></math> id="icf0002" file="A2008100167480002C2.tif" wi="34" he="5" top= "72" left = "31" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>l＝0，1，...，L-1可由hi，jm(n-1)得到新的粒子集hi，jm(n)；3)对于每一个粒子，应用似然函数计算粒子权值，<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><msubsup> <mi>lik</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow> <mo>(</mo> <mi>y</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>|</mo> <msubsup><mi>h</mi><mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi></mrow><mi>m</mi> </msubsup> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>1</mn> <mrow><msub> <mi>&sigma;</mi> <mi>n</mi></msub><msqrt> <mn>2</mn> <mi>&pi;</mi></msqrt> </mrow></mfrac><msup> <mi>e</mi> <mrow><mo>-</mo><mfrac> <mrow><mo>(</mo><mi>y</mi><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>-</mo><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mover> <mi>x</mi> <mo>^</mo></mover><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mrow><mn>2</mn><msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn></msubsup> </mrow></mfrac> </mrow></msup> </mrow>]]></math></maths>4)计算归一化权值wi，jm(n)；5)重采样具有大权值的样本，重复采样；同时小权值样本，尽量抛弃，如此获得M个随机样本，其权值设为<maths id="math0004" num="0004" ><math><![CDATA[ <mrow><msubsup> <mover><mi>w</mi><mo>~</mo> </mover> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>/</mo><mi>M</mi><mo>,</mo> </mrow>]]></math> id="icf0004" file="A2008100167480002C4.tif" wi="26" he="4" top= "123" left = "79" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>样本分布渐趋于后验分布，<maths id="math0005" num="0005" ><math><![CDATA[ <mrow><mi>p</mi><mrow> <mo>(</mo> <msubsup><mi>h</mi><mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi></mrow><mi>m</mi> </msubsup> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>|</mo> <mi>y</mi> <mrow><mo>(</mo><mn>1</mn><mo>:</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>&ap;</mo><msubsup> <mover><mi>w</mi><mo>~</mo> </mover> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><msubsup> <mi>&Sigma;</mi> <mrow><mi>m</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></msubsup><msub> <mi>&delta;</mi> <mrow><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow> </mrow></msub><mrow> <mo>(</mo> <mi>dh</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>1</mn> <mi>M</mi></mfrac><msubsup> <mi>&Sigma;</mi> <mrow><mi>m</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></msubsup><msub> <mi>&delta;</mi> <mrow><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow> </mrow></msub><mrow> <mo>(</mo> <mi>dh</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow>]]></math></maths>6)由后验分布得到信道估计值为<maths id="math0006" num="0006" ><math><![CDATA[ <mrow><mover> <mi>h</mi> <mo>^</mo></mover><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>=</mo><mi>E</mi><mrow> <mo>(</mo> <mi>h</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mo>&Integral;</mo><mi>h</mi><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mi>p</mi><mrow> <mo>(</mo> <mi>h</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow> <mo>|</mo> <mi>y</mi> <mrow><mo>(</mo><mn>1</mn><mo>:</mo><mi>n</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mi>dh</mi><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>&ap;</mo><mfrac> <mn>1</mn> <mi>M</mi></mfrac><msubsup> <mi>&Sigma;</mi> <mrow><mi>m</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></msubsup><msubsup> <mi>h</mi> <mrow><mi>i</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo></mrow><mo>.</mo> </mrow>]]></math></maths>
全文摘要
本发明涉及一种基于粒子滤波器和神经网络的双选择性信道跟踪方法，其特征在于，采用基于粒子滤波和神经网络的方法对无线通信中的时间-频率双选择性信道进行良好的信道跟踪，跟踪过程不再需要插入训练序列或者导频符号，可提高系统的频谱利用率，得到高的误码率性能。本发明具有重要的实际应用价值。
文档编号H04L1/06GK101340266SQ200810016748
公开日2009年1月7日 申请日期2008年6月6日 优先权日2008年6月6日
发明者波 刘, 刘彦君, 刘玉龙, 丽 彭, 李海燕, 江铭炎, 程永明 申请人:山东大学