基于最小二乘法的系统参数辨识吴令红，熊晓燕，张涛太原理工大学机械电子研究所，太原 (030024)E-mail lhwu0818@摘要：系统辨识是自动控制学科的一个重要分支，由于其特殊作用，已经广泛应用于各种领域，尤其是复杂系统或参数不容易确定的系统的建模。
过去，系统辨识主要用于线性系统的建模，经过多年的研究，已经形成成熟的理论。
但随着社会、科学的发展，非线性系统越来越受到人们的关注，其控制与模型之间的矛盾越来越明显，因而非线性系统的辨识问题也越来越受到重视，其辨识理论不断发展和完善本。
文重点介绍了系统参数辨识中最小二乘法的基本原理，并通过悬臂梁模型的辨识实例，具体说明了基于最小二乘法参数辨识在Matlab 中的实现方法。
结果表明基于最小二乘法具有算法简单、精度较高等优点。
关键词：系统辨识；参数辨识；滑动平均模型(ARX)；最小二乘法；Matlab中图分类号：TH-91. 引言所谓辨识就是通过测取研究对象在人为输入作用下的输出响应，或正常运行时的输入输出数据记录，加以必要的数据处理和数学计算，估计出对象的数学模型。
这是因为对象的动态特性被认为必然表现在它的变化着的输入输出数据之中，辨识只不过是利用数学的方法从数据序列中提炼出对象的数学模型而已[1]。
最小二乘法是系统参数辨识中最基本最常用的方法。
最小二乘法因其算法简单、理论成熟和通用性强而广泛应用于系统参数辨识中。
本文基于悬臂梁的实测数据，介绍了最小二乘法的参数辨识在Matlab中的实现。
2. 系统辨识一般而言，建立系统的数学模型有两种方法：激励分析法和系统辨识法。
前者是按照系统所遵循的物化（或社会、经济等）规律分析推导出模型。
后者则是从实际系统运行和实验数据处理获得模型。
如图1所示，系统辨识就是从系统的输入输出数据测算系统数学模型的理论和方法。
更进一步的定义是L.A.Zadeh曾经与1962年给出的，即“系统辨识是在输入和输出的基础上，从系统的一类系统范围内，确立一个与所实验系统等价的系统”。
另外，系统辨识还应该具有3个基本要素，即模型类、数据和准则[5]。
被辨识系统模型根据模型形式可分为参数模型和非参数模型两大类。
所谓参数模型是指微分方程、差分方程、状态方程等形式的数学模型；而非参数模型是指频率响应、脉冲响应、传递函数等隐含参数的数学模型。
在辨识工程中，模型的确定主要根据经验对实际对象的特性进行一定程度上的假设，如对象的模型是线性的还是非线性的、是参数模型还是非参数模型等。
在模型确定之后，就可以根据对象的输入输出数据，按照一定的辨识算法确定模型的参数[4]。
y图1 被研究的动态系统3. 参数辨识模型结构及参数估计方法3.1 参数辨识模型结构在工程技术领域，描述离散系统的动态特性的一般模型是n 阶差分方程。
对于一个所谓的“黑盒子”系统，其差分方程的最一般的表达式可定义为：()()()()()()()()k e q d q C nk k u q F q B k y qA 11111−−−−−+−= (1)其中：()()f d c b a x F D CB A X q x q x x q X nx nx ，，，，，，，，==+++=−−−;1101Λ；1−q 为单位滞后算子；()nf nd nc nb na nx ，，，，分别为多项式F D C B A ，，，，的阶次，nk为系统的延迟；()k e 为系统的噪声。
根据式(1)中的不同参数取值，可将参数辨识模型分为以下几种常见模型：自回归模型(AR模型)、滑动平均模型(ARX 模型)和自回归滑动平均模型(AMAX)等[3]。
在本文中选用ARX模型结构，即滑动平均模型：()()()()()k e nk t u q B k y q A +−=−−11。
因为此模型结构具有如下优点：在噪声较小时，ARX 模型具有鲁棒性、算法简单等优点；在噪声较大时，可适当提高模型的阶次，以补偿噪声对辨识精度的影响。
3.2 参数估计方法对于参数模型辨识结构，系统辨识的任务是参数估计，即利用输入输出数据估计这些参数，建立系统的数学模型。
在参数估计中最常用的是最小二乘法(LS)、误差预测估计法、辅助变量法(IV 估计)、神经网络法等[4]。
由于最小二乘法容易理解和掌握，利用最小二乘法原理所拟定的辨识算法在实施上比较简单，并且不需要数理统计的知识，使得最小二乘法广泛应用于系统辨识领域，但它也存在着一定的局限和不足，当系统噪声为有色噪声时，最小二乘法不能给出无偏一致估计[5]。
本文应用最小二乘法对悬臂梁数学模型进行了参数辨识。
3.3 最小二乘法(LS)设一个SISO(单输入/单输出)过程的ARX 模型：()()()()()k e nk t u q B k y q A +−= (2)其中：()k u 、()k y 为系统的输入输出量，()k e 为噪声。
式中()()nb nb nana q b q b q B q a q a q A −−−−++=+++=ΛΛ11111，1,1−=t nk Λ模型(2)可化为最小二乘格式：()()()k e k h k y +=θτ (3)记[]T nb na b b a a a b 2121ΛΛ，，，，=θ为待估计的参数。
()()()()()[]τnb k u k u na k y k y k h −−−−−−=,,1,,1ΛΛ，对于L k Λ2,1=(L 为数据长度)，方程(3)构成一个线性方程组，写成()()K e k H y L L L +=θ()()()()()()()()()⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=L e e e e L h h h H L y y y y L L L ΜΜΜ212121，，τττ 根据最小二乘一次完成算法，其参数估计值为：()L L L L Z H H H LS ττθ1−∧= (4) 4. 实例分析本文以悬臂梁为研究对象，图2为悬臂梁振动实验的接线原理图，图3为实验采集的某一段输出信号。
本文利用系统辨识理论知识对悬臂梁建立模型，采用Matalb 最小二乘法[2](LS)进行参数估计。
图2 悬臂梁振动实验原理图图3 输出信号在Matlab环境下，编辑M函数文件，Load shiyanshuju.mat; %导入用于模型辨识的实验数据。
shuju=iddata(double(Data1_AI_16),double(Data1_AI_0),0.001);%把待辨识的数据由单精度转换为双精度，生成一个包含这些输入/输出频道的iddata对象。
ze=shuju(1000:2000);%取采样数据中的一段数据m2=arx(ze,[10 10 1]); %%基于最小二乘法估计的ARX模型辨识。
Present(m2); %%显示模型的信息。
figure,compare(ze,m2,1,100:300); %预测输出与实际输出比较。
利用最小二乘估计的辨识方法，辨识结果如下：>> Discrete-time IDPOL Y model: A(q)y(t) = B(q)u(t) + e(t)A(q) = 1 - 1.415 q^-1 + 0.7886 q^-2 - 0.137 q^-3 + 0.1522 q^ -4 - 0.283 q^-5 + 0.4221 q^-6 -0.2915 q^-7 - 0.04392 q^-8 + 0.5723 q^-9 - 0.4517 q^-10 B(q) = -48.99 q^-1 + 80.47 q^-2 - 0.4402 q^-3 - 97.77 q^-4 - 126.5 q^-5 + 213.4 q^-6 + 129.3 q^-7 - 41.96 q^-8 - 139.1 q^-9 + 31.36 q^-10Estimated using ARX from data set zeLoss function 0.291794 and FPE 0.297686Sampling interval: 0.001图4 基于最小二乘法的系统辨识模型同实际输出比较结果从参数辨识的结果可得出结论：从系统参数、损失函数、最终预测误差和辨识模型同实际输出结果的拟合度可以看出，最小二乘法辨识的ARX模型参数估计精度比较高，证明了最小二乘法的有效性。
4. 结束语1、ARX模型是时间序列中应用较多的一种模型，通过建立合理的模型并对参数进行精确的估计，可以得到较满意的结果，其模型的定阶和参数的估计以及实验的精度还需要进一步研究。
2、系统辨识的方法有很多,最小二乘法是其中最基本、应用最广泛的一种系统辨识方法。
最小二乘法可用于动态系统，也可以用于静态系统；可用于线性系统，也可以用于非线性系统；可用于离线估计，也可以用于在线估计，这种方法不仅在今天，而且在未来都会有非常广泛的应用前景。
参考文献[1] 李言俊，张科. 系统辨理论及应用[M]. 北京：国防工业出版社，2006[2] 刘叔军，盖晓华，樊京，崔世林等.MA TLAB7.0控制系统应用与实例[M]. 北京：机械工业出版社，2005[3] 石贤良，吴成富. 基于MTALAB的最小二乘法参数辨识与仿真[J]. 微机处理，2005:44-46[4] 严晓久，周爱国，林建平，吴晔. 基于辅助变量法的系统参数辨识[J]. 机床与液压，2006:180-184[5] 王秀峰，卢桂章. 系统建模与辨识[M]. 北京：电子工业出版社，2004Parametric System Identification Based on Least-squaresMethodWu Ling Hong，Xiong Xiao Yan，Zhang TaoTai Yuan University of Technology, Department of Mechatronical Engineering , Tai yuan(030024)AbstractAs one of the main branches of the automatic control science, systemidentification has been applied in many fields. In the past, system identification was mostly applied in linear system modeling, and the flawless theory for linear system had come into being for many years. However, with the development of society and science, nonlinear system is more and more important. The conflict between control and model is getting more and more evident. This fact results in the development of nonlinear system identification theory. The fundamental principle of least-squares method was introduced in this paper, and amodeling of cantilever beam experiment was given to show how to implement system identification based on Least-squares Method in Matlab. The results of modeling indicate that system identification with least-squares method is more feasible and accurate.Keywords: System identification; Parametric identification; ARX; Least square method; Matlab作者简介：吴令红，男，1981年出生，硕士研究生，主要研究方向是机电系统动态测试与故障诊断。
