数字滤波器，基于matlab的数字滤波器程序，包扩低通、高通、带通、带阻滤波器，IIR，切比雪夫I型。 可提供FFT分析程序，进行波形频谱分析。 图1为低通滤波器及相关波形频谱分析，图2为高通滤波器及相关波形频谱分析，图3为带通滤波器及相关波形频谱分析。

数字滤波器的魔法在信号处理领域持续发威，咱们今天用Matlab实操一把。先看低通滤波器怎么拦截高频噪音——假设采样率8kHz，想保留500Hz以下信号。核心代码就这几行：

Fs = 8000;
Fc = 500;
N = 6; % 切比雪夫I型阶数
Rp = 1; % 通带波纹(dB)
[b,a] = cheby1(N, Rp, Fc/(Fs/2), 'low');
freqz(b,a) % 看频响曲线

这里的cheby1函数参数暗藏玄机：归一化截止频率必须用Fc/(Fs/2)处理，相当于把数字频率转换到0-1区间。滤波器系数b和a直接决定系统函数，用fvtool可视化更直观。

当处理实际信号时，先做个混频实验：

t = 0:1/Fs:1;
signal = sin(2*pi*300*t) + 0.5*sin(2*pi*2000*t);
filtered = filter(b,a,signal);

原始信号里300Hz和2000Hz成分叠加，滤波后2000Hz应该被压制。FFT验证效果：

L = length(signal);
Y = fft(signal);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1) % 原始频谱

高通滤波器只需把'low'改为'high'，但要注意相位畸变问题。带通更有趣，比如想提取1000-2000Hz频段：

Wn = [1000 2000]/(Fs/2);
[b,a] = cheby1(5, 0.5, Wn, 'bandpass');

这里Wn参数变成二元数组，相当于设置通带边界。当信号包含800Hz、1500Hz、3000Hz三个成分时，滤波后1500Hz成分会被完美捕获。

数字滤波器，基于matlab的数字滤波器程序，包扩低通、高通、带通、带阻滤波器，IIR，切比雪夫I型。 可提供FFT分析程序，进行波形频谱分析。 图1为低通滤波器及相关波形频谱分析，图2为高通滤波器及相关波形频谱分析，图3为带通滤波器及相关波形频谱分析。

切比雪夫I型的波纹特性在阻带边缘尤为明显，对比巴特沃斯滤波器，它的过渡带更陡峭但代价是通带波动。实际调试中发现，当阶数超过8阶时，数值稳定性开始变差，这时候需要改用二阶节实现。

频谱分析时有个坑要注意：做FFT前必须加窗！矩形窗泄露严重，建议用汉宁窗：

window = hann(length(signal));
Y_windowed = fft(signal .* window');

最后给个实用技巧：用fdatool交互设计滤波器后，直接导出MATLAB代码，比自己写参数方便得多。滤波后的时域波形可能看起来有点"糊"，这时候别急着调参，先看频谱确认实际效果——眼见不一定为实，频谱才不说谎。