数字信号处理课程设计

1 数字信号处理课程设计 课程总结报告题 目 运用 matlab 的数字音效处理系统设计学 院 电子信息工程学院专 业 通信工程姓 名 赵亚学 号 P21314067队 员 吴双胜 程勇20 15 年 12 月 21 日安徽大学电子学院通信工程系2目录一、 设计的目的和意义 -31.1 设计的目的和意义 -31.2 关于 MATLAB 软件： -31.3 课题要求及分工 -3二、 实现方案 -32.1 总体方案 -32.2 设计流程图 -4三、设计原理 -53.1 回声部分设计 -53.2 混响部分设计 -63.3 和声部分设计 -6四、设计过程 -64.1 设计方案 -64.2 实验程序 -74.2.1 原始信号分析 -74.2.2 回声信号的产生与分析 -94.2.3 混响的产生与分析 -134.2.4 和声的产生与分析 -164.4.5 简易的界面设计 -184.5 实验结果分析 -183五、设计总结 -185.1 课程设计的可取之处： -185.2 课程设计的不足 -185.3 编程中遇到的问题 -18六、收获与感受 -19七、参考文献 -1941、设计的目的和意义1.1 设计的目的和意义现在随着时代的发展，人类已开始进入了信息化时代，用现代手段研究语音信号，使人们能更加有效地产生、传输、存储、获取和应用语音信息，这对于促进社会的发展具有十分重要的意义。学习了数字信号处理以及相关课程，就想自己动手做一段音频处理，来锻炼自己的动手能力，把知识学以致用。一般的录音是单声道的，放音时也是由一个扬声器发出声音，感觉不到立体感，因此就需要我们对声音进行处理，使得声音信号在平时听起来也可以产生立体感。 综合所学的数字信号处理知识，利用 MATLAB 语言编程进行单回声、多回声、混音和和声的设计，并对处理后的信号进行频谱分析，实现的单回声、多回声、混音及和声的产生。同时通过处理音乐，我们也能更加的了解编码，熟悉运用 MATLAB 软件来处理音乐、声音，并熟悉数字信号处理方面的理论知识。1.2 关于 MATLAB 软件：1.MATLAB 语言具备高效的数值计算及符号计算功能，能使用户从繁杂的 数学运算分析中解脱出来；2. 具有完备的图形处理功能，实现计算结果和编程的可视化；3. 友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；4. 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) 5，为用户提供了大量方便实用的处理工具。1.3 课题要求1.课题要求(1).录制采集一段自己的语音信号；(2).设计滤波器，实现单回声、多回声、混响效果和和声效果；2、实现方案2.1 总体方案首先，对原始的音频信号进行频谱分析；然后通过数字信号处理对原声信号进行延时，衰减，产生单回声信号及多回声信号，通过图形与原声信号进行比较。混响效果主要是用于增加音源的融合感。自然音源的延时声阵列非常密集、复杂，所以模拟混响效果也复杂多变。声音在一个密闭的控件到达听众时，包含直达声、早期的反射和混响几部分。早期的反射有几个空间上临近的直达声的基本延迟和衰减组成，而混响则有密集的回声组成。当振幅有少量变换并声音之间有很少的时间差别时，会产生和声效果。延时可用如下函数实现： )(1(2)( 0nCOSRn注意：在时刻 n 处 的值通常是一个非整数值，在实际视线中，)(n输出样本值需要用到内插法进行计算。62.2 设计流程图1.通过计算机中的 y=fft(x)傅里叶变换，对原声信号和回声信号进行频谱分析。2.通过设计 FIR、IIR 滤波器对原声信号进行衰减、延时产生一次、打开界面选择功能原信号分析输出原信号分析波形原信号取样设计一个 FIR滤波器，进行延时及衰减处理，得到回声信号。频谱分析，绘出信号的频谱图以及时域波形单回声 设计一个 IIR滤波器，进行延时及衰减处理，得到回声信号。将每一级回声相对于上一级回声一次衰减和延时R 个周期频谱分析，绘出信号的频谱图以及时域波形多回声 设计全通滤波器，进行延时和叠加，达到混响效果。频谱分析，绘出信号的频谱图以及时域波形混响 设计滤波器，把原信号进行不同延迟和衰减，叠加达到和声效果。频谱分析，绘出信号的频谱图以及时域波形和声7多次（次数可通过 N 改变，此处以三次为例）。三、设计原理3.1 回声部分设计原理回声的定义：当声音投射到距离声源有一段距离的大面积上时，声能的一部分被吸收，而另一部分声能会被反射会俩，如果听者听到由声源直接发出的声和反射回来的声的时间间隔超过 1/10 秒，就能分辨这是两个声音，将返回的声音成为“回声”。1、单回声滤波器：回声可以由简单的延时单元产生。直达声和在 R 抽样周期后出现的一种单个回音可以用 FIR 滤波器产生。微分方程为：y(n)=x(n)+x(n-R)该差分方程的系统函数可表示为：H（z）=1+ Rz传递函数的幅频响应形状像梳子，这种滤波器又叫做梳状滤波器 此滤波器为 FIR 滤波器，其结构图如下：图 2.2FIR 滤波器结构2、多重回声滤波器：y(n) x(n) Rz8为了产生以间隔 R 个抽样周期分开的具有指数衰减振幅的多重回声，可用一个以下形式传输函数的 IIR 滤波器： RNRNRR zzzzzH 1)1(121)( 无限个振幅以指数衰减间隔为 R 个抽样周期的多重回声可用以下形式传输函数的 IIR 滤波器生成：1,1)(RzzHIIR 滤波器结构图:图 2.3IIR 滤波器结构图3.2 混响部分设计原理混响则有密集的回声组成。上述多重回声滤波器不能提供自然声音混响。有其幅频特性可知其幅度响应对于所有频率不是常数，不能达到应产生的效果。为了产生更真实的混响效果，采用全通结构的混响器器，其传输函数为： 1,1)(RzzH3.3 和声部分设计原理x(n) y(n)Rz9和声是对原始信号进行不同程度的衰减和延时，最后这些信号和原信号进行叠加会产生和声效果。延时可用如下函数实现： ncos12)(0Rn四、设计过程4.1 设计方案1.读取语音信号与回放一段合适的原声信号，用“x,fs=audioread(d:v4.mp3);”来读取录制好的语音信号，这是以采样频率 fs 来读取语音，并将采样信号值放入 x 中。用 sound(x,fs)对声音信号进行回放。2.傅里叶变换 FFT傅里叶变换调用格式：y=fft(x)其中，x 是序列，y 是序列的 FFT,x 可以为一向量或矩阵，若x 为一向量，y 是 x 的 FFT。且和 x 相同长度。若 x 为一矩阵，则 y是对矩阵的每一列向量进行 FFT。如果 x 长度是 2 的幂次方，函数fft 执行高速基 -2FFT 算法；否则 fft 执行一种混合基离散傅里叶变换算法，计算速度大较慢。3.函数 filter函数 filter 的调用格式为：y=filter(b,a,x)该格式采用数字滤波器对数据进行滤波，既可以用于 IIR 滤波器，也可以用于 FIR 滤波器。其中向量 b 和 a 分别表示系统函数的分子、分母多项式的系数，若 a=1,此时表示 FIR 滤波器，否则就是10IIR 滤波器。该函数是利用给出的向量 b 和 a,对 x 中的数据进行滤波，结果放入向量 y。4.函数 zearos函数的一种调用格式为：B=zeros(n)返回一个 n*n 的矩阵，其元素都为 0，如果参数 n 不是一个标量就会报错。函数的另一种调用格式为：B=zeros(m,n)或 B=zeros(m n)返回一个 m*n 的零矩阵。函数的另一种调用格式为：B=zeros(m,n,p,.)或 B=zeros(m n p .)返回一个 m*n*p*.的零矩阵。注意：参数 m,n,p 都是非负整数，负整数都会看成 0。4.2 实验程序4.2.1