还剩12页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
语音信号处理与滤波班级XXXXXX学号XXXXX姓名XXX指导老师XXX成绩________________________________一五年五月二十七O比较滤波先后语音信号的波形及频谱
3.6利用切比雪夫滤波器,程序如下程序结果如图
3.6wp=1000*2*pi;ws=1200*2*pi;Rp=
0.3;Rs=40;[N,Wn]=cheb1ordwp,ws,Rp,Rs/s;[z,p,kj=cheb1apN,Rp;[A,B,C,D]=zp2ssz,p,k;[AT,BT,CT,DT]=lp21pA,B,C,D,Wn;[b,a]=ss2tfAT,BT,CT,DT;[bz,az]=impinvarb,a,fs;[H,W]=freqzbz,az;X=filterbz,az,x%soundX,fs,bits;figure2Xl=fftXJ024subplot221;处理后的信号波形上pk tXsubplot222;处理后的信号频谱;plotn1:512,absX11:512;titleo处理后的佰号渡形处理后的循号频谱40246823x10图
3.6
四、心得体味在本次课程设计的过程中,主要遇到的问题是原始语音信号的采集有些语音信号在MATLAB中运行时不能达到理想中的状态,后来经过多次实验和对照,发现,
一、录音时声音要足够的大,才干采集好的语音信号;
二、要将录音后的文件后缀wma转换成wav的格式其次就是在滤波器参数的选取上,开始时不清晰怎么选取,后来通过查找书本相关知识和与同学讨论才有效的解决这个问题因此学习的过程是相互讨论共同进步的,多多讨论课题中遇到的问题,可以巩固我们的知识掌握能力,增加熟练运用度从本次课程设计的中心来看,课题是希翼将数字信号处理技术应用于某一实际领域,这里就是指对语音信号的处理作为存储于计算机中的语音信号,其本身就是离散化了的向量,我们只需将这些离散的量提取出来,就可以对其进行处理了在这里,用到了处理数字信号的强有力工具MATLAB,通过MATLAB里几个命令函数的调用,很轻易的在实际化语音与数字信号的理论之间搭起了一座桥梁课题的特色在于它将采集到的语音看做了一个向量,于是语音数字化了,则可以彻底利用数字信号处理的知识来解决我们可以像给普通信号做频谱分析一样,来给语音信号做频谱分析,也可以较容易的用数字滤波器来对语音进行滤波处理同时改变相关参数,理论结合实际,可以分析出各参数对图形的影响,从而加深对各个参数的理解在完成这次课程设计过程中我学到了许多东西,进一步理解了滤波器设计的方法和各参数意义,通过分析信号时域和频域的关系等,加深了对滤波性能的理解,而且学会了使用Matlab一些基本函数,增加了进一步学习Matlab软件的兴趣同时,通过本次课程设计,使我懂得了理论与实际相结合是很重要的,惟独理论知识是远远不够的,惟独把所学的理论知识与实践相结合起来,从理论中得出结论,才干真正学有所用,从而立锻炼了我的动手能力,和提高了我分析问题,解决问题的能力
五、源程序代码fs=8000,bits=8,T=l,Ts=1/fs;N=T/Ts;;x=x:,lfigure1subplot22l;plotx;%soundx,fs,bits;时域波形图titleoy=fftx,1024;n=fs/1024*[l:1024];subplot222;plotnl:512,absyl:512;频谱图’;title odt=
0.01*sin2*pi*7000*l:sizex/fs;%dtl=fftdt,1024;%plotnl:512,absdtll:512;xl=x+dt;%soundxl,fs,bits;yl=fftxl,1024;subplot223;plotxl;污染信号波形;titleosubplot224;plotnl:1024,absy11:1024;axis[0,12000,0,
0.4J;污染信号频谱;titleowp=1000*2*pi;ws=1200*2*pi;Rp=
0.3;Rs=40;[N,Wn]=cheblordwp,ws,Rp,Rs,s;[z,p k]=cheb1apN,Rp;9[A,B,C,D]=zp2ssz,p,k;[AT,BT,CT,DT]=lp21pA,B,C,D,Wn;[b,a]=ss2tfAT,BT,CT,DT;[bz,az]=impinvarb,a,fs;[H,W]=freqzbz,az;X=filterbz,az,x%soundX,fs,bits;figure2Xl=fftX,1024subplot221;,处理后的信号波形;plotX;titlesubplot222;处理后的信号频谱;%主程序plotnl:512,absXll:512;titled wp=1000*2*pi;ws=1200*2*pi;fs=2500;Rp=
0.3;Rs=40;[N,WnJ=cheb1ordwp,ws,Rp,Rs/s;[z,p,k]=cheblapN,Rp;[A,B,C,D]=zp2ssz,p,k;[AT,BT,CT,DT]=lp21PA,B,C,D,Wn;[b,a]=ss2tfAT,BT,CT,DT;[bz,az]=impinvarb,a,fs;[H,W]=freqzbz,az;plotW*fs/2*pi,absH;切比雪夫滤波器幅频特性grid;titleoxlabelCTrequency/Hz1;ylabelmagnitude;%切比雪夫滤波器
六、参考文献
1.张建生等.数字信号处理原理与应用.东南大学出版社.
2022.
2.胡航.语音信号处理.哈尔滨工业大学出版社.
2022.
3.张威.MATLAB基础与编程入门.西安电子科技大学出版社.
2022.目录
一、设计要求1
二、设计步骤
12.1理论依据
12.2信号采集
12.3构造受干扰信号并对其进行FFT频谱分析
12.4数字滤波器设计
12.5信号处理2
三、课程设计实现
23.1语音信号的采集
23.2语音信号的FFT频谱分析
23.3构造受干扰信号并对其进行FFT频谱分析
43.4设计数字滤波器
63.5用滤波器对加噪语音信号进行滤波
73.6比较滤波先后语音信号的波形及频谱7
四、心得体味9
五、参考文献错误!未定义书签
六、源程序代码11
一、设计要求本次课程设计要求利用MATLAB对语音信号进行数字信号处理和分析,要求采集语音信号后,在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号
二、设计步骤理论依据
1.1根据设计要求分析系统功能,掌握设计中所需理论(采样频率、采样位数的概念,采样定理;时域信号的FFT分析;数字滤波器设计原理和方法,各种不同类型滤波器的性能比较),阐明设计原理信号采集
1.2采集语音信号,并对其进行频谱分析,画出信号时域波形图和频谱图FFT构造受干扰信号并对其进行频谱分析
1.3FFT对所采集的语音信号加入干扰噪声,对语音信号进行回放,感觉加噪先后声音的变化,分析原因,得出结论并对其进行FFT频谱分析,比较加噪先后语音信号的波形及频谱,对所得结果进行分析,阐明原因,得出结论数字滤波器设计
1.4根据待处理信号特点,设计合适数字滤波器,绘制所设计滤波器的幅频和相频特性信号处理
2.5用所设计的滤波器对含噪语音信号进行滤波对滤波后的语音信号进行FFT频谱分析画出处理过程中所得各种波形及频谱图对语音信号进行回放,感觉滤波先后声音的变化比较滤波先后语音信号的波形及频谱,对所得结果和滤波器性能进行频谱分析,阐明原因,得出结论
三、课程设计实现语音信号的采集
3.1利用PC机上的声卡和WINDOWS操作系统可以进行数字信号的采集将话筒输入计算机的语音输入插口上,启动录音机按下录音按钮,接着对话筒说话,说完后住手录音,屏幕左侧将显示所录声音的长度点击放音按钮,可以实现所录音的重现以文件名“
12022.wav”保存入中,将文件以wav,格式保存语音信号的频谱分析
3.1FFT通过用windows录音之后将录音的文件导入到MATLAB中,利用wavread函数对语音信号进行采样,并使用MATLAB绘出采样后的语音信号的时域波形和频谱图Wavread函数调用格式%采样值放在向量中,表示采样频率[y,fs,nbits]=wavreadfile yfs读取所规定的文件,返回采样值放在向量中y=wavread file%file wavy表示采样位数hz,nbits%读取钱点的采样值放在向量中y=wavread file,N Ny读取从到点的采样值放在向量中对语音信号y=wavreadfile,[N1,N2]%N1N2y speech,进行采样其程序如下wav[y,fs,nbits]=wavered speech1;%巴语音信号进行加载入Mat lab仿真软件平台中fs=8000nbits=8然后画出语音信号的时域波形,然后对语音信号进行频谱分析在matlab中利用ff对信号进行快速傅里叶变换,得到信号的频谱特性其程序如下fs=8000,bits=8,T=1,Ts=1/fs;N=T/Ts;x=x:,l;subplot32l;plotx;soundx,fs,bits;时域波形图,titly=fftx,1024;n=fs/1024*[l:1024];subplot322;plotnl:512,absyl:512;频憎图,;titl程序结果如图
3.2附辘腿加豳图
3.2构造受干扰信号并对其进行频谱分析
3.2FFT在MATLAB软件平台下,给原始的语音信号叠加之噪声噪声类型分为如下几种1白噪声;2单频噪声正弦干扰;3多频噪声多正弦干扰;4其它干扰,可设置为低频、高频、带限噪声,或者Chirp干扰、冲激干扰这里我们叠加单频正弦干扰噪声其程序如下dt=
0.01*sin2*pi*7000*1:sizex/fs;%dtl=fftdt,1024;%plotnl:512,absdt11:512;二xl x+dt;%sound xl,f s,bi ts;yl=fftxl,1024;subplot223;plot xl;污染信号波形,;titleo subplot224;plot nl:1024,abs yll:1024;污染信号频谱,;axis[0,12000,0,
0.4];titleo程序结果如图
3.3:之后对语音信号进行回放,感觉加噪先后声音的变化,并对其进行FFT频谱分析,比较加噪先后语音信号的波形及频谱,并对所得结果进行分析设计数字滤波器
3.3在设计一个数字滤波器之前,必须首先根据实际的需要来确定滤波器的技术指标,设计合适的数字滤波器,并画出滤波器的幅频特性这里我们设计了一个切比雪夫滤波器其程序如下:wp=1000*2*pi;ws=1200*2*pi;fs=2500;Rp=
0.3;Rs=40;[N,Wn]=cheblordwp,ws,Rp,Rs/sf;[z,p,k]=cheblapN,Rp;[A,B,C,D]=zp2ssz,p,k;[AT,BT,CT,DT]=lp21pA,B,C,D,Wn;[b,a]=ss2tfAT,BT,CT,DT;[bz,az]=impinvarb,a,fs;[H,W]=freqzbz,az;plotW*fs/2*pi,absH;切比雪夫滤波器幅频特性,grid;titleC xlabeKTrequency/Hz*;ylabelmagnitude;%切比雪夫滤波器切比雪夫滤波器的幅频特性如图
3.
41.
41.
20.
40.20图
3.4用滤波器对加噪语音信号进行滤波
3.4用自己设计的切比雪夫滤波器对加噪的语音信号进行滤波,在Mat lab中,HR滤波器利用函数filt对信号进行滤波切比雪夫滤波器X=filterbz,az,x。