基于MATLAB的声音的延迟与混响

1、课程设计作业书标题：基于 MATLAB 的声音延迟和混响初始状态：MATLAB软件数字信号处理和声音处理基础需要的主要任务：对常用的数据分析处理原理和方法有较全面的了解；能够使用相关软件进行仿真分析；掌握文献检索和文献阅读的基本方法；提高正确写论文的基本能力。参考书目：1.全，数字信号处理原理与实现（第二版），电子工业，2009；2. Majestic，DSP集成开发与应用实例，第1版，电子工业，2002；3. 洪涛等，数字信号处理，第一版，华中科技大学，20074. 齐聪等，DSP技术的发展与应用，第二版，高等教育，20075. 舒斌等，数字信号处理实战教程，第一版，华中科技大学，2007日

2、程17 周模拟设计18周箭竹15楼防守讲师签名：_年月日系主任（或负责教师）签名：_年月日目录TOC o 1-3 h u HYPERLINK l _Toc356504307 摘要PAGEREF _Toc356504307 h3 HYPERLINK l _Toc356504308 1 简介PAGEREF _Toc356504308 h5 HYPERLINK l _Toc356504309 1.1 设计任务PAGEREF _Toc356504309 h5 HYPERLINK l _Toc356504310 1.2 设计要求PAGEREF _Toc356504310 h6 HYPERLINK l _

3、Toc356504311 2 总体方案及设计原则PAGEREF _Toc356504311 h7 HYPERLINK l _Toc356504312 2.1 设计总体方案 图PAGEREF _Toc356504312 h7 HYPERLINK l _Toc356504313 2.2 程序流程图PAGEREF _Toc356504313 h8 HYPERLINK l _Toc356504314 2.3 详细设计原则PAGEREF _Toc356504314 h9 HYPERLINK l _Toc356504315 2.3.1 信号采样PAGEREF _Toc356504315 h9 HYPERL

4、INK l _Toc356504316 2.3.2 混响和延迟PAGEREF _Toc356504316 h9 HYPERLINK l _Toc356504317 2.3.3 离散傅里叶变换PAGEREF _Toc356504317 h10 HYPERLINK l _Toc356504318 2.3.4 滤波器设计PAGEREF _Toc356504318 h11 HYPERLINK l _Toc356504319 3 编程的步骤和过程PAGEREF _Toc356504319 h13 HYPERLINK l _Toc356504320 3.1 原始语音信号的时域波形和频谱PAGEREF _T

5、oc356504320 h13 HYPERLINK l _Toc356504321 3.2 采样语音信号的时域波形和频谱PAGEREF _Toc356504321 h13 HYPERLINK l _Toc356504322 3.3 延迟采样信号PAGEREF _Toc356504322 h14 HYPERLINK l _Toc356504323 3.4 混响采样信号PAGEREF _Toc356504323 h14 HYPERLINK l _Toc356504324 3.5 单回波滤波器编程PAGEREF _Toc356504324 h15 HYPERLINK l _Toc356504325

6、3.6 多重回波滤波器编程PAGEREF _Toc356504325 h15 HYPERLINK l _Toc356504326 3.7 无限回波滤波器编程PAGEREF _Toc356504326 h16 HYPERLINK l _Toc356504327 3.8 全通结构化滤波器编程PAGEREF _Toc356504327 h17 HYPERLINK l _Toc356504328 4 结果分析与经验PAGEREF _Toc356504328 h18 HYPERLINK l _Toc356504329 5 设计经验PAGEREF _Toc356504329 h25 HYPERLINK l

7、 _Toc356504330 6 主要参考文献PAGEREF _Toc356504330 h26 HYPERLINK l _Toc356504331 附录：通用仿真程序PAGEREF _Toc356504331 h27摘要数字信号处理技术自诞生以来，发展迅速。它主要研究用数字或符号序列对信号的表示和处理，广泛应用于各个领域。本设计采用MATLAB语言对语音信号进行采样分析，设计数字滤波器对信号进行滤波，比较滤波前后信号特性的变化。使用 MATLAB 开发环境设计图形用户界面可简化布局编程。1 绪论数字信号处理是随着计算机技术的发展而迅速发展起来的一门新兴而古老的学科。它发挥着重要的作用，其方法

