版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数字信号处理课程设计报告专 业:班 级:学 号:姓 名:指导教师:二0 一二年六月十四日一:课程设计目的 4二:课程设计的题目描述和要求 41:设计题目的描述 42:课程设计描述 43:课程设计要求 4三:课程设计报告内容 51:总体设计 51.1 :采集语音信号 51.2 :声音信号的时域和频域分析 51.3 :设计一个余弦噪声,并产生噪声的时域和频域图,并分析51.4 :对声音信号加噪声61.5 :设计bf低通滤波器,并对其进行频谱分析 61.6 :对加噪声的信号滤波,并回放,并分析频谱 72:软件仿真调试结果分析 72.1 :运行语音信号程序 82.2 :运行噪声信号 82.3 :运行加
2、噪声后的语音信号 92.4 运行滤波器程序 102.5 :滤波后的语音信号 11四:设计总结 12源程序代码: 1318一:课程设计目的1 :通过对课程设计任务的完成,进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法;2:进一步扩展我们掌握的基本理论和分析方法方面的知识,能有效地将理论和实际紧密结合;3:同时增强自己的软件编程实现能力和解决实际问题的能力。4: 练习能够熟练地用 matlab 语言编程实现iir 数字滤波器和 fir 数字滤 波器,进一步明确数字信号处理的工程应用。二:课程设计的题目描述和要求1: 设计题目的描述基于 matlab 的 iir 的 butterwort
3、h 滤波器设计目标: ( 1)语音采集( 2)设计滤波器,将采集的语音信号进行iir 滤波 .( 3)回放语音信号2:课程设计描述要求录制一段自己的语音信号后,格式为 wav在matla啾件中采集语音 信号、回放语音信号。画出语音信号的时域波形和频谱图, 对所采集的信号加入干扰噪声, 要求噪声在3000hz 以上,对加入噪声进行播放,并进行时域和频谱分析,对比加噪声前后的时域图和频谱图, 对比加噪声音信号和通过低通滤波器处理的音频信号不同的时域和频域波形。3: 课程设计要求通过本次课程设计加深我对所学知识的理解和认识。并通过的实验来验证仿真自己的理解。理解了总体滤波器的设计。运用 matlab
4、 观察信号的处理操作过程。 在加深自己 matlab 动手实验能力的同时,独立运用自己的数字信号处理知识来处理音频信号。通过本次课程设计基本掌握matlab 软件的编程方法。并进行相应的拓展。练习matlab gui的编写,通过按钮来调用程序,进行操作三:课程设计报告内容1 :总体设计1.1 :采集语音信号通过录音机,录一段自己的声音,存名为 xx.wav保存到桌面上1.2 :声音信号的时域和频域分析通过wavread函数读取录制的声音文件,进彳t fft转换。显示其时域和频域的 波形习 editor - lkitihed1file edh text cell took debug deskt
5、op window help口安。. 咯6小 昌晶九层由cd曰日叱y; f s, nbi t s j -wavread ( c i aadmi ni at ratoree ski opqj 3.2 ound(y, f/ nbit/;拓播放语音信号3 n-length -y),用计算采样点的长.至4 y=fft (y, n). 它n变换5 subplot (2f 1, l);6 riot(y);7 tit is c原始信号波形);xlabelc时间),了1。子(幅值).8 subplot (2; l 2);9 plot(abs(y):10 titlec 原始借号频谱);xlabelc 频率) ;
6、ylabelc11 -bl 苏 |wti /1.3:设计一个余弦噪声,并产生噪声的时域和频域图,并分析1.4:对声音信号加噪声对声音信号加一个频率在3000hz以上的余弦信号,然后对加噪声以后的信号进行时域和频域分析1.5:设计bf低通滤波器,并对其进行频谱分析1.6:对加噪声的信号滤波,并回放,并分析频谱2:软件仿真调试结果分析2.1:运行语音信号程序结果分析:观察上图可以看到,在采集的原始信号进行付氏变换后所得的频域波形就是原语音信号的频谱。2.2 :运行噪声信号结果分析:可以看到噪声的频谱是对称的冲击2.3 :运行加噪声后的语音信号结果分析:观察上图可以看见加噪的音频时域频域波形均出现了
7、很大的变化, 通 过matla歆件所播放的加噪声音也是相应的出现了很大的噪声。 加噪程序成功。2.4 :运行滤波器程序结果分析:观察上图可以清晰的看到器幅频特性为低通滤波器,还可以分析上图得我们所设置的采样频率、通带截止频率、阻带截止频率和阶数。2.5 :滤波后的语音信号结果分析:通过观察两个不同信号的频谱我们可以清晰的看到在截止频率周围器 频谱的变化。因为是iir滤波器的缘故,其频谱是相应周期延托。故而低通滤波 器处理信号结果是如同上右下图的。试验成功。四:设计总结通过这次数字信号处理的课程设计, 我懂得了理论与实际相结合的重要性, 从理 论中得出结论,从实践中提高自己。才能真正的锻炼自己,
8、虽然这次的课程设计 经历了半个月,但是我还是学会了很多很多的东西,在巩固以前所学过的知识同 时,而且学到了很多在书本上所没有学到过的知识。通过这次设计,进一步加深了真正培养了自己的独立思考能力与动手能力。这次试验也让我认识到了自身的很多不足。 其一就是对于数字信号处理课程知识 的不透彻,让我在一些基础东西的处理上显得不是很游刃有余。特别是在编写程 序时,发现自己还是有很大的提升空间的。 实验虽然结束了,但我们学习还在继 续。我将更加努力。虽然也受了很多的挫折,不过也正是这些的挫折让我收获了更多,不只是信心的 提升。同样强化了我以后面对挫折的不懈心态。这次课程设计收获很多。源程序代码:原始语音信
9、号y,fs,nbits=wavread (c:usersadministratordesktopqjs.wav);% 读取语音信号的数据,赋给变量ysound(y,fs,nbits); %播放语音信号n=length (y) ;%计算采样点的长度y=fft(y,n); %fft 变换subplot(2,1,1);plot(y);title(原始信号波形);xlabel(时间 n);ylabel(幅值);subplot(2,1,2);plot(abs(y);title(原始信号频谱);xlabel(频率hz);ylabel(幅值);噪声信号y,fs,nbits=wavread (c:usersa
10、dministratordesktopqjs.wav);% 读取语音信号的数据,赋给变量yn=length (y) ;%计算采样点的长度t=0:1/44100:(size(y)-1)/44100; %将所加噪声信号的点数调整到与原始信号相同au=0.08d=au*cos(2*pi*3500*t); % 噪声为 3.5khz 的余弦信号y1=d;%将噪声进行转置,以便使噪声与信号的行列相同,进行相加sound(y1,fs,nbits);y1=fft(d,3800); %fft 变换subplot(2,1,1);plot(d);title( 噪声的波形);subplot(2,1,2);plot(a
11、bs(y1); title( 噪声的频谱);加噪声之后的语音信号y,fs,nbits=wavread (c:usersadministratordesktopqjs.wav);% 读取语音信号的 数据,赋给变量y%sound(y,fs,nbits); %播放语音信号 n= length (y) ; %计算采样点的长度 y=fft(y,n); %fft 变换 subplot(2,2,1);plot(y);title(原始信号波形);subplot(2,2,2);plot(abs(y); title(原始信号频谱);t=0:1/44100:(size(y)-1)/44100; %将所加噪声信号的点
12、数调整到与原始信号相同au=0.08d=au*cos(2*pi*3500*t); % 噪声为 3.5 khz 的余弦信号y1=d;%将噪声进行转置,以便使噪声与信号的行列相同,进行相加x1=y+y1;sound(x1,fs,nbits);subplot(2,2,3);plot(x1);title( 加噪语音信号的时域波形 );s=fft(x1);%傅里叶变换subplot(2,2,4);plot(abs(s); title(加噪语音信号的频域波形);滤波器的程序wp=0.08*pi;ws=0.1*pi;% 滤波器的指标rp=1;%通带内的最大衰减rs=15;%阻带内的最小衰减fs=44100;
13、ts=1/fs;wp1=2/ts*tan(wp/2); %将数字指标转换成模拟指标ws1=2/ts*tan(ws/2);n,wn=buttord(wp1,ws1,rp,rs,s);% 选择滤波器的最小阶数z,p,k=buttap(n); %创建 butterworth 模拟滤波器bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn); % 低通到低通bz,az=bilinear(b,a,fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换h,w=freqz(bz,az); %绘制频率响应曲线figure(1)plot(w*fs/(2*pi),abs(h)gr
14、id;% 网格xlabel(频率/ hz);ylabel(频率响应幅度)title(butterworth);滤波的程序y,fs,nbits=wavread (c:usersadministratordesktopqjs.wav);% 读取语音信号的数据,赋给变量yn=length (y) ; %计算采样点的长度t=0:1/44100:(size(y)-1)/44100; %将所加噪声信号的点数调整到与原始信号相同au=0.08d=au*cos(2*pi*3500*t); % 噪声为 3.5 khz 的余弦信号y1=d;%将噪声进行转置,以便使噪声与信号的行列相同,进行相加x1=y+y1;wp
15、=0.08*pi;ws=0.1*pi;rp=1;rs=15;fs=44100;ts=1/fs;wp1=2/ts*tan(wp/2); %将数字指标转换成模拟指标ws1=2/ts*tan(ws/2);n,wn=buttord(wp1,ws1,rp,rs,s);% 选择滤波器的最小阶数z,p,k=buttap(n); %创建 butterworth 模拟滤波器bap,aap=zp2tf(z,p,k);b,a=lp2lp(bap,aap,wn);bz,az=bilinear(b,a,fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换h,w=freqz(bz,az); %绘制频率响应曲线subplot(2,3,1)plot(w*fs/(2*pi),abs(h)gridxlabel(频率/ hz)ylabel(频率响应幅度)title(butterworth)f1=filter(bz,az,x1);sound(f1,fs,nbits);subplot(2,3,2)plot(t,x1) %画出滤波前的时域图title(滤波前的时域波形);subplot(2,3,3)plot(t,f1); %画出滤波后
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二五年度集装箱运输企业信用评价与风险管理合同3篇
- 二零二五年环保节能型监控设备采购与技术支持合同2篇
- 二零二五版房屋租赁及转让合同全方位权益创新协议2篇
- 二零二五版文化创意产业园区使用权转让合同3篇
- 二零二五年度国际公路运输代理合同2篇
- 二零二五版城市绿化苗木租赁合同3篇
- 二零二五版环保设备质押贷款合同模板3篇
- 二零二五年度高级管理人员出差责任免除服务合同范本2篇
- 二零二五版体育行业劳动合同管理规范及运动员权益保障协议3篇
- 二零二五年度节水减排供水合同范本3篇
- 2023年山东省青岛市中考化学试题(含答案解析)
- 商业计划书(BP)产品与服务的撰写秘籍
- 安徽华塑股份有限公司年产 4万吨氯化石蜡项目环境影响报告书
- 公司章程(二个股东模板)
- 世界奥林匹克数学竞赛6年级试题
- 药用植物学-课件
- 文化差异与跨文化交际课件(完整版)
- 国货彩瞳美妆化消费趋势洞察报告
- 云南省就业创业失业登记申请表
- UL_标准(1026)家用电器中文版本
- 国网三个项目部标准化手册(课堂PPT)
评论
0/150
提交评论