系统辨识最小二乘法解剖.docx

1、系统辨识最小二乘法解剖最小二乘法参数辨识系统辨识是根据系统的输入输出时间函数来确定描述系统行为的数学模型。 现代控制理论中的一个分支。通过辨识建立数学模型的U的是估计表征系统行为 的重要参数，建立一个能模仿真实系统行为的模型，用当前可测量的系统的输入 和输出预测系统输出的未来演变，以及设讣控制器。对系统进行分析的主要问题 是根据输入时间函数和系统的特性来确定输出信号。对系统进行控制的主要问题 是根据系统的特性设计控制输入，使输出满足预先规定的要求。而系统辨识所研 究的问题恰好是这些问题的逆问题。通常，预先给定一个模型类u=M（即给 定一类已知结构的模型），一类输入信号u和等价准则J=L（y,

2、yM）（般情况下, J是误差函数，是过程输出y和模型输出yM的一个泛函）；然后选择使误差函数 J达到最小的模型，作为辨识所要求的结果。系统辨识包括两个方面：结构辨识 和参数估计。在实际的辨识过程中，随着使用的方法不同，结构辨识和参数估计 这两个方面并不是截然分开的，而是可以交织在一起进行的。2系统辨识的目的在提出和解决一个辨识问题时，明确最终使用模型的的是至关重要的。 它对模型类（模型结构）、输入信号和等价准则的选择都有很大的影响。通过辨 识建立数学模型通常有四个目的。1估计具有特定物理意义的参数有些表征系统行为的重要参数是难 以直接测量的，例如在生理、生态、环境、经济等系统中就常有这种情况。

3、这就 需要通过能观测到的输入输出数据，用辨识的方法去估计那些参数。2仿真仿真的核心是要建立一个能模仿真实系统行为的模型。用于 系统分析的仿真模型要求能真实反映系统的特性。用于系统设讣的仿真，则强调 设计参数能正确地符合它本身的物理意义。3预测这是辨识的一个重要应用方面，其目的是用迄今为止系统的 可测量的输入和输岀去预测系统输出的未来的演变。例如最常见的气象预报，洪 水预报，其他如太阳黑子预报，市场价格的预测，河流污染物含量的预测等。预 测模型辨识的等价准则主要是使预测误差平方和最小。只要预测误差小就是好的 预测模型，对模型的结构及参数则很少再有其他要求。这时辨识的准则和模型应 用的目的是一致的

4、，因此可以得到较好的预测模型。4控制 为了设计控制系统就需要知道描述系统动态特性的数学模 型，建立这些模型的U的在于设计控制器。建立什么样的模型合适，取决于设计 的方法和准备采用的控制策略。3系统辨识的方法经典方法：经典的系统辨识方法的发展已经比较成熟和完善，他包括阶跃响应法、脉冲 响应法、频率响应法、相关分析法、谱分析法、最小二乘法和极大似然法等。其 中最小二乘法(LS)是一种经典的和最基本的，也是应用最广泛的方法。但是，最 小二乘估计是非一致的，是有偏差的，所以为了克服他的缺陷，而形成了一些以 最小二乘法为基础的系统辨识方法：广义最小二乘法(GI S)、辅助变量法(IV). 增广最小二乘法

5、(El, S)和广义最小二乘法(GIS),以及将一般的最小二乘法与其 他方法相结合的方法，有最小二乘两步法(COR-I S)和随机逼近算法等。经典的系统辨识方法还存在着一定的不足：(1)利用最小二乘法的系统辨识 法一般要求输入信号已知，并且必须具有较丰富的变化，然而，这一点在某些动 态系统中，系统的输入常常无法保证；(2)极大似然法计算耗费大，可能得到的 是损失函数的局部极小值；(3)经典的辨识方法对于某些复杂系统在一些情况下 无能为力。现代方法：随着系统的复杂化和对模型精确度要求的提高，系统辨识方法在不断发展， 特别是非线性系统辨识方法。主要有：1、 集员系统辨识法在1979年集员辨识首先岀

6、现于Fogel撰写的文献中，1982年Fogel和Huang 乂对其做了进一步的改进。集员辨识是假设在噪声或噪声功率未知但有界 UBB(Unknown But Bounded)的情况下，利用数据提供的信息给参数或传递函数 确定一个总是包含真参数或传递函数的成员集(例如椭球体、多面体、平行六边 体等)。不同的实际应用对象，集员成员集的定义也不同。集员辨识理论已广泛 应用到多传感器信息融合处理、软测量技术、通讯、信号处理、鲁棒控制及故障 检测等方面。2、 多层递阶系统辨识法多层递阶方法的主要思想为：以时变参数模型的辨识方法作为基础，在输入 输出等价的意义下，把一大类非线性模型化为多层线性模型，为非

