复杂重尾噪声影响下多变量Hammerstein模型的辨识方法技术

针对工业过程中广泛存在的重尾噪声影响下的非线性系统的辨识，目前还没有被过多的研究。本发明专利技术公开一种复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，利用高斯混合分布来改进布谷鸟算法(CS)，所形成新算法为GMDA，RBF神经网络可以近似复杂非线性系统，以及所改进布谷鸟算法GMDA强大的全局搜索能力。利用GMDA训练RBF神经网络的结构，以此来近似Hammerstein模型的静态非线性部分，同时利用GMDA算法来辨识线性部分的参数值。提供了一种解决复杂重尾噪声影响下多变量Hammerstein模型辨识问题的有效途径。

【技术实现步骤摘要】
复杂重尾噪声影响下多变量Hammerstein模型的辨识方法
本专利技术公开了一种复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，属于工业过程模型辨识领域。
技术介绍
非线性系统在工业过程中是广泛存在的，由于非线性本身特有的复杂性，导致对于非线性系统辨识研究方面有很大的不足，因此对于以系统辨识为基础的现代控制论来说，有很大的发展限制。Hammerstein模型是一种由一个静态非线性模块和一个动态线性模块串联构成的典型的模块化非线性系统模型，由于其结构简单，且能表征一大类实际非线性系统，如PH中和反应、挖掘臂伺服系统、增压锅炉等工程应用领域，引起了很多学者的关注。然而对于Hammerstein模型辨识方面的研究，目前多数集中在单变量系统，且噪声模型为白噪声或由白噪声产生的有色噪声的研究上。对于多变量Hammerstein模型，以及近年来普遍报道的重尾噪声(即一类包含有离群点的噪声模型)影响下的模型研究，还比较少见。尤其是目前多数的方法所对应的系统均是在白噪声假设下的系统辨识问题。因此，对于在重尾噪声影响下的多变量Hammerstein模型的辨识问题的研究，刻不容缓。
技术实现思路
为了解决上述现有技术中的不足，本专利技术研究的对象为重尾噪声影响下的多变量Hammerstein模型，其数学表达为：xr(k)＝Fr(u(k)),r＝1,2,…,m其中u(k)＝[u1(k),u2(k),…,um(k)]Ty(k)＝[y1(k),y2(k),…,yn(k)]Tx(k)＝[x1(k),x2(k),…,xm(k)]Tε(k)＝[ε1(k),ε2(k),…,εm(k)]Tu(k)和y(k)分别为系统在k时刻的输入和输出，x(k)为非线性部分的输出向量，Fr(u(k))(r＝1,2,…,m)是Hammerstein模型的静态非线性部分，Ai∈Rn×n和Bj∈Rn×m是Hammerstein模型的静态线性部分的参数矩阵。ε(k)表示可由确定分布函数描述的重尾噪声向量。目的是辨识线性部分的参数矩阵Ai∈Rn×n和Bj∈Rn×m，同时训练RBF神经网络。本专利技术提出了针对复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，包括以下步骤：S1、产生辨识系统所用的测试激励信号，并将所产生的测试激励信号施加在待辨识的工业过程中；S2、采集工业过程对施加的测试激励信号产生的响应信号，也即系统的输出；S3、根据输入与输出数据之间的关系，待辨识工业过程的模型结构，利用GMDA-RBF方法，辨识多变量非线性系统的模型参数；S4、步骤S3中获得的模型参数，代入到Hammerstein模型结构中，得到工业过程模型；S5、向S4得到的工业过程模型施加测试信号，测试所得模型是否符合模型辨识要求。所述S1中提到的激励信号，在这里采用的是M序列。所述S3的具体过程为：S301、替换标准布谷鸟算法(CS)的种群更新方式，列维飞行，利用混合高斯分布产生随机行走，替代列维飞行，产生新的种群，具体操作如下：标准布谷鸟算法(CS)的种群更新方式如下：其中是当前时刻的个体值，是下一时刻的个体值，θ＞0是步长，代表矩阵元素乘法，列维飞行实际上是一种随机行走，它的步长是由列维分布产生：Lévy～u＝t-λ,(1＜λ≤3)利用二项高斯混合分布：其中是高斯分布，其均值和方差分别为μi＝0和表示脉冲组成部分，均值和方差分别为μi和其中ki＞＞1，脉冲发生的概率为αi，替换列维分布更新种群，如下：S302、利用混合高斯分布产生的随机序列，初始化布谷鸟算法的种群、对超出边界的值进行优化、调节步长，具体做法如下：首先通过二项高斯混合分布，产生随机数，长度与算法的种群大小相同为N，序列可以被描述为{GMD(N)}，这里假设GMDmax:＝max{GMD(N)}，GMDmin:＝min{GMD(N)}，利用归一化方法，建立从[GMDmin，GMDmax]到[0,1]的线性映射，得到归一化后的序列{GMD0,1(N)}。接下来就可以顺序提取序列{GMD0,1(N)}中相关的元素，对应到GMD(0,1)中，此处GMD(0,1)可直接用于接下来的改进；S303、利用GMD序列，初始化布谷鸟算法的种群初值，在标准布谷鸟算法中，随机化初始值产生方式如下：nest＝Lb+(Ub-Lb)×rand(0,1)在我们提出的GMDA算法中，设置初始值的产生方式如下：nest＝Lb+(Ub-Lb)×GMD(0,1)其中Ub和Lb分别是搜索区间的上下界，rand(0,1)是遵循0和1之间的正态分布的随机数，GMD(0,1)为遵循0和1之间二项高斯混合分布的随机数。S304、利用GMD序列优化边界调整过程：在标准布谷鸟算法中，一旦个体超出搜索区间，个体值就会被立即重置为边界值。这个过程可以被表达如下：在我们提出的GMDA算法中，我们使用以下方法来处理超出边界的个体：其中是当前时刻的候选解，Lb,Ub和GMD(0,1)与S303含义一致。S305、使用GMD序列调节步长：标准布谷鸟算法中步长是一个恒定值，一般来说取1，在GMDA算法中，我们使用以下方式来调节步长：θ＝θmax-GMD(0,1)*(θmax-θmin)其中θmax和θmin分别是步长的上下界，GMD(0,1)为遵循0和1之间二项高斯混合分布的随机数。S306、RBF神经网络部分的结构如图1以及需要训练的参数如下：本方法所采用的RBF神经网络的径向基函数是高斯函数，因此神经网络的输入和输出可由下式表示：其中i,j∈N+分别代表神经网络隐含层和输出层神经元的个数。u(k)＝[u1(k),u2(k),…,um(k)]是输入向量，是第j个输出，wij是第i个隐含层神经元到第j个输出的权值，‖·‖是欧几里得范数，mi和σi是第i个隐含层神经元的中心向量和宽度，参数l∈N+是隐含层神经元的个数。在这里，我们取mi,σi,和wij为需要GMDA优化的参数。S307、利用GMDA进行参数优化时，所需的目标函数设置如下，利用MSE指标来表示：其中，k是训练数据的个数，i多变量Hammerstein模型的通道数，yi(k)为系统真实输出，为辨识出的模型输出。搜索目标使MSE指标最小。另一方面，本专利技术提供复杂重尾噪声影响下多变量Hammerstein模型的辨识系统，包括依次连接的测试信号产生单元，系统数据输出单元，算法运行单元和辨识结果验证分析单元。测试信号产生单元：产生用于辨识的系统激励信号。系统数据输出单元：收集系统在接受激励信号刺激后的响应信号。算法运行单元：根据系统的输入输出信号，设立响应的辨识目标函数，运行辨识算法获得辨识结果。结果验证分析单元：将辨识结果代入系统的结构模型中，检验其是否达到辨识要求。附图说明图1本专利技术实施实例复杂重尾噪声影响下多变量Hammerstein模型的辨识方法流程。图2为RBF神经网络结构图。图3为多变量Hammerstein模型辨识的原理图。图4为蒸汽-水热交换器结构图。图5为辨识结果与原蒸汽-水热交换器系统对比图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下介绍一种蒸汽-水热交换器，结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用本文档来自技高网...

