基于最小二乘法的系统辨识问题研究综述概要Word下载.docx
《基于最小二乘法的系统辨识问题研究综述概要Word下载.docx》由会员分享,可在线阅读,更多相关《基于最小二乘法的系统辨识问题研究综述概要Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
2基于最小二乘法的系统辨识的理论基础及应用
2.1系统辨识的理论基础
从字面上讲,系统辨识(SystemIdentification)就是识别一个系统、辨识一个系统[3]。
系统通常是由表征系统输入输出关系的数学模型描述的,这个模型有其特定的结构和参数。
因此,系统辨识包含系统结构辨识(SystemStructureIdentification)和参数估计(ParameterEstimation).
系统结构(或模型结构)就是系统数学表达式的形式。
对单输入单输出线性系统而言,模型结构就是系统的阶次(Order);
对多变量线性系统而言,模型结构就是系统的能控性结构指数(ControllabilityStructureIndex)或能观测性结构指数(ObservabilityStructureIndex),系统阶次等于系统的能控性结构指数或能观测性结构指数之和。
对传递函数而言,系统参数就是传递函数分子分母多项式的系数(Coeffi-cient),系统阶次就是传递函数分母多项式的次数(Degree);
对状态空间模型而言,系统参数就是状态空间模型的A,B,C,D矩阵,就是状态向量的维数或矩阵的维数,它等于系统的能控性结构指数系统阶次或能观测性结构指数之和。
求解系统辨识问题实质上就是找到合适的数学方法来判断系统的结构以及得到系统参数。
2.2最小二乘法的理论基础
最小二乘法作为一种传统的参数估计方法,早已经被大家所了解。
然而大多同学对最小二乘法的认识都比较模糊,仅仅把最小二乘法理解为简单的线性参数估计。
事实上,最小二乘法在参数估计、系统辨识以及预测、预报等众多领域都有着广泛的应用。
特别是针对动态系统辨识的方法有很多,但其中应用最广泛,辨识效果良好的就是最小二乘辨识方法,研究最小二乘法在系统辨识中的应用具有现实的、广泛的意义。
因此要用最小二乘法解决实际的辨识问题,首先要对最小二乘法有深刻理解。
下面是一般的最小二乘法问题:
实系数线性方程组
a11x1+a12x2+…+a1nxn-b1=0
a21x1+a22x2+…+a2nxn-b2=0
…………(2.1)
am1x1+am2x2+…+amnxn-bm=0
方程组可能无解。
即可能不存在一组实数x1,x2,……,xs使
(2.2)
恒成立。
因此我们转而求其次,设法找到实数组x1,x2,…,xs使误差的平方和最小,这样的x1,x2,…,xs称为方程组的最小二乘解,这样问题就叫最小二乘法问题。
2.3最小二乘法处理辨识问题的应用举例
考虑如下线性系统:
(1)
其中,u(k)为系统激励信号,y(k)为系统输出,e(k)为模型噪声。
其系统模型如图1所示:
图1SISO的系统模型结构图
其中G(z-1)是系统函数模型,N(z-1)为有色噪声系统模型,e(k)为白噪声v(k)经过系统函数为N(z-1)的系统后的输出。
通常
(2)
式中:
(3)
(4)
则系统可表示为:
(5)
设样本和参数集为:
(6)
h(k)为可观测的量,差分方程可写为最小二乘形式
(7)
如何在系统噪声e(k)存在的情况下从该方程中正确的解出
,即是系统辨识的任务。
为了求出
,我们面临三大问题:
一是输入信号的选择,二是判决准则的选取,三是辨识算法的选择,下面一一探讨。
一.选择输入
为了准确辨识系统参数,我们对输入信号有两大要求,一是信号要能持续的激励系统所有状态,二是信号频带能覆盖系统的频带宽度。
除此之外还要求信号有可重复性,不能是不可重复的随机噪声,因此我们通常选择M序列或逆M序列作为输入。
二.准则函数
因为本文主要探讨最小二乘类辨识方法,在此选取准则函数
(8)
使准则函数
的
估计值记做
,称作参数
的最小二乘估计值。
在式(7)中,令k=1,2,3,……L,可构成线性方程组:
(9)
式中
(10)
准则函数相应变为:
(11)
极小化
求得参数
的估计值,将使模型更好的预报系统的输出。
三.最小二乘辨识算法
设
使得
,则有
(12)
展开上式,并根据以下两个向量微分公式:
(13)
得正则方程:
(14)
当
为正则阵时,有
(15)
且有
,所以满足式(15)的
唯一使得
,这种通过极小化式(11)计算
的方法称作最小二乘法。
而且可以证明,当噪声e(k)是均值为0的高斯白噪声时,可实现无偏估计
3最小二乘法在辨识时存在的问题及其改进型的系统辨识应用举例
3.1传统最小二乘存在的问题
最小二乘法存在一些缺陷制约着最小二乘法在系统辨识中的应用,在处理日益复杂的系统辨识问题中,最小二乘法在系统辨识中存在的缺陷逐渐显现出来。
如传统的最小二乘法不适合在动态辨识系统中使用,而且其参数估计存在偏差,耗时较长的问题。
3.2递推最小二乘算法
为了减少计算量,减少数据在计算机中占用的内存,并实时辨识出系统动态特性,我们常利用最小二乘法的递推形式。
下面我们来推导递推最小二乘算法的原理。
首先,将式(11)的最小二乘一次完成算法写为
(16)
定义
(17)
(18)
式中,h(i)是一个列向量,也就是HL的第i行的倒置,P(k)是一个方阵,它的维数取决于未知参数的个数,假设未知参数的个数是n,则P(k)的维数是n×
n.。
由式17可得P(k)的递推关系为:
(19)
(20)
则
(21)
由此可得:
(22)
由式19和22可得
(23)
引进增益矩阵K(k),定义
(24)
式23可以进一步写为
(25)
接下来可以进一步把式20写为
(26)
利用矩阵反演公式
将式(26)演变成
(27)
将上式代入式24,整理后可得
(28)
综合式25、27和28可得最小二乘递推参数估计算法RLS
3.3广义最小二乘法
设SISO系统采用如下模型:
(30)
其中A(z-1),B(z-1)和C(z-1)的定义见式?
假定模型阶次na,nb和nc已知,用广义最小二乘法可以得到无偏一致估计。
令
(31)
及
(32)
将模型化为最小二乘格式:
(33)
由于v(k)是白噪声,所以用最小二乘可以获得参数θ的无偏估计,由于噪声模型C(z-1)未知,还需要用迭代的方法来求得C(z-1)。
(34)
置
(35)
这样就把噪声模型也转变为最小二乘格式:
(36)
由于上式中的噪声已为白噪声,所以用最小二乘也可获得参数θe的无偏估计,但是数据向量中依然含有不可测的噪声量
,可用相应的估计值来代替,置
,其中k<
0时,e(k)=0;
k>
0时,按照
(37)
计算,式中
(38)
综上所述,广义最小二乘法可归纳为
4系统辨识的发展及对策
最小二乘法(LS)是一种经典的和最基本的,也是应用最广泛的方法。
但是,最小二乘估计是非一致的,是有偏差的,所以为了克服他的缺陷,而形成了一些以最小二乘法为基础的系统辨识方法:
广义最小二乘法、增量最小二乘法、渐消记忆的最小二乘法以及将一般的最小二乘法与其他方法相结合的方法,有最小二乘两步法(COR-LS)和随机逼近算法等。
广义最小二乘法的基本思想是引入一个白化滤波器,把相关噪声转换为白噪声,基于对观测数据先进行一次滤波处理,然后利用普通最小二乘法对滤波后的数据进行辨识。
递推算法的基本思想是用新估计值对老的估计值进行修正,应用的数据是实时采集的系统输入输出数据,应用递推算法对参数估计值进行不断修正,以取得更为准确的参数估计值,而且此方法占据计算机存储量小,因此在辨识系统中得到了广泛的应用。
【最小二乘参数估计的递推算法及其C语言实现】
渐消记忆的最小二乘法在对系统模型参数进行辨识时强调新数据的作用,贬低老数据的作用,将遗忘因子λ引入系统中,在实际应用中遗传因子λ的大小对参数估计的精度以及参数估计值跟踪真值的变化的能力都有很大的影响,所以选取合适的遗传因子会显著提高系统的辨识能力。
5结论
系统辨识的方法有很多,但是最小二乘法是其中最重要、应用最广泛的一种系统辨识方法。
针对经典的最小二乘法存在的一些不足,广义最小二乘法、递推式最小二乘法以及渐消记忆的最小二乘法等是一些在科学研究中经常采用的系统辨识方法,这些方法不仅在今天,而且在未来都会有非常广泛的应用前景,此外新的改进型最小二乘法也会在其他学科的发展下涌现,使基于最小二乘法的系统能适应更多的场合,得到更加强大的系统辨识能力。
参考文献
[1]ZadehLA.FromCircuitTheorytoSystemTheory[J].Proc.IRE,1962,50(5):
856-865.
[2]维基百科
[3]系统辨识
(1):
辨识导引
[1]佚名.最小二乘法的论述与证明[J].刊名缺失,出版年缺失,卷缺失(期缺失):
页码范围缺失.