数字信号处理实验-音频信号处理.pdfVIP

实验课程数字信号处理

开课学期2022-2023学年第一学期

实验题目实验五：音频信号分析及处理

学生专业王大珩未来技术学院创新实验班

学生姓名班级学号

实验时间11.23实验台号线上实验

报告成绩

报告评语

一、实验目的

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

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

3.掌握数字音频信号合成的方法；

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

5.掌握运用Matlab实现对加噪音频信号进行去噪滤波的方法。锻炼学生运用所学知

识独立分析问题解决问题的能力，培养学生创新能力。

二、实验性质

设计性实验

三、实验任务分析

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

通过对任务具体内容的分析，可以建立出任务一程序框图如下，之后将对编程思

想及思路进行介绍：

图1任务一程序流程图

1、音频信号采集

在获取windows系统中的“ding”声并另存为”ding.wav”后,先用audioread函

数读取音频信号得到采样序列signal及对应采样频率fs,由于获取的音频信号是双声

道，只取第一个声道进行处理,接着使用sound函数以fs频率进行音频回放。

2、音频信号频域分析

以采样间隔T划分时域并绘制出signal信号的时域波形；调用fft函数,对signal

进行快速傅里叶变换，用abs函数取傅里叶变换后结果的幅值进行幅频分析，绘制出

频谱图。在绘制频谱图时由于考虑到快速傅里叶变换的对称性，只取序列的前半部分

进行观察分析。

3、音频信号分解

为了实现音频信号的分解及合成，先对原信号的频谱图进行观察分析，发现原信

号的主要能量集中在三个主要频率上，于是考虑用这三频率的正弦信号合成原信号。

为了求得这三个频率，先调用findpeaks函数找到频谱图上的各个局部极大值peak及

其对应的位置locs,然后用sort对峰值点进行排序，找到最大的三个值，接着用find

函数找到这三个最大值在locs中的位置，也就知道了对应的频率。这里有一个问题就

是最小的峰值频率并不是在sort排序后的第三位而是在第四位，需要有一个调整；确

定了主要谱线后，使用text函数进行峰值标注；

4、音频信号合成

接着将这三个谱线还原回时域正弦信号，幅度的比例等于对应频率上的幅度比例

然后然后叠加，得到合成后的信号，绘制出时域波形，与原信号波形进行比较,接着对

两个正弦信号进行fft,绘制出他们的频谱，然后对合成的信号进行fft,做出频谱图和

原信号的频谱图进行比较.

5、音频信号回放

用sound函数进行原信号和合成信号的回放，比较差异。

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

通过对任务具体内容的分析，可以建立出任务二程序框图如下，之后将对编程思想

及思路进行介绍：

图2任务二程序流程图

1、音频信号采样

自己录音频并另存为”ding.wav”后,先用audioread函数读取音频信号得到采样

序列signal及对应采样频率fs,由于获取的音频信号是双声道，只取第一个声道进行

处理。

2、时域采样

使用