7、线性系统的建模给出了一个十分有效的途 径。3、 神经网络系统辨识法山于人工神经网络具有良好的非线性映射能力、自学习适应能力和并行信息 处理能力，为解决未知不确定非线性系统的辨识问题提供了一条新的思路。与传统的基于算法的辨识方法相比较，人工神经网络用于系统辨识具有以下 优点：(1)不要求建立实际系统的辨识格式，可以省去对系统建模这一步骤；(2) 可以对本质非线性系统进行辨识；(3)辨识的收敛速度仅与神经网络的本身及所 采用的学习算法有关:(4)通过调节神经元之间的连接权即可使网络的输出来逼 近系统的输出；(5)神经网络也是系统的一个物理实现，可以用在在线控制。4、 模糊逻辑系统辨识法模糊逻辑理论

8、用模糊集合理论，从系统输入和输出的量测值来辨识系统的模 糊模型，也是系统辨识的一个新的和有效的方法，在非线性系统辨识领域中有十分广泛的应用。模糊逻辑辨识 具有独特的优越性：能够有效地辨识复杂和病态结构的系统：能够有效地辨识具 有大时延、时变、多输入单输出的非线性复杂系统；可以辨识性能优越的人类控 制器；可以得到被控对象的定性与定量相结合的模型。模糊逻辑建模方法的主要 内容可分为两个层次：一是模型结构的辨识，另一个是模型参数的佔计。典型的 模糊结构辨识方法有：模糊网格法、自适应模糊网格法、模糊聚类法及模糊搜索 树法等。5、小波网络系统辨识法小波网络是在小波分解的基础上提出的一种前馈神经网络口，使

9、用小波网 络进行动态系统辨识，成为神经网络辨识的一种新的方法。小波分析在理论上保 证了小波网络在非线性函数逼近中所具有的快速性、准确性和全局收敛性等优 点。小波理论在系统辨识中，尤其在非线性系统辨识中的应用潜力越来越大，为 不确定的复杂的非线性系统辨识提供了一种新的有效途径，其具有良好的应用前 景。1.设计内容设SISO系统的差分方程为：Z伙)+ qz伙-1) + d忆伙 一 2) = buk -1)4- b：u(k -2)4- u(k) 式(IT)参数取真值为：少T1H6 0.483 0.57 0.42,利用NIATLAB的M语言辨识系统中 的未知参数、S、仇。2.设计过程2.1 建立系统设

10、SIS0系统的差分方程为：z伙)+ dZ伙1) + n2z(k 2) = hxii(k 1) + 伙2) + u(k) 式(2-1)参数取真值为：少=1376 0.483 0.57 0.42,利用陋巩屈的M语言辨识系统中 的未知参数坷、幻、9、b2 o要求：用参数的真值利用差分方程求出z伙)作为测量值，u伙)是均值为0, 方差为0.1、0.5和0.01的不相关随机序列。选取一种最小二乘算法辨识。2. 2 最小二乘简介2. 2.1 最小二乘法的概念与应用对工程实践中测得的数据进行理论分析，用恰当的函数去模拟数据原型是一 类十分重要的问题，最常用的逼近原则是让实测数据和估计数据之间的距离平方 和最

11、小，这即是最小二乘法。最小二乘法是一种经典的数据处理方法。在系统辨 识领域中，最小二乘法是一种得到广泛应用的估计方法，可用于动态系统，静 态系统,线性系统，非线性系统。可用于离线估计，也可用于在线佔计。这种 辨识方法主要用于在线辨识。在随机的环境下，利用最小二乘法时，并不要求观 测数据提供其概率统计方面的信息，而其佔计结果，却有相出好的统计特性。MATLAB是一套高性能数字计算和可视化软件，它集成概念设计，算法开 发，建模仿真，实时实现于一体，构成了一个使用方便、界面友好的用户环境, 其强大的扩展功能为各领域的应用提供了基础。对于比较复杂的生产过程，山于 过程的输入输出信号一般总是可以测量的，

12、而且过程的动态特性必然表现在这 些输入输出数据中,那么就可以利用输入输出数据所提供的信息来建立过程的 数学模型。这种建模方法就称为系统辨识。把辨识建模称作“黑箱遂模2. 2. 2 最小二乘法系统辨识结构:本文把待辨识的过程看作“黑箱S只考虑过程的输入输出特性，而不强调 过程的内部机理。图中，输入u(k)和输出z(k)是可以观测的；G (z-1)是系统模型，用来描述系 统的输入输出特性；N (z-1)是噪声模型，v(k)是白噪声，e(k)是有色噪声，根 据表示定理：可以表示为e(k) =N (z-1) v(k)G(zr = 3(=)7V(z)=(r)A(z ) = 1 + ez i + + +

13、n,V耐)詁宀+jC(z_1) = 1 + qz-1 + c2z-2 +- - -+ cnn zTHa(z_) = z+2广2 + + dnt znh2. 2.3 准则函数设一个随机序列&伙)，(12,厶)的均值是参数0的线性函数： Ez伙)=hT(k)0f其中爪幻是可测的数据向量，那么利用随机序列的一个实现, 使准则函数：&）=习粼）一只&e （式 2-2）达到极小的参数估计值&称作&的最小二乘估计。最小二乘格式：Z伙）= /r（w+e，&为模型参数向量，讹）为零均值 随机噪声。2. 3 广义最小二乘法2. 3.1 广义最小二乘数学模型川以比伙）=B（Q加伙）+ 呛）式中，u（k）和Z伙）表