【技术保护点】
一种复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，其特征在于，包括以下步骤，S1、产生辨识系统所用的测试激励信号，并将所产生的测试激励信号施加在待辨识的工业过程中；S2、采集工业过程对施加的测试激励信号产生的响应信号，也即系统的输出；S3、根据输入与输出数据之间的关系，待辨识工业过程的模型结构，利用GMDA‑RBF方法辨识多变量非线性系统的模型参数；S4、将步骤S3中获得的模型参数，代入到Hammerstein模型结构中，得到工业过程模型；S5、向S4得到的工业过程模型施加测试信号，验证所得模型是否符合模型辨识要求。

【技术特征摘要】
1.一种复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，其特征在于，包括以下步骤，S1、产生辨识系统所用的测试激励信号，并将所产生的测试激励信号施加在待辨识的工业过程中；S2、采集工业过程对施加的测试激励信号产生的响应信号，也即系统的输出；S3、根据输入与输出数据之间的关系，待辨识工业过程的模型结构，利用GMDA-RBF方法辨识多变量非线性系统的模型参数；S4、将步骤S3中获得的模型参数，代入到Hammerstein模型结构中，得到工业过程模型；S5、向S4得到的工业过程模型施加测试信号，验证所得模型是否符合模型辨识要求。2.根据权利要求1所述的复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，其特征在于，所述S1中所提到的激励信号，在这里采用的是M序列。3.根据权利要求1所述的复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，其特征在于，所述的重尾噪声影响下的多变量Hammerstein模型，其数学表达为：xr(k)＝Fr(u(k))，r＝1，2，...，m其中u(k)＝[ui(k)，u2(k)，...，um(k)]Ty(k)＝[y1(k)，y2(k)，...，yn(k)]Tx(k)＝[x1(k)，x2(k)，...，xm(k)]Tε(k)＝[ε1(k)，ε2(k)，...，εm(k)]Tu(k)和y(k)分别为系统在k时刻的输入和输出，x(k)为非线性部分的输出向量Fr(u(k))(r＝1，2，...，m)是Hammerstein模型的静态非线性部分，Ai∈Rn×n和Bj∈Rn×m是Hammerstein模型的静态线性部分的参数矩阵。ε(k)表示可由确定分布函数描述的重尾噪声向量。目的是辨识线性部分的参数矩阵Ai∈Rn×n和Bj∈Rn×m，同时训练RBF神经网络。4.根据权利要求1所述的复杂重尾噪声影响下多变量Hammerstein模型的辨识方法，其特征在于，S3的步骤为：S301、替换标准布谷鸟算法(CS)的种群更新方式，列维飞行，利用混合高斯分布产生的产生随机行走，替代列维飞行，产生新的种群，具体操作如下：标准布谷鸟算法(CS)的种群更新方式如下：其中，是当前时刻的个体值，是下一时刻的个体值，θ＞0是步长，代表矩阵元素乘法，列维飞行实际上是一种随机行走，它的步长是由列维分布产生：Lévy～u＝t-λ，(1＜λ≤3)利用二项高斯混合分布：其中是高斯分布，其均值和方差分别为μi＝0和表示脉冲组成部分，均值和方差分别为μi和其中ki＞＞1，脉冲发生的概率为αi，替换列维分布更新种群，如下：S302、利用混合高斯分布产生的随机序列，初始化布谷鸟算法的种群、对超出边界的值进行优化、调节步长，具体做法如下：首先通过二项高斯混合分布，产生随机数，长度与算法的种群大小相同为N，序列可以...

【专利技术属性】
技术研发人员：靳其兵，王贺贺，蒋北艳，周星，
申请(专利权)人：北京化工大学，北京国控天成科技有限公司，
类型：发明
国别省市：北京,11