简介：这份文档面向学习数字信号处理、需要在MATLAB中实现滤波器的学生与工程人员，围绕长度为N=45、阻带衰减AS=60dB的FIR带阻滤波器展开设计讲解。内容以凯塞-贝塞尔窗函数法为主线，说明参数beta如何影响主瓣宽度、旁瓣大小与过渡带宽度，并给出Beta=0.1102*(As-8.7)的计算依据，同时提供freqz.m、ideal_lp、bsfilter.m等源程序，演示频率响应、相位响应与群延时的求解过程，以及理想冲激响应、窗函数和实际冲激响应的绘图对比。资源包为单个doc文件，约60KB，结构紧凑，便于直接查阅与复现。目前已有1376人学习下载，适合希望掌握窗函数选型、滤波器性能分析与MATLAB实现思路的读者参考。
1. 从一段能跑的 FIR 带阻代码说起
很多人在 MATLAB 里做数字滤波器设计，第一步就卡在“窗函数参数到底怎么定”。这份matlab设计FIR带阻滤波器.doc给了一个很具体的答案：N=45、阻带衰减 As=60dB、用凯塞—贝塞尔窗（Kaiser 窗）来逼近。它没有停在公式推导，而是直接给了bsfilter.m、ideal_lp.m、freqz.m三个可运行片段，把理想冲激响应、窗函数、实际冲激响应和幅度响应四张图一次性画出来。
带阻滤波器要干的事很明确：把某一段频率（这里是 π/3 到 2π/3）压下去，其余频段尽量原样保留。FIR 的好处是线性相位、天然稳定，代价是要用足够长的阶数去换过渡带陡度。这份文档的价值在于，它把“As 决定 beta、beta 决定窗形、窗形决定过渡带”这条链路用代码串了起来，适合正在做课程设计、信号处理大作业，或者第一次用 MATLAB 手写 FIR 的人照着复现。
2. Kaiser 窗参数与带阻滤波器的理论映射
2.1 为什么带阻要用“低通相减”来构造
FIR 带阻没有现成的闭式公式，常见做法是用频域拼接：一个截止在 wc1 的低通，加上一个截止在 π 的全通（其实就是 δ 函数），再减去一个截止在 wc2 的低通。文档里这一行是关键：
bd = ideal_lp(wc1,N) + ideal_lp(pi,N) - ideal_lp(wc2,N);逻辑上，ideal_lp(pi,N)在 n=0 处为 1、其余为 0，相当于全通；ideal_lp(wc1,N)保留 0~wc1；减去ideal_lp(wc2,N)就把 wc1~wc2 这段挖掉了。这样得到的bd是理想带阻的冲激响应，但它无限长且非因果，必须加窗截断。
2.2 Kaiser 窗的 beta 与 As 的定量关系
Kaiser 窗的可调性全在 beta 上。文档给出的经验式是：
Beta = 0.1102*(As-8.7);这个式子只在 As > 50dB 时成立。As 在 21~50dB 区间要用0.5842*(As-21)^0.4 + 0.07886*(As-21)，As < 21dB 时 beta 取 0。很多人直接套 0.1102 那一行，结果 As=40 时过渡带偏宽，就是没分段。
| 参数 | 含义 | 本例取值 | 影响 |
|---|---|---|---|
| N | 滤波器长度（阶数+1） | 45 | 越大过渡带越窄，计算量越大 |
| As | 阻带最小衰减 | 60 dB | 决定 beta，进而决定旁瓣 |
| beta | Kaiser 窗形状参数 | 0.1102*(60-8.7)≈5.65 | 越大旁瓣越低、主瓣越宽 |
| wc1 | 下通带截止 | π/3 | 阻带左边界 |
| wc2 | 上通带截止 | 2π/3 | 阻带右边界 |
提示：N=45 是奇数，
alpha=(N-1)/2=22为整数，ideal_lp里m=n-alpha+eps加 eps 是为了避免 n=alpha 时除零。如果换成偶数 N，群延迟会落在半整数点，线性相位仍然成立，但绘图时 stem 的对称中心会偏移半个点。
2.3 过渡带宽度与 N 的取舍
Kaiser 窗的过渡带近似满足Δw ≈ (As-8)/(2.285*N)。把 As=60、N=45 代进去，Δw≈0.507 rad，换算成归一化频率约 0.08。也就是说 wc1 和 wc2 之间至少要留出这么宽，否则阻带边缘会翘起来。文档选 π/3 和 2π/3，间隔 π/3≈1.047 rad，远大于 0.507，所以过渡带是够用的。如果想把阻带压得更窄，要么加大 N，要么降低 As 要求，二者不能同时满足。
3. 三个子程序的 MATLAB 实现与逐行拆解
3.1 ideal_lp：理想低通的 sinc 生成
function hd = ideal_lp(wc, N) % 生成截止频率为 wc、长度为 N 的理想低通冲激响应 alpha = (N-1)/2; % 滤波器中心，保证线性相位 n = 0:N-1; m = n - alpha + eps; % 加 eps 防止 m=0 时除零 hd = sin(wc*m) ./ (pi*m); % sinc 形式，pi*m 是归一化 end这里wc用弧度制，sin(wc*m)/(pi*m)就是理想低通的标准形式。eps是 MATLAB 的最小浮点数，加在分母上让 n=alpha 那一点不报 NaN，实际值趋近 wc/pi。调用时ideal_lp(pi,N)会得到 n=alpha 处为 1、其余接近 0 的序列，正好当全通用。
3.2 freqz.m：频率响应与群延迟的封装
function [db,mag,pha,grd,w] = freqz_m(b,a) % db 相对振幅(dB)，mag 绝对振幅，pha 相位，grd 群延迟，w 频率采样点 [H,w] = freqz(b,a,1000,'whole'); % 取整个单位圆 1000 点 H = H(1:1:501)'; % 只取 0~pi 半圈 w = w(1:1:501)'; mag = abs(H); db = 20*log10((mag+eps)/max(mag)); % 归一化到 0dB 峰值 pha = angle(H); grd = grpdelay(b,a,w); % 群延迟，衡量相位线性度 endfreqz的'whole'参数让采样覆盖 0~2π，取前 501 点就是 0~π。db里除以max(mag)是为了把通带峰值归一到 0dB，看阻带衰减时直接读负值即可。grpdelay对 FIR 来说理论上应恒等于 (N-1)/2=22，如果画出来波动大，说明系数算错了。
3.3 bsfilter.m：主程序串起全流程
N = 45; As = 60; n = 0:N-1; beta = 0.1102*(As-8.7); w_kai = kaiser(N, beta)'; % 注意转置成列向量 wc1 = pi/3; wc2 = 2*pi/3; bd = ideal_lp(wc1,N) + ideal_lp(pi,N) - ideal_lp(wc2,N); h = bd .* w_kai; % 加窗截断 [db,mag,pha,grd,w] = freqz_m(h,[1]);kaiser(N,beta)返回列向量，ideal_lp返回行向量，所以w_kai要转置，否则.*会因维度不匹配报错——这是新手最常见的坑。freqz_m(h,[1])里 a=1 表示 FIR，分母只有常数项。
画图部分用subplot(2,2,1)到(2,2,4)排四张图，stem画离散冲激响应，plot画连续幅度曲线。axis([0,1,-80,10])把横轴限制在归一化频率 0~1，纵轴 -80~10dB，方便看 60dB 的阻带深度是否达标。
4. 复现时的排错与性能验证
4.1 常见报错与对应修法
| 现象 | 原因 | 修法 |
|---|---|---|
Matrix dimensions must agree | w_kai与bd行列不一致 | 给kaiser结果加'转置 |
| 阻带衰减只有 40dB | beta 用了 0.1102 公式但 As<50 | 换分段公式算 beta |
| 群延迟曲线波动 | 系数非对称或 N 取偶数未调整 | 检查ideal_lp的 alpha |
freqz_m未定义 | 文件名与函数名不一致 | 存为freqz_m.m，调用同名 |
4.2 用数据验证 60dB 是否真的达到
画图只能看个大概，要确认指标得直接读db数组：
% 找阻带区间对应的索引 idx = find(w/pi >= 1/3 & w/pi <= 2/3); max_stop = max(db(idx)); % 阻带内最大 dB 值 fprintf('阻带最大衰减: %.2f dB\n', max_stop);如果max_stop在 -60 以下，说明设计达标；如果在 -55 左右，说明 N 偏小或 beta 偏小，把 N 加到 51 或 61 再试。注意db已经归一化，通带峰值是 0dB，所以阻带值直接就是衰减量。
4.3 过渡带与阶数的边界
前面算过 Δw≈0.507，实际过渡带会略宽于理论值。如果发现 wc1 附近通带就开始掉，说明过渡带侵入了通带，需要把 wc1 往左挪或加大 N。FIR 带阻没有免费午餐：As、过渡带宽度、N 三者固定两个，第三个必然被约束。课程设计里常见的要求是 As=60dB、过渡带 0.1π，那 N 至少要 100 以上，45 是压不住的。
5. 从脚本到可复用函数：参数化与批量验证
把bsfilter.m改造成函数，才能快速扫参数：
function [h,db,w] = design_bs(N, As, wc1, wc2) % 参数化 FIR 带阻设计，返回系数与幅度响应 if As > 50 beta = 0.1102*(As-8.7); elseif As >= 21 beta = 0.5842*(As-21)^0.4 + 0.07886*(As-21); else beta = 0; end n = 0:N-1; w_kai = kaiser(N, beta)'; bd = ideal_lp(wc1,N) + ideal_lp(pi,N) - ideal_lp(wc2,N); h = bd .* w_kai; [db,~,~,~,w] = freqz_m(h,[1]); end调用时一行就能换参数：[h,db,w]=design_bs(61,60,pi/3,2*pi/3);。想批量看 N 对阻带的影响，套个循环把max(db(idx))打出来，比一张张画图快得多。实际工程里还会把h导出成系数文件给 FPGA 或 DSP 用，注意 MATLAB 的freqz是归一化频率，移植时要按采样率换算成实际 Hz。这套脚本跑通一次，后面改指标就是改几个数字的事。
本文还有配套的精品资源，点击获取
