控制相关专业研究生选修课程系统建模方法马宏军东北大学信息学院控制理论与导航技术研究所2013年3月逸夫楼203第9章系统建模灰箱方法1、问题的提出热敏电阻的测量值t)(C203251738895R)(76582687394210101032例:表中是在不同温度下测量同一热敏电阻的阻值,根据测量值确定该电阻的数学模型。vbtaR1、问题的提出辨识目的:根据过程所提供的测量信息,在某种准则意义下,估计模型的未知参数。ProcessInputOutputvbtaRt(℃)203251738895R(Ω)76582687394210101032a,b1、问题的提出辨识目的:根据过程所提供的测量信息,在某种准则意义下,估计模型的未知参数。ProcessInputOutput工程实践目的模型结构参数辨识模型校验模型确定m次独立试验的数据),(11zt),(22zt),(mmzt)()()()(22110thathathaatynn)(kG)(kt)(ky)(kv)(kz)(kG)(kt)(ky1、问题的提出)()()(kvkykzm次独立试验的数据),(11zt),(22zt),(mmzt)()()()()(22110kvkhakhakhaakznnzt)(tf)(kG)(kt)(ky)(kv)(kz1、问题的提出未知量的最可能值是使各项实际观测值和计算值之间差的平方乘以其精确度的数值以后的和为最小。1795年,高斯提出的最小二乘的基本原理是1、问题的提出Gauss(1777-1855))()()(kvkykz使最小mkkykzkw12|)()(|)(2、最小二乘辨识方法的基本概念通过试验确定热敏电阻阻值和温度间的关系t)(C1t2t1NtNtR)(1R2R1NRNR•当测量没有任何误差时,仅需2个测量值。•每次测量总是存在随机误差。btaRiiivRy或iivbtayiiiiiibtayvRyv=或2.1利用最小二乘法求模型参数根据最小二乘的准则有NiiiNiibtaRvJ1212min)]([根据求极值的方法,对上式求导NiiiibbNiiiaatbtaRbJbtaRaJ1ˆ1ˆ0)(20)(2NiiiibbNiiiaatbtaRbJbtaRaJ1ˆ1ˆ0)(20)(2NiNiiiNiiiNiNiiitRtbtaRtbaN111211ˆˆˆˆ2112111211211121ˆˆNiiNiiNiiNiiNiiiNiiNiiNiiNiiiNiiNiittNtRtRNbttNttRtRa例:表1中是在不同温度下测量同一热敏电阻的阻值,根据测量值确定该电阻的数学模型,并求出当温度在C70时的电阻值。表1热敏电阻的测量值t)(C20.52632.740516173808895.7R)(765790826850873910942980101010322.1利用最小二乘法求模型参数表1热敏电阻的测量值t)(C20.52632.740516173808895.7R)(76579082685087391094298010101032762.702ˆa4344.3ˆbCt70168.943R2112111211211121ˆˆNiiNiiNiiNiiNiiiNiiNiiNiiNiiiNiiNiittNtRtRNbttNttRtRabtaR2.2一般最小二乘法原理及算法)(kG)(ku)(ky)(kv)(kz图3.4SISO系统的“灰箱”结构nnnnzazazazbzbzbzuzyzG221122111)()()(niiniiikubikyaky11)()()(2.2一般最小二乘法原理及算法)(zG)(ku)(ky)(kv)(kz图SISO系统的“灰箱”结构若考虑被辨识系统或观测信息中含有噪声)()()()(11kvikubikyakzniinii)(kz为系统输出量的第k次观测值;)(ky为系统输出量的第k次真值;)(ku为系统的第k个输入值;)(kv是均值为0的随机噪声。)()()()(11kvikubikyakzniinii)](,),2(),1(),(,),2(),1([)(nkukukunkykykykh如果定义Tnnbbbaaa],,,,,,,[2121)()()(kvkhkz式中为待估参数。2.2一般最小二乘法原理及算法)()()(kvkhkz令mk,,2,1,则有)()2()1(mzzzZm)()1()()1()2()1()2()1()1()0()1()0()()2()1(nmumunmymynuunyynuunyymhhhHmTnnbbaa11TmmvvvV)()2()1(2.2一般最小二乘法原理及算法mmmVHZ)(zG)(ku)(ky)(kv)(kz图SISO系统的“灰箱”结构最小二乘的思想就是寻找一个的估计值ˆ,使得各次测量的),1(miZi与由估计ˆ确定的量测估计ˆˆiiHZ之差的平方和最小,即2.2一般最小二乘法原理及算法min)ˆ()ˆ()ˆ(mmTmmHZHZJ0)ˆ(2ˆmmTmHZHJmTmmTmZHHHˆ如果mH的行数大于等于列数,即nm2,mTmHH满秩,即nHHmTm2)(rank,则1)(mTmHH存在。则的最小二乘估计为2.2一般最小二乘法原理及算法mTmmTmZHHH1)(ˆ最小二乘估计虽然不能满足式mmmVHZ中的每一个方程,使每个方程都有偏差,但它使所有方程偏差的平方和达到最小,兼顾了所有方程的近似程度,使整体误差达到最小,这对抑制测量误差),,1)((miiv是有益的。2.2一般最小二乘法原理及算法mTmmTmZHHH1)(ˆ最小二乘估计虽然不能满足式mmmVHZ中的每一个方程,使每个方程都有偏差,但它使所有方程偏差的平方和达到最小,兼顾了所有方程的近似程度,使整体误差达到最小,这对抑制测量误差),,1)((miiv是有益的。mmmVHZzt)(tf2.2一般最小二乘法原理及算法最小二乘法的几何解释mmmVHZ)()2()1(mzzzZm)()1()()1()2()1()2()1()1()0()1()0()()2()1(nmumunmymynuunyynuunyymhhhHmTnnbbaa11的线性组合维空间中基向量是)}(),2(),1({mhhhmZm的近似是在最小二乘意义下对mmZHˆ的张成的空间的投影。在应该等于)}(),2(),1({ˆmhhhZHmm2.2一般最小二乘法原理及算法最小二乘法的几何解释222VHZ2Z)1(h)2(hVˆˆ2H)2()1(2zzZTnnbbaa11)2()1(2hhH当系统的量测噪声mV是均值为0,方差为R的随机向量,则最小二乘估计有如下性质。2.2一般最小二乘法原理及算法(1)最小二乘估计是无偏估计(无偏性))ˆ(E或0)~(E(2)最小二乘估计是有效估计(有效性)11)()()~~(mTmmTmmTmTHHRHHHHE(3)最小二乘估计是一致估计(一致性)0)|ˆ(|limmmp最小(1)最小二乘估计是无偏估计,即)ˆ(E或0)~(E证明:mTmmTmZHHH1)(ˆ~])([)ˆ()~(1mTmmTmVHHHEEE)()(1mmTmmTmZHHHHmTmmTmVHHH1)(0mTmmTmmTmmTmZHHHHHHH11)()()()()(1mTmmTmVEHHH2.2一般最小二乘法原理及算法如果由测量噪声及模型误差等引起的误差V的均值为0,且V与输入矢量Hm是统计独立,最小二乘的估计值是无偏的。(2)最小二乘估计为有效估计。11)()()~~(mTmmTmmTmTHHRHHHHE证明:11)()()()~~(mTmmTmmTmmTmTHHHVVEHHHE11)()(mTmmTmmTmHHRHHHH根据第(1)式的证明,显然有2.2一般最小二乘法原理及算法112)()(mTmmTmmTmHHHHHH12)(mTmHH独立随机变量同分布、零均值、中的各个量是mVIR2(3)最小二乘估计是一致估计,其估计值依概率收敛于真值,即证明:2.2一般最小二乘法原理及算法0)|ˆ(|limmmp最小二乘估计ˆ依概率收敛于。等价于随着m,0)~~(TE12121)(limlimmTmmmTmmHHmmHH0)(12limmTmmHH奇异常数阵是非mTmHH例3.2用2台仪器对未知标量各直接测量一次,量测量分别为1z和2z,仪器的测量误差均值为0,方差分别为r和r4的随机量,求的最小二乘估计,并计算估计的均方误差。2.2一般最小二乘法原理及算法解:由题意得量测方程222VHZ212zzZ112HrrR400)(21111111ˆ21211zzzzrrrrET45111111400111111)~~(112.2一般最小二乘法原理及算法222VHZ212zzZ112HrrR4002.3加权最小二乘法原理及算法一般最小二乘估计精度不高的原因之一是对测量数据同等对待各次测量数据很难在相同的条件下获得的有的测量值置信度高,有的测量值置信度低对不同置信度的测量值采用加权的办法分别对待置信度高的,权重取得大些;置信度低的,权重取的小些min)ˆ()ˆ()ˆ(mmTmmHZHZJmin)ˆ()ˆ()ˆ(mmmTmmHZWHZJ])(,),2(),1(diag[mmin)ˆ()ˆ()ˆ(mmmTmmHZWHZJ2.3加权最小二乘法原理及算法0)ˆ(2ˆmmmTmHZWHJmmTmmmTmZWHHWH1)(ˆ当系统的量测噪声mV是均值为0,方差为R的随机向量,则加权最小二乘估计有如下性质。2.2加权最小二乘法原理及算法(1)加权最小二乘估计是无偏估计(无偏性))ˆ(E或0)~(E(2)加权最小二乘估计是有效估计(有效性)(3)最小二乘估计是一致估计(一致性)0)|ˆ(|limmmp1)()~~(mTmTRHHE当系统的量测噪声mV是均值为0,方差为R的随机向量,则加权最小二乘估计有如下性质。(1)加权最小二乘估计是无偏估计,即)ˆ(E或0)~(E(2)加权最小二乘