RBF算法

  • 格式:doc
  • 大小:241.00 KB
  • 文档页数:4

1引言

作为一种单隐层前馈神经网络,径向基函数(RBF)网络已经成功地应用于模式识别、函数逼近信号处理、系统建模和控制等领域【】。RBF网络的广泛应用,是与其具有的网络结构简单、非线性逼近能力强、收敛速度快以及全局收敛等优点密不可分的[3]。

对于RBF神经网络的学习算法,关键问题是隐层神经元中心参数的合理确定。在已有的常用学习算法中,中心参数要么是从给定的训练集样本中按照某种方法直接选取,要么采用聚类的方法进行确定。实际应用表明,这些学习算法均有不足之处,使之应用范围受到限制。如正交优选法【】,其缺点是隐层神经元中心的取值是训练集样本中的数据,这在多数情况下难以反映系统的真正映射关系,且在中心点的优选过程中会出现病态现象,导致训练失败。再如Moody和Darken算法[],其缺点是无法合理地确定隐层神经元的数目,所得到的中心值也未必是合理的。

本文拟采用最近邻聚类和可变速率的最小均方(LMS)算法相结合的方法来给出RBF网络的学习算法。其中用最近邻聚类来确定径向基函数的中心,

2 RBF神经网络结构

最基本的RBF神经网络的构成包括三层,分别为输入层、隐层(中间层)和输出层。其中输入层由一些源点(感知单元)组成,它们将网络与外部环境连接起来,仅起到数据信息的传递作用,对输入信息不进行任何变换;隐层神经元的核函数(或称作用函数)取为径向基函数,对输入信息到隐层空间之间进行非线性变换,通常具有较高的维数;输出层是线性的,为输入层的激活模式提供响应。

设隐层、输出层上的神经元数分别为M,Q,输入模式记为X,12[,,,]TRXxxx,输出记为Y,12[,,,]TQYyyy。本文取径向基函数为Gauss函数,隐单元输出则为

2expjjjXCz 1,2,,jM 1

式中:jz为隐层第j个神经元的输出值;jC为隐层第j个神经元的中心,由隐层第j个神经元对应于输入层所有神经元的中心分量构成,12,,,TjjjjRCccc;j为隐层第j个神经元的宽度,与jC相对应;为欧氏范数。

输出层神经元的输入输出关系表达式是:

1Mkkjjjywz 1,2,,kQ 2

式中:ky为输出层第j个神经元的输出值;kjw为输出层第k个神经元与隐层第j个神经元间的权值。

RBF神经网络的参数在此主要是指网络的中心、宽度、和调节权重。

图1 RBF网络的拓朴结构

3学习算法

RBF网络设计是在于网络的学习算法设计。RBF网络的学习算法一般包括两个不同的阶段:一是隐层径向基函数中心的确定阶段;二是径向基函数权值学习调整阶段。因此,RBF网络所用的非线性函数形式并不影响网络的性能,关键在于基中心、宽度和权值调整算法的确定

3.1、网络设计步聚

神经网络设计大致需要经过四个步骤:

(1)网络的生成。主要是根据输入样本数据和应用要求来确定网络的结构、层数、输入层、输出层和隐层的神经元的数目、传递函数及训练算法。

(2)网络的初始化。在网络生成的同时需根据不同的要求对网络各层的权值和阀值进行初始化。

(3)网络的训练。根据提供的样本数据对“输入矢量——输出矢量”和训练算法对网络进行训练。

(4)网络的仿真。根据测试样本数据对已训练好的网络进行仿真计算,也可对训练后的网络实际输出与仿真输出进行误差比较。

3.2 基于最近邻聚类的中心选取

给定高斯函数宽度,作为所分配隐单元的初始宽度,并把其作为衡量输入数据与已存在隐单元中心距离的标准,的大小决定了网络的最终聚类数目。

设RBF网络已存在个M隐单元,对于一个新的样本,如果iXC成立,1,2,,iM,则生成一个新的隐单元。否则,将该样本归到最小距离所属的类中,中心保持不变。

具体步骤如下:

(1)选择一个适当的高斯函数宽度,定义一个矢量Al用于存放属于各类的输出向量之和,定义一个计数器Bl用于统计属于各类的样本个数,其中l为类别数;

(2)从第一个样本111121,,,RXxxx开始,建立第一个聚类中心,令11CX,11B。这样建立的RBF网络,只有一个隐单元,该单元的中心是1C,该隐单元到输出层的权为11/1WAB;

(3)考虑第2个样本数据221222,,,RXxxx,求出2X到1C这个聚类中心的距离21XC。如果21XC,则1C为2X的最近邻聚类,且令12B;如果21XC,则将2X作为一个新的聚类中心,并令22CX,21B。在上述建立的RBF网络中再添加一个隐单元,该隐单元到输出层的权为22/2WAB;

(4)假设我们考虑第k个样本时,已存在M个聚类中心,其中心分别为12,,,MCCC,则在上述建立的RBF网络中已有M个隐单元。这时再分别求出kX到这M个聚类中心的距离kiXC,1,2,,iM。

设kjXC为这些距离中的最小距离,即jC为kX的最近邻聚类:此时,如果kjXC,则将kX作为一个新的聚类中心,令1MkCX,1MM,1kAMY,1BM,并保持Bi的值不变,1,2,,iM,在上述建立的RBF网络中再添加第1M个隐单元,该隐单元到输出层的权矢量为11/1MWAMBM;如果kjXC,作如下计算:kAjAjY,1BjBj。当ij时,1,2,,iM,保持Ai、Bi的值不变。隐单元到输出层的权/iWAiBi。

宽度的大小决定动态自适应RBF网络的复杂程度。越小,所得到的聚类数目就越多,计算量就越大,精度也越高;越大,所得到的聚类数目就越少,计算量就相对较小,从而导致精度也较低。

3.3权值更新如下 对所有输入样本,定义代价函数瞬时值:2112Njje

其中N是用于学习的训练样本数目,je是误差信号,定义如下:

2*1expMjijjjjiiiXCedFXdw则

权值更新如下:

21expNjijjiiXCnenwn

1002111expiiijjjjjjijinwnwnwnwnwnnzndnwnznnnXCzn

中心更新:

在LMS算法中,学习率参数在计算过程中保持不变,从而导致收敛速度较慢,并且对输入结构特征结构的变化反应较灵敏。LMS算法一般需要输入空间维数十倍的迭代次数才能达到稳定状态,当输入空间维数较高时缓慢的收敛速度会变得特别严重。因此我们提出了可变速率的学习:01nn, 0和是用户选择的常数。在自适应的早期阶段,即迭代次数相对搜寻时间常数较小时,学习率参数近似等于,算法运行实际上也是与“标准”的LMS算法一样的。因此,通过在允许范围内选择一个较大,我们希望对可调权值找到一组找到一组较好的值并在其中上下浮动。然后,当迭代次数比搜寻时间常数大时,学习率参数近似为。此时,该算法以一个传统的随机逼近算法运行,且权值收敛到它们的最优值。这样搜寻后的收敛进度是把标准LMS算法的期望特征和传统随机逼近理论结合起来。

1jjjjjwnwnnzndnwnzn

01nn

2expjijiXCzn

其中0和是用户选择的常数