Adaptive Filter Learing Notes 自适应滤波学习笔记08 Levinson-Durbin算法

Levinson-Durbin算法是用来解线性预测误差滤波的系数(prediction-error filter coefficients)和预测误差功率(prediction-error power)的递归算法,优点是减少计算量和存储空间

Levinson-Durbin算法的推导

先参考上一篇笔记Adaptive Filter Learing Notes 自适应滤波学习笔记07 线性预测的记号, R m + 1 a m = [ P m 0 ] , R m + 1 a m B ∗ = [ 0 P m ] . \boldsymbol{R}_{m+1}\boldsymbol{a}_{m}=\left[\begin{matrix} P_m\\ \boldsymbol{0} \end{matrix}\right],\quad \boldsymbol{R}_{m+1}\boldsymbol{a}_{m}^{B*}=\left[\begin{matrix} \boldsymbol{0}\\ P_m \end{matrix}\right]. Rm+1am=[Pm0],Rm+1amB=[0Pm]. R m + 1 [ a m − 1 0 ] = [ R m r m B ∗ r m B T r ( 0 ) ] [ a m − 1 0 ] = [ R m a m − 1 r m B T a m − 1 ] = : [ P m − 1 0 Δ m − 1 ] . \begin{aligned} \boldsymbol{R}_{m+1}\left[\begin{matrix} \boldsymbol{a}_{m-1}\\ 0 \end{matrix}\right]=&\left[\begin{matrix} \boldsymbol{R}_{m} & \boldsymbol{r}_{m}^{B*} \\ \boldsymbol{r}_m^{BT} & r(0) \end{matrix}\right]\left[\begin{matrix} \boldsymbol{a}_{m-1}\\ 0 \end{matrix}\right]\\ =&\left[\begin{matrix} \boldsymbol{R}_m\boldsymbol{a}_{m-1}\\ \boldsymbol{r}_{m}^{BT}\boldsymbol{a}_{m-1} \end{matrix}\right]\\ =:&\left[\begin{matrix} P_{m-1}\\ \boldsymbol{0}\\ \Delta_{m-1} \end{matrix}\right]. \end{aligned} Rm+1[am10]===:[RmrmBTrmBr(0)][am10][Rmam1rmBTam1]Pm10Δm1. Δ m − 1 = r m B T a m − 1 = a m − 1 B T r m \Delta_{m-1}=\boldsymbol{r}_{m}^{BT}\boldsymbol{a}_{m-1}=\boldsymbol{a}_{m-1}^{BT}\boldsymbol{r}_m Δm1=rmBTam1=am1BTrm. R m + 1 [ 0 a m − 1 B ∗ ] = [ r ( 0 ) r m H r m R m ] [ 0 a m − 1 B ∗ ] = [ r m H a m − 1 B ∗ R m a m − 1 B ∗ ] = [ Δ m − 1 ∗ 0 P m − 1 ] . \begin{aligned} \boldsymbol{R}_{m+1}\left[\begin{matrix} 0\\ \boldsymbol{a}_{m-1}^{B*} \end{matrix}\right]=&\left[\begin{matrix} r(0) & \boldsymbol{r}_{m}^{H} \\ \boldsymbol{r}_m & \boldsymbol{R}_m \end{matrix}\right]\left[\begin{matrix} 0\\ \boldsymbol{a}_{m-1}^{B*} \end{matrix}\right]\\ =& \left[\begin{matrix} \boldsymbol{r}_m^{H}\boldsymbol{a}_{m-1}^{B*}\\ \boldsymbol{R}_{m}\boldsymbol{a}_{m-1}^{B*} \end{matrix}\right]\\ =&\left[\begin{matrix} \Delta_{m-1}^*\\ \boldsymbol{0}\\ P_{m-1} \end{matrix}\right]. \end{aligned} Rm+1[0am1B]===[r(0)rmrmHRm][0am1B][rmHam1B

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值