主要内容

firpmord

帕克斯-麦克莱伦最优FIR滤波器阶估计

描述

naow) = firpmord (f一个dev返回近似顺序n,归一化频带边缘,频带振幅ao,重量w满足输入规范f一个,dev

例子

___) = firpmord (___fs指定采样频率fsfs默认为2hz,意味着Nyquist频率为1hz。您可以指定按特定应用程序的采样率缩放的频带边缘。您可以将它与前面的任何输入语法一起使用。

例子

c= firpmord (___“细胞”返回单元格数组c谁的元素是参数firpm

例子

全部折叠

设计一个最小阶低通滤波器,通带截止频率为500hz,阻带截止频率为600hz。指定采样频率为2000hz。要求阻带至少40 dB的衰减,通带纹波小于3 dB。

rp = 3;通频带纹波百分比,单位为dBrs = 40;阻带纹波百分比,单位为dBfs = 2000;%采样频率F = [500 600];%截止频率A = [1 0];%的振幅

将偏差转换为线性单位。设计滤波器并可视化其幅值和相位响应。

dev = [(10 ^ (rp / 20) 1) / (10 ^ (rp / 20) + 1) 10 ^ (- r / 20)];[n, fo ao w] = firpmord (f, a, dev, fs);b = firpm (n, fo ao w);freqz (b, 1024 fs)标题(“按规格设计的低通滤波器”

图中包含2个轴。标题为“设计到规格的低通滤波器”的轴1包含一个类型为line的对象。Axes 2包含一个类型为line的对象。

滤波器略低于满足阻带衰减和通带纹波规格。使用n + 1而不是n在号召中firpm达到所需的振幅特性。

设计一个通带截止频率为1500hz、阻带截止频率为2000hz的低通滤波器。请指定采样频率为8000hz。要求最大阻带振幅为0.1,最大通带误差(纹波)为0.01。

[n, fo ao w] = firpmord (2000 [1500], [1 0], [0.01 - 0.1], 8000);b = firpm (n, fo ao w);

通过…得到相等的结果firpmord生成单元格数组。可视化滤波器的频率响应。

C = firpmord([1500 2000],[1 0],[0.01 0.1],8000,“细胞”);B = firpm (c {:});freqz (B, 1024, 8000)

图中包含2个轴。Axes 1包含一个类型为line的对象。Axes 2包含一个类型为line的对象。

输入参数

全部折叠

频带边缘,指定为实值向量。参数必须在范围[0,F年代/ 2)F年代为奈奎斯特频率。向量中的元素数总是2的倍数。频率必须是递增的。

所包含的点的期望振幅f,指定为向量。f一个必须满足条件长度(f) = 2长度(一个) 2.期望的函数是分段常数。

最大允许偏差,指定为向量。dev有相同的尺寸一个.它指定每个频带的频率响应和输出滤波器的期望振幅之间的最大允许偏差或波纹。

采样率,指定为实标量。

输出参数

全部折叠

筛选顺序,作为正整数返回。

归一化频率点,指定为实值向量。参数必须在[0,1]范围内,其中1对应于Nyquist频率。向量中的元素数总是2的倍数。频率必须是递增的。

振幅响应,返回实值向量。

权值用于调整每个频带的适合度,指定为实值向量。的长度w长度的一半是f一个,所以每条带子只有一个重量。

FIR滤波器参数,作为单元格数组返回。

算法

firpmord使用中建议的算法[1].该函数对接近0或Nyquist频率的频带边缘产生不准确的结果,fs / 2

请注意

在某些情况下,firpmord低估或高估订单n.如果过滤器不符合规格,请尝试更高的顺序,如n + 1n + 2

参考文献

[1] Rabiner, Lawrence R.和Otto Herrmann。某些最佳有限脉冲响应数字滤波器的可预测性IEEE®电路理论汇刊.第20卷第4期,1973年,401-408页。

[2] Rabiner, Lawrence R.和Bernard Gold。数字信号处理理论与应用。Englewood Cliffs, NJ: Prentice-Hall, 1975, pp. 156-157。

扩展功能

之前介绍过的R2006a