基于输入输出数据非线性系统建模.docVIP

基于输入输出数据非线性系统建模 　　【摘 要】随着机电系统的复杂程度和控制要求不断提高，需要运用控制效果更高的控制算法设计控制器，这就必然需要知道系统模型。针对工程对象普遍含有的非线性，对模块化非线性系统建模进行研究。利用正交基辨识法对输入非线性模型（Hammerstein-型）进行建模，对于正交基函数的获取进行改进，对输入输出数据的特殊处理，利用输出信号恢复中间变量，再利用最小二乘法得到模型参数，仿真结果表明了方法有效性。 　　【关键词】系统模型；非线性；Hammerstein-型 　　0 引言 　　线性系统建模方法虽然对于线性系统具有很好的效果。随着控制过程要求的不断提高，相对于一个线性模型在当前工作点的局部近似，非线性模型能更好地描述过程在整个运行过程的整体特性，在实际应用中具有很好的效果。模块化非线性模型是一种非线性模块与线性子系统的串联型的模型，由于其结构清晰，并能描述常见的非线性系统，所以得到了广泛地关注与应用[1-2]。 　　近年来，对于模块化非线性Hammerstein型的辨识文献相对较多。Narendrad等[3-4]提出了迭代方法，将参数化为线性模块参数和非线性模块参数两个集合，计算一个参数集最优估计值时固定另一个，两参数集估计轮换计算，但两个参数集之间的链接矩阵是一个输入变量的函数矩阵，由此导致迭代最小二乘法的协方差矩阵计算量大。Chang等[5]提出了过参数辨识方法，是把非线性展开为一些基函数的和，经参数化后得到一个过参数化模型，然而得到的模型参数向量包含了原始的静态非线性模块与线性模块参数的乘积，使得参数向量维数大大增加，导致算法计算量相应增加。Pawlak[6]给出随机方法，利用白噪声性质分离线性部分与非线性部分。Bai等[7-8]给出了分离最小二乘法、盲辨识和频域辨识法。 　　为此，本文利用正交基辨识法对输入非线性模型（Hammerstein-型）进行建模，其优点是避免了迭代算法以及参数向量维数增大所带来的计算量，对于正交基函数的获取进行改进，对输入输出数据做特殊处理，仅利用利用输出信号恢复中间变量，最终利用最小二乘法得到模型参数，仿真结果表明了方法有效性。 　　1 输入型非线性系统建模 　　1.1 输入型非线性系统描述 　　Hammerstein模型是一种输入端具有非线性的串联型非线性系统模型图1所示，被应用于许多工程问题中。考虑离散的Hammerstein型系统，建模的目的是仅基于输入数据u输出数据y，估计线性部分的传递函数G（z）以及非线性函数f（u），其中间变量x是不可测量的，预先设定线性部分的模型阶数为n。 　　图1 Hammerstein模型 　　输入端的非线性模块，通常以泰勒展开多项式的形式进行描述，即 　　其中r是非线性种多项式的个数，线性子系统可以通过输出信号与中间变量的离散传递函数描述为 　　1.2 线性传递函数分母参数确定 　　首先对系统的输入输出信号进行特殊处理，根据系统的特点消除中间变量，利用最小二乘法估计出传递函数分母参数，具体方法描述如下： 　　输入信号的采样间隔为T，不在采样点上的数据设置为零，对于输出使用间隔h=T/（n+1）来进行采样，如图2、图3所示。 　　图2 输入数据采样 　　图3 输出数据采样 　　将式（2）的线性模型转换为时域表达式 　　1.3 线性传递函数分子与非线性参数的确定 　　传统的正交基建模的方法，对于正交基函数的获取是通过先验知识或者相关的经验，本文在求得线性部分传递函数分母的基础上，可以求得传递函数的极点值，然后采用正交基辨识思想来利用极点构造标准正交基函数，将线性子系统表示为基函数的形式，结合非线性模块表达形式构造出关系矩阵，最后利用奇异值分解的方法获得各个模块的参数，从而减少了对于先验知识的依赖，并具有更高的准确性，具体描述如下所示： 　　参考图1，对线性部分采用基函数的表达形式： 　　其中，p=n+1，Bl（z）是线性部分的正交基函数。根据正交基系统辨识的方法在这里对于线性模块的基函数取做如下形式[9] 　　ξi是线性部分的极点，可以通过已经求的分母构造多项式进行求解得到。将tk时刻的输入表示为uk，将线性子系统与非线性模型相结合通过式（2）和（14）可以将tk时刻输出yk用如下的公式表示 　　观察上面的式子，通过最小二乘法或其他算法进进行估计参数的时候得到的是bc的形式，这也是所有的串联模块化非线性辨识的通病，结果经过分离后得到的参数往往是[αc，α-1b]的形式，其中α为一常数，因为在分离两者时没有一个统一的标准，即前面有可能多乘了一个常数，后面少乘了个常数，但是整体的效果是一样。这就要求有一个标准使得到的参数保持唯一性，大多方法都是使分子参数的二范数为1，这种参数的不唯一性对于