fir1：基于窗函数法的 FIR 滤波器设计

b = fir1(n,Wn) b = fir1(n,Wn,ftype)

b = fir1(___,window) b = fir1(___,scaleopt)

b = fir1(n,Wn) 使用 Hamming 窗设计具有线性相位的 n 阶低通，带通或多频带 FIR 滤波器。滤波器的类型取决于 Wn 的元素数。举例

b = fir1(n,Wn,ftype) 根据 ftype 的值和 Wn 的元素数，设计一个低通，高通，带通，带阻或多频带滤波器。举例

b = fir1(___,window) 使用 window 中指定的向量和上述语法中的任何参数来设计滤波器。举例

b = fir1(___,scaleopt) 额外指定滤波器的幅度响应是否已归一化。

注意： 使用 fir2 设计具有任意频率响应的加窗滤波器。

设计一个通带为 0.35 π ⩽ ω ⩽ \pi\leqslant \omega \leqslant π⩽ω⩽ 0.65 π \pi π rad / sample 的 48 阶 FIR 带通滤波器。将其幅度和相位响应可视化。

加载 chirp.mat。 该文件包含一个信号 y，其大部分功率都高于 Fs/4 或奈奎斯特频率的一半。 采样率为 8192 Hz。

设计一个 34 阶 FIR 高通滤波器，以衰减低于 Fs/4 的信号分量。使用 0.48 的截止频率和 30 dB 纹波的切比雪夫窗。

对信号进行滤波。显示原始和高通滤波后的信号。两图都使用相同的 y 轴比例。

设计具有相同规格的低通滤波器。对信号进行滤波，并将结果与原始信号进行比较。两图都使用相同的 y 轴比例。

设计一个46阶FIR滤波器，该滤波器可以衰减低于 0.4 π \pi π rad / sample 和介于 0.6 π \pi π 和 0.9 π \pi πrad / sample之间的归一化频率。将滤波器称为 bM。

重新设计 bM，交换其通带和阻带。将新的滤波器称为 bW。 使用 fvtool 显示滤波器的频率响应。

使用 Hann 窗重新设计 bM。 （DC-0 是可选的。）比较 Hamming 窗和 Hann 窗设计的幅度响应。

使用 Tukey 窗重新设计 bw。 比较 Hamming 窗和 Tukey 窗设计的幅度响应。

滤波器阶数，指定为整数标量。

对于高通和带阻滤波器，fir1 始终使用偶数阶滤波器。此时阶数必须是偶数，因为奇数阶对称 FIR 滤波器在奈奎斯特频率下必须具有零增益。如果为高通或带阻滤波器指定奇数 n，则 fir1 将 n 增加 1。

数据类型：double

频率约束，指定为标量，二元向量或多元向量。

Wn 的所有元素必须严格大于 0 且严格小于 1，其中 1 对应于奈奎斯特频率：0 < Wn