文档

yulewalk

递归数字滤波器设计

语法

[b] = yulewalk (n、f m)

描述

yulewalk设计递归的IIR数字滤波器使用最小二乘适合一个特定的频率响应。

[b] = yulewalk (n、f m)返回行向量,b一个,包含n + 1阶系数n其频率幅值特性近似匹配矢量给出的IIR滤波器f

  • f是频率点的矢量,取值范围在0到1之间,其中1对应采样频率(奈奎斯特频率)的一半。第一点f必须是0,最后一点是1。所有中间点的顺序必须递增。允许重复频率点,对应于频率响应中的阶跃。

  • 一个向量是否包含在指定的点上所需的幅值响应f

  • f必须是相同的长度。

  • 情节(f,米)显示筛选器形状。

输出滤波器系数按的降序幂排列z

B z ) 一个 z ) b 1 ) + b 2 ) z 1 + + b n + 1 ) z n 一个 1 ) + 一个 2 ) z 1 + + 一个 n + 1 ) z n

当指定频率响应时,要避免从通带到阻带的过渡过于尖锐。您可能需要试验过渡区域的斜率,以获得最佳的滤波器设计。

例子

全部折叠

设计一个归一化截止频率为0.6的8阶低通滤波器。绘制其频率响应并覆盖相应理想滤波器的响应。

F = [0 0.6 0.6 1];M = [1 1 0 0];[b] = yulewalk (8 f m);[h, w] = freqz (b, a, 128);情节(w /π,abs (h)、f m,“——”)包含的弧度频率ω(\ / \π)', ylabel传奇(“Yule-Walker”“理想”),传说boxoff

算法

yulewalk在时域内进行最小二乘拟合。它使用修正的Yule-Walker方程计算分母系数,相关系数由特定频率响应的傅里叶反变换计算。为了计算分子,yulewalk采取以下步骤:

  1. 计算对应于工频响应加性分解的分子多项式。

  2. 计算分子和分母多项式对应的完整频率响应。

  3. 利用频谱分解技术得到滤波器的脉冲响应。

  4. 通过最小二乘拟合此脉冲响应得到分子多项式。

参考文献

[1] Friedlander, B.和Boaz Porat。ARMA光谱估计的改进Yule-Walker方法IEEE®航空航天电子系统汇刊。1984年第2期,第158-173页。

扩展功能

之前介绍过的R2006a

这个话题有用吗?