系统辨识
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
系统辨识
1 最小二乘参数估计法
（1）批处理最小二乘法 CAR 模型：
()()()
()()11A z y k B z u k d k ξ--=-+
(1.1)
式中，()k ξ为白噪声，结构参数a n 、b n 和d 均已知，且：
()()
11111
011a
a b b n n n n A z a z a z B z b b z b z ------⎧=+++⎪
⎨=+++⎪⎩
参数估计的任务是根据可测量的输入和输出，确定如下1a b n n ++个参数：
1201,,,,,,,a b n n a a a b b b
将式（1.1）写成最小二乘形式：
()()()()()()
()()
101 a b n a n b T y k a y k a y k n b u k d b u k d n k k k ξϕθξ=----+-++--+=+ (1.2)
式中，()k ϕ为数据向量，θ为待估计参数向量，且：
()()()()()(
)()11
11
101,,,,
,,,,,,a b a b a b T
n n a b T
n n n n k y k y k n u k d u k d n R
a a
b b R ϕθ++⨯++⨯⎧⎡⎤=-------∈⎣⎦⎪
⎨⎡⎤⎪=∈⎣⎦⎩
LS 法：
对于L 次观测，取如下性能指标：
()()()()
2
2
1
1
ˆˆˆˆ2L
L
T
T T T T T k k J k y k k Y Y Y
ε
ϕθθ
θθ==⎡⎤=
=-=-Φ+ΦΦ⎣⎦
∑
∑ (1.3)
参数的最小二乘估计就是是目标函数取极小值的参数ˆθ。
现有L 组输入输出观测数据()(){},,1,,y k u k k L =，利用批处理法得到系统参数的最小二乘估计ˆθ为：
()
1
ˆT T Y θ
-=ΦΦΦ (1.4)
式中
()()()()()()()1
11122,a b T T L n n L T y y Y R R y L L ⨯++⨯⎡⎤Φ⎡⎤
⎢⎥⎢⎥
Φ⎢⎥⎢⎥=∈Φ=∈⎢⎥⎢⎥
⎢⎥⎢⎥
⎢⎥⎢⎥Φ⎣⎦
（2）递推最小二乘参数估计
解决处理数据量大，占用内存大的问题 递推公式：
()()()()()()()()()()()()()()()()ˆˆˆ1111111T T
T k k K k y k k k P
k k K k k P k k P k K k k P k θθϕθϕϕϕϕ⎧⎡⎤
=-+--⎣⎦⎪
⎪-⎪
=⎨+-⎪
⎪
⎡⎤=--⎪⎣
⎦⎩ (1.5)
初值()0P 、()ˆ0θ
有如下两种方法： ①若已取得L 组数据（1a b L n n >++），利用批处理最小二乘参数估算法：
()
()()()
1
1
T L L T
T L L L L
P L L Y θ--⎧=ΦΦ
⎪⎨
⎪=ΦΦΦ⎩
②直接令
()()0ˆ
0P I
αθε⎧=⎪⎨=⎪
（3）遗忘因子递推最小二乘法 针对参数时变系统 性能指标：
()()2
1
ˆL
L k T
k J y k k λϕθ-=⎡⎤=
-⎣⎦∑
(1.6)
递推公式：
()()()()()()()()()()()()()()()()ˆˆˆ11111
11T T
T k k K k y k k k P k k K k k P k k P k K k k P k θθϕθϕλϕϕϕλ⎧
⎡⎤=-+--⎪⎣⎦⎪
-⎪
=⎨+-⎪
⎪⎡⎤⎪=--⎣⎦⎩
(1.7)
（4）递推增广最小二乘法 增加白噪声
()k ξ不可测，用其估计值代替： ()()()()()ˆˆˆT k y k y
k y k k ξϕθ=-=- (1.8)
式中：
()()()()()()()()()1111101ˆˆ1,,,,,,1,,ˆˆ,,,,,,,,a b c a b c a b c T
n n n a b c T n n n n n n k y k y k n u k d u k d n k k n R a a b b c c R ϕξξθ+++⨯+++⨯⎧⎡⎤=---------∈⎪⎣⎦
⎨⎪⎡⎤=∈⎣⎦⎩
递推公式：
()()()()()()()()()()()()()()()()ˆˆˆˆ11ˆ1ˆˆ11ˆ11T T
T k k K k y k k k P k k K k k P k k P k K k k P k θθϕθ
ϕϕϕϕ⎧⎡⎤=-+--⎣⎦⎪
⎪-⎪
=⎨+-⎪
⎪
⎡⎤=--⎪
(1.9)
2 梯度矫正参数估计法
相比最小二乘法，计算量显著减小，但是收敛变慢
基本思想：沿着准则函数（目标函数）的负梯度方向，逐步修正模型参数估计值，直至准则函数达到最小值。其简单易懂、实时计算量小（没有标量运算），但是收敛速度比较慢。
（1）确定性系统的梯度矫正参数估计法（适用于不含噪声的系统）
()()()()()
()()()ˆˆˆ11T T k k k y k k k c k k αϕθθϕθϕϕ⎡⎤=-+--⎣⎦+ (1.10) 式中，0c >，02α<<。