数字信号处理实验报告.docVIP

长春理工大学

电工电子实验教学中心

学生实验报告

2014——2015学年第一学期

实验课程数字信号处理实验

实验地点东1教学楼414实验室

学院电子信息工程学院

专业通信工程

学号120421101

姓名杨杰

实验工程

音频信号采集、分析及处理

实验时间

第十四周星期四9、10节

实验台号

操作成绩

报告成绩

一、实验目的

以音频信号为例，熟悉模拟信号数字处理过程，进一步理解数字信号处理概念。

掌握运用Matlab实现对音频信号的时频分析方法；

初步掌握数字音频信号合成的方法。

掌握运用Matlab设计IIR和FIR滤波系统的方法；

掌握运用Matlab实现对加噪音频信号进行去噪滤波的方法。锻炼学生运用所学知识独立分析问题解决问题的能力，培养学生创新能力。

二、实验仪器

计算机、Matlab软件

三、实验性质

综合分析、设计性实验

四、实验任务

实验内容一：windows系统中的“ding”音频信号的采集、分析、合成

音频信号的采集

编写Matlab程序，采集windows系统中的“ding”声，得到*.wav音频文件，而后实现音频信号回放。

音频信号的频谱分析

运用Matlab软件实现对音频信号的时域分析和频域分析，并打印相应的图形，完成在实验报告中。

注意：此音频信号的频谱包含两条主要谱线，在进行频谱分析时，注意频谱的完整性，利用MATLAB实现对两条主要谱线的定位并计算谱线所对应的模拟频率。

音频信号的分解和合成

运用Matlab软件实现音频信号的分解与合成，将音频信号的频谱中两局部频谱成分进行分解，分别绘制出分解后的两个信号的频谱图；然后将分解后的两个信号再合成为一个新的信号，将合成后的新信号的时域、频域图与原来的信号时域、频域图相比拟，绘制出比照效果图。

音频信号的回放

运用Matlab软件实现音频信号的回放，将合成后的新信号和原音频信号分别进行回放，比照两个信号的声音效果。

音频信号分段傅里叶分析〔选作〕

分析对一般音频.wav信号进行一次性傅里叶分析时存在的主要问题，利用分段傅里叶变换对该音频信号重新分析并合成。比照一次傅里叶分析结果并进行总结。

实验内容二：任意音频信号的时域和频域分析及数字滤波器设计

音频信号的采集

音频信号的采集可以通过Windows自带的录音机也可以用专用的录制软件录制一段音频信号〔尽量保证无噪音、干扰小〕，也可以直接复制一段音频信号〔时间为1s〕，但必须保证音频信号保存为.wav的文件。

音频信号的时域、频域分析

运用Matlab软件实现对音频信号的翻开操作以及时域分析和频域分析，并画出相应的图形，打印在实验报告中。

引入干扰信号

在原有的音频信号上，叠加一个频率为100KHz的正弦波干扰信号〔幅度自定，可根据音频信号幅度情况而定〕。

数字滤波器设计

运用Matlab实现IIR或FIR滤波系统，要求参加干扰的音频信号经过滤波系统后，能够滤除100KHz的干扰信号，同时保存原有的音频信号，要求绘制出滤波系统的频谱图，同时绘制出滤波前后音频信号的时域、频域比照图。

实验内容三：编写FFT算法程序〔此局部为选作内容〕

分析FFT算法实现过程，运用Matlab软件实现FFT函数功能，要求运用自行编写的Matlab程序实现对固定信号〔自定义信号〕的频谱分析，并绘制出频谱图；然后将自己设计的FFT算法程序与FFT函数实现的频谱图相比拟，分析两个程序结果的差异及原因。

五、实验任务分解及思路

实验内容一：windows系统中的“ding”音频信号的采集、分析、合成

通过对任务内容要求的分析，建立程序框图如图一〔任务一程序框图〕所示，详细程序及结果见附录。编程思想及思路如下：

首先用audioread读取音频信号得到signal序列及采样频率fs，由于音频信号是双声道，这里只取第一个声道进行处理，由于原信号本身就较短，这里对长度就不作处理。接着按实验要求进行音频回放，并绘制出音频信号的时域波形。

调用fft函数，对signal进行快速傅里叶变换，并绘制出频谱图。在绘制频谱图时考虑到快速傅里叶变换的对称性，这里只取序列的前半局部进行观察分析。

为了实现音频信号的分解及合成，先对原信号的频谱图进行观察分析，发现原信号的主要能量集中在两个频率上，于是考虑用这两个频率的正弦信号合成该信号。这样首先要找到这两个频率。为了求得这两个频率，我先调用findpeaks函数找到频谱图上的各个峰值点SIGNAL_pks及其对应的位置locs，然后用sort对峰值点进行排序，找到最大的两个值，接着用find函数找到这两个峰