主要内容

yulewalk

递归数字滤波器设计

折叠所有页面

语法

[b,a] = yulewalk(n,f,m)

描述

例子

［b，一个= yulewalk(n，f，米）返回an的传递函数系数n阶IIR滤波器，其频幅响应近似匹配式中给出的值f而且米．

例子

全部折叠

打开实时脚本

设计一个归一化截止频率为0.6的8阶低通滤波器。画出它的频率响应，并叠加相应的理想滤波器的响应。

F = [0 0.6 0.6 1];M = [1 1 0 0];[b,a] = yulewalk(8,f,m);[h,w] = freqz(b,a,128);Plot (w/pi,mag2db(abs(h))) yl = ylim;持有在情节(f(2:3)、黄“——”)包含(“ω\ / \π”) ylabel (“级”网格)

图中包含一个轴对象。axis对象包含2个line类型的对象。

通过指定更宽的过渡带来增加阻带衰减。

F = [0 0.55 0.6 0.65 1];M = [1 1 0.5 0 0];[b,a] = yulewalk(8,f,m);H = freqz(b,a,128);持有在情节(w /π,mag2db (abs (h)))从ylim(黄)

图中包含一个轴对象。axis对象包含3个line类型的对象。

输入参数

全部折叠

筛选器顺序，指定为正整数标量。

数据类型:单|双

频率点，指定为0到1之间范围内的点的向量，其中1对应奈奎斯特频率，或采样率的一半。第一点f必须是0，最后一点是1。所有中间点必须按递增顺序排列。f可以有重复的频率点对应的阶跃频率响应。

例子:[0 0.25 0.4 0.5 0.5 0.7 1]指定不规则抽样的奈奎斯特范围。

数据类型:单|双

大小响应，指定为一个矢量，其中包含在所指定的点上所需的响应f．米长度必须和f．

例子:[0 1 1 1 0 0 0]指定带通幅度响应。

数据类型:单|双

输出参数

全部折叠

过滤系数，作为行向量返回。输出滤波器系数按的降幂排序z：

B （ z ） 一个 （ z ） ＝ b （ 1 ） + b （ 2 ） z − 1 + ⋯ + b （ n + 1 ） z − n 一个 （ 1 ） + 一个 （ 2 ） z − 1 + ⋯ + 一个 （ n + 1 ） z − n ．

提示

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

算法

yulewalk设计递归IIR数字滤波器使用最小二乘拟合到指定的频率响应。函数在时域内进行拟合。

  • 为了计算分母系数，yulewalk采用修正的Yule-Walker方程，相关系数由指定频率响应的傅里叶反变换计算。

  • 为了计算分子，yulewalk以下步骤:

    1. 计算与功率频率响应的加性分解相对应的分子多项式。

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

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

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

参考文献

[1]弗里德兰德，B.和波阿斯·波拉特。ARMA谱估计的改进Yule-Walker方法。IEEE®航空航天电子系统汇刊．卷AES-20，第2期，1984，第158-173页。

扩展功能

版本历史

R2006a之前介绍

另请参阅

|||||||

Baidu
map