巴特沃斯数字滤波器设计与仿真系统

资 源 简 介

本项目专注于在MATLAB计算环境中实现基于巴特沃斯（Butterworth）逼近理论的数字滤波器设计全流程。核心功能在于利用巴特沃斯滤波器在通带内具有最大平坦幅度响应（无纹波）的特性，为信号处理任务提供稳定的滤波解决方案。系统首先根据用户输入的设计指标（包括采样频率、通带截止频率、阻带截止频率、通带最大衰减及阻带最小衰减），通过`buttord`函数自动计算满足条件所需的滤波器最小阶数和3dB截止频率。随后，利用`butter`函数生成滤波器的传递函数系数（分子核分母多项式），支持低通、高通、带通及带阻四种基本滤波类型的构建。项目集成了详细的分析模块，利用`freqz`和`zplane`函数绘制滤波器的幅频特性曲线、相频特性曲线及零极点分布图，以验证设计的准确性和稳定性。此外，系统包含模拟仿真测试功能，能够生成含噪信号（如混叠了高频噪声的正弦波），通过`filter`函数将设计好的滤波器应用于该信号，并对比输出滤波前后的时域波形与频谱图，直观展示巴特沃斯滤波器在滤除特定频段噪声同时保持通带信号平滑性的实际效果。

详 情 说 明

基于MATLAB的巴特沃斯数字滤波器设计与仿真系统

项目简介

本项目是一个基于MATLAB环境开发的数字信号处理仿真系统，旨在演示巴特沃斯（Butterworth）数字滤波器的完整设计与应用流程。系统通过编写的MATLAB脚本，实现了从信号生成、滤波器指标设定、阶数自动计算、系数生成、频率特性分析到时域滤波与频域分析的全过程。

项目的核心在于利用巴特沃斯滤波器在通带内“最大平坦”（无纹波）的幅度响应特性，设计一个低通滤波器，从混叠了高频干扰和随机噪声的混合信号中精确提取出有用的低频信号。

功能特性

• 混合信号仿真生成：能够生成包含基波低频信号、多重高频正弦干扰以及高斯白噪声的复杂仿真信号。
• 自动化滤波器设计：通过设定通带/阻带截止频率及衰减指标，算法自动计算满足条件所需的最小滤波器阶数和3dB截止频率。
• 全面的特性分析：自动绘制滤波器的幅频响应、相频响应曲线以及零极点分布图，用于评估滤波器的稳定性与频率特性。
• 时域/频域对比验证：对比滤波前后的时域波形（直观展示去噪效果）以及FFT频谱（展示特定频段噪声的滤除情况）。

系统要求

• MATLAB R2016a 或更高版本
• Signal Processing Toolbox（信号处理工具箱）

使用方法

1. 确保MATLAB环境已安装信号处理工具箱。
2. 直接运行主程序脚本。
3. 程序将自动清理工作区，执行计算，在控制台输出滤波器的设计参数（阶数N、截止频率Wn）。
4. 程序运行结束后，将弹出四个图形窗口，分别展示：

* 滤波器频率响应特性（幅频与相频）。 * 零极点分布状况。 * 时域滤波效果对比（局部放大视图）。 * 频域滤波效果对比（单边频谱）。

核心算法与代码实现逻辑

本项目的主程序（main.m）严格按照数字信号处理的流水线进行设计，具体实现逻辑如下：

1. 环境初始化与信号源构建

程序首先清理工作环境。随后定义了采样频率为 2000 Hz。

• 有用信号：构建了一个频率为 50 Hz 的纯净正弦波。
• 噪声模拟：合成了一个复合噪声信号，包含 350 Hz 和 600 Hz 的高频正弦干扰，并叠加了随机高斯白噪声。
• 输入信号：将有用信号与噪声叠加，形成待处理的时间序列数据。

2. 滤波器设计指标设定

代码中硬编码了低通滤波器的设计指标，旨在保留50Hz信号并滤除350Hz以上的干扰：

• 通带截止频率 (Wp)：120 Hz
• 阻带截止频率 (Ws)：300 Hz
• 通带最大衰减 (Rp)：1 dB
• 阻带最小衰减 (Rs)：40 dB

程序会将频率参数归一化为相对于奈奎斯特频率（Fs/2）的值，以适配MATLAB函数要求。

3. 巴特沃斯滤波器核心计算

这是代码的核心算法部分：

• 利用 buttord 函数，根据归一化的频率指标和衰减指标，自动计算出满足设计要求所需的最小滤波器阶数 N 和 3dB 截止频率 Wn。此步骤确保了滤波器既能满足性能指标，又不会因阶数过高而增加计算负担。
• 利用 butter 函数，基于计算出的 N 和 Wn，生成低通滤波器的传递函数系数（分子向量 b 和分母向量 a）。

4. 频率响应与稳定性分析

设计完成后，程序对滤波器本身的特性进行分析：

• 使用 freqz 函数计算滤波器的复频响应，并绘制幅频响应（dB）和相频响应曲线。图中会使用参考线标出通带和阻带的边界，便于验证设计是否达标。
• 使用 zplane 函数绘制零极点图，通过观察极点是否位于单位圆内来验证系统的稳定性。

5. 信号滤波仿真

• 调用 filter 函数，利用生成的系数 (b, a) 对含噪输入信号进行时域差分方程滤波。该函数模拟了实时处理过程（注：会引入自然的相位延迟）。

6. 频谱分析 (FFT)

为了从频域验证效果，代码执行了快速傅里叶变换（FFT）：

• 计算原始含噪信号和滤波后输出信号的单边频谱。
• 对频谱幅度进行归一化处理，以便于直观对比。

7. 结果可视化

程序最后生成两组对比图：

• 时域对比：截取前0.1秒（200个点）的数据，在一个窗口中分三栏展示理想信号、含噪信号和滤波输出信号。通过对比可以看出毛刺和高频振荡被平滑滤除。
• 频域对比：展示滤波前后的频谱。可以清晰观察到350Hz和600Hz的干扰谱峰在滤波后被显著抑制，而50Hz的主信号得到保留。

关键函数说明

• buttord(Wp, Ws, Rp, Rs)

* 功能：巴特沃斯滤波器阶数选择函数。 * 作用：根据通带/阻带频率和波纹/衰减指标，返回最低阶数 N 和对应的截止频率 Wn。

• butter(N, Wn, 'low')

* 功能：巴特沃斯滤波器设计函数。 * 作用：根据阶数和截止频率，生成数字滤波器的传递函数系数（b, a）。

• freqz(b, a, ...)

* 功能：计算数字滤波器的频率响应。 * 作用：用于分析并绘制滤波器的幅度和相位特性。

• filter(b, a, input_signal)

* 功能：一维数字滤波器。 * 作用：直接利用差分方程对数据进行滤波处理。

• zplane(b, a)

* 功能：零极点图绘制。 * 作用：直观展示系统零点和极点在复平面上的分布，判断系统稳定性（极点需在单位圆内）。