基于最小二乘法的系统辨识问题研究综述

基于最小二乘法的系统辨识问题研究综述

摘要：对基于最小二乘法的系统辨识方法进行了介绍。首先对系统辨识概念以及最小二乘法原理进行了介绍，然后根据例子来说明怎样运用最小二乘法来解决实际辨识问题。而且本文针对最小二乘存在的缺陷进一步阐述了一些改进型最小二乘法在系统辨识中的应用，最后对系统辨识的发展趋势做了预测。

关键字：系统辨识最小二乘法改进型最小二乘法发展趋势

1引言

系统辨识归根到底是一种数学建模的过程，而建模过程中运用的方法并不唯一，最小二乘法是较早被应用于系统辨识中的一类方法。1962年，L. A. Zadeh 最先提出了系统辨识的定义[1]：“辨识就是在输入和输出数据的基础上，从一组给定的模型类中，确定一个与所测系统等价的模型。”简单的来说，就是在现有数据的基础上，按照一个准则在一组模型类中选择一个与提供的数据拟合得最好的模型。而根据最小二乘法的定义[2]：“最小二乘法是一种数学优化技术，它通过最小化误差的平方和寻找数据的最佳函数匹配。”其基本思想就是让实测数据和估计数据之间的平方和最小，这恰恰是系统辨识所需要解决的问题，所以最小二乘法很早就被用来求解辨识中需要的拟合数学模型。

本文在阐述最小二乘法理论的基础上对于其在系统辨识中的应

用做了介绍，并指出实际应用中存在的不足，列举了几种改进型的最小二乘算法来改进系统辨识能力，最后给出了系统辨识的发展趋势。

2 基于最小二乘法的系统辨识的理论基础及应用

2.1 系统辨识的理论基础

从字面上讲，系统辨识( System Identification) 就是识别一个系统、辨识一个系统[3]。系统通常是由表征系统输入输出关系的数学模型描述的，这个模型有其特定的结构和参数。因此，系统辨识包含系统结构辨识( System Structure Identification) 和参数估计( Parameter Estimation) ．

系统结构(或模型结构) 就是系统数学表达式的形式。对单输入单输出线性系统而言，模型结构就是系统的阶次(Order) ；对多变量线性系统而言，模型结构就是系统的能控性结构指数(Controllability Structure Index) 或能观测性结构指数(Observability Structure Index) ，系统阶次等于系统的能控性结构指数或能观测性结构指数之和。对传递函数而言，系统参数就是传递函数分子分母多项式的系数(Coeffi-cient) ，系统阶次就是传递函数分母多项式的次数(Degree) ; 对状态空间模型而言，系统参数就是状态空间模型的A，B，C，D 矩阵，就是状态向量的维数或矩阵的维数，它等于系统的能控性结构指数系统阶次或能观测性结构指数之和。

求解系统辨识问题实质上就是找到合适的数学方法来判断系统

的结构以及得到系统参数。 2.2 最小二乘法的理论基础

最小二乘法作为一种传统的参数估计方法，早已经被大家所了解。 然而大多同学对最小二乘法的认识都比较模糊，仅仅把最小二乘法理解为简单的线性参数估计。 事实上，最小二乘法在参数估计、系统辨识以及预测、预报等众多领域都有着广泛的应用。特别是针对动态系统辨识的方法有很多，但其中应用最广泛，辨识效果良好的就是最小二乘辨识方法，研究最小二乘法在系统辨识中的应用具有现实的、广泛的意义。因此要用最小二乘法解决实际的辨识问题，首先要对最小二乘法有深刻理解。

下面是一般的最小二乘法问题：实系数线性方程组 a 11x 1+a 12x 2+…+a 1n x n -b 1=0 a 21x 1+a 22x 2+…+a 2n x n -b 2=0

………… (2.1) a m1x 1+a m2x 2+…+a mn x n -b m = 0

方程组可能无解。 即可能不存在一组实数 x 1,x2,……,x s 使

211

21

20()i i in n i m

i a

x a x a x b =++⋯+-=∑ (2.2)

恒成立。因此我们转而求其次，设法找到实数组 x 1，x 2，…，x s 使误差的平方和最小，这样的 x 1，x 2，…，x s 称为方程组的最小二乘解，这样问题就叫最小二乘法问题。

2.3 最小二乘法处理辨识问题的应用举例

考虑如下线性系统：

()()()()()()1111a b n a n b z k a z k a z k n bu k b u k n e k +-+

+-=-++-+ （1）

其中，u(k)为系统激励信号，y(k)为系统输出，e(k)为模型噪声。

其系统模型如图1所示：

图1 SISO 的系统模型结构图

其中G(z -1)是系统函数模型，N(z -1)为有色噪声系统模型，e(k)为白噪声v(k)经过系统函数为N(z -1)的系统后的输出。通常

()()

()

()()

()

111

1

1

1

, B z D z G z

N z

A

z C

z ------== （2）

式中：

()()11212112

121a

a

b b n n n n A z a z a z a z B z b z b z b z

--------⎧=++++⎪⎨=+++⎪⎩ （3） ()()11212112

121c

c

d d n n n n C z c z c z c z D z d z d z b z

--------⎧=++++⎪⎨=+++⎪⎩ （4） 则系统可表示为：()()()

()()

()

()1111B z D z z k u k v k A z C z ----=+

（5）

设样本和参数集为：

1212()[-(-1) , -(-2), ...... -(-), (-1),(-2), ......, (-)][,,......,,,,......,]

T

T

n n h k z k z k z k n u k u k u k n a a a b b b θ⎧=⎨=⎩ （6） h(k)为可观测的量， 差分方程可写为最小二乘形式

()()()T z k h k e k θ=+ （7）

如何在系统噪声e(k)存在的情况下从该方程中正确的解出θ，即是系统辨识的任务。 为了求出θ，我们面临三大问题：一是输入信号的选择，二是判决准则的选取，三是辨