主要内容

yulewalk

递归数字滤波器设计

所有的页面崩溃

语法

[b] = yulewalk (n、f m)

描述

例子

(b,一个)= yulewalk (n,f,米)返回的传递函数系数n级阶IIR滤波器的频率响应近似匹配给定的值f和米。

例子

全部折叠

打开生活的脚本

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

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

图包含一个坐标轴对象。坐标轴对象包含ω/π,ylabel级包含2线类型的对象。

增加阻带衰减通过指定一个更广泛的过渡带。

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

图包含一个坐标轴对象。坐标轴对象包含ω/π,ylabel级包含3线类型的对象。

输入参数

全部折叠

过滤器订单,指定为一个正整数标量。

数据类型:单|双

频率点,指定为一个向量点的范围在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。波阿斯,波拉特。“修改Yule-Walker ARMA谱估计的方法。”IEEE®在航空电子系统交易。AES-20卷,第二,1984年,页158 - 173。

扩展功能

版本历史

之前介绍过的R2006a

另请参阅

|||||||