14、示系统的输入输出；v（k）是均值为零的不相关的随机序列；且4（昇）=1 + 昇+厂+ % B（d+b2厂C（zl） = l + clZl + c2-2 + + c.厂2. 3. 2 广义最小二乘递推算法如下0伙）=6伙1） + K, （k）Z/ 伙）-h；伙）0伙-1）K,伙）=P 伙l）h 伙）h；伙比伙l）hy伙）+ lpP 伙）= I K ）h：伙）P 伙1）Q 伙）= 伙1） + K伙）歟）-h：伙）0伙1） K） = 伙l）h,伙）h：伙）P,伙一 l）h0.5,u(i)=-l; %如果M序列的值为1,辨识的输入信号取“-1” else u(i)=l; %如果M序列的值为0,辨识的输

15、入信号取“1” end %小循环结束yl=xl; y2=x2; y3=x3: y4=x4; %为F次的输入信号做准备end %大循环结束，产生输入信号u figure(l); %第一个图形stem(u),grid on %显示岀输入信号M序列径线图并给图形加上网格v=normrnd(0, sqrt(O.Ol), 1, 300);%均值为零的,方差为 0. 01 或 0. 5 或 0.1 不相 关的随机噪声ze(2)=0;ze(l)=0;for k=3:301;ze(k) =0*ze(k-1) +0*ze(k-2) +v(k-1) ;%C(zl)=l,即取 cl=0, c2=0endz(2)=0

16、;z(l)=0; %设2的前两个初始值为零for k=3:301; %循环变量从3到301z(k)=-l. 376*z(k-l)-0. 483*z(k-2)+57*u(k-l)+0. 42*u(k-2)+ze(k-l); %输出采 样信号(测量值)end%RGLS广义最小二乘辨识C0二0.0001 0. 0001 0. 0001 0.0001; %直接给岀被辨识参数的初始值，即一个充分 小的实向量pfO二10飞*eye (4, 4); %直接给出初始状态P0,即一个充分大的实数单位矩阵 ce0=0. 001 0. 0011，;pe0=eye (2, 2);C二cO, zeros (4, 299

17、); %被辨识参数矩阵的初始值及大小 ce=ce0, zeros (2, 299);e二zeros(4, 300); %相对误差的初始值及大小ee=zeros(2,300);s二0;%广义最小二乘递推算法的计算步骤for k=3:300;zf (k)=z (k) +ce(l, k-2)*z (k-l)+ce(2, k-2)*z(k-2);uf (k)=u(k) +ce(l, k-2)*u (k-l)+ce(2, k-2)*u(k-2);hfl=-zf(k-l), -zf(k-2), uf(k-1), uf(k-2);x=hf r *pf0*hf 1+1; xl=inv (x); %开始求 K(

18、k) kl=pfO*hfl*xl;%求出 K 的值dl=zf(k)-hfr*cO; cl=c0+kl*dl; %求被辨识参数 cel=cl-c0;鸟求参数当前值与上一次的值的差值e2二el./cO; %求参数的相对变化e(：,k)二e2; %把当前相对变化的列向量加入误差矩阵的最后一列c0二cl; %新获得的参数作为下一次递推的旧参数c(：,k)=cl; %把辨识参数c列向量加入辨识参数矩阵的最后一列 pfi=pfo-ki*hfr*pfo； %求出 p(k)的值pfo=pf 1; %给下次用hl=一z(k-l), -z(k-2), u(kl), u(k-2);s=s+(z(k)-hr *1.

19、642 0. 715 0. 39 0. 35) 2%求准则函数 ee(k) =z (k)-hl *cl;hel=-ee(kT),-ee(k-2);x=heT *peO*hel+l; xl=inv (x);kl=peO*hel*xl;dl=ee(k) -heT *ce0;cel=ceO+kl*dl;pel=peO-kl*her *peO;ce0=cel;ce(:, k) =cel;peO=pel;end $大循环结朿c%辨识参数变化矩阵%显示被辨识参数及其误差(收敛)情况%分离参数al=c(1,1:300); a2=c(2, 1:300) ;bl=c(3,1:300) ;b2=c(4, 1:300);cl=ce(1, 1:300);c2=ce(2,1:300);eal=e(l, 1:300); ea2=e (2, 1:300): ebl=e (3, 1:300); eb2=e(4, 1:300);figure(2); %第二个图形i二1:300; %横坐标从1