i =1 i =0 i =1
n1
计算不同n1、n2、n3时的AIC值,取最小的AIC值对应的n1、 n2、n3值为系统的阶次。
8.3 按残差白色定阶
定阶原理: 定阶原理:若阶次n设计合适,则残差近似为白噪声。因此可 利用计算残差e(k)的自相关函数来检查白色性。 自相关函数的计算如下:
1 ˆ R (i ) = N
H( l,k ) = 0
但存在噪声则无上述结论,因此定义指标:
D= H (l , k )行列式的平均值 H (l + 1, k )行列式的平均值
当D达到极大时L值即为系统阶次n。
另一种求D 另一种求D的方法
计算脉冲响应序列的自相关值:
N −i 1 R g (i ) = ∑ g k g k −i N − i + 1 k =0
ε(k)为服从正态分布的白噪声, 经推导,得:
ˆ AIC = N ln σ ε2 + 2 ( n1 + n 2 + n 3 )
式中:
1 ˆ σε = N
2
∑ εˆ
k =1
N
2
(k )
n2 n3
ˆ ˆ ˆ ˆˆ ε (k) = y(k) + ∑ai y(k − i) − ∑bi u(k − i) − ∑ciε (k − i)
由
∂ lnL ˆ = 0 ⇒ θ = ( Φ T Φ ) −1 Φ T Y ∂θ ∂ lnL 1 T 2 ˆ ˆ ˆ = 0 ⇒ σe = e e 2 N ∂ σe
⇒ lnL( | θ) = − Y
N 2 lnσe + const 2
N 2 Y lnL( | θ) = − lnσe + const 2 由 p = n1+ n2 + 2