带限滤波器matlab代码

### 如何用Matlab编写带限滤波器代码

为了设计和实现一个带限滤波器，在 MATLAB 中可以采用多种方式。下面提供了一种基于 FIR 滤波器的设计方案，该方案利用 `fir1` 函数来创建低通、高通或带通滤波器。

#### 使用 fir1 设计带限滤波器

% 参数设置
Fs = 8000; % 采样频率 (Hz)
Fc_low = 500; % 带通下截止频率 (Hz)
Fc_high = 2000; % 带通上截止频率 (Hz)

% 计算归一化截止频率
Wn_low = Fc_low / (Fs/2);
Wn_high = Fc_high / (Fs/2);

% 创建窗函数
windowSize = 127;
win = hamming(windowSize); 

% 设计算法选择 'bandpass' 表明这是一个带通滤波器
b = fir1(windowSize, [Wn_low Wn_high], 'bandpass', win);

% 绘制频响图
fvtool(b, 1);

这段代码展示了如何定义所需的参数，并调用 `fir1` 来生成具有指定特性的 FIR 滤波器系数[^1]。接着使用 `fvtool` 可视化工具查看所设计的滤波器特性。

对于更复杂的场景，比如当需要处理实时数据流时，则可以通过如下所示的方式应用此滤波器：

% 初始化状态向量
z = zeros(length(b)-1, 1);

while true
    % 获取新样本 x_new （此处应替换为实际的数据获取逻辑）
    
    % 应用 filtfilt 或 filter 对新的输入执行零相位失真过滤
    y_filtered = filter(b, 1, x_new, z);
    
    % 更新状态向量用于下一个循环迭代
    z = y_filtered(end-length(z)+1:end);
end

上述代码片段说明了怎样在一个持续运行的过程中不断更新输入信号并对其进行滤波处理。