维纳滤波的线性预测

1.维纳滤波线性预测概念

线性预测(LP)

{u⃗(n)=[u(n−1)u(n−2)⋯u(n−M)]d(n)=u(n) \begin{align} \left\{ \begin{aligned} &\vec{u}(n)=\begin{bmatrix}u(n-1)& u(n-2)& \cdots& u(n-M)\end{bmatrix}\\ &d(n)=u(n) \end{aligned} \right. \end{align} {u(n)=[u(n1)u(n2)u(nM)]d(n)=u(n)

前后向线性预测(FBLP)

(1)前向线性预测(FLP){u⃗(n)=[u(n−1)u(n−2)⋯u(n−M)]d(n)=u(n)(2)后向线性预测(BLP){u⃗(n)=[u(n−M+1)u(n−M+2)⋯u(n)]d(n)=u(n−M) \begin{align} &(1)前向线性预测(FLP)\\ &\left\{ \begin{aligned} &\vec{u}(n)=\begin{bmatrix}u(n-1)& u(n-2)& \cdots& u(n-M)\end{bmatrix}\\ &d(n)=u(n) \end{aligned} \right.\\ &(2)后向线性预测(BLP)\\ &\left\{ \begin{aligned} &\vec{u}(n)=\begin{bmatrix}u(n-M+1)& u(n-M+2)& \cdots& u(n)\end{bmatrix}\\ &d(n)=u(n-M) \end{aligned} \right. \end{align} (1)前向线性预测(FLP){u(n)=[u(n1)u(n2)u(nM)]d(n)=u(n)(2)后向线性预测(BLP){u(n)=[u(nM+1)u(nM+2)u(n)]d(n)=u(nM)

2.线性预测与AR模型

线性预测与AR模型互为逆系统(本质是线性预测误差作为输出时与AR模型互为逆系统)

(1)输入输出关系

{v(n)=u(n)−w⃗oHu⃗(n)(v(n)=e(n)=u(n)−d^(n))u(n)=−∑k=1Maku(n−k)+v(n) \begin{align} \left\{ \begin{aligned} &v(n)=u(n)-\vec{w}_o^H\vec{u}(n)(v(n)=e(n)=u(n)-\hat{d}(n))\\ &u(n)=-\sum\limits_{k=1}^{M}a_ku(n-k)+v(n) \end{aligned} \right. \end{align} v(n)=u(n)woHu(n)(v(n)=e(n)=u(n)d^(n))u(n)=k=1Maku(nk)+v(n)

(2)系统函数关系

HLP(z)HAR(z)=1{HLP(z)=1+a1z−1+a2z−1+⋯+aMz−MHAR(z)=11+a1z−1+a2z−1+⋯+aMz−M \begin{align} &H_{LP}(z)H_{AR}(z)=1\\ &\left\{ \begin{aligned} &H_{LP}(z)=1+a_1z^{-1}+a_2z^{-1}+\cdots+a_Mz^{-M}\\ &H_{AR}(z)=\frac{1}{1+a_1z^{-1}+a_2z^{-1}+\cdots+a_Mz^{-M}}\\ \end{aligned} \right. \end{align} HLP(z)HAR(z)=1HLP(z)=1+a1z1+a2z1++aMzMHAR(z)=1+a1z1+a2z1++aMzM1

(3)滤波器权向量与模型参数关系

ai=−wi−1∗ a_i=-w_{i-1}^* ai=wi1

线性滤波器是白化滤波器

若线性滤波器输入时AR(M)过程,则其输出为均值为零,方差为Jmin=σv2J_{min}=\sigma_v^2Jmin=σv2的白噪声
![[655FD1CD304F9837A0C31E79DC70322C.png]]

3.前后向线性预测

(1)前向预测(FLP))和后向预测(BLP)合并使用,且其滤波器权重满足wfo=wbow_{fo}=w_{bo}wfo=wbo
(2)前后向线性预测可由格型滤波器实现,需要确定预测误差emf(n),emb(n)e^f_m(n),e^b_m(n)emf(n),emb(n)以及反射系数κm\kappa_mκm
(3)Burg算法

a.计算预测误差emf(n),emb(n)e^f_m(n),e^b_m(n)emf(n),emb(n)以及反射系数κm\kappa_mκm
b.计算LP(M)LP(M)LP(M)
c.计算AR(M)AR(M)AR(M)
d.计算AR功率谱估计结果

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值