学 号：
课 程 设 计
题 目
信息处理课群综合训练设计
学 院
信息工程学院
专 业
通信工程
班 级
通信
姓 名
指导教师
2015
年
月
日
课程设计务书
学生姓名： 专业班级： 通信
指导教师： 工作单位： 信息工程学院
题 目声音延时混响仿真
初始条件：Matlab应软件 Windows带录音功
延时混响相关知识
求完成务 （包括课程设计工作量技术求说明书撰写等具体求）
(1)利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
(2)语音信号频谱分析画出采样语音信号时域波形频谱图
(3)信号加入延时混响分析频谱原始信号频谱进行较
(4)设计种特殊类型滤波器：单回声滤波器重回声滤波器限回声滤波器全通结构混响器画出滤波器频域响应
(5)设计滤波器采集语音信号进行滤波
(6)分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
(7)回放语音信号
时间安排：
1理讲解老师布置课程设计题目学生根选题开始查找资料
2课程设计时间2周
（1）理解相关技术原理确定技术方案 时间2天
（2）选择仿真工具进行仿真设计分析时间6天
（3）总结结果完成课程设计报告时间2天
指导教师签名： 2015 年 6 月 10 日
系（责教师）签名： 2015年 6月 10 日
目录
摘 4
Abstract 5
1 绪 6
11 简介 6
12 设计务求 6
13研究背景意义 7
2 课题描述 8
21课题背景意义 8
22课题目标 8
3系统框架实现 9
31总体方案图 9
32程序流程图 10
4设计原理 11
41混响延时基原理 11
411混响 11
412延时 12
42离散傅里叶变换 13
43滤波器设计 13
431特殊滤波器系统函数 14
432滤波函数filter 14
433滤波器设计 14
5设计步骤 16
51windows系统中读取1s语音信号 16
52采样语音信号时域波形频谱 16
53采样信号延时混响 16
54信号进行混响 17
55单回声滤波器 17
56限回声滤波器 17
57重回声滤波器 18
58全通结构混响器 18
6运行结果 20
7结体会 26
参考文献 27
附录 28
摘
数字信号处理信号数字方式表示处理理技术MATLAB款功强集数学运算仿真体软件次课设利MATLAB软件声音信号进行数字处理程然次运软件强计算功尤软件带较复杂函数调方说：快速傅里叶变换时利软件仿真功设计种滤波器读入声音信号进行滤波处理总体讲通课设学数字信号处理相关理知识熟悉软件滤波器设计解语音信号特点
关键词：数字信号处理快速傅里叶变换延时混响滤波器
Abstract
Digital signal processing is will signal to digital way said the theory and technique of the and processing and MATLAB is a powerful mathematical operation and simulation of the set for one of the software This class set is the use of the software MATLAB voice signal digital processing process Of course the main using software powerful computing functions especially software bringing some complex function call for example fast Fourier transformation At the same time also using the software simulation function design of filter to read in a voice signal filtering processing In general through the class set study the digital signal processing related theory knowledge be familiar with the use of the software and the design of filter understand the characteristics of the speech signal
Keywords digital signal processing fast Fourier transform delay and reverb filter
1 绪
11 简介
数字信号处理研究数字符号序列表示处理信号处理目削弱信号中余容滤混杂噪声干扰者信号变换容易分析识形式便估计选择特征参数
语音信号分析语音信号处理前提基础分析出表示语音信号质特征参数利参数进行高效语音通信语音合成语音识等处理语音合成音质坏语音识率高低取决语音信号分桥准确性精确性语音信号分析语音信号处理应中具举足轻重位根分析出参数性质语音信号分析分时域分析频域分析倒领域分析等时域分析方法具简单计算量物理意义明确等优点语音信号重感知特性反映功率谱中相位变化起着作相时域分析说频域分析更重文简介绍时域分析频域分析
声音信号维连续信号计算机处理离散信号离散信号原连续信号根采样定理确定采样频率值wav文件种数字声音文件格式课程设计基Matlab分析wav声音文件频谱声音关系通采集段声音进行频谱分析等处理然设计数字滤波器处理原始声音wav文件较滤波输出声音信号原声音信号异
12 设计务求
设计务
(1)利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
(2)语音信号频谱分析画出采样语音信号时域波形频谱图
(3)信号加入延时混响分析频谱原始信号频谱进行较
(4)设计种特殊类型滤波器：单回声滤波器重回声滤波器限回声滤波器全通结构混响器画出滤波器频域响应
(5)设计滤波器采集语音信号进行滤波
(6)分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
(7)回放语音信号
设计求
(1)熟悉离散信号系统时域特性
(2)熟悉语音信号特点
(3)掌握数字信号处理基概念基理基方法
(4)掌握序列快速傅里叶变换方法
(5)学会MATLAB掌握MATLAB程序设计方法
(6)掌握MATLAB设计种数字滤波器方法信号进行滤波方法
13研究背景意义
音频信号类社会中发挥着举足轻重作古类始终声音信息传递信息感情现着计算机通信微电子网络技术迅猛发展已改变传统利声音传递声音方式摇旗呐喊线通信类利声音方式始终时代高科技水联系起声音类重交流工具
混响生活中种常见声学现象会议室里开会音乐厅里欣赏音乐会明显感觉混响声音存类型声音声学环境求某固定声学设计法满足现代功厅担负佳混响效果求传统采房间布置吸音系数反射板升降吊顶改变厅堂容积调节厅混响时间方法显然仅麻烦难调佳效果着电子技术发展数字信号处理器性断提高述原通建筑声学解决声学效果问题现利电声学知识数字信号处理手段较方便解决说电子技术发展数字信号处理技术完善声学带次前未重变革
滤波信号处理中直占十分重位现需信号处理领域数字信号处理技术精确快速复杂数进行处理数字滤波数字信号处理技术基方法信号滤检测预测等广泛滤波器数字滤波器分成两种形态：限长脉响应滤波器（FIR）限长脉响应滤波器（IIR）种滤波器设计需量计算手工计算已现实必须助计算机完成中MATLAB款十分优秀计算仿真软件MATLAB种面科学工程计算语言集数值分析矩阵运算信号处理图形显示体具编程效率高调试手段丰富扩充力强等特点MATLAB信号处理工具箱具强函数功仅设计数字滤波器设计达优化数字滤波器设计强力工具
2 课题描述
21课题背景意义
语音信号处理门较实电子工程专业课程语音类获取信息重源利信息重手段通语言相互传递信息类重基功语言类特功创造记载千年类文明史根手段没语言没天类文明语音语言声学表现相互传递信息重手段类重效常方便交换信息形式
年着计算机规模数字集成电路迅速发展语音数字信号处理相应发展语音信号分析模拟语音合成语音识等研究已较成熟种声码器声控器语声识系统语声合成器等已逐渐商品出现
22课题目标
录制wav语音文件进行采样信号加入延时混响信号通种特殊滤波器信号分画出应时域波形频谱图进行分析
3系统框架实现
31总体方案图
利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
语音信号频谱分析画出采样语音信号时域波形频谱图
信号加入延时混响分析频谱原始信号频谱进行较
设计种特殊类型滤波器：限回声滤波器全通结构混响器画出滤波器频域响应设计滤波器采集语音信号进行滤波
分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
回放语音信号
图31 程序总体框图
32程序流程图
wavread电脑导入Windows1s语音信号
画出采样语音信号时域波形频谱图N2001调fft函数快速傅里叶变换
采样语音信号进行延时30调整采样信号长度两者相加混响
单回声滤波器
a05
yy1x+z*05
YY1fft(yy12001)
限回声滤波器
Bz[00000000001]
Az[1000000000a]
yy2filter(BzAzx)
YY2fft(yy22001)
重回声滤波器
N5
Bz1[1000000000a^N] Az1[1000000000a]
yy3filter(Bz1Az1x)
YY3fft(yy32001)
全通结构混响器
Bz1[a0000000001]
Az1[1000000000a]
yy4filter(Bz1Az1x)
YY4fft(yy420001)
sound(yy1fsbits)
sound(yy2fsbits)
sound(yy3fsbits)
sound(yy3fsbit)
图32 程序流程框图
4设计原理
41混响延时基原理
411混响
声音通媒质传入听觉器官媒质声音传播途径声音载体日常生活中司空见惯例子说媒质湖中传递波浪湖水静湖面投块石子石子击起水波波纹越越水波越传越远湖水传播水波媒质声音仿佛某声源发音空气中声波震荡会改变周围空气静止状态形成时压缩时稀疏前进波声波样断散发开传播声音媒质空气
闭室中声源发出声音信号(例脉声信号)时位室点听者接收声音中应包含两部分成份部分声源直接传听者声音称直达声部分声波传墙壁障碍物处反射出声音中听者接收外会传墙壁处次发生反射反射声音中听者接收类似持续听者接收部分声音统称反射声时间反射声较直达声落(称延迟声)量角度次反射会部分量吸收反射声系列量逐渐衰减延迟声
数字混响通工创作回音加入原始信号里隔音录音室里录制声音转换种然声音形式回音简单延迟单元产生
混响效果增加音源融合感然音源延时声阵列非常密集复杂模拟混响效果程序复杂变常见参数种：
（1）混响时间：逼真模拟然混响数码混响器套复杂程序然技术参数调然技术参数调整会原效果更然尤混响时间
（2）高频滚降：项参数模拟然混响中空气高频吸收效应产生较然混响效果般高频混降调范围01~10值较高时混响效果较接然混响值较低时混响效果较清澈
（3）扩散度：项参数调整混响声阵密度增长速度调范围0~10值较高时混响效果较丰厚温暖值较低时混响效果较空旷冷僻
（4）预延时：然混响声阵建立会延迟段时间预延时模拟次效应设置
（5）声阵密度：项参数调整声阵密度值较高时混响效果较温暖明显声染色值较低时混响效果较深邃切声染色较弱
（6）频率调制：项技术性参数电子混响声阵密度然混响稀疏混响声音较滑连贯需混响声阵列延时时间进行调制项技术效消延时声阵列段裂声增加混响声柔感
（7）调治深度：指述调频电路调治深度
412延时
延时音源延迟段时间欲播放效果处理延迟时间分产生合唱镶边回音等效果
延迟时间3~35ms间时耳感觉滞音存原音源叠加会相位干涉产生梳状滤波效应镶边效果果延迟时间50ms时延迟音清晰辨时处理效果回音回音处理般产生简单混响效果
延时合唱镶边回音等效果调参数差具体项：
延时时间（Dly）延时电路延时时间调整
反馈增益（FB Gain）延时反馈增益控制
反馈高频（Hi Ratio）反馈回路高频衰减控制
调制频率（Freq）指延时调频周期
调制深度（Depth）指述调频电路调制深度
高频增益（HF）指高频均衡控制
预延时（Ini Dly）指延时电路预延时时间调整
均衡频率（EQ F）里频率均衡音色调整均衡中点频率选择
延时产生效果较复杂变果效果处理专家建议设备提供预置参数预置参数出处理效果般较
42离散傅里叶变换
MATLAB信号处理工具箱中函数FFTIFFT快速傅立叶变换逆变换面介绍函数
函数FFT序列快速傅立叶变换
函数种调格式 yfft(x) （式41）
中x序列y序列FFTx量矩阵x量yxFFTx相长度x矩阵y矩阵列量进行FFT
果x长度2幂次方函数fft执行高速基－2FFT算法否fft执行种混合基离散傅立叶变换算法计算速度较慢
函数FFT种调格式 yfft(xN) （式42）
式中xy意义前N正整数
函数执行N点FFTx量长度N函数x补零长度N量x长度N函数截短x长度Nx 矩阵相方法x进行处理
函数fft求序列y般复序列通常求幅值相位MATLAB提供求复数幅值相位函数：absangle函数般FFT时
MATLAB工具箱函数fft进行频谱分析时需注意：
（1） 函数fft返回值y数结构称性
（2） 频率计算
（3） 作FFT分析时幅值FFT选择点数关影响分析结果
43滤波器设计
数字滤波器作利离散时间系统特性输入信号波形(频谱)进行加工处理者说利数字方法预定求信号进行变换
数字滤波器理解计算程序算法代表输入信号数字时间序列转化代表输出信号数字时间序列转化程中信号预定形式变化
431特殊滤波器系统函数
1 单回声滤波器系统函数：
H(z) <1 （式43）
2 限回声滤波器系统函数
H(z) <1 （式44）
3 重回声滤波器系统函数 ：
H(z) <1 （式45）
4 全通结构混响器系统函数
H(z) <1 （式46）
432滤波函数filter
系统离散时间LTI系统系统输入信号响应实质输入信号频域角度连续时间LTI频谱进行选择处理程程称滤波MATLAB信息处理工具箱中提供维滤波器函数filter二维滤波器函数filter2
函数filter调格式 yfilter(bax) （ 式47）
该格式采数字滤波器数进行滤波IIR滤波器FIR滤波器中量ba分表示系统函数分子分母项式系数a＝1时表示FIR滤波器否IIR滤波器该函数利出量bax中数进行滤波结果放入量y
433滤波器设计
先设计中种滤波器进行简介
（1）单回声滤波器
回声简单延时单元产生直达声m抽样周期出现种单回声FIR滤波器产生差分方程：
y[n]x[n]+αx[nm]|α|<1
单回声滤波器传输函数：
传输函数幅频响应形状象梳子种滤波器梳状滤波器
（2）重回声滤波器
产生间隔m抽样周期分开具指数衰减振幅重回声形式传输函数FIR滤波器：
限振幅指数衰减间隔m抽样周期重回声形式传输函数IIR滤波器生成：
（3）限回声滤波器
限回声滤波器传输函数：
述重回声滤波器提供然声音混响幅频特性知幅度响应频率常数收听效果令满意次秒回声数量太少会引起合成声颤动需秒约1000回声生成没颤动反射声
（4）全通结构混响器
生成更真实混响提出种全通结构混响器传输函数：
外函数filter语音信号进行滤波调格式 ：
yfilter(bax)
该格式采数字滤波器数进行滤波IIR滤波器FIR滤波器中量ba分表示系统函数分子分母项式系数a＝1时表示FIR滤波器否IIR滤波器该函数利出量bax中数进行滤波结果放入量y
5设计步骤
51windows系统中读取1s语音信号
Matlab语句[xfsbits]wavread(＇makewav＇[500 50000])读取语音采样值放量x中fs表示采样频率(Hz)bits表示采样位数[500 50000]表示读取501点50000点值（N点表示读取前N点采样值）
52采样语音信号时域波形频谱
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
sound(5*xfs) 声音回放
n102000
Nsize(x1)
figure（2）subplot(211)plot(x)title(＇采样语音信号时域波形＇)
Yfft(x2001) subplot(212) 信号做2001点FFT变换
plot(n1(11000)Y(11000)) title(＇采样语音信号频谱图＇) 信号幅度
中[xfsbits]wavread(＇makewav＇[10000 40000])读取语音采样值放量x中fs表示采样频率(Hz)bits表示采样位数
Sound函数声音回放X代表信号fft函数求解傅里叶变换出采样信号
53采样信号延时混响
z[zeros(2001)x] figure(3) 信号延时
subplot(211)plot(z) title(＇延时时域图＇) 画出延时信号时域图
Zfft(z2001)
subplot(212)plot(n1(11000)Z(11000)) 延时信号频谱图
title(＇延时频谱图 ＇)
中zeros函数信号延时zeros生成全0矩阵zeros(MN) 表示M行*N列0矩阵
54信号进行混响
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
y1x+z 信号混响
figure(4)subplot(211)plot(y1) title(＇混响时域图＇) 混响时域图
Y1fft(y12001) 混响信号2001点FFT变换
subplot(212)plot(n1(11000)Y1(11000)) 混响频谱图
title(＇混响频谱图＇)
语音信号进行延时原语音信号长度变换延时信号长时两信号相加便混响信号
55单回声滤波器
a05 a取等1
yy1x+z*05 信号单回声滤波
figure（5）subplot(211)plot(yy1)title(＇单回声滤波器时域图＇)
YY1fft(yy12001) 单回声信号做2001点FFT变换
subplot(212)plot(n1(11000)Y2(11000)) 单回声信号频谱图
title(＇单回声滤波器频谱图＇)
中函数yy1x+z*05信号进行单回声滤波
56限回声滤波器
a05 a取等1
Bz[00000000001] 分子系数
Az[1000000000a] 分母系数
yy2filter(BzAzx) 滤波器进行滤波
figure(6)subplot(211)plot(yy2) 限回声滤波器时域波形
title(＇限回声滤波器时域波形＇)
YY2fft(yy22001) 限回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 限回声滤波器频谱图
title(＇限回声滤波器频谱图 ＇)
中编写出应限回声滤波器函数然filter函数调滤波器
57重回声滤波器
a05 a取等1
N5
Bz1[100000000005^N] 分子系数
Az1[100000000005] 分母系数
yy3filter(Bz1Az1x) figure(7) 滤波器进行滤波
subplot(211)plot(yy3) title(＇重回声滤波器时域图＇)重回声滤波器混响器时域波形
YY3fft(yy32001) 重回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY3(11000)) 重回声滤波器频谱图
title(＇重回声频谱图＇)
58全通结构混响器
a05 a取等1
Bz1[a0000000001] 分子系数
Az1[1000000000a] 分母系数
yy3filter(Bz1Az1x) figure(6) 滤波器进行滤波
subplot(211)plot(yy2) title(＇全通应时域图＇)全通结构混响器时域波形
YY2fft(yy22001) 全通结构混响器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 全通结构混响器频谱图
title(＇全通应频谱图＇)
理出全通结构滤波器函数filter函数调滤波器fft函数进行傅里叶变换出应时域波形频谱结构
6运行结果
图61 原始信号时域波形
图62 采样时域波形图频谱图
图63 延时时域波形频谱图
分析：1）时域延时波形右移动
2）频谱图延时前曲线延时曲线幅值变
图64混响时域波形图频谱图
分析：1)时域混响时域波形发生定变化24000HZ前幅度定变化24000HZ波形明显变化
2)频谱图混响频谱跳动幅值定改变采样点200点处混响较明显
混响波信原始波形相已明显差异时域图混响波形前前部分没变化波形方幅度增强方幅度减程序中设计混响效果相：原始语音信号延时信号（充回音）叠加引起波形幅度变化
图65单回声滤波器时域波形图频谱图
分析：1) 时域滤波波形定变化特24000HZ波形变化较明显
2)频谱滤波频谱微跳动
图66限回声滤波器时域波形图频谱图
分析：1) 时域波形滤波前波形体相幅度微变化
2)频谱图幅度定变化采样点100点200点间幅值变化明显曲线相滤波前总体变化趋势
图67重回声滤波器时域波形图频谱图
分析：1)时域波形变化20000HZ观察波形形状明显滤波前波形30000HZ便没重回声滤波器时域波形延长80000HZ
2)频谱图幅值总体减
图68全通结构混响时域波形图频谱图
分析：1） 时域波形滤波前波形体相幅度微变化
2）频谱图幅值略微减曲线总体趋势
7结体会
说实话Matlab编程强项时更C编程次设计程中先基原理概念致弄清楚然找例程句句分析发现该matlab程序核心实际数学表达式转变成Matlab语言然Matlab软件便识执行试着种滤波器系统函数表达式写Matlab程序断修改调试终获成功触类旁通接设计工作利
通次信息处理课群综合训练设计数字信号处理更进步熟悉亲手设计程中加深课知识理解次课题基MATLAB语音信号特技处理语音信号进行延时混响滤波等操作通网搜索资料查阅课课外书籍动手设计滤波器采集语音语音分析等工作加强MATLAB程序编写力数字信号处理相关知识理解
课程设计培养学生综合运学知识发现分析解决实际问题锻炼实践力重环节学生实际工作力具体训练考察程通课程设计够较系统解理知识理实践相结合生活中做设计程中总会出现种问题种情况会努力寻求佳路径解决问题形间提高动手动脑力学间相互探讨问题研究解决方案增加家团队意识
通课程设计知道时学知识果加实践话等纸谈兵课程设计理知识延伸目设计中发现问题找解决问题方案形成种独立意识设计中检验学理知识底少巩固已学会断学忽略新知识期提升专业技
参考文献
[1] 邹彦DSP原理应第版电子工业出版社2006
[2] 周开利MATLAB基础应教程北京学出版社
[3] 张雄伟DSP集成开发应实例第版电子工业出版社2002
[4] 刘幺宋庭新语音识控制技术科学出版社
[5] 彭启琮等DSP技术发展应第二版高等教育出版社2007
[6] 杨述斌等数字信号处理实践教程第版华中科技学出版社2007
附录
程序：
1 原始信号
x1wavread(＇makewav＇) 读取信号
figure(1)plot(x1)title(＇原始信号＇)
sound(5*x140000)
2 采样信号
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
sound(5*xfs) 声音回放
n102000
Nsize(x1)
figure(2)subplot(211)plot(x)title(＇采样语音信号时域波形＇)
Yfft(x2001) subplot(212) 信号做2001点FFT变换
plot(n1(11000)Y(11000)) title(＇采样语音信号频谱图＇) 信号幅度
3 延时信号
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
n102000
z[zeros(2001)x] figure(3) 信号延时
subplot(211)plot(z) title(＇延时时域图＇) 画出延时信号时域图
Zfft(z2001)
subplot(212)plot(n1(11000)Z(11000)) 延时信号频谱图
title(＇延时频谱图 ＇)
sound(5*zfs)
4 混响信号
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
y1x+z 信号混响
figure(4)subplot(211)plot(y1) title(＇混响时域图＇) 混响时域图
Y1fft(y12001) 混响信号2001点FFT变换
subplot(212)plot(n1(11000)Y1(11000)) 混响频谱图
title(＇混响频谱图＇)
sound(5*y1fs) 回放混响信号
5 单回声滤波器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
yy1x+z*05 信号单回声滤波
figure(5)subplot(211)plot(yy1)title(＇单回声滤波器时域图＇)
YY1fft(yy12001) 单回声信号做2001点FFT变换
subplot(212)plot(n1(11000)YY1(11000)) 单回声信号频谱图
title(＇单回声滤波器频谱图＇)
sound(5*yy1fsbits) 回放单回声滤波器语音信号
6 限回声滤波器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
Bz[00000000001] 分子系数
Az[1000000000a] 分母系数
yy2filter(BzAzx) 滤波器进行滤波
figure(6)subplot(211)plot(yy2) 限回声滤波器时域波形
title(＇限回声滤波器时域波形＇)
YY2fft(yy22001) 限回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 限回声滤波器频谱图
title(＇限回声滤波器频谱图 ＇)
sound(5*yy2fsbits) 回放限回声滤波器语音信号
pause(2)
7 重回声滤波器
[xfsbits]wavread(＇makewav＇)
xx(1) 取单声道
n102000
a05
N5
Bz1[100000000005^N] 分子系数
Az1[100000000005] 分母系数
yy3filter(Bz1Az1x) figure(7) 滤波器进行滤波
subplot(211)plot(yy3) title(＇重回声滤波器时域图＇)重回声滤波器混响器时域波形
YY3fft(yy32001) 重回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY3(11000)) 重回声滤波器频谱图
title(＇重回声频谱图＇)
sound(5*yy3fsbits) 回放限回声滤波器语音信号
8 全通结构混响器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
Bz1[a0000000001] 分子系数
Az1[1000000000a] 分母系数
yy4filter(Bz1Az1x) figure(8) 滤波器进行滤波
subplot(211)plot(yy4) title(＇全通结构混响器时域图＇)全通结构混响器时域波形
YY4fft(yy42001) 全通结构混响器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY4(11000)) 全通结构混响器频谱图
title(＇全通结构混响器频谱图＇)
sound(5*yy4fsbits) 回放全通结构混响器语音信号
附件6：
科生课程设计成绩评定表
姓 名
性
专业班级
课程设计题目：
课程设计答辩质疑记录：
课程设计评分标准
序 号
容
分值
1
完成相关算法设计
25分
2
通编写调试程序完成基功
35分
3
设计报告完成
20分
4
课程设计答辩
20分
终评定成绩（优良中格格评定）
指导教师签字：
2015年 7 月 10 日
武han理工学
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档
学 号：
课 程 设 计
题 目
信息处理课群综合训练设计
学 院
信息工程学院
专 业
通信工程
班 级
通信
姓 名
指导教师
2015
年
月
日
课程设计务书
学生姓名： 专业班级： 通信
指导教师： 工作单位： 信息工程学院
题 目声音延时混响仿真
初始条件：Matlab应软件 Windows带录音功
延时混响相关知识
求完成务 （包括课程设计工作量技术求说明书撰写等具体求）
(1)利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
(2)语音信号频谱分析画出采样语音信号时域波形频谱图
(3)信号加入延时混响分析频谱原始信号频谱进行较
(4)设计种特殊类型滤波器：单回声滤波器重回声滤波器限回声滤波器全通结构混响器画出滤波器频域响应
(5)设计滤波器采集语音信号进行滤波
(6)分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
(7)回放语音信号
时间安排：
1理讲解老师布置课程设计题目学生根选题开始查找资料
2课程设计时间2周
（1）理解相关技术原理确定技术方案 时间2天
（2）选择仿真工具进行仿真设计分析时间6天
（3）总结结果完成课程设计报告时间2天
指导教师签名： 2015 年 6 月 10 日
系（责教师）签名： 2015年 6月 10 日
目录
摘 4
Abstract 5
1 绪 6
11 简介 6
12 设计务求 6
13研究背景意义 7
2 课题描述 8
21课题背景意义 8
22课题目标 8
3系统框架实现 9
31总体方案图 9
32程序流程图 10
4设计原理 11
41混响延时基原理 11
411混响 11
412延时 12
42离散傅里叶变换 13
43滤波器设计 13
431特殊滤波器系统函数 14
432滤波函数filter 14
433滤波器设计 14
5设计步骤 16
51windows系统中读取1s语音信号 16
52采样语音信号时域波形频谱 16
53采样信号延时混响 16
54信号进行混响 17
55单回声滤波器 17
56限回声滤波器 17
57重回声滤波器 18
58全通结构混响器 18
6运行结果 20
7结体会 26
参考文献 27
附录 28
摘
数字信号处理信号数字方式表示处理理技术MATLAB款功强集数学运算仿真体软件次课设利MATLAB软件声音信号进行数字处理程然次运软件强计算功尤软件带较复杂函数调方说：快速傅里叶变换时利软件仿真功设计种滤波器读入声音信号进行滤波处理总体讲通课设学数字信号处理相关理知识熟悉软件滤波器设计解语音信号特点
关键词：数字信号处理快速傅里叶变换延时混响滤波器
Abstract
Digital signal processing is will signal to digital way said the theory and technique of the and processing and MATLAB is a powerful mathematical operation and simulation of the set for one of the software This class set is the use of the software MATLAB voice signal digital processing process Of course the main using software powerful computing functions especially software bringing some complex function call for example fast Fourier transformation At the same time also using the software simulation function design of filter to read in a voice signal filtering processing In general through the class set study the digital signal processing related theory knowledge be familiar with the use of the software and the design of filter understand the characteristics of the speech signal
Keywords digital signal processing fast Fourier transform delay and reverb filter
1 绪
11 简介
数字信号处理研究数字符号序列表示处理信号处理目削弱信号中余容滤混杂噪声干扰者信号变换容易分析识形式便估计选择特征参数
语音信号分析语音信号处理前提基础分析出表示语音信号质特征参数利参数进行高效语音通信语音合成语音识等处理语音合成音质坏语音识率高低取决语音信号分桥准确性精确性语音信号分析语音信号处理应中具举足轻重位根分析出参数性质语音信号分析分时域分析频域分析倒领域分析等时域分析方法具简单计算量物理意义明确等优点语音信号重感知特性反映功率谱中相位变化起着作相时域分析说频域分析更重文简介绍时域分析频域分析
声音信号维连续信号计算机处理离散信号离散信号原连续信号根采样定理确定采样频率值wav文件种数字声音文件格式课程设计基Matlab分析wav声音文件频谱声音关系通采集段声音进行频谱分析等处理然设计数字滤波器处理原始声音wav文件较滤波输出声音信号原声音信号异
12 设计务求
设计务
(1)利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
(2)语音信号频谱分析画出采样语音信号时域波形频谱图
(3)信号加入延时混响分析频谱原始信号频谱进行较
(4)设计种特殊类型滤波器：单回声滤波器重回声滤波器限回声滤波器全通结构混响器画出滤波器频域响应
(5)设计滤波器采集语音信号进行滤波
(6)分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
(7)回放语音信号
设计求
(1)熟悉离散信号系统时域特性
(2)熟悉语音信号特点
(3)掌握数字信号处理基概念基理基方法
(4)掌握序列快速傅里叶变换方法
(5)学会MATLAB掌握MATLAB程序设计方法
(6)掌握MATLAB设计种数字滤波器方法信号进行滤波方法
13研究背景意义
音频信号类社会中发挥着举足轻重作古类始终声音信息传递信息感情现着计算机通信微电子网络技术迅猛发展已改变传统利声音传递声音方式摇旗呐喊线通信类利声音方式始终时代高科技水联系起声音类重交流工具
混响生活中种常见声学现象会议室里开会音乐厅里欣赏音乐会明显感觉混响声音存类型声音声学环境求某固定声学设计法满足现代功厅担负佳混响效果求传统采房间布置吸音系数反射板升降吊顶改变厅堂容积调节厅混响时间方法显然仅麻烦难调佳效果着电子技术发展数字信号处理器性断提高述原通建筑声学解决声学效果问题现利电声学知识数字信号处理手段较方便解决说电子技术发展数字信号处理技术完善声学带次前未重变革
滤波信号处理中直占十分重位现需信号处理领域数字信号处理技术精确快速复杂数进行处理数字滤波数字信号处理技术基方法信号滤检测预测等广泛滤波器数字滤波器分成两种形态：限长脉响应滤波器（FIR）限长脉响应滤波器（IIR）种滤波器设计需量计算手工计算已现实必须助计算机完成中MATLAB款十分优秀计算仿真软件MATLAB种面科学工程计算语言集数值分析矩阵运算信号处理图形显示体具编程效率高调试手段丰富扩充力强等特点MATLAB信号处理工具箱具强函数功仅设计数字滤波器设计达优化数字滤波器设计强力工具
2 课题描述
21课题背景意义
语音信号处理门较实电子工程专业课程语音类获取信息重源利信息重手段通语言相互传递信息类重基功语言类特功创造记载千年类文明史根手段没语言没天类文明语音语言声学表现相互传递信息重手段类重效常方便交换信息形式
年着计算机规模数字集成电路迅速发展语音数字信号处理相应发展语音信号分析模拟语音合成语音识等研究已较成熟种声码器声控器语声识系统语声合成器等已逐渐商品出现
22课题目标
录制wav语音文件进行采样信号加入延时混响信号通种特殊滤波器信号分画出应时域波形频谱图进行分析
3系统框架实现
31总体方案图
利Windows录音机软件录制段语音信号时间控制1s左右录制信号进行采样
语音信号频谱分析画出采样语音信号时域波形频谱图
信号加入延时混响分析频谱原始信号频谱进行较
设计种特殊类型滤波器：限回声滤波器全通结构混响器画出滤波器频域响应设计滤波器采集语音信号进行滤波
分析信号频谱画出滤波信号时域波形频谱滤波前信号进行分析信号变化
回放语音信号
图31 程序总体框图
32程序流程图
wavread电脑导入Windows1s语音信号
画出采样语音信号时域波形频谱图N2001调fft函数快速傅里叶变换
采样语音信号进行延时30调整采样信号长度两者相加混响
单回声滤波器
a05
yy1x+z*05
YY1fft(yy12001)
限回声滤波器
Bz[00000000001]
Az[1000000000a]
yy2filter(BzAzx)
YY2fft(yy22001)
重回声滤波器
N5
Bz1[1000000000a^N] Az1[1000000000a]
yy3filter(Bz1Az1x)
YY3fft(yy32001)
全通结构混响器
Bz1[a0000000001]
Az1[1000000000a]
yy4filter(Bz1Az1x)
YY4fft(yy420001)
sound(yy1fsbits)
sound(yy2fsbits)
sound(yy3fsbits)
sound(yy3fsbit)
图32 程序流程框图
4设计原理
41混响延时基原理
411混响
声音通媒质传入听觉器官媒质声音传播途径声音载体日常生活中司空见惯例子说媒质湖中传递波浪湖水静湖面投块石子石子击起水波波纹越越水波越传越远湖水传播水波媒质声音仿佛某声源发音空气中声波震荡会改变周围空气静止状态形成时压缩时稀疏前进波声波样断散发开传播声音媒质空气
闭室中声源发出声音信号(例脉声信号)时位室点听者接收声音中应包含两部分成份部分声源直接传听者声音称直达声部分声波传墙壁障碍物处反射出声音中听者接收外会传墙壁处次发生反射反射声音中听者接收类似持续听者接收部分声音统称反射声时间反射声较直达声落(称延迟声)量角度次反射会部分量吸收反射声系列量逐渐衰减延迟声
数字混响通工创作回音加入原始信号里隔音录音室里录制声音转换种然声音形式回音简单延迟单元产生
混响效果增加音源融合感然音源延时声阵列非常密集复杂模拟混响效果程序复杂变常见参数种：
（1）混响时间：逼真模拟然混响数码混响器套复杂程序然技术参数调然技术参数调整会原效果更然尤混响时间
（2）高频滚降：项参数模拟然混响中空气高频吸收效应产生较然混响效果般高频混降调范围01~10值较高时混响效果较接然混响值较低时混响效果较清澈
（3）扩散度：项参数调整混响声阵密度增长速度调范围0~10值较高时混响效果较丰厚温暖值较低时混响效果较空旷冷僻
（4）预延时：然混响声阵建立会延迟段时间预延时模拟次效应设置
（5）声阵密度：项参数调整声阵密度值较高时混响效果较温暖明显声染色值较低时混响效果较深邃切声染色较弱
（6）频率调制：项技术性参数电子混响声阵密度然混响稀疏混响声音较滑连贯需混响声阵列延时时间进行调制项技术效消延时声阵列段裂声增加混响声柔感
（7）调治深度：指述调频电路调治深度
412延时
延时音源延迟段时间欲播放效果处理延迟时间分产生合唱镶边回音等效果
延迟时间3~35ms间时耳感觉滞音存原音源叠加会相位干涉产生梳状滤波效应镶边效果果延迟时间50ms时延迟音清晰辨时处理效果回音回音处理般产生简单混响效果
延时合唱镶边回音等效果调参数差具体项：
延时时间（Dly）延时电路延时时间调整
反馈增益（FB Gain）延时反馈增益控制
反馈高频（Hi Ratio）反馈回路高频衰减控制
调制频率（Freq）指延时调频周期
调制深度（Depth）指述调频电路调制深度
高频增益（HF）指高频均衡控制
预延时（Ini Dly）指延时电路预延时时间调整
均衡频率（EQ F）里频率均衡音色调整均衡中点频率选择
延时产生效果较复杂变果效果处理专家建议设备提供预置参数预置参数出处理效果般较
42离散傅里叶变换
MATLAB信号处理工具箱中函数FFTIFFT快速傅立叶变换逆变换面介绍函数
函数FFT序列快速傅立叶变换
函数种调格式 yfft(x) （式41）
中x序列y序列FFTx量矩阵x量yxFFTx相长度x矩阵y矩阵列量进行FFT
果x长度2幂次方函数fft执行高速基－2FFT算法否fft执行种混合基离散傅立叶变换算法计算速度较慢
函数FFT种调格式 yfft(xN) （式42）
式中xy意义前N正整数
函数执行N点FFTx量长度N函数x补零长度N量x长度N函数截短x长度Nx 矩阵相方法x进行处理
函数fft求序列y般复序列通常求幅值相位MATLAB提供求复数幅值相位函数：absangle函数般FFT时
MATLAB工具箱函数fft进行频谱分析时需注意：
（1） 函数fft返回值y数结构称性
（2） 频率计算
（3） 作FFT分析时幅值FFT选择点数关影响分析结果
43滤波器设计
数字滤波器作利离散时间系统特性输入信号波形(频谱)进行加工处理者说利数字方法预定求信号进行变换
数字滤波器理解计算程序算法代表输入信号数字时间序列转化代表输出信号数字时间序列转化程中信号预定形式变化
431特殊滤波器系统函数
1 单回声滤波器系统函数：
H(z) <1 （式43）
2 限回声滤波器系统函数
H(z) <1 （式44）
3 重回声滤波器系统函数 ：
H(z) <1 （式45）
4 全通结构混响器系统函数
H(z) <1 （式46）
432滤波函数filter
系统离散时间LTI系统系统输入信号响应实质输入信号频域角度连续时间LTI频谱进行选择处理程程称滤波MATLAB信息处理工具箱中提供维滤波器函数filter二维滤波器函数filter2
函数filter调格式 yfilter(bax) （ 式47）
该格式采数字滤波器数进行滤波IIR滤波器FIR滤波器中量ba分表示系统函数分子分母项式系数a＝1时表示FIR滤波器否IIR滤波器该函数利出量bax中数进行滤波结果放入量y
433滤波器设计
先设计中种滤波器进行简介
（1）单回声滤波器
回声简单延时单元产生直达声m抽样周期出现种单回声FIR滤波器产生差分方程：
y[n]x[n]+αx[nm]|α|<1
单回声滤波器传输函数：
传输函数幅频响应形状象梳子种滤波器梳状滤波器
（2）重回声滤波器
产生间隔m抽样周期分开具指数衰减振幅重回声形式传输函数FIR滤波器：
限振幅指数衰减间隔m抽样周期重回声形式传输函数IIR滤波器生成：
（3）限回声滤波器
限回声滤波器传输函数：
述重回声滤波器提供然声音混响幅频特性知幅度响应频率常数收听效果令满意次秒回声数量太少会引起合成声颤动需秒约1000回声生成没颤动反射声
（4）全通结构混响器
生成更真实混响提出种全通结构混响器传输函数：
外函数filter语音信号进行滤波调格式 ：
yfilter(bax)
该格式采数字滤波器数进行滤波IIR滤波器FIR滤波器中量ba分表示系统函数分子分母项式系数a＝1时表示FIR滤波器否IIR滤波器该函数利出量bax中数进行滤波结果放入量y
5设计步骤
51windows系统中读取1s语音信号
Matlab语句[xfsbits]wavread(＇makewav＇[500 50000])读取语音采样值放量x中fs表示采样频率(Hz)bits表示采样位数[500 50000]表示读取501点50000点值（N点表示读取前N点采样值）
52采样语音信号时域波形频谱
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
sound(5*xfs) 声音回放
n102000
Nsize(x1)
figure（2）subplot(211)plot(x)title(＇采样语音信号时域波形＇)
Yfft(x2001) subplot(212) 信号做2001点FFT变换
plot(n1(11000)Y(11000)) title(＇采样语音信号频谱图＇) 信号幅度
中[xfsbits]wavread(＇makewav＇[10000 40000])读取语音采样值放量x中fs表示采样频率(Hz)bits表示采样位数
Sound函数声音回放X代表信号fft函数求解傅里叶变换出采样信号
53采样信号延时混响
z[zeros(2001)x] figure(3) 信号延时
subplot(211)plot(z) title(＇延时时域图＇) 画出延时信号时域图
Zfft(z2001)
subplot(212)plot(n1(11000)Z(11000)) 延时信号频谱图
title(＇延时频谱图 ＇)
中zeros函数信号延时zeros生成全0矩阵zeros(MN) 表示M行*N列0矩阵
54信号进行混响
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
y1x+z 信号混响
figure(4)subplot(211)plot(y1) title(＇混响时域图＇) 混响时域图
Y1fft(y12001) 混响信号2001点FFT变换
subplot(212)plot(n1(11000)Y1(11000)) 混响频谱图
title(＇混响频谱图＇)
语音信号进行延时原语音信号长度变换延时信号长时两信号相加便混响信号
55单回声滤波器
a05 a取等1
yy1x+z*05 信号单回声滤波
figure（5）subplot(211)plot(yy1)title(＇单回声滤波器时域图＇)
YY1fft(yy12001) 单回声信号做2001点FFT变换
subplot(212)plot(n1(11000)Y2(11000)) 单回声信号频谱图
title(＇单回声滤波器频谱图＇)
中函数yy1x+z*05信号进行单回声滤波
56限回声滤波器
a05 a取等1
Bz[00000000001] 分子系数
Az[1000000000a] 分母系数
yy2filter(BzAzx) 滤波器进行滤波
figure(6)subplot(211)plot(yy2) 限回声滤波器时域波形
title(＇限回声滤波器时域波形＇)
YY2fft(yy22001) 限回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 限回声滤波器频谱图
title(＇限回声滤波器频谱图 ＇)
中编写出应限回声滤波器函数然filter函数调滤波器
57重回声滤波器
a05 a取等1
N5
Bz1[100000000005^N] 分子系数
Az1[100000000005] 分母系数
yy3filter(Bz1Az1x) figure(7) 滤波器进行滤波
subplot(211)plot(yy3) title(＇重回声滤波器时域图＇)重回声滤波器混响器时域波形
YY3fft(yy32001) 重回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY3(11000)) 重回声滤波器频谱图
title(＇重回声频谱图＇)
58全通结构混响器
a05 a取等1
Bz1[a0000000001] 分子系数
Az1[1000000000a] 分母系数
yy3filter(Bz1Az1x) figure(6) 滤波器进行滤波
subplot(211)plot(yy2) title(＇全通应时域图＇)全通结构混响器时域波形
YY2fft(yy22001) 全通结构混响器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 全通结构混响器频谱图
title(＇全通应频谱图＇)
理出全通结构滤波器函数filter函数调滤波器fft函数进行傅里叶变换出应时域波形频谱结构
6运行结果
图61 原始信号时域波形
图62 采样时域波形图频谱图
图63 延时时域波形频谱图
分析：1）时域延时波形右移动
2）频谱图延时前曲线延时曲线幅值变
图64混响时域波形图频谱图
分析：1)时域混响时域波形发生定变化24000HZ前幅度定变化24000HZ波形明显变化
2)频谱图混响频谱跳动幅值定改变采样点200点处混响较明显
混响波信原始波形相已明显差异时域图混响波形前前部分没变化波形方幅度增强方幅度减程序中设计混响效果相：原始语音信号延时信号（充回音）叠加引起波形幅度变化
图65单回声滤波器时域波形图频谱图
分析：1) 时域滤波波形定变化特24000HZ波形变化较明显
2)频谱滤波频谱微跳动
图66限回声滤波器时域波形图频谱图
分析：1) 时域波形滤波前波形体相幅度微变化
2)频谱图幅度定变化采样点100点200点间幅值变化明显曲线相滤波前总体变化趋势
图67重回声滤波器时域波形图频谱图
分析：1)时域波形变化20000HZ观察波形形状明显滤波前波形30000HZ便没重回声滤波器时域波形延长80000HZ
2)频谱图幅值总体减
图68全通结构混响时域波形图频谱图
分析：1） 时域波形滤波前波形体相幅度微变化
2）频谱图幅值略微减曲线总体趋势
7结体会
说实话Matlab编程强项时更C编程次设计程中先基原理概念致弄清楚然找例程句句分析发现该matlab程序核心实际数学表达式转变成Matlab语言然Matlab软件便识执行试着种滤波器系统函数表达式写Matlab程序断修改调试终获成功触类旁通接设计工作利
通次信息处理课群综合训练设计数字信号处理更进步熟悉亲手设计程中加深课知识理解次课题基MATLAB语音信号特技处理语音信号进行延时混响滤波等操作通网搜索资料查阅课课外书籍动手设计滤波器采集语音语音分析等工作加强MATLAB程序编写力数字信号处理相关知识理解
课程设计培养学生综合运学知识发现分析解决实际问题锻炼实践力重环节学生实际工作力具体训练考察程通课程设计够较系统解理知识理实践相结合生活中做设计程中总会出现种问题种情况会努力寻求佳路径解决问题形间提高动手动脑力学间相互探讨问题研究解决方案增加家团队意识
通课程设计知道时学知识果加实践话等纸谈兵课程设计理知识延伸目设计中发现问题找解决问题方案形成种独立意识设计中检验学理知识底少巩固已学会断学忽略新知识期提升专业技
参考文献
[1] 邹彦DSP原理应第版电子工业出版社2006
[2] 周开利MATLAB基础应教程北京学出版社
[3] 张雄伟DSP集成开发应实例第版电子工业出版社2002
[4] 刘幺宋庭新语音识控制技术科学出版社
[5] 彭启琮等DSP技术发展应第二版高等教育出版社2007
[6] 杨述斌等数字信号处理实践教程第版华中科技学出版社2007
附录
程序：
1 原始信号
x1wavread(＇makewav＇) 读取信号
figure(1)plot(x1)title(＇原始信号＇)
sound(5*x140000)
2 采样信号
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
sound(5*xfs) 声音回放
n102000
Nsize(x1)
figure(2)subplot(211)plot(x)title(＇采样语音信号时域波形＇)
Yfft(x2001) subplot(212) 信号做2001点FFT变换
plot(n1(11000)Y(11000)) title(＇采样语音信号频谱图＇) 信号幅度
3 延时信号
[xfsbits]wavread(＇makewav＇[10000 40000])
xx(1) 取单声道
n102000
z[zeros(2001)x] figure(3) 信号延时
subplot(211)plot(z) title(＇延时时域图＇) 画出延时信号时域图
Zfft(z2001)
subplot(212)plot(n1(11000)Z(11000)) 延时信号频谱图
title(＇延时频谱图 ＇)
sound(5*zfs)
4 混响信号
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
y1x+z 信号混响
figure(4)subplot(211)plot(y1) title(＇混响时域图＇) 混响时域图
Y1fft(y12001) 混响信号2001点FFT变换
subplot(212)plot(n1(11000)Y1(11000)) 混响频谱图
title(＇混响频谱图＇)
sound(5*y1fs) 回放混响信号
5 单回声滤波器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
z[zeros(2001)x] 语音信号进行延时
x[xzeros(2001)] 语音信号延时信号等长度
yy1x+z*05 信号单回声滤波
figure(5)subplot(211)plot(yy1)title(＇单回声滤波器时域图＇)
YY1fft(yy12001) 单回声信号做2001点FFT变换
subplot(212)plot(n1(11000)YY1(11000)) 单回声信号频谱图
title(＇单回声滤波器频谱图＇)
sound(5*yy1fsbits) 回放单回声滤波器语音信号
6 限回声滤波器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
Bz[00000000001] 分子系数
Az[1000000000a] 分母系数
yy2filter(BzAzx) 滤波器进行滤波
figure(6)subplot(211)plot(yy2) 限回声滤波器时域波形
title(＇限回声滤波器时域波形＇)
YY2fft(yy22001) 限回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY2(11000)) 限回声滤波器频谱图
title(＇限回声滤波器频谱图 ＇)
sound(5*yy2fsbits) 回放限回声滤波器语音信号
pause(2)
7 重回声滤波器
[xfsbits]wavread(＇makewav＇)
xx(1) 取单声道
n102000
a05
N5
Bz1[100000000005^N] 分子系数
Az1[100000000005] 分母系数
yy3filter(Bz1Az1x) figure(7) 滤波器进行滤波
subplot(211)plot(yy3) title(＇重回声滤波器时域图＇)重回声滤波器混响器时域波形
YY3fft(yy32001) 重回声滤波器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY3(11000)) 重回声滤波器频谱图
title(＇重回声频谱图＇)
sound(5*yy3fsbits) 回放限回声滤波器语音信号
8 全通结构混响器
[xfsbits]wavread(＇makewav＇[10000 40000]) 读取语音信号
xx(1) 取单声道
n102000
a05 a取等1
Bz1[a0000000001] 分子系数
Az1[1000000000a] 分母系数
yy4filter(Bz1Az1x) figure(8) 滤波器进行滤波
subplot(211)plot(yy4) title(＇全通结构混响器时域图＇)全通结构混响器时域波形
YY4fft(yy42001) 全通结构混响器信号做2001点FFT变换
subplot(212)plot(n1(11000)YY4(11000)) 全通结构混响器频谱图
title(＇全通结构混响器频谱图＇)
sound(5*yy4fsbits) 回放全通结构混响器语音信号
附件6：
科生课程设计成绩评定表
姓 名
性
专业班级
课程设计题目：
课程设计答辩质疑记录：
课程设计评分标准
序 号
容
分值
1
完成相关算法设计
25分
2
通编写调试程序完成基功
35分
3
设计报告完成
20分
4
课程设计答辩
20分
终评定成绩（优良中格格评定）
指导教师签字：
2015年 7 月 10 日
武han理工学
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档