8、多样，应用众多，源远流长；也是一门理论性强、实践性强、理论技术发展较快的前沿学科。随着数字时代的到来，科学技术的进步和生产开发需求的不断增加，推动了数字信号处理学科的发展，产生了各种巧妙的信号处理算法；强大的生命力。数字信号处理主要是研究用数字或符号序列来表示和处理信号。处理的目的可以是衰减信号中的过剩容量，滤除混合的噪声和干扰，或者将信号变换为易于分析和识别的形式，便于对其特征参数的估计和选择。例如，通过分析计算，可以估计出EEG或ECG中的某些特征参数，帮助医生找到原因，分析病情，确定合理的治疗方案；另一个例子，当信号传输时，它会受到各种干扰，包括失真。 、衰落和混合背景噪声、信号处理以消

9、除这些干扰。声音信号是一维连续信号，而计算机只能处理离散信号。为了从离散信号中恢复出连续信号，根据采样定理，可以确定采样频率的最小值。 wav 文件是一种数字声音文件格式。本课程的设计是基于Matlab来分析wav声音文件的频谱与声音的关系。通过采集一段个人声音进行频谱分析等处理，然后设计一个数字滤波器对原声的wav文件进行处理，比较滤波后输出的声音信号与原声信号的异同。1.1 设计任务(1)使用录音机或Windows下的其他软件录制一段自己的语音信号，时间控制在1s左右，对录制的信号进行采样。(2)语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图。(3) 对信号添加延迟和混响，然后分

10、析其频谱，并与原始信号频谱进行比较。(4)设计几种特殊类型的滤波器：单回波滤波器、多回波滤波器、无限回波滤波器、全通结构的混响器，并画出滤波器的频域响应。(5) 用自己设计的滤波器对采集到的语音信号进行滤波。(6)对得到的信号进行频谱分析，画出滤波后信号的时域波形和频谱，比较滤波前后的信号，分析信号的变化。(7) 播放语音信号。1.2 设计要求(1) 熟悉离散信号和系统的时域特性。(2)熟悉语音信号的特点。（3）掌握数字信号处理的基本概念、基本理论和基本方法。(4)掌握序列快速傅里叶变换方法。(5)学习使用MATLAB，掌握MATLAB编程方法。（6）掌握MATLAB中设计各种数字滤波器的方法

11、和滤波信号的方法2 总体方案及设计原则2.1 设计总体方案语音信号延迟与混响设计方案框图12.2 程序流程图语音信号延迟与混响设计程序流程图22.3 详细设计原则2.3.1 信号采样(1) 采样频率采样频率是指计算机每秒收集多少声音样本。采样频率越高，采样间隔越短，计算机在单位时间内可以获得的声音样本数据越多，声音波形的表示就越准确。只有当采样频率高于声音信号最高频率的两倍时，才能将数字信号所代表的声音还原为原始声音。(2) 采样位即采样值或采样值，用于衡量声音波动的参数，是指声卡在采集和播放声音文件时所使用的数字声音信号的二进制位。声卡的位数客观地反映了数字声音信号描述输入声音信号的准确程度

12、。声卡的主要功能之一是记录和播放声音信息。在这个过程中，比特数和采样频率决定了声音采集的质量。2.3.2 混响和延迟(1)混响效果主要用于增加声源的融合感。自然声源的延迟声阵列非常密集和复杂，因此模拟混响效果的程序也复杂多样。常用参数如下：混响时间：可以真实模拟自然混响的数字混响具有一组复杂的程序。虽然有很多技术参数可以调整，但这些技术参数的调整不会比原来更有效。自然，尤其是混响时间。高频滚降：该参数用于模拟自然混响，空气对高频的吸收效果，从而产生更自然的混响效果。一般高频混频的可调范围为0.11.0。当该值较高时，混响效果更接近自然混响；当此值较低时，混响效果更清晰。Diffusion：这个

13、参数可以调整混响声阵列密度的增长速度。可调范围为010。数值高时，混响效果更丰富、更温暖；值越低，混响效果越强烈。空虚而寂寞。预延迟：自然混响声阵的建立会延迟一段时间，预延迟设置为模拟二次效果。Sound Array Density：此参数可以调整声音阵列的密度。数值高时，混响效果偏暖，但声音染色明显；值越低，混响效果越深，音色越明显。虚弱的。调频：这是一个技术参数，因为电子混响的声音阵列密度比自然混响要稀疏。为了使混响的声音更加流畅和连贯，需要对混响声阵列的延迟时间进行调制。该技术可以有效消除延迟声阵列的分裂声，增加混响声的柔和度。调整深度：指上述调频电路的调整深度。(2)延迟是音源延迟一段

14、时间后要播放的效果处理。根据延迟时间，可以单独产生合唱、镶边和回声等效果。当延迟时间在335ms之间时，人耳无法感知滞后声的存在，与原声源叠加后，会因其相位干扰而产生“梳状滤波器”效应，即翻边效果。如果延迟时间超过50ms，则延迟声音清晰可辨，此时的处理效果就是回声。回声处理通常用于产生简单的混响效果。延迟、合唱、镶边、回声等效果的可调参数类似，具体项目如下：*延迟时间（Dly），即主延迟电路延迟时间的调整。*反馈增益（FB Gain），即延迟反馈的增益控制。*反馈高频比（Hi Ratio），即对反馈回路的高频衰减控制。*调制频率（Freq），是指主延迟的频率调制周期。*调制深度（Depth）

15、是指上述FM电路的调制深度。*高频增益（HF），是指高频均衡控制。*预延迟（Ini Dly）是指调整主延迟电路的预延迟时间。*均衡频率（EQ F），这里的频率均衡是用来调节音色的，也就是均衡的中点频率选择。由于延时产生的效果比较复杂多变，如果你不是效果处理专家，建议使用设备提供的预设参数，因为这些预设参数给出的处理效果一般比较好。2.3.3 离散傅里叶变换FFT 和 IFFT 函数用于 MATLAB 信号处理工具箱中的快速傅里叶变换和逆变换。这些功能如下所述。函数 FFT 用于序列的快速傅里叶变换。调用函数的一种形式是 y=fft(x)其中，x为序列，y为序列的FFT，x可为向量或矩阵，若x为

16、向量，y为x的FFT。和 x 的长度相同。如果 x 是矩阵，则 y 是矩阵的每个列向量的 FFT。如果 x 的长度是 2 的幂，则函数 fft 执行高速 radix-2FFT 算法；否则，fft 执行混合基数离散傅里叶变换算法，计算速度较慢。FFT函数的另一种调用格式是y=fft(x,N)式中x和y的含义同上，N为正整数。该函数执行 N 个点的 FFT。如果 x 是长度小于 N 的向量，则函数将 x 补零到长度 N。如果向量 x 的长度大于 N，则函数截断 x 以使其长度为 N。如果 x 是矩阵，以同样的方式处理 x。函数fft得到的序列y一般是一个复数序列，通常需要它的幅值和相位。 MATL

17、AB提供了复杂的幅度和相位函数：abs、angle，这些函数一般和FFT同时使用。函数 abs(x) 用于计算复矢量 x 的大小，函数 angle(x) 用于计算复矢量的相位角，弧度之间和弧度之间。函数 unwrap(p) 用于展开以弧度为单位的相位角 p，当相位角的绝对变化超过时，函数将其展开为。使用 MATLAB 工具箱函数 fft 进行频谱分析时请注意：(1)函数fft的返回值y的数据结构对称性。(2) 频率计算。(3)在做FFT分析时，幅度的大小与FFT选择的点数有关，但不影响分析结果。2.3.4 滤波器设计本设计使用 IIR 滤波器对信号进行滤波，函数名称为 filter函数filt

18、er的调用格式为y=filter(b,a,x)这种格式使用数字滤波器来过滤数据，可用于 IIR 和 FIR 滤波器。向量b和a分别表示系统函数的分子和分母多项式的系数。如果a=1，则代表FIR滤波器，否则为IIR滤波器。作用是使用给定的向量 b 和 a 对 x 中的数据进行过滤，并将结果放入向量 y 中。(1) 单回波滤波器回声可以由简单的延迟单元产生。在 R 采样周期之后出现的直达声和单个回声可以通过 FIR 滤波器产生，其微分方程为：yn=xn+xnR |1传递函数为：H(z)=1+z -R传递函数的幅频响应呈梳状，这种滤波器也称为梳状滤波器。(2) 多重回波滤波器为了生成多个回波，其幅度

19、呈指数衰减，间隔 R 个采样周期，FIR 滤波器具有以下形式的传递函数：可以使用具有以下形式的传递函数的 IIR 滤波器生成无限数量的具有 R 个采样周期的指数衰减间隔的多重回波：IIR 多回声滤波器的基频 F R =F s /R 通常被锁定为伴随声音设备的基频，例如基鼓节拍。(3) 无限回波滤波器的系统函数：, 1(4)全通结构混响器系统功能：全通滤波器的传递函数公式为：H(Z)=y(Z)/X(Z)=(-K+Z(-m)/(1-K*Z(-m)其中 m 是回波延迟采样，k 是反馈系数。以直接 1 形式表示此传递函数为：y(n)=k*x(n)+x(nm)+k*y(nm)可以看出其实就是一个简单的I

20、IR滤波器，n时刻的输出是根据n时刻的输入和m点之前的输入输出计算出来的。由于该 IIR 滤波器的频率响应是一条水平线，因此称为全通滤波器。3 程序设计的步骤和过程3.1 原始语音信号的时域波形和频谱y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ); %用于读取语音，样本值放在向量y中，fs代表采样频率（Hz），bits代表采样位数。 50000 150000表示读取50000点到150000点的值（如果只有一个N点，表示读取前N点的采样值）。y=y(:,1);% 只取单通道声音（y，fs，位）； %播放原始信号Y=

21、fft(y); % 对原始信号进行傅里叶变换subplot(2,2,1:2);plot(y);title( 原始信号时域波形 );subplot(2,2,3);plot(abs(Y);title( 原始信号幅度和频率 );subplot(2,2,4);plot(angle(Y);title( 原始信号相位频率 );3.2 语音采样信号的时域波形和频谱y,fs,bits=wavread( C:UsersAdministratorDesktopaudio signal.wav ,50000 150000); % 只取原始信号的一部分，即时间间隔50000 150000y=y(:,1);Y=fft(

22、y,6001);% 提取 6001 个点subplot(2,2,1:2);plot(y);title( 采样信号时域波形 );subplot(2,2,3);plot(abs(Y);title( 采样信号幅度和频率 );subplot(2,2,4);plot(angle(Y);title( 采样信号相位频率 );声音（y，fs，位）； %回放采样信号3.3 延迟采样信号y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y; % 延迟 3000Z=fft

23、(z,6001);subplot(2,2,1:2);plot(z);title( 延迟后的时域波形 );subplot(2,2,3);plot(abs(Z);title( 延迟后的幅频 );subplot(2,2,4);plot(angle(Z);title( 延迟后的相位频率 );声音（y，fs，位）； % 播放延迟信号3.4 混响采样信号y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1); % 补码 3000

24、y1=y+z； % 将原始信号和延迟信号相加生成混响信号Y1=fft(y1,6001);subplot(2,2,1:2);plot(y1);title( 混响时域波形 );subplot(2,2,3);plot(abs(Y1);title( 混响幅度频率 );subplot(2,2,4);plot(angle(Y1);title( 混响的相位频率 );声音（y，fs，位）； % 播放混响信号3.5 单回波滤波器编程y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3

25、000,1);y;y=y;zeros(3000,1);a=0.5；R=5000;%过滤阶数设置，数值越高，回声越明显Bz1=1,zeros(1,R-1),a; % 单回波滤波器系统函数分子az1=(1); % 单回波滤波器系统函数分母y2=过滤器（Bz1，Az1，y）； % 单回波过滤系统功能Y2=fft(y2,6001);h,w=freqz(Bz1,Az1); % 找到设计的滤波器光谱子图（3,2,1）；情节（abs（h）；title( 单回波滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 单回波滤波器相位频率响应 );subplot(3,2,3:4);plot(y2

26、);title( 单回波滤波器时域图 );subplot(3,2,5);plot(abs(Y2);title( 单回波滤波器幅频 );subplot(3,2,6);plot(angle(Y2);title( 单回波滤波器相位频率 );声音（y2，fs，位）；3.6 多个回波滤波器的编程y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1);a=0.5；N=5；R=9000;%过滤阶数设置，数值越高，回声越明显Bz1=

27、1,zeros(1,N*R-1),-aN; %多重回波滤波器系统函数分子Az1=1,zeros(1,R-1),-a; % 多重回波滤波器系统函数的分母y2=过滤器（Bz1，Az1，y）；Y2=fft(y2,6001);h,w=freqz(Bz1,Az1);子图（3,2,1）；情节（abs（h）；title( 多重回波滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 多重回波滤波器相位频率响应 );subplot(3,2,3:4);plot(y2);title( 多重回波滤波器时域图 );subplot(3,2,5);plot(abs(Y2);title( 多重回波滤波器幅

28、度和频率 );subplot(3,2,6);plot(angle(Y2);title( 多重回波滤波器相位频率 );声音（y2，fs，位）；3.7 无限回波滤波器编程y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1);a=0.5；R=5000;%过滤阶数设置，数值越高，回声越明显Bz=0,zeros(1,R-1),1; % 无限回波滤波器系统函数分子Az=1,zeros(1,R-1),-a; % 无限回波滤波器系

29、统函数分母y2=过滤器（Bz，Az，y）；Y2=fft(y2,6001);h,w=频率z(Bz,Az);子图（3,2,1）；情节（abs（h）；title( 无限回波滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 无限回波滤波器相位频率响应 );subplot(3,2,3:4);plot(y2);title( 无限回波滤波器时域图 );subplot(3,2,5);plot(abs(Y2);title( 无限回波滤波器信号幅度 );subplot(3,2,6);plot(angle(Y2);title( 无限回波滤波器信号相位频率 );声音（y2，fs，位）；3.8 全

30、通结构滤波器编程y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1);a=0.5；R=5000;%过滤阶数设置，数值越高，回声越明显Bz1=a,zeros(1,R-1),1; % 全通滤波器系统函数分子Az1=1,zeros(1,R-1),a; % 全通滤波器系统函数分母h,w=freqz(Bz1,Az1);yy2=过滤器（Bz1，Az1，y）；YY2=fft(yy2,6001);子图（3,2,1）；情节（abs（

31、h）；title( 全通滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 全通滤波器相位频率响应 );subplot(3,2,3:4);plot(yy2);title( 全通结构时域图 );subplot(3,2,5);plot(abs(YY2);title( 全通结构幅频 );subplot(3,2,6);plot(angle(YY2);title( 全通结构相位频率 );声音（yy2，fs，位）；4 结果分析与经验运行结果：对信号进行采样后，信号的时域波形呈离散分布，不像原始图像那样分布密集。频谱图也以同样的方式进行变换，幅度和频率根据采样间隔而变化。时域图的延迟是把

32、原来的波形向右移动了3000，频域的幅度发生了变化。混响：在时域图中，前面几乎没有变化，后面的幅度有一定的变化。谱图幅度变化较大，中间上下相互抵消，两侧跳动较大。在时域中，形状大致相同，幅度变化很小。在频谱中，滤波前后两侧的拍频变换也较大，幅度变化较大。时域滤波前幅度大于滤波后值，滤波后形状连续，滤波前不连续光谱上的变化是在光谱两侧上下跳跃的动态变化，中间的变化相互抵消在时域，波形变得稀疏，幅度变化比较大。光谱上的变化是在光谱两侧上下跳跃的动态变化，中间的变化相互抵消全通结构的时域大致相同，滤波前的幅度比滤波后的幅度变化了6000倍。频域变化不明显，幅度发生了变化。5 设计经验课程设计是培养我

33、们综合运用所学、发现、提出、分析和解决实际问题，锻炼实践能力的重要环节。这是对我们实际工作能力的具体培训和检验过程。理论知识，理论与实践相结合，并在生活中应用。在设计的过程中，总会遇到各种各样的问题。在这种情况下，我们会努力寻找解决问题的最佳方法，这无形中提高了我们的动手能力和大脑能力，学生们也可以互相讨论问题，研究解决方案，增强大家的团队意识。本课程设计结合了数字信号处理的理论知识和实践经验。期间复习了相关资料，结合之前学过的数字信号处理课程和实验，如期完成了设计任务要求，加深了自己的认识。了解并提高他们的实践能力。当我接到课程设计任务时，我不知道从哪里开始。选好题目后，查了网上的相关资源，

34、去图书馆取材，确定了基本的设计方案和步骤，开始用软件实现。问了很多问题，上网查资料，在专业论坛看，下载MATLAB视频，慢慢开始有了更好的了解。在之前的实验中，没有接触到GUI界面的问题，也不知道从何下手。查阅了资料，看了视频，有了一定的了解。但是，由于是第一次接触，出现问题时仍然很难解决。经过和同学们的讨论，最终解决了这些问题。通过本课程设计，我对数字信号处理的理论知识有了更深入的了解，对使用MATLAB语言对信号进行采样和滤波有了更深入的了解。第一次接触到MATLAB提供的用户界面设计，对一些基本问题有所了解，为以后的实践打下了一定的基础。相信在以后的实践过程中会有更好的改进。另外，这门课

35、的设计让我更加意识到，平时学到的知识，如果不付诸实践，就一文不值。课程设计主要是我们理论知识的延伸。其目的是在设计中发现问题，并能找到解决问题的办法，形成独立的意识。我们还可以测试我们从设计中学到了多少理论知识，巩固所学，继续学习我们错过的新知识，使这门课程扎实。6 主要参考文献1.全，数字信号处理原理与实现（第二版），电子工业，2009；2. Majestic，DSP集成开发与应用实例，第1版，电子工业，2002；3. 洪涛等，数字信号处理，第一版，华中科技大学，20074. 齐聪等，DSP技术的发展与应用，第二版，高等教育，20075. 舒斌等，数字信号处理实战教程，第一版，华中科技大学，

36、2007附录：向导界面下的通用仿真程序% 在按钮 1 中按下按钮时执行。函数pushbutton1_Callback(hObject, eventdata, 句柄)% hObject 到 pushbutton1 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使用句柄和用户数据处理结构（参见 GUIDATA）y,fs,bits=wavread( C:UsersAdministratorDesktopaudio signal.wav ); %读取原始音频信号y=y(:,1); % 只占用一个通道声音（y，fs，位）； %播放原始信号Y=

37、fft(y); % 对原始信号进行傅里叶变换subplot(2,2,1:2);plot(y);title( 原始信号时域波形 );subplot(2,2,3);plot(abs(Y);title( 原始信号幅度和频率 );subplot(2,2,4);plot(angle(Y);title( 原始信号相位频率 );% 在按钮 2 中按下按钮时执行。函数pushbutton2_Callback(hObject, eventdata, 句柄)% hObject 到 pushbutton2 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使

38、用句柄和用户数据处理结构（参见 GUIDATA）y,fs,bits=wavread( C:UsersAdministratorDesktopaudio signal.wav ,50000 150000); % 只取原始信号的一部分，即时间间隔50000 150000y=y(:,1);Y=fft(y,6001); % 提取 6001 点subplot(2,2,1:2);plot(y);title( 采样信号时域波形 );subplot(2,2,3);plot(abs(Y);title( 采样信号幅度和频率 );subplot(2,2,4);plot(angle(Y)频);title( 相位信号相

39、 );声音（y，fs，位）； % 反信号标志% 在按钮 3 中按下按钮时执行。函数pushbutton3_Callback(hObject, eventdata, 句柄)% hObject 到 pushbutton3 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使用句柄和用户数据处理结构（参见 GUIDATA）y,fs,bits=wavread( C:UsersAdministratorDesktop音频信号.wav ,50000 150000);y=y(:,1);z=零(3000,1);y; %延迟3000Z=fft(z,600

40、1);subplot(2,2,1:2);plot(z);title( 延迟后的时域波形 );subplot(2,2,3);plot(abs(Z);title( 延迟后的幅频 );subplot(2,2,4);plot(angle(Z);title( 延迟后的相位频率 );声音（y，fs，位）； % 播放延迟信号% 在按钮 4 中按下按钮时执行。函数pushbutton4_Callback（hObject，事件数据，句柄）% hObject 到 pushbutton4 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使用句柄和用户数据处

41、理结构（参见 GUIDATA）y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1); % 补码 3000y1=y+z； % 将原始信号和延迟信号相加生成混响信号Y1=fft(y1,6001);subplot(2,2,1:2);plot(y1);title( 混响时域波形 );subplot(2,2,3);plot(abs(Y1);title( 混响幅度频率 );subplot(2,2,4);plot(angle(Y

42、1);title( 混响的相位频率 );声音（y，fs，位）； % 播放混响信号% 在按钮 5 中按下按钮时执行。函数pushbutton5_Callback（hObject，事件数据，句柄）% hObject 到 pushbutton5 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使用句柄和用户数据处理结构（参见 GUIDATA）y,fs,bits=wavread( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y

43、=y;zeros(3000,1);a=0.5；R=5000； %过滤顺序设置，数值越高，回声越明显Bz1=1,zeros(1,R-1),a; % 单回波滤波器系统函数分子az1=(1); % 单回波滤波器系统函数分母y2=过滤器（Bz1，Az1，y）； % 单回波过滤系统功能Y2=fft(y2,6001);h,w=freqz(Bz1,Az1); % 找到设计的滤波器光谱子图（3,2,1）；情节（abs（h）；title( 单回波滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 单回波滤波器相位频率响应 );subplot(3,2,3:4);plot(y2);title(

44、单回波滤波器时域图 );subplot(3,2,5);plot(abs(Y2);title( 单回波滤波器幅频 );subplot(3,2,6);plot(angle(Y2);title( 单回波滤波器相位频率 );声音（y2，fs，位）；% 在按钮 6 中按下按钮时执行。函数pushbutton6_Callback(hObject, eventdata, 句柄)% hObject 到 pushbutton6 的句柄（见 GCBO）% eventdata reserved - 将在未来版本的 MATLAB 中定义% 使用句柄和用户数据处理结构（参见 GUIDATA）y,fs,bits=wavr

45、ead( C:UsersAdministratorDesktopAudioSignal.wav ,50000 150000);y=y(:,1);z=零(3000,1);y;y=y;zeros(3000,1);a=0.5；N=5；R=9000； %过滤顺序设置，数值越高，回声越明显Bz1=1,zeros(1,N*R-1),-aN; %多重回波滤波器系统函数分子Az1=1,zeros(1,R-1),-a; % 多重回波滤波器系统函数的分母y2=过滤器（Bz1，Az1，y）；Y2=fft(y2,6001);h,w=freqz(Bz1,Az1);子图（3,2,1）；情节（abs（h）；title( 多

46、重回波滤波器幅频响应 );子图（3,2,2）；情节（角度（h）；title( 多重回波滤波器相位频率响应 );subplot(3,2,3:4);plot(y2);title( 多重回波滤波器时域图 );subplot(3,2,5);plot(abs(Y2);title( 多重回波滤波器幅度和频率 );subplot(3,2,6);plot(angle(Y2);title( 多重回声滤波器相频 );声音（y2，fs，位）；% 在按钮 7 中按下按钮时执行。函数pushbutton7_Callback(hObject, eventdata, 句柄)% hObject 到 pushbutton7 的句柄