基于MATLAB GUI的数字信号处理音频FIR去噪滤波器设计

biyezuopinvip

已于 2024-05-13 10:45:33 修改

阅读量3.3k

点赞数 3

于 2022-03-21 19:01:24 首次发布

本文链接：https://blog.csdn.net/newlw/article/details/123643134

版权

MATLAB GUI 数字信号处理音频 FIR去噪滤波器源码报告

资源下载地址：https://download.csdn.net/download/sheziqiong/85663076
资源下载地址：https://download.csdn.net/download/sheziqiong/85663076

在数字信号处理领域，音频去噪是一个非常重要的应用。通过使用MATLAB GUI，我们可以方便地设计和实现FIR（有限冲激响应）滤波
器，以对音频信号进行去噪处理。以下是关于基于MATLAB GUI的数字信号处理音频FIR去噪滤波器设计的详细介绍：

1. 内容简介

本设计利用MATLAB GUI平台，采用窗函数法设计FIR数字滤波器。通过对含有噪声的音频信号进行数字滤波处理，我们可以得到降噪后的声音信号，并对其进行时域和频域分析。此外，还可以比较不同窗函数对滤波效果的影响。用户只需将相关文件解压至一个目录下，并运行相应的m文件，即可使用该系统。

2. 函数使用

在设计过程中，我们使用了以下MATLAB函数：

audioread()：用于读取.wav格式的音频文件。
sound()：播放音频信号。
clear sound：停止播放音频。
audiowrite()：将处理后的音频信号写入.wav文件。
加入噪声的方法：通过特定的数学表达式向音频信号中添加白噪声或特定频率的噪声。
频谱分析：使用fft()和fftshift()函数对信号进行频域分析。
fir1()：使用窗函数法设计FIR滤波器。

3. 实现功能

该系统实现了以下功能：

打开音频文件：选择并打开.wav格式的音频文件，自动生成音频的原始波形和频谱图。
加入噪声：可以选择加入白噪声或特定频率的噪声，并自动绘制加入噪声后的波形和频谱图。
滤波处理：输入滤波器的通/阻带频率，选择窗函数和滤波器类型，生成滤波后的波形和频谱图。
音频播放/停止：允许用户播放或停止播放原始、加噪或滤波后的音频。
图片导出：将波形和频谱图导出为jpg、png、bmp或eps格式。
保存文件：将加噪或滤波后的音频保存为.wav文件。

3、实现功能

实现的功能有：

打开文件：选择路径打开wav格式的音频文件，自动生成音频的原始波形与频谱。
加入噪声：有两种噪声可以选择加入，一种是白噪声，其频率蔓延整个频谱；一种是特定频率的噪声，可通过输入频率加入单一频率的噪声。加入噪声后自动绘制加入噪声后的波形与频谱。
滤波处理：首先输入滤波器通/阻带的开始频率与截止频率（若为低/高通类型滤波，则只需输入开始频率；若为带通/阻类型，则开始与截止都要输入；输入频率值为真实频率值，可根据频谱图进行判断），之后选取窗函数和滤波类型，将会生成滤波处理后的波形与频谱。www.biyezuopin.cc
音频播放/停止：可随时播放/停止原始、加噪、滤波处理后的音频。
图片导出：将波形、频谱图片一张张导出保存，可选的格式有jpg、png、bmp、eps。
保存文件：将加躁/滤波后的音频导出保存。

4. 操作实例

以“卢本伟语音包”为例，首先将其转换为wav格式并导入系统。然后，添加白噪声并观察频谱变化。接着，使用低通滤波器进行滤波处理，以去除噪声。如果噪声是特定频率的，可以设计带阻滤波器来去除该频率范围内的噪声。

添加白噪声后，此时播放音频能听到显著杂音。而从原始信号的频谱来看，初始音频的频率主要集中在0-1000Hz，因此我们可以选用低通滤波器，阻带开始频率设为1000Hz，选用矩形窗进行滤波，得到结果如下：
由于白噪声遍布于整个频谱，对于噪声频谱于音频频谱的重叠部分，我们无法通过FIR滤波器进行滤除，依然会有小部分杂音存在。若噪声为特定单一频率的噪声，我们可以较好地将其去除。对于该音频添加5000Hz的特定频率，通过设计带阻滤波器，阻带范围为4500-5500Hz对其进行滤波，如下图所示：

5. 窗函数对比

通过向音频中添加特定频率的噪声，并使用不同的窗函数设计带阻滤波器，我们可以比较不同窗函数对滤波效果的影响。实验结果表明，凯塞窗的滤波效果最好，其次是矩形窗，然后是汉宁窗，而海明窗、三角窗和布莱克曼窗的滤波效果相对较差。

通过上述设计，我们可以有效地对音频信号进行去噪处理，提高音频质量。此外，通过比较不同窗函数的滤波效果，我们可以为特定应用场景选择最合适的窗函数。这种基于MATLAB GUI的音频FIR去噪滤波器设计方法，不仅易于实现，而且具有很强的实用性和灵活性。
选用带阻滤波器，设置阻带范围4000-6000Hz，观察各窗函数对其滤波的效果。(每行从左到右分别是：矩形窗，三角窗，海明窗，汉宁窗，布莱克曼窗，凯塞窗，下同)