中科院研究生院信息工程学院课件系统辨识第六讲

合集下载
相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《系统辨识》第6讲要点

第4章 线性动态模型参数辨识-最小二乘法

4.10.3 限定记忆法

限定记忆法依赖于有限长度的数据,每增加一个新的数据信息,就要去掉一个老数据的信息,数据长度始终保持不变。这种方法的参数估计递推算法如下:

[]

[]

⎪⎪⎪⎪⎩⎪⎪

⎪⎪⎨⎧-+++=++-++++-+=+-++-+++-+=++++=+++-+=+++-++-+=++--.

)1,()](),([),()

()1,()(1)()1,(),()]1,(ˆ

)()()[,()1,(ˆ),(ˆ),()](),1([),1()(),()(1)(),(),1()],(ˆ)()()[,1(),(ˆ),1(ˆ11L k k L k L k k L k k L k L k k L k L k L k k L k k L k k L k L k z L k k L k k L k k L k k k L k k L k k k L k k k k L k k L k k L k k k k z L k k L k k L k k P h K P h P h h P K h K P h K P h P h h P K h K ττ

ττττ

θθθθθθ-I +I 算法前三个式子用于去掉老数据的信息,后三个式子用来增加新数据的信

息,初始值取

P (,)

(,),00002==⎧⎨⎩

a I ,

θε 其中a 为充分大实数,ε 为充分小实向量.相应的准则函数递推计算式为:]

,)

()1,()(1)

(~)

(),()(1)

(~)1,(),1(2221L k L k k L k L k z k L k k k k z L k k J L k k J +-++++++--

-+=++h P h h P h τ

τ

其中

⎪⎩⎪⎨⎧-++-+=++-=)1,(ˆ)()()(~)

,(ˆ)()()(~21L k k L k L k z L k z L k k k k z k z θθττ

h h 4.10.4 折息法

折息法把加权最小二乘法和遗忘因子法融合起来,形成如下算法:

⎪⎪

⎪⎪

⎪⎨⎧--=⎥⎦⎤⎢⎣⎡+--=--+-=-)1()]()([)(1)()()()()1()()()1()()]

1(ˆ

)()()[()1(ˆ)(ˆ1k k k k k k k k k k k k k k k k z k k k P h K P h P h h P K h K τττθθθI μμΛ 折息因子与加权因子和遗忘因子之间的关系为ΓΛ(,)()()k i i j j i k

==+∏μ1

当遗忘因子取常数时,折息因子又可表示成ΓΛ(,)()k i i k i =-μ。折息法同时

具备加权最小二乘法和遗忘因子法的作用,既可获得系统的平均特性,又具有时变跟踪能力。

4.10.5 协方差重调最小二乘法

在辨识递推计算过程中,协方差矩阵P (k )衰减很快,此时算法的增益矩阵K (k )也急剧衰减。这种现象的出现,促使人们去考虑一种修正的方案,即在指定的时刻重新调整协方差矩阵P (k ),使算法始终保持较快的收敛速度。这种协方差重调的最小二乘算法描述如下:

[]

⎪⎪⎩

⎪⎪⎨

⎧--=+--=--+-=-)1()]()([)(1)()1()()()1()()]

1(ˆ)()()[()1(ˆ)(ˆ1T

k k k k k k k k k k k k k z k k k P h K P h P h h P K h K ττ

θθθI 当k k k k l ∉{,,,}12 时,P (k )按上式算法计算;当k k k k k i l =∈{,,,}12 时,把P (k )重调为P ()k a i i =I , 0

4.10.6 协方差修正最小二乘法

对时变系统辨识来说,为了防止矩阵P (k )趋于零,当参数估计值超过某阀值时,矩阵P (k )自动加上附加项Q , 具体算法如下:

[]

⎪⎩⎪⎪⎨

⎧≥+=--=+--=--+-=-0

I Q Q P P P h K P h P h h P K h K ,)()()1()]()([)(1)()1()()()1()()]1(ˆ)()()[()1(ˆ)(ˆ1

k k k k k k k k k k k k k k k z k k k τ

ττ

θθθ

引言

最小二乘法是一种最基本的辨识方法,但如果模型的噪声不是白噪声,最小二乘法不一定能给出无偏、一致估计。以下着重讨论模型噪声是有色噪声时的各种最小二乘辨识方法。

4.11 增广最小二乘法 4.11.1 增广最小二乘原理 考虑如下模型

A z z k

B z u k N z v k ()()()()()()---=+111

式中u (k )和z (k ) 分别为模型输入和输出变量;v (k ) 是均值为零、方差为σv 2的不相关随机噪声或称白噪声;N z ()-1为噪声模型;A z ()-1 和B z ()-1为迟延算子多项式,记作

A z a z a z a z

B z b z b z b z n n

n n a a

b b

()()--------=++++=+++⎧⎨⎪⎩

⎪11122111221 其中n a 和n b 为模型阶次。为了运用最小二乘原理来辨识这种模型的参数,

需要把模型(4-1)式写成最小二乘格式

)()()(k v k k z +=θτh

这样就必须把噪声模型的参数包含在参数向量θ 中,从而引出增广概念,用来构造参数向量θ 和数据向量h(k ),具体的构成形式会因噪声模型的结构不同而不同。下面是三种不同噪声模型的向量构成方法: ① 若N z D z d z d z d z n n d d ()()==++++----111221 ,可按下式构成参数向量和数据向量:

⎪⎩⎪⎪⎨⎧--------+==--------=∑∑∑===)(ˆ

)1(ˆ)()1(ˆ)()1(ˆ)()(ˆ]

,,,,,,,,[)](ˆ,),1(ˆ),(,),1(),(,),1([)(111111i k v k d i k u k b i k z k a k z k v d d b b a a n k v k v n k u k u n k z k z k d b a d b a n i i n i i n i i n n n d b

a ττθ h ② 若N z C z c z c z c z n n c c

()()=

=

++++----11111122 ,参数向量和数据向量的构成形式为:

相关文档
最新文档