华科电信数字信号处理实验报告.doc_第1页
华科电信数字信号处理实验报告.doc_第2页
华科电信数字信号处理实验报告.doc_第3页
华科电信数字信号处理实验报告.doc_第4页
华科电信数字信号处理实验报告.doc_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数字信号处理实验报告姓 名:_ _学 号:_ _专业班级:_ _联系电话:_ _2016年11月 24 日一实验目的1、掌握使用 FFT 进行信号谱分析的方法2、设计数字滤波器对指定的语音信号进行滤波处理二实验内容1、使用 Matlab 的 fft 函数对语音信号进行频谱分析,找出干扰信号的频谱;2、设计数字滤波器滤除语音信号中的干扰分量,并进行播放对比。三实验环境Window xpMatlab 2008b四实验内容1. 观察分析原始语音信号特征：原始信号时域图：（横坐标为以1/fs为单位的时间，纵坐标为信号幅值）原始信号频域图：（横坐标为频率，纵坐标为信号fft幅值）分析：从频域图可以看出噪声在0Hz，1575Hz，3150Hz，4725Hz。用hn=1,A,1进行滤波，|H(ejw)|=|A+2cosw|，在噪声频率处令|H(ejw)|=0算出A，用hn和原始信号卷积即可滤除噪声。2.时域和频域滤波效果滤波后信号时域图：（横坐标为以1/fs为单位的时间，纵坐标为信号幅值）滤波后信号频域图：（横坐标为频率，纵坐标为信号fft幅值）五实验结果分析从滤波后信号时域图可以看出原来幅值反常噪声已经滤除，滤波后信号频域图原来幅值很大的8个冲激已经被滤除了。六MATLAB源代码clear; %读取声音信号xn,fs = wavread(SunshineSquare.wav); n = 0:length(xn)-1;f = n*fs/length(xn);XK = fft(xn); figure(1); subplot(2,2,1);plot(xn);title(原始信号时域图);ylabel(xn);subplot(2,2,2);plot(f, abs(XK);title(原始信号频域图);xlabel(f/Hz);ylabel(XK); %4次滤波yn = band(xn, 0, fs);yn = band(yn, 1575, fs);yn = band(yn, 3150, fs);yn = band(yn, 4725, fs); yn(85440:85455) = xn(1:16);yn(118500:length(yn)=xn(1:length(yn)-118500+1); n1 = 0:length(yn)-1;f1 = n1*fs/length(yn);YK = fft(yn);subplot(2,2,3);plot(yn);title(滤波后信号时域图);ylabel(yn);subplot(2,2,4);plot(f1, abs(YK);title(滤波后信号频域图);xlabel(f/Hz);ylabel(YK); wavwrite(yn,fs,SunshineSquare_New.wav);function y = band(xn, f, fs)%设计滤
人人文库> 全部分类> 应用文书 > 年终总结

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论