《2022年基于MATLAB的声音的延迟与混响课程设计任务书.docx》由会员分享，可在线阅读，更多相关《2022年基于MATLAB的声音的延迟与混响课程设计任务书.docx（40页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精品学习资源课程设计任务书题目: 基于 MATLAB地声音推迟与混响初始条件 : MATLAB软件 数字信号处理与声音处理基础学问要求完成地主要任务 : 较全面明白常用地数据分析与处理原理及方法； 能够运用相关软件进行模拟分析； 把握基本地文献检索和文献阅读地方法； 提高正确地撰写论文地基本才能.参考书目 :1. 刘泉，数字信号处理原理与实现（其次版），电子工业出版社，2021；2. 张雄壮， DSP 集成开发与应用实例，第1 版，电子工业出版社，2002；3. 张洪涛等，数字信号处理，第一版，华中科技高校出版社，20074. 彭启琮等， DSP 技术地进展与应用，其次版，高等训练出版社，20
2、075. 杨述斌等，数字信号处理实践教程，第一版，华中科技高校出版社，2007时间支配17 周 仿真设计18 周 鉴主 15 楼答辩指导老师签名： 年 月日系主任（或责任老师）签名： 年月 日欢迎下载精品学习资源目录摘要 .31 绪论 .41.1 设计任务51.2 设计要求52 设计地总体方案与原理62.1 设计地总体方案图62.2 程序流程图72.3 设计地具体原理82.3.1 信号采样82.3.2 混响与延时82.3.3 离散傅立叶变换92.3.4 滤波器设计103 程序设计地步骤和过程113.1 原始语音信号时域波形和频谱113.2 采样后语音信号地时域波形和频谱123.3 对采样后地信
3、号延时123.4 对采样后地信号混响123.5 单回声滤波器程序设计133.6 多重回声滤波器程序设计143.7 无限个回声滤波器程序设计153.8 全通结构滤波器程序设计154 结果分析与体会.165 设计心得体会 .246 主要参考文献 .25附录：仿真总程序.25欢迎下载精品学习资源摘要数字信号处理技术自产生以来，有了快速地进展，主要是争论用数字或符号序列表示和处理信号，被广泛应用于各个领域.本次设计是用MATLAB语言对语音信号进行采样分 析，并设计数字滤波器对信号进行滤波，比较滤波前后信号特性地变化.用 MATLAB开发环境设计用户图形界面使布局编程简化.ABSTRACTDigita
4、l signal processing technology since its birth had rapid development. the technology isto欢迎下载精品学习资源study the sequence of numbers or symbols used to represent and signal processing, it is widely used in various field. the design is to use theMATLAB language sampling and analysis of voice signal, and
5、digital design for signal filtering, comparing before and after filteringthe signal properties change .GUIinterface providedwith MATLABdevelopment environment for designing the layout of graphical user interface makes programming simpler.1 绪论数字信号处理是随着运算机技术地进展而快速进展起来地一门新兴而古老地学科，它在新地领域如生物医学工程、声学、雷达、地震
6、不、语音通信、数据通信、核科学等学科发挥着重要地作用，而它所采纳地各种方法及众多应用已有悠久地历史；同时也是一门具有很强地理论性与实践性，且理论和技术进展都特别快速地前沿性学科.随着数字化时代地来临，科学技术地进步而生产进展需求地与日俱增，促进了数字信号处理学科地进展，产生了各种奇妙地信号处理算法；特殊是运算机技术地飞速进展，为数字信号处理增加了庞大地生命力 .数字信号处理主要是争论用数字或符号序列表示和处理信号.处理地目地可以是减弱信号中地余外内容，滤除混杂地噪声和干扰，或者是将信号变换为简洁分析和识别地势式，便于估量和挑选它地特点参数.例如通过分析和运算，可以估量脑电图或心电图中地某种特点
7、参数，帮忙医生查找病因和分析病情，确定合理地治疗方案；又如，信号在传输时，要欢迎下载精品学习资源受到各种干扰，包括失真、衰落和混入地背景噪声，信号处理要排除这些干扰.声音信号是一维连续信号,而运算机只能处理离散信号.为了从离散信号仍原连续信号, 依据采样定理 ,可以确定采样频率地最小值.wav 文件是一种数字声音文件格式,本课程设计基于 Matlab 分析了 wav 声音文件频谱与声音地关系.通过采集个人地一段声音进行频谱分析等处理，然后设计数字滤波器处理这个原始声音地wav 文件 ,并比较滤波以后输出声音信号与原声音信号地异同.1.1 设计任务(1) 利用 Windows 下地录音机或其他软
8、件，录制一段自己地语音信号，时间掌握在1s 左右，并对录制地信号进行采样.(2) 语音信号地频谱分析，画出采样后语音信号地时域波形和频谱图.(3) 将信号加入延时和混响，再分析其频谱，并与原始信号频谱进行比较.(4) 设计几种特殊类型地滤波器：单回声滤波器，多重回声滤波器，无限个回声滤波器， 全通结构地混响器，并画出滤波器地频域响应.(5) 用自己设计地滤波器对采集地语音信号进行滤波.(6) 分析得到信号地频谱，画出滤波后信号地时域波形和频谱，并对滤波前后地信号进行对比，分析信号地变化 .(7) 回放语音信号 .1.2 设计要求(1) 熟识离散信号和系统地时域特性. 2熟识语音信号地特点 .(
9、3) 把握数字信号处理地基本概念，基本理论和基本方法.(4) 把握序列快速傅里叶变换方法.(5) 学会 MATLAB地使用，把握 MA TLAB地程序设计方法 .(6) 把握 MATLAB设计各种数字滤波器地方法和对信号进行滤波地方法欢迎下载精品学习资源2 设计地总体方案与原理2.1 设计地总体方案图语音信号地延时和混响设计方案框图1欢迎下载精品学习资源2.2 程序流程图语音信号地延时和混响设计程序流程图2欢迎下载精品学习资源2.3 设计地具体原理2.3.1 信号采样（1） 采样频率采样频率是指运算机每秒钟采集多少个声音样本，采样频率越高，即采样地间隔时间越短，就在单位时间内运算机得到地声音样
10、本数据就越多，对声音波形地表示也越精确.只有采样频率高于声音信号最高频率地两倍时，才能把数字信号表示地声音仍原成为原先地声音.（2） 采样位数即采样值或取样值，用来衡量声音波动变化地参数，是指声卡在采集和播放声音文件时所使用数字声音信号地二进制位数.声卡位置客观地反映了数字声音信号对输入声音信号描述地精确程度 .声卡地主要地作用之一是对声音信息进行录制与回放，在这个过程中采样位置数和采样地频率打算了声音采集地质量.2.3.2 混响与延时(1) 混响成效主要是用于增加音源地融合感 .自然音源地延时声阵列特别密集、复杂，所以模拟混响成效地程序也复杂多变 .常见参数有以下几种：混响时间：能逼真地模拟
11、自然混响地数码混响器上都有一套复杂地程序，其上虽然有许多技术参数可调，然而对这些技术参数地调整都不会比原有地成效更为自然，特殊是混响时间 .高频滚降：此项参数用于模拟自然混响当中，空气对高频地吸取效应，以产生较为自然地混响成效 .一般高频混降地可调范畴为0.11.0. 此值较高时，混响成效也较接近自然混响；此值较低时，混响成效就较清亮.扩散度：此项参数可调整混响声阵密度地增长速度，其可调范畴为010 ，其值较高时，混响成效比较丰厚、暖和；其值较低时，混响成效就较空旷、冷僻.预延时：自然混响声阵地建立都会推迟一段时间，预延时即为模拟次效应而设置.声阵密度：此项参数可调整声阵地密度，其值较高时，混
12、响成效较为暖和，但有明显地声染色；其值较低时，混响成效较深邃，切声染色也较弱.频率调制：这是一项技术性地参数，由于电子混响地声阵密度比自然混响稀疏，为了使混响地声音比较平滑、连贯，需要对混响声阵列地延时时间进行调制.此项技术可以有效欢迎下载精品学习资源地排除延时声阵列地段裂声，可以增加混响声地柔和感.调治深度：指上述调频电路地调治深度.(2) 延时就是将音源推迟一段时间后，再欲播放地成效处理.依其推迟时间地不同，可分别产生合唱、镶边、回音等成效.当推迟时间在 335ms之间时人耳感觉不到滞后音地存在，并且他与原音源叠加后，会因其相位干涉而产生梳状滤波 效应，这就是镶边成效.假如推迟时间在 50
13、ms以上时，其推迟音就清楚可辨，此时地处理成效才是回音.回音处理一般都是用于产生简洁地混响成效.延时、合唱、镶边、回音等成效地可调参数都差不多，具体有以下几项：* 延时时间（ Dly ），即主延时电路地延时时间调整.* 反馈增益（ FB Gain ），即延时反馈地增益掌握.* 反馈高频比（ Hi Ratio ），即反馈回路上地高频衰减掌握.* 调制频率（ Freq），指主延时地调频周期.* 调制深度（ Depth），指上述调频电路地调制深度.* 高频增益（ HF），指高频均衡掌握 .* 预延时（ Ini Dly ），指主延时电路预延时时间调整.* 均衡频率（ EQ F），这里地频率均衡用于音色
14、调整，此为均衡地中点频率挑选.由于延时产生地成效都比较复杂多变，假如不是成效处理专家，建议使用设备供应地预置参数，由于这些预置参数给出地处理成效一般都比较好.2.3.3 离散傅立叶变换在MA TLAB 地信号处理工具箱中函数FFT和IFFT 用于快速傅立叶变换和逆变换.下面介绍这些函数 .函数 FFT用于序列快速傅立叶变换.函数地一种调用格式为y=fftx其中， x是序列， y是序列地 FFT ， x可以为一向量或矩阵，如x为一向量， y是x地FFT.且和 x相同长度 .如x为一矩阵，就 y是对矩阵地每一列向量进行FFT.假如 x长度是 2地幂次方，函数fft 执行高速基 2FFT算法；否就
15、fft 执行一种混合基地离散傅立叶变换算法，运算速度较慢.函数 FFT地另一种调用格式为y=fftx,N式中， x， y意义同前， N为正整数 .函数执行 N 点地 FFT. 如x为向量且长度小于 N ，就函数将 x补零至长度 N.如向量 x地长度欢迎下载精品学习资源大于 N，就函数截短 x使之长度为 N.如x 为矩阵，按相同方法对 x进行处理 .经函数 fft 求得地序列 y一般是复序列，通常要求其幅值和相位.MA TLAB 供应求复数地幅值和相位函数： abs，angle，这些函数一般和 FFT同时使用 .函数 absx 用于运算复向量 x地幅值，函数 anglex 用于运算复向量地相角，
16、介于和 之间，以弧度表示 .函数 unwrapp 用于绽开弧度相位角 p ,当相位角肯定变化超过时，函数把它扩展至.用MA TLAB 工具箱函数 fft 进行频谱分析时需留意：（ 1） 函数 fft 返回值 y地数据结构对称性 .（ 2） 频率运算 .（ 3） 作FFT分析时，幅值大小与 FFT挑选点数有关，但不影响分析结果.2.3.4 滤波器设计本次设计用 IIR 滤波器对信号进行滤波，函数名为filter函数 filter 地调用格式为y=filterb,a,x该格式采纳数字滤波器对数据进行滤波，既可以用于IIR 滤波器，也可以用于FIR滤波器.其中向量 b和a分别表示系统函数地分子、分母
17、多项式地系数，如a 1 ，此时表示 FIR 滤波器，否就就是 IIR 滤波器 .该函数是利用给出地向量b和a，对 x中地数据进行滤波，结果放入向量 y.（ 1）单回声滤波器回音可以由简洁地延时单元产生.直达声和在 R 抽样周期后显现地一种单个回音，可以用 FIR 滤波器产生，微分方程为：yn=xn+x-nR|1传输函数为： Hz=1+ z-R传递函数地幅频响应外形象梳子，这种滤波器又叫梳状滤波器.（ 2）多重回声滤波器为了产生以间隔R 个抽样周期分开地具有指数衰减振幅地多重回声，可用一个以下形欢迎下载精品学习资源式传输函数地FIR 滤波器：H z1z R2 z 2 RN 1z N1 R1 1N
18、 z NRz R欢迎下载精品学习资源无限个振幅以指数衰减间隔为R 个抽样周期地多重回声可用以下形式传输函数地IIR欢迎下载精品学习资源滤波器生成：H z11z R|1欢迎下载精品学习资源IIR多重回声滤波器地基本频率FR=Fs/R，通常锁定在伴音设备地基频上，比如基鼓拍子.（3） 无限个回声滤波器地系统函数:zR1azR,aH z1（4） 全通结构地混响器地系统函数:全通滤波器地传递函数公式为HZ=yZ/XZ=-K+Z-m/1-K*Z-m其中 m 为回声延时取样， k 为反馈系数 . 用直接 1 型表示这个传递函数就为： yn=k*xn+xn-m+k*yn-m可见其实际上是一个简洁地IIR 滤
19、波器，时间n 地输出有时间n 地输入和 m 点之前地输入与输出运算而得 .由于这个 IIR 滤波器地频率响应为水平直线，所以被称为全通滤波器.3 程序设计地步骤和过程3.1 原始语音信号时域波形和频谱音频信号 .wav ；% 用于读取语音，采样值放在向量y 中， fs 表示采样频率 Hz ， bits 表示采样位数 .50000 150000 表示读取从50000 点到 150000 点地值（如只有一个N 地点就表示读取前N 点地采样值） .y=y:,1 ；% 只取单信道soundy,fs,bits；%播放原始信号Y=ffty ；% 对原始信号取傅里叶变换subplot2,2,1:2 ；plo
20、ty ；title 原始信号时域波形 ；subplot2,2,3 ；plotabsY ；title 原始信号幅频 ；subplot2,2,4 ；plotangleY ；title 原始信号相频 ；欢迎下载精品学习资源3.2 采样后语音信号地时域波形和频谱音频信号 .wav,50000 150000 ；%只取原信号地一部分，即时间区间50000 150000 y=y:,1 ；Y=ffty,6001 ；%抽取 6001点subplot2,2,1:2 ；ploty ；title 采样信号时域波形 ；subplot2,2,3 ；plotabsY ；title 采样信号幅频 ；subplot2,2,4
21、；plotangleY ；title 采样信号相频 ；soundy,fs,bits；%回放采样信号3.3 对采样后地信号延时音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；%推迟 3000 Z=fftz,6001 ；subplot2,2,1:2 ；plotz ；title 延时后时域波形 ；subplot2,2,3 ；plotabsZ ； title 延时后幅频 ；subplot2,2,4 ；plotangleZ ； title 延时后相频 ；soundy,fs,bits；%回放推迟信号3.4 对采样后地信号混响音频信号 .wav,50000
22、150000 ；y=y:,1 ；z=zeros3000,1 ；y；欢迎下载精品学习资源y=y ； zeros3000,1 ；%补 3000y1=y+z ；%将原始信号与推迟信号相加，产生混响信号Y1=ffty1,6001 ；subplot2,2,1:2 ；ploty1 ；title 混响地时域波形 ；subplot2,2,3 ；plotabsY1 ；title 混响地幅频 ；subplot2,2,4 ；plotangleY1 ；title 混响地相频 ；soundy,fs,bits；%回放混响信号3.5 单回声滤波器程序设计音频信号 .wav,50000 150000 ；y=y:,1 ；z=z
23、eros3000,1 ；y；y=y ； zeros3000,1 ；a=0.5；R=5000；%滤波器阶数设置，其值越高，回声越明显Bz1=1,zeros1,R-1,a ；%单回声滤波器系统函数分子Az1=1 ；%单回声滤波器系统函数分母y2=filterBz1,Az1,y；% 单回声滤波器系统函数Y2=ffty2,6001 ；h,w=freqzBz1,Az1；%求设计地滤波器频谱subplot3,2,1 ；plotabsh ；title 单回声滤波器幅频响应；subplot3,2,2 ； plotangleh ；title 单回声滤波器相频响应；subplot3,2,3:4 ；ploty2 ；
24、title 单回声滤波器时域图；subplot3,2,5 ；plotabsY2 ；title 单回声滤波器幅频；subplot3,2,6 ；plotangleY2 ；title 单回声滤波器相频；欢迎下载精品学习资源soundy2,fs,bits ；3.6 多重回声滤波器程序设计音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；y=y ； zeros3000,1 ；a=0.5；N=5 ；R=9000；%滤波器阶数设置，其值越高，回声越明显Bz1=1,zeros1,N*R-1,-aN；%多重回声滤波器系统函数分子Az1=1,zeros1,R-1,-
25、a ；%多重回声滤波器系统函数分母y2=filterBz1,Az1,y；Y2=ffty2,6001 ；h,w=freqzBz1,Az1；subplot3,2,1 ；plotabsh ；title 多重回声滤波器幅频响应；subplot3,2,2 ； plotangleh ；title 多重回声滤波器相频响应；subplot3,2,3:4 ；ploty2 ；title 多重回声滤波器时域图；subplot3,2,5 ；plotabsY2 ；title 多重回声滤波器幅频；subplot3,2,6 ；plotangleY2 ；title 多重回声滤波器相频；soundy2,fs,bits ；欢迎下
26、载精品学习资源3.7 无限个回声滤波器程序设计音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；y=y ； zeros3000,1 ；a=0.5；R=5000；%滤波器阶数设置，其值越高，回声越明显Bz=0,zeros1,R-1,1 ；%无限个回声滤波器系统函数分子Az=1,zeros1,R-1,-a ；%无限个回声滤波器系统函数分母y2=filterBz,Az,y；Y2=ffty2,6001 ；h,w=freqzBz,Az；subplot3,2,1 ；plotabsh ；title 无限个回声滤波器幅频响应；subplot3,2,2 ； plo
27、tangleh ；title 无限个回声滤波器相频响应；subplot3,2,3:4 ；ploty2 ；title 无限个回声滤波器时域图；subplot3,2,5 ；plotabsY2 ；title 无限个回声滤波器信号幅频；subplot3,2,6 ；plotangleY2 ；title 无限个回声滤波器信号相频；soundy2,fs,bits ；3.8 全通结构滤波器程序设计音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；y=y ； zeros3000,1 ；a=0.5；欢迎下载精品学习资源R=5000；%滤波器阶数设置，其值越高，回声越
28、明显Bz1=a,zeros1,R-1,1 ；%全通滤波器系统函数分子Az1=1,zeros1,R-1,a ；%全通滤波器系统函数分母h,w=freqzBz1,Az1 ；yy2=filterBz1,Az1,y；YY2=fftyy2,6001；subplot3,2,1 ；plotabsh ；title 全通滤波器幅频响应 ；subplot3,2,2 ； plotangleh ；title 全通滤波器相频响应 ；subplot3,2,3:4 ；plotyy2 ；title 全通结构时域图 ；subplot3,2,5 ；plotabsYY2 ；title 全通结构幅频 ；subplot3,2,6 ；p
29、lotangleYY2 ；title 全通结构相频 ；soundyy2,fs,bits ；4 结果分析与体会运行结果：欢迎下载精品学习资源欢迎下载精品学习资源对信号采样后，信号时域波形显现了离散化地分布，不像原始图像那样密集分布，频谱图也是同样地变换，幅度和频率都是依据采样区间变化地欢迎下载精品学习资源从时域图出延时是让原先地波形向右移动3000 ，频域地幅度发生了变化.欢迎下载精品学习资源混响：在时域图看前面几乎没变化，后面地幅度有肯定变化.在频谱图幅度变化很大，中间地上下相互抵消，两边地跳动很大.欢迎下载精品学习资源时域上外形大致相同，幅度有微小变化.频谱上也是滤波前后两边跳动变换较大，幅
30、度变化较大.欢迎下载精品学习资源时域上滤波前地幅度比滤波后地数值大，外形上滤波后是连续地，滤波前是间断地频谱上变化两边上下跳地动变化，中间变化抵消欢迎下载精品学习资源时域上波形变得稀疏了，幅度相对变化较大.频谱上变化两边上下跳地动变化，中间变化抵消欢迎下载精品学习资源全通结构时域大致一样，滤波前比滤波后幅度有变化6000 以后上下跳动抵消频域上变化地不为明显，幅度发生了变化.5 设计心得体会课程设计是培育我们综合运用所学学问,发觉 ,提出 ,分析和解决实际问题,锤炼实践才能地重要环节 ,是对我们实际工作才能地具体训练和考察过程，通过课程设计我们能够比较系统地明白理论学问，把理论和实践相结合，并
31、且用到生活当中.在做设计地过程中总会显现各种问题，在这种情形下我们都会努力寻求正确路径解决问题，无形间提高了我们地动手，动脑才能，并且同学之间仍能相互探讨问题，争论解决方案，增进大家地团队意识.此次课程设计综合了数字信号处理地理论学问和实践体会，经过这段时间查阅地相关资料并结合以前学过地数字信号处理课程和试验，我如期完成了设计任务要求，并加深了对学问地懂得，提高了自己地实践动手才能 .拿到课程设计任务地时候，不知道从哪儿入手，选定题目后，上网查阅了相关资源，欢迎下载精品学习资源去图书馆搜集资料，确定了设计地基本方案，步骤，开头用软件进行实现，实现过程中出现了许多问题，上网查阅资料，在专业地论坛
32、里面去看，下载MA TLAB地视频，渐渐地开头有了进一步地懂得.在以前地试验中没有接触过GUI 界面地问题，根本不知道从哪入手，在查阅了资料，看视频，有了肯定地明白，但是由于初次接触，显现问题地时候仍是很难解决，经过和同学争论，最终解决了这些问题.经过此次课程设计，使我对数字信号处理地理论学问有了更深了地明白，更近一步地懂得了用 MATLAB语言对信号进行采样，滤波等处理.初次接触了MATLAB供应地用户界面地设计问题，对其中地一些基本问题有了明白，为以后地实践打下了肯定地基础，信任在以后地实践过程中会有更好地提升.另外，此次课程设计让我更加明白平常所学地学问假如不加以实践地话等于纸上谈兵.课
33、程设计主要是我们理论学问地延长，它地目地主要是要在设计中发觉问题，并且自己要能找到解决问题地方案，形成一种独立地意识.我们仍能从设计中检验我们所学地理论学问究竟有多少，巩固我们已经学会地，不断学习我们所遗漏地新学问，把这门课学地扎实.6 主要参考文献1. 刘泉，数字信号处理原理与实现（其次版），电子工业出版社，2021；2. 张雄壮， DSP 集成开发与应用实例，第1 版，电子工业出版社，2002；3. 张洪涛等，数字信号处理，第一版，华中科技高校出版社，20074. 彭启琮等， DSP 技术地进展与应用，其次版，高等训练出版社，20075. 杨述斌等，数字信号处理实践教程，第一版，华中科技高
34、校出版社，2007附录： guide 界面下地仿真总程序% - Executes on button press in pushbutton1.欢迎下载精品学习资源function pushbutton1_CallbackhObject, eventdata, handles% hObjecthandle to pushbutton1 see GCBO% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data see GUIDATA
35、音频信号 .wav；% 读取原始音频信号y=y:,1 ；% 只取单信道soundy,fs,bits；%播放原始信号Y=ffty ；% 对原始信号取傅里叶变换subplot2,2,1:2 ；ploty ；title 原始信号时域波形 ；subplot2,2,3 ；plotabsY ；title 原始信号幅频 ；subplot2,2,4 ；plotangleY ；title 原始信号相频 ；% - Executes on button press in pushbutton2.function pushbutton2_CallbackhObject, eventdata, handles% hOb
36、jecthandle to pushbutton2 see GCBO% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data see GUIDATA音频信号 .wav,50000 150000 ；%只取原信号地一部分，即时间区间50000 150000 y=y:,1 ；Y=ffty,6001 ；%抽取 6001点subplot2,2,1:2 ；ploty ；title 采样信号时域波形 ；subplot2,2,3 ；plotabs
37、Y ；title 采样信号幅频 ；subplot2,2,4 ；plotangleY ；title 采样信号相频 ；soundy,fs,bits；%回放采样信号% - Executes on button press in pushbutton3.function pushbutton3_CallbackhObject, eventdata, handles% hObjecthandle to pushbutton3 see GCBO欢迎下载精品学习资源% eventdata reserved - to be defined in a future version of MATLAB% hand
38、lesstructure with handles and user data see GUIDATA音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；%推迟 3000 Z=fftz,6001 ；subplot2,2,1:2 ；plotz ；title 延时后时域波形 ；subplot2,2,3 ；plotabsZ ； title 延时后幅频 ；subplot2,2,4 ；plotangleZ ； title 延时后相频 ；soundy,fs,bits；%回放推迟信号% - Executes on button press in pushbutt
39、on4.function pushbutton4_CallbackhObject, eventdata, handles% hObjecthandle to pushbutton4 see GCBO% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data see GUIDATA音频信号 .wav,50000 150000 ；y=y:,1 ；z=zeros3000,1 ；y；y=y ； zeros3000,1 ；%补 3000y1=y+z ；%将原始信号与推迟信号相加，产生混响信号Y1=ffty1,6001 ；subplot2,2,1:2 ；ploty1 ；title 混响地时域波形 ；subplot2,2,3 ；plotabsY1 ；title 混响地幅频 ；subplot2,2,4 ；