基于matlab的回声信号产生器

《基于matlab的回声信号产生器》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、WORD完美格式电子信息工程学院《DSP技术及应用》课程设计报告题目：基于MATLAB的回声信号产生器专业班级：通信工程专业　　　　　二〇一三年六月十六日目录一、设计目的…………………………………………………………1技术资料专业整理WORD完美格式二、设计要求…………………………………………………………1三、设计原理及方案…………………………………………………1四、设计流程…………………………………………………………5五、调试分析…………………………………………………………5六、GUI界面…………………………………………………………17七、设计总结.………

2、………………………………………………..20技术资料专业整理WORD完美格式一、设计目的综合运用我们掌握的数字信号处理知识，利用matlab语言编程进行回声信号产生器的设计再对其进行频谱分析，实现回声的产生，消除和隐藏。1、首先熟练掌握对声音采集、处理、衰减、延时等过程的理论知识。2、根据设计思路，学习资料的搜集与整理。3、通过理论与实际的结合，在理论的基础结合实际情况完成课程的设计。关于MATLAB软件：MATLAB语言具备高效、可及推理能力强等特点,是目前工程界流行最广泛的科学计算语言。MATLAB强大的运算和图形显示功能，可使信号与系统上机实验效率大

3、大提高。特别是它的频谱分析和滤波器分析与设计功能很强，使数字信号处理工作变得十分简单，直观。关于GUI界面：图形用户界面GUI(GraphicsUserInterface)是由各种图形对象,如图形窗口、图轴、菜单、按钮、文本框等构建的用户界面,是人机交流信息的工具和方法。GUI设计即可以基本的MATLAB程序设计为主,也可以鼠标为主利用GUIDE工具进行设计。二、设计要求首先，自己录取一个合适的原声信号，并对其进行频谱分析；然后通过数字信号处理对已录制好的原声信号进行延时，衰减，产生出回声信号，通过图形与原声信号进行比较；再把信息转换成二进制信号将其叠加在

4、回声中，通过数字信号处理实现信息的隐藏；最后，通过设计GUI界面将结果展示出来。1、通过计算机中的y=fft（x）傅里叶变换，对原声信号和回声信号进行频谱分析。2、通过设计FIR、IIR滤波器对原声信号进行衰减、延时产生一次、三次、六次无限次回声。3、学习并掌握MATLAB的使用方法，用MATLAB编程和会使用MATLAB中的GUI界面设计。三、设计原理及方案技术资料专业整理WORD完美格式1、设计原理（1）回声的定义：当声投射到距离声源有一段距离的大面积上时，声能的一部分被吸收，而另一部分声能要反射回来，如果听者听到由声源直接发来的声和由反射回来的声的时

5、间间隔超过十分之一秒，它就能分辨出两个声音这种反射回来的声叫“回声”。（2）回声信号产生：根据回声的定义，回声信号可以是原始信号通过延时、衰减来生成。（3）回声信号消除：数字滤波器的作用是滤除信号中某一部分频率分量。可以将回声信号消除。（4）回声信号隐藏：回声隐藏利用了人耳的听觉掩蔽效应，是一种有效的音频信息隐藏方法。（5）处理单元原理：假设x(n)为原始信号，y(n)为经过延时、衰减迭加后产生的“回声信号”，声音信号和经延迟R个周期，衰减系数α（α<1）的单个回声的差分方程有如下表示:y(n)=x(n)+α·x(n-R)上述差分方程可用系统函数表示为：H

6、(z)=1+α·z-R此滤波器为FIR滤波器，其结构图如下：在实际生活中，往往不止存在一重回声，因此，需要实现多重回声的设计，关于均匀间隔的多重回声的实现，是将每一级回声相对于上一级回声依次衰减和延时R个周期，此时不能再用FIR滤波器来实现，应该用IIR滤波器设计。IIR滤波器设计的系统函数如下：技术资料专业整理WORD完美格式当回声次数很多时多重回声的系统函数如下：IIR滤波器结构图：IIR滤波器结构隐藏目的是以添加回声的方式在原有音频信息中嵌入新信息，实现信息隐藏。它的数学模型可以表示为：典型的回声核脉冲嵌入回声的声音可表示为和的卷积，和，分别为原始声

7、音信号和回声核的单位脉冲响应。回声信号由引入到原始声音当中，其中，d为延迟时间，为衰减系数。嵌入回声后的声音信号可以表示为技术资料专业整理WORD完美格式2、设计方案（1）读取语音信号与回放录制一段合适的原声信号，用“[x,fs,bits]=wavread('y',[140000])”来读取录制好的语音信号，这是以采样频率fs、采样位数bits从[140000]来读取语音，并将采样值放在x中。用sound(5*x,fs,bits)对声音信号进行回放，乘以5以加强信号，便于听取。（2）傅里叶变换FFT傅里叶变换调用格式:y=-fft(x)其中，x是序列，y是

8、序列的FFT，,x可以为一向量或矩阵，若x为一向量，y是x的FFT