带阻滤波器matlab代码

带阻滤波器是一种电子电路或信号处理技术，用于抑制特定频率范围内的信号而允许其他频段通过。在MATLAB中，可以使用`fir1`函数设计一个有限 impulse response (FIR) 带阻滤波器。以下是一个简单的例子，展示如何创建一个低通部分和高通部分串联而成的带阻滤波器：

% 设定所需频率特性
fc = [50 100]; % 阻止频率范围
fs = 1000;      % 总采样率
order = 64;     % 滤波器阶数

% 定义带阻区域的中心频率和衰减点
[numerator, denominator] = fir1(order, [fc(1)/fs - 1, fc(2)/fs + 1], 'bandstop');

% 创建带阻滤波器系统函数
H = tf(numerator, denominator);

% 打印滤波器参数
disp('Filter Parameters:');
disp(['Order: ', num2str(order)]);
disp(['Passband edge frequencies: ', num2str(fc)]);
disp(['Stopband attenuation: High-pass filter at DC and low-pass filter at ', num2str(fc(2))]);

% 测试滤波器效果（假设输入信号为单位冲激）
x = ones(1000, 1);
y = filter(H, x);

在这个例子中，`fir1`函数的第一个参数指定滤波器的阶数，第二个参数是一个向量，表示了所需的截止频率及其两侧衰减90分贝的边界。`'bandstop'`参数指示我们需要的是一个带阻滤波器。