一、随机信号的参数建模法
随机信号的参数建模法为随机信号建立参数模型是研究随机信号的一种基本方法。
对于平稳随机信号,主要有三种常用的线性模型:AR(Auto-Regression,自回归)模型、MA(Moving Average,滑动平均)模型和ARMA(Auto-Regression-Moving Average,自回归滑动平均)模型。
这里我们选定 AR 模型,用适当的算法估计模型参数 ( a k , b k , p , q ) (a_k,b_k,p,q) (ak,bk,p,q),以便用模型对随机信号进行预测。
二、AR模型参数和自相关函数的关系
由上图参数模型及系统函数得到: x ( n ) = w ( n ) − ∑ k = 1 p a k x ( n − k ) x(n)=w(n)-\sum_{k=1}^{p} a_{k} x(n-k) x(n)=w(n)−k=1∑pakx(n−k)
两边同时乘以 x ( n − m ) x(n−m) x(n−m) 并求均值,得: E [ x ( n ) x ( n − m ) ] = E [ w ( n ) x ( n − m ) − ∑ k = 1 p a k x ( n − k ) x ( n − m ) ] E[x(n) x(n-m)]=E\left[w(n) x(n-m)-\sum_{k=1}^{p} a_{k} x(n-k) x(n-m)\right] E[x(n)x(n−m)]=E[w(n)x(n−m)−k=1∑pakx(n−k)x(n−m)]
根据自相关函数的偶对称特性,可进一步化简为: R x x ( m ) = R x w ( m ) − ∑ k = 1 p a k R x x ( m − k ) R_{x x}(m)=R_{x w}(m)-\sum_{k=1}^{p} a_{k} R_{x x}(m-k) Rxx(m)=Rxw(m)−k=1∑pakRxx(m−k)
根据 x ( n ) = w ( n ) ∗ h ( n ) = ∑ k = 0 ∞ h ( k ) w ( n − k ) x(n)=w(n) * h(n)=\sum_{k=0}^{\infty} h(k) w(n-k) x(n)=w(n)∗h(n)=k=0∑∞h(k)w(n−k)
有: R x w ( m ) = E [ x ( n ) w ( n + m ) ] = E [ ∑ k = 0 ∞ h ( k ) w ( n − k ) w ( n + m ) ] = ∑ k = 0 ∞ h ( k ) E [ w ( n − k ) w ( n + m ) ] = ∑ k = 0 ∞ h ( k ) R w w ( m + k ) = ∑ k = 0 ∞ h ( k ) σ w 2 δ ( m + k ) = σ w 2 h ( − m ) \begin{aligned} R_{x w}(m) &=E[x(n) w(n+m)] \\ &=E\left[\sum_{k=0}^{\infty} h(k) w(n-k) w(n+m)\right] \\ &=\sum_{k=0}^{\infty} h(k) E[w(n-k) w(n+m)] \\ &=\sum_{k=0}^{\infty} h(k) R_{w w}(m+k) \\ &=\sum_{k=0}^{\infty} h(k) \sigma_{w}^{2} \delta(m+k) \\ &=\sigma_{w}^{2} h(-m) \end{aligned} Rxw(m)=E[x(n)w(n+m)]=E[k=0∑∞h(k)w(n−k)w(n+m)]=k=0∑∞h(k)E[w(n−k)w(n+m)]=k=0∑∞h(k)Rww(m+k)=k=0∑∞h(k)σw2δ(m+k)=σw2h(−m)
由于系统的单位脉冲响应 h ( n ) h(n) h(n) 因果,故 R x w ( m ) = { σ w 2 h ( − m ) , m ≤ 0 0 , m > 0 R_{x w}(m)=\left\{\begin{array}{cc} \sigma_{w}^{2} h(-m) & , m \leq 0 \\ 0 & , m>0 \end{array}\right. Rxw(m)={ σw2h(−m)0,m≤0,m>0
带入得: R x x ( m ) = { R x x ( − m ) , m < 0 − ∑ k = 1 p a k R x x ( m − k ) + h ( 0 ) σ w 2 , m = 0 − ∑ k = 1 p a k R x x ( m − k ) , m > 0 R_{x x}(m)=\left\{\begin{array}{ll} R_{x x}(-m) & , m<0 \\ -\sum_{k=1}^{p} a_{k} R_{x x}(m-k)+h(0) \sigma_{w}^{2} & , m=0 \\ -\sum_{k=1}^{p} a_{k} R_{x x}(m-k) & , m>0 \end{array}\right. Rxx(m)=⎩⎨⎧R