- 1、本文档共65页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理语音处理课程设实验报告数字信号处理语音处理课程设计实验报告数字信号处理语音处理课程设计实验报告数字信号处理语音处理课程设计实验报告
第十章
一、实验目的
综合运用所学数字信号处理理论知识进行语音信号的采集、频谱分析、延时、滤波处理等,通过理论推导得出相应结论并利用MATLAB作为编程工具进行计算机实现,从而对数字信号处理有一个系统的掌握和直观理解。
二
1、语音信号采集与处理
2、信号频谱分析
3、数字滤波器分析与设计
三实验
四实验指导
通过录音软件录制一段语音,说出自己的名字和学号,命名为“学号+姓名”,时长大约1到2秒,在MATLAB中,通过使用wavread函数,对语音进行采样:
[y,fs,nbits]=wavread(xxxxxxx); %语音信号的采集
采样值放在向量y中,采样频率为fs,采样位数为nbits。
2、语音信号的频谱分析[y,fs,nbits]=wavread(xxxxxxxx); %语音信号的采集
sound(y,fs,nbits); %语音信号的播放
n=length(y) ; %计算语音信号的长度
Y=fft(y,n); %快速傅里叶变换
figure;
subplot(2,1,1); %绘出时域波形
3、原始信号延时混响和回声的处理
(2) 设计如下几种特殊类型滤波器,对原始信号分别进行如下几种回声滤波或混响处理,并画出滤波后信号的时域波形和频谱图,然后分析其变化,最后对声音进行回放,对感觉到的滤波前后的声音变化作出说明。
回声或者混响滤波器类型:
单回声滤波器系统函数(可以用FIR滤波器实现只需要
回声就是原始信号与他自己延时衰减后信号的声音。
方程:y(n)=x(n)+ax(n-r),如上。
实现
③无限回声滤波器系统函数:
④全通结构混响器系统函数:
4、原始信号的滤波处理
(1) 分别用FIR和IIR滤波器对采集的语音信号进行滤波。
(2) 对滤波后的信号进行傅里叶变换,分别得到FIR和IIR滤波后信号的频谱特性,同时画出其时域波形,与原始信号进行比较,分析信号的变化。
(3) 对FIR滤波和IIR滤波后的信号进行比较并作出说明。
(4) 对滤波后的声音进行回放并对比分析变化。
根据频谱特征设计FIR和IIR滤波器。在Matlab中,可以利用函数fir1设计FIR滤波器,利用函数butter,cheby1设计IIR滤波器,利用Matlab中的函数freqz画出各步滤波器的频率响应。
①低通滤波器的性能指标:fp=1000Hz,fc=1200Hz,As=100db ,Ap=1dB
②高通滤波器的性能指标:fp=3500Hz,fc=4000Hz,As=100dB,Ap=1dB;
③带通滤波器的性能指标:fp1=1200Hz,fp2=3000hZ,fc1=1000Hz,fc2=3200Hz,As=100dB,Ap=1dB
在MATLAB中,利用[N,wc]=butter(N,wc,Rp,As,s)设计并计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wc;[B,A]=cheby1(N,Rp,wpo,ftypr)设计切比雪夫I型滤波器。
五实验设计报告要求
通过理论分析和回放等说明滤波前后的信号变化
答:实验,能发现明显变化,时域波形和频域波形都有改变,不同滤波器的不同具体分析在的实验结果中。比如带通滤波器频谱信号几乎没有,通过全通系统则没有变化,经过低通滤波器变化很小。
答:
3、通过使用randn函数产生随机噪声信号,并加到语音信号中得到被污染的语音信号叠加到原始的语音信号中,分析叠加前后信号频谱的变化,设计一个合适的滤波器,能够把噪声部分滤除。
参考文献:
杨述斌,李永全 主编,《数字信号处理实践教程》,华中科技大学出版社,2007.1
实验报告
(1)语音采样和观察clear,clc;
[y,fs]=audioread(E:\大学课程\大三上\数字信号处理\201400121184吴蔓.mp3); %语音信号的采集,把采样值放在y中subplot(3,1,1)
plot(y);
title(时域波形);
sound(y,fs); %语音信号的播放
n=length(y) %计算语音信号的长度
Y=fft(y) ; %快速傅里叶变换
subplot(3,1,2)
plot(abs(Y)); %绘出频域波形
title(幅频特性);
subplot(3,1,
您可能关注的文档
- 数据结构 实验五数据结 实验五数据结构 实验五数据结构 实验五.doc
- 数据从业者的职业规划与大数数据从业者的职业规划与大数据数据从业者的职业规划与大数据数据从业者的职业规划与大数据.ppt
- 数据泛化 data cub computation and data generalization数据泛化 data cube computation and data generalization数据泛化 data cube computation and data generalization数据泛化 data cube computation and data generalization.ppt
- 数据分层汇总交叉报表sql句实现方法数据分层汇总交叉报表sql语句实现方法数据分层汇总交叉报表sql语句实现方法数据分层汇总交叉报表sql语句实现方法.doc
- 数据采集实习报告数据采集实报告习报告.doc
- 数据结构经典算法试题数据结经典算法试题数据结构经典算法试题数据结构经典算法试题.doc
- 数据结构版答案.ppt
- 数据结构 - 递归算法数据构 - 递归算法数据结构 - 递归算法数据结构 - 递归算法.ppt
- 数据结构课程设计说明书规范式戴成秋数据结构课程设计说明书规范格式戴成秋数据结构课程设计说明书规范格式戴成秋数据结构课程设计说明书规范格式戴成秋.doc
- 输液接头的临床应用输液接头临床应用的临床应用.ppt
文档评论(0)