【MATLAB第17期】基于MATLAB的音频数字处理系统设计，含GUI和报告（数字信号课程作业）
一． 研究背景
在数字信号处理上已经学了很多数字信号处理的模型和算法，例如基本的 DFT，FFT，IIR 和 FIR
等方法的学习，将这些基本的处理方法和滤波器实现应用在音频的处理上更有助于我们对这些方法的理解。因此，此次制作一个音频数字处理的系统，可以实现基本的时域频域分析和噪声的一些添加处理，以及
IIR 和 FIR 滤波器的实现。
二． 实验环境
- 操作系统：Window 10； 设计软件：MATLAB R2020a；
三． 设计原理
1.理论知识
1.1.DFT 的基本原理
散傅里叶变换(Discrete Fourier Transform，DFT)实质是有限长序列傅里叶变换的有限点离散采样，从而实现了频域离散化，使数字信号处理可以在频域采用数值运算的方法进行，这样就大大增加了数字信号处理的灵活性。DFT 有多种快速算法，统称为快速傅里叶变换(Fast Fourier Transform，FFT)
- )
FFT(Fast Fourier Transform)是离散傅立叶变换的快速算法，可以将一个信号变换到频域。有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易看出特征了。这就是很多信号分析采用 FFT 变换的原因。另外，FFT 可以将一个信号的频谱提取出来，这在频谱分析方面也是经常用的。
FFT 能够更高效的计算离散傅里叶变换，因此常用于计算机编程的计算，用来节省复杂信号的数字处理，应用较广泛的就是数学软件 MATLAB 对 FFT 算法的设计。
1.2.IIR 滤波器
一般将网络结构分成两类，一类称为有限长单位脉冲响应网络，简称 FIR
（Finite Impulse Response）网络，另一类称为无限长单位脉冲响应网络，简称 IIR
（Infinite Impulse Response）网络。 IIR 网络结构存在输出对输入的反馈支路，也就是说，信号流图中存在反馈环路。这类网络的单位脉冲响应是无限长的。其 N 阶差分方程：
模拟滤波器的理论和设计方法也已经相当成熟，有多种典型的模拟滤波器，例如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Ellipse)滤波器、贝塞尔(Bessel)滤波器等。
1.3.FIR 滤波器
FIR 网络中一般不存在输出对输入的反馈支路，因此差分方程描述为：
：
- FIR 滤波器的设计主要包括线性相位 FIR 数字滤波器设计，利用窗函数法设计，频率采样法设计，等波纹最佳逼近法设计。其中窗函数的设计方式较为常见，基本窗函数特性如图 2 所示。
2.算法原型
DFT（或 FFT）MATLAB 提供了用快速傅里叶变换算法 FFT 计算 DFT 的函数 fft，其调用格式为: Xk = fft (xn, N)；Ifft 函数计算 IDFT，其调用格式与 fft 函数相同，可参考 help 文件。
Matlab 中利用巴特沃斯滤波器和双线性变换法设计实现 IIR 滤波器的基本实现，具体包含的函数 buttord（），buttap（）和 bilinear（）函数。 FIR 滤波器的实现利用窗函数法设计，并选用哈明窗设计，具体函数包括hamming（），fir1（）函数等。
四． 系统方案
1.音频处理系统界面搭建
利用 MATLAB 的 GUI 图形化界面搭建工具和函数搭建基本的音频分析界面， 便于后续的音频处理的选择和参数设定，来方便观察不同参数的影响。工具栏主要包括文件——打开声音文件，录音，关闭三个选项，用于导入音频文件以便后续分析；叠加噪声——叠加白噪声，叠加 单频噪声和多频噪声三个选项；滤波器设计——IIR 滤波器设计和 FIR 滤波器设计。
2.音频的时域频域分析
导入音频文件或录音文件后，对音频信息进行时域分析和频域分析，具体核心代码如下：左边是时域分析，右边是频域分析。
3.音频信号加噪声分析
- 对音频基本的时域和频域观察后，增加高斯白噪声核心代码如下：
- 添加单频噪声和多频噪声核心代码如下：（左边单频噪声，右边多频噪声）
4.音频信号经 IIR 和 FIR 滤波器分析
- 利用巴特沃斯滤波器和双线性变换法设计实现 IIR 滤波器，代码如下：
利用窗函数法设计，并选用哈明窗设计 FIR 滤波器，代码如下：
五． 实验测试与结果分析
1.音频时域频域分析结果
导入一段长约 9 秒的音频到设计的音频数字处理系统中，获得的时域频域图如图 3 所示。
通过分析可以看出，该段音频的的低频和高频成分相对较多，中频成分占比较少。
2.音频加噪声分析结果
在源声音文件上添加高斯白噪声后的得到的结果如图 4 所示。以看到高斯白噪声的添加使得音频频谱在各个频段普遍叠加上了一段信号，该信号即使高斯白噪声的影响。
在源声音文件上添加单频噪声得到的结果如图 5 所示.
在源声音文件上添加多频噪声得到的结果如图 6所示.
3.音频经 IIR 和 FIR 滤波器分析结果
IIR 滤波器设置参数（通带频率 Wp=0.2kHz；阻带截至 Ws=0.25kHz；通带衰减 Rp=1dB；阻带衰减 Rs=15dB）后得到的滤波器特性如图 7所示，经过该滤波器得到的音频特性如图 8所示。可以看到，经过 IIR 滤波器后的滤去了很多幅值较低的中频信号，基本保留了低通信号
————————————————
版权声明：本文为CSDN博主「随风飘摇的土木狗」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
。
Figure 7：IIR 滤波器幅值相位特性曲线
Figure 8：经过 IIR 滤波器得到的音频特性
FIR 滤波器设置参数（通带频率 Wp=0.2kHz；阻带截至 Ws=0.25kHz；阻带衰减 Rs=50dB）后得到的滤波器特性如图 9 所示，经过该滤波器得到的音频特性如图 10 所示。
Figure 9：FIR 滤波器幅值相位特性曲线
Figure 10：经过 FIR 滤波器得到的音频特性参考文献
[1]曹亮.基于 DSP 的音频信号处理系统设计[J].电子技术与软件工程,2019(10):95.
[2]詹颖珩.基于 MATLAB 的椭圆滤波器对音频信号的分离[J].通讯世界,2018(06):273-274. [3]刘海波,汤群芳.基于 MATLAB 的 IIR 滤波器的设计研究[J].软件导刊,2008(11):187-188.
[4]周红鸥.基于 DSP 的音频信号采集处理系统设计[J].西南民族大学学报(自然科学版),2011,37(S 1):100-103.
[5]冯志鸿.DSP 设计的数字音频信号处理[J].电子技术与软件工程,2019(10):93.
代码获取
