课程设计说明书-语音信号的数字滤波-应用MATLAB仿真.doc

课程设计说明书 no.1 语音信号的数字滤波 -应用matlab仿真一、课程设计的目的课程设计的目的，是使学生融会贯通本课程所学专业理论知识，完成一个较完整的设计过程，以加深学生对所学理论的理解与应用，培养学生综合运用基础理论知识和专业知识去解决实际工程设计问题的能力。课程设计是工科学不可少的一个综合性实践环节。 课程设计是课程学习后的一个综合性实践环节，是对课程理论和课程实验的综合和补充，是学完一门课后应用本课知识及以前的知识积累而进行的综合性、开放性的训练，是培养学生工程意识和创新能力的重要环节，培养学生综合运用已学过的理论和技能去分析和解决问题的能力，对加深课程理论的理解和应用具有重要意义。二、课程设计方案与论证线性相位fir滤波器通常采用窗函数法设计。窗函数法设计fir滤波器的基本思想是：根据给定的滤波器技术指标，选择滤波器长度n和窗函数（n），使其具有最窄宽度的主瓣和最小的旁瓣。其核心是从给定的频率特性，通过加窗确定有限长单位脉冲响应序列h(n)。工程中常用的窗函数共有6种，即矩形窗、巴特利特（bartlett）窗、汉宁（hanning）窗、汉明（hamming）窗、布莱克曼（blackman）窗和凯塞（kaiser）窗。设计一个线性相位带通fir滤波器指标如下：fn=1000,1375,3625,4000;a=0,1,0;dev=0.0005,0.05,0.0005;阻带最小衰减60db，再设采样频率fs=10khz，则根据阻带最小衰减来选择凯塞（kaiser）窗，利用matlab的求阶函数和fir滤波器的设计函数，可以快速地设计出所需的数字滤波器。这两个设计函数如下： 沈 阳 大 学 课程设计说明书 no.2 n,wn,beta,ftype=kaiserord(fn,a,dev,fs)b=fir1(n,wn,ftype,kaiser(n+1,beta)；最后，利用的滤波器分析函数freqz分析所设计出的滤波器的幅频特性和相频特性，并用图形显示函数plot将它们显示出来，如图1所示。由图1可见，设计结果满足指标要求。(1)根据任务,确定性能指标:在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:带通滤波器的阻带边界频率关于中心频率p0几何对称，因此ws1=wp0- (ws2-wp0)=0.3通带截止频率wc1=0.4,wc2=0.6;阻带截止频率wr1=0.3，wr2=0.7；阻带最小衰减s=3db和通带最大衰减p=15db;(2)用=2/t*tan(w/2)对带通数字滤波器h(z)的数字边界频率预畸变，得到带通模拟滤波器h(s)的边界频率主要是通带截止频率p1,p2;阻带截止频率s1，s2的转换。为了计算简便，对双线性变换法一般t=2s 沈 阳 大 学 课程设计说明书 no.3 通带截止频率wc1=(2/t)*tan(wp1/2)=tan(0.4/2)=0.7265wc2=(2/t)*tan(wp2/2)=tan(0.6/2)=1.3764阻带截止频率wr1=(2/t)*tan(ws1/2)=tan(0.3/2)=0.5095wr2=(2/t)*tan(ws2/2)=tan(0.7/2)=1.9626阻带最小衰减s=3db和通带最大衰减p=15db;(3)运用低通到带通频率变换公式=(2)-(02)/(b*)将模拟带通滤波器指标转换为模拟低通滤波器指标。b=wc2-wc1=0.6499normwr1=(wr12)-(w02)/(b*wr1)=2.236normwr2=(wr22)-(w02)/(b*wr2)=2.236normwc1=(wc12)-(w02)/(b*wc1)=1normwc2=(wc22)-(w02)/(b*wc2)=1得出，normwc=1，normwr=2.236模拟低通滤波器指标：normwc=1，normwr=2.236，p=3db，s=15db(4)设计模拟低通原型滤波器。用模拟低通滤波器设计方法得到模拟低通滤波器的传输函数ha(s);借助巴特沃斯(butterworth)滤波器、切比雪夫(chebyshev)滤波器、椭圆(cauer)滤波器、贝塞尔(bessel)滤波器等。(5)调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。(6)利用窗函数变换法将模拟带通滤波器ha(s)转换成数字带通滤波器h(z).语音信号的频谱分析：语音文件“zjh.wav”频谱分析程序：z=wavread(zjh.wav);y1=z(1:16384);y1=fft(y1);n=0:16383;plot(n,y1);grid; 沈 阳 大 学 课程设计说明书 no.4 沈 阳 大 学 课程设计说明书 no.5三、课程设计结果分析程序1:用matlab对原始语音信号进行分析，画出它的时域波形和频谱fs=22050; %语音信号采样频率为22050x1=wavread(d:数字信号原音.wav ); %读取语音信号的数据，赋给变量x1sound(x1,22050); %播放语音信号y1=fft(x1,1024); %对信号做1024点fft变换f=fs*(0:511)/1024;figure(1)plot(x1) %做原始语音信号的时域图形title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1) %绘制原始语音信号的频率响应图title(频率响应图)figure(3)subplot(2,1,1);plot(abs(y1(1:512) %做原始语音信号的fft频谱图 title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1) %绘制原始语音信号的频率响应图title(频率响应图)figure(3)subplot(2,1,1);plot(abs(y1(1:512) %做原始语音信号的fft频谱图title(原始语音信号fft频谱)subplot(2,1,2); 沈 阳 大 学 课程设计说明书 no.6 plot(f,abs(y1(1:512);title(原始语音信号频谱)xlabel(hz);ylabel(fuzhi); 沈 阳 大 学 课程设计说明书 no.7 沈 阳 大 学 课程设计说明书 no.8 程序2:给原始的语音信号加上一个高频余弦噪声，频率为5khz。画出加噪后的语音信号时域和频谱图，与原始信号对比，可以很明显的看出区别。fs=22050; x1=wavread( d:数字信号原音.wav );f=fs*(0:511)/1024;t=0:1/22050size(x1)-1)/22050;%将所加噪声信号的点数调整到与原始信号相同au=0.03;d=au*cos(2*pi*5000*t);%噪声为5khz的余弦信号x2=x1+d;sound(x2,22050);%播放加噪声后的语音信号y2=fft(x2,1024);figure(1) plot(t,x2)title(加噪后的信号);xlabel(time n);ylabel(fuzhi n);figure(2)subplot(2,1,1);plot(f,abs(y1(1:512);title(原始语音信号频谱);xlabel(hz);ylabel(fuzhi); 沈 阳 大 学 课程设计说明书 no.9 subplot(2,1,2); plot(f,abs(y2(1:512); title(加噪后的信号频谱); xlabel(hz); ylabel(fuzi)程序3:窗函数法设计滤波器：fs=22050;x1=wavread(d:数字信号原音.wav );t=0:1/22050size(x1)-1)/22050;au=0.03; d=au*cos(2*pi*5000*t); x2=x1+d; wp=0.25*pi; ws=0.3*pi; wdelta=ws-wp; n=ceil(6.6*pi/wdelta); wn=(0.2+0.3)*pi/2;b=fir1(n,wn/pi,hamming(n+1); figure(1)freqz(b,1,512)f2=filter(bz,az,x2)figure(2) 沈 阳 大 学 课程设计说明书 no.10 subplot(2,1,1)plot(t,x2) title(滤波前的时域波形); subplot(2,1,2) plot(t,f2); title(滤波后的时域波形); sound(f2,22050); %播放滤波后的语音信号 f0=fft(f2,1024); f=fs*(0:511)/1024; figure(3) y2=fft(x2,1024); subplot(2,1,1); plot(f,abs(y2(1:512); title(滤波前的频谱) xlabel(hz); ylabel(fuzhi); subplot(2,1,2) f2=plot(f,abs(f0(1:512); title(滤波后的频谱) xlabel(hz); ylabel(fuzhi); 沈 阳 大 学 课程设计说明书 no.11 四、设计体会通过这次课程设计，不但了解了数字滤波器的设计过程，同时也学会怎么使用matlab软件。为以后能更好的运用matlab这个工具打下良好的基础，也为我以后的工作生活带来一定的帮助。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说是困难重重，同时在设计的过程中发现了自己的不足之处，对以前所学过的知识理解得不够深刻，掌握得不够牢固。总的来说，这次设计的语音信号的数字滤波器还是比较成功的，在设计中遇到了很多问题，最后在老师的指导和同学的帮助下，终于游逆而解，觉得平时所学的知识 有了实用的价值，达到了理论与实际相结合的目的，不仅学到了不少知识，而且锻炼了自己的能力，使自己对以后的路有了更加清楚的认识。 五、参考文献1.sanjit k. mitral .数字信号处理基于计算机的方法m电子工业出版社.20022.方勇 .数字信号处理原理与实践m.清华大学出版社.200