MATLAB结课论文设计 (3)

数字信号处理课程设计报告

图3.3 滤波前后波形

由图1.2中滤波前后波形比较可看出，经过滤波后的波形比原波形的振幅有所减小，去除了很多由于噪声所产生的干扰；从滤波前后的频谱比较可以看出经过滤波后除了原本的声音外，中间由于噪声产生的频谱波形已经滤除；经过MATLAB仿真，听滤波前后的声音，可以听出有明显的滤波效果。

7

数字信号处理课程设计报告

第 4 章 总结

语音信号处理是语音学与数字信号处理技术相结合的交叉学科，课题在这里不讨论语音学，而是将语音当做一种特殊的信号，即一种“复杂向量”来看待。也就是说，课题更多的还是体现了数字信号处理技术。

从课题的中心来看，课题“基于MATLAB的有噪声语音信号处理”是希望将数字信号处理技术应用于某一实际领域，这里就是指对语音及加噪处理。作为存储于计算机中的语音信号，其本身就是离散化了的向量，我们只需将这些离散的量提取出来，就可以对其进行处理了。这一过程的实现，用到了处理数字信号的强有力工具MATLAB。通过MATLAB里几个命令函数的调用，很轻易的在实际语音与数字信号的理论之间搭了一座桥。

课题的特色在于它将语音信号看作一个向量，于是就把语音数字化了。那么，就可以完全利用数字信号处理的知识来解决语音及加噪处理问题。我们可以像给一般信号做频谱分析一样，来对语音信号做频谱分析，也可以较容易的用数字滤波器来对语音进行滤波处理。通过比较加噪前后，语音的频谱和语音回放，能明显的感觉到加入噪声后回放的声音与原始的语音信号有很大的不同,前者随较尖锐的干扰啸叫声。从含噪语音信号的频谱图中可以看出含噪声的语音信号频谱,在整个频域范围内分是布均匀。其实，这正是干扰所造成的。通过滤波前后的对比，低通滤波后效果最好，高通滤波后的效果最差。由此可见，语音信号主要分布在低频段，而噪声主要分布在高频段。

本人签字：

8

数字信号处理课程设计报告

参考文献

1.高西全，丁玉美．数字信号处理．第3版．北京：西安电子科技大学出版社，2008 2.刘泉，阙大顺．数字信号处理原理与实现．北京：电子工业出版社，2005 3 .张磊，毕靖，郭莲英．MATLAB实用教程．北京：人民邮电出版社，2008 4 .张威．MATLAB基础与编程入门．西安：西安电子科技大学出版社，2006 5 .周利清，苏菲．数字信号处理基础．北京：北京邮电大学出版社，2005 6.蒋瑞艳 振动，噪声处理系统[J]. 大连理工大学学院学报(自然科学版) 2002（1） 7.崔晓娟 基于MATLAB的舰船辐射噪声信号小波消噪处理[J]. 《舰船科学技术》 2006（06）

9

数字信号处理课程设计报告

附录

源代码

%对语言信号做原始的时域波形分析和频谱分析 [y,fs,bits]=wavread('abc.wav'); % sound(y,fs) % 回放语音信号 n=length(y) %选取变换的点数

y_p=fft(y,n); %对n点进行傅里叶变换到频域 f=fs*(0:n/2-1)/n; % 对应点的频率 figure(1)

subplot(2,1,1);

plot(y); %语音信号的时域波形图 title('原始语音信号采样后时域波形'); xlabel('时间轴') ylabel('幅值 A') subplot(2,1,2);

plot(f,abs(y_p(1:n/2))); %语音信号的频谱图 title('原始语音信号采样后频谱图'); xlabel('频率Hz'); ylabel('频率幅值');

%对加噪的语音信号进行去噪程序如下：

fp=1500;fc=1700;As=100;Ap=1; %(以上为低通滤波器的性能指标）

wc=2*pi*fc/fs; wp=2*pi*fp/fs; wdel=wc-wp;

beta=0.112*(As-8.7);

N=ceil((As-8)/2.285/wdel); wn= kaiser(N+1,beta); ws=(wp+wc)/2/pi; b=fir1(N,ws,wn); figure(2);

freqz(b,1);

%（此前为低通滤波器设计阶段）——接下来为去除噪声信号的程序—— x=fftfilt(b,y); X=fft(x); figure(3);

10

数字信号处理课程设计报告

subplot(2,2,1);plot(f,abs(y_p(1:n/2))); title('滤波前信号的频谱');

subplot(2,2,2);plot(f,abs(X(1:n/2))); title('滤波后信号的频谱'); subplot(2,2,3);plot(y);

title('滤波前信号的时域波形'); subplot(2,2,4);plot(x);

title('滤波后信号的时域波形'); sound(x,fs,bits) den=2*ifft(y_p); wavwrite('2.wav')11

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新综合文库MATLAB结课论文设计 (3)全文阅读和word下载服务。