FIR滤波器的设计报告

-.z---.--总结资料?数字信号处理课程设计?课程设计报告专业：电子信息工程班级：****学号：*****姓名：**指导教师：********2011年6月28日目录1、课程设计目的22、课程设计题目描述和要求22.1题目描述22.2设计要求23、课程设计报告内容33.1总体设计33.1.1采集声音信号3声音信号的时域与频谱分析33.1.3对声音信号加噪43.1.4滤波器设计53.1.5加噪信号的滤波与回放63.2软件仿真调试结果分析63.2.1原始信号的读取与分析63.2.2加噪信号的时域与频域分析73.2.3滤波器的设计及分析83.2.4加噪信号的滤波与图形分析93.3设计中出现的问题与解决方案114、设计总结11参考文献12附录：〔源程序代码〕121、课程设计目的通过对课程设计任务的完成，使学生进一步稳固数字信号处理的根本概念、理论、分析方法和实现方法；使学生掌握的根本理论和分析方法方面的知识得到进一步扩展；使学生能有效地将理论和实际严密结合；增强学生软件编程实现能力和解决实际问题的能力。要求学生能够熟练地用Matlab语言编程实现IIR数字滤波器和FIR数字滤波器，进一步明确数字信号处理的工程应用。2、课程设计题目描述和要求2.1题目描述FIR数字滤波器：FIR数字滤波器的冲激响应的Z变换为：其中滤波器最重要的两个特性为线性相位特性和幅度特性。本次课程设计要求录制一段自己的语音信号后，格式为WAV。在MATLAB软件中采集语音信号、回放语音信号。画出语音信号的时域波形和频谱图，对所采集的信号参加干扰噪声，对参加噪声进展播放，并进展时域和频谱分析；比照早前后的时域图和频谱图即比照加噪声音信号和通过低通滤波器处理的音频信号不同的时域和频域波形。2.2设计要求理解课题教学的理论内容→掌握和熟悉matlab软件的编程方法和仿真过程→编写课程设计所需的程序→对程序进展仿真测试，并进展分析→写出课程设计报告3、课程设计报告内容3.1总体设计3.1.1采集声音信号利用PC机上的声卡和WINDOWS操作系统自带的录音机，来进展声音的采集，将录的音以文件名“a〞保存在F：\中，文件存储的后缀默认为.wav。3.1.2声音信号的时域与频谱分析在MATLAB软件平台下，利用MATLAB中的“wavread〞命令来读入〔采集〕语音信号，将它赋值给*一向量。再对其进展采样，记住采样频率和采样点数。通过使用wavread函数，可以深入理解采样频率、采样位数等概念，下面介绍Wavread函数几种调用格式。〔1〕y=wavread〔file〕；读取file所规定的wav文件，返回采样值放在向量y中。〔2〕[y,fs,nbits]=wavread(file)；采样值放在向量y中，fs表示采样频率〔hz〕，nbits表示采样位数。〔3〕y=wavread〔file，N〕；读取钱N点的采样值放在向量y中。〔4〕y=wavread〔file，[N1,N2]〕；读取从N1到N2点的采样值放在向量y中。接下来，对语音信号a.wav进展采样，采样程序如下：[y,fs,nbits]=wavread('F\a');即把语音信号加载在Matlab仿真软件平台上。然后，画出语音信号的时域波形，再对语音信号进展频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：*k=fft(*n,N)参数*n为被变换的时域序列向量，N是DFT变换区间长度，当N大于*n的长度时，fft函数自动在*n后面补零。，当N小于*n的长度时，fft函数计算*n的前N个元素，忽略其后面的元素。在本次设计中，我们利用fft对语音信号进展快速傅里叶变换，就可以得到信号的频谱特性。截图如下：3.1.3对声音信号加噪先对噪声信号进展设置，将噪声信号和原信号进展叠加。然后对叠加后的信号进展时域与频域的仿真及分析。如下：3.1.4滤波器设计用汉明窗〔hamming〕函数对滤波器进展设计，选择适宜的参数，然后分析其幅频特性。如下：3.1.5加噪信号的滤波与回放对叠加后的信号进展滤波，观察滤波前后波形在时域与频域的变化：3.2软件仿真调试结果分析3.2.1原始信号的读取与分析对3.1.2中的程序进展仿真，所采集的信号的时域及频谱图形如下:3.2.2加噪信号的时域与频域分析叠加信号和原信号的时域波形仿真如下：叠加信号和原信号的频域特性如下：3.2.3滤波器的设计及分析由频谱图可以看出信号能量集中在低频局部，为了减少无用的高频分量我设计了以下的FIR低通滤波器将其滤除。该设计采用窗函数法〔汉明窗〕来完成以下指标的低通滤波器。通带截频fp=3000Hz，阻带截频fs=3500Hz，通带最大衰减0.3dB，阻带最小衰减-50dB，采样频率fs为信号源采样频率22050Hz。分析：在仿真图上可以看出来所设计的滤波器的技术指标：通带截止频率fp、阻带截止频率fs、通带最大衰减0.3db、阻带最小衰减-50db、采样频率为22050hz。3.2.4加噪信号的滤波与图形分析三种信号的时域比较：分析：从仿真的结果上可以看出来，原信号与经过滤波后的加噪信号在时域波形产生了变化，但总体上还是与原信号近似；原因是所设计的滤波器比较简单，没有充分考虑所有的因素，从而造成两个信号的差异性。结果证明：加噪及简单的滤波器的设计良好。三种信号的频域比较：3.3设计中出现的问题与解决方案在设计的过程中，刚开场的时候对于用MATLAB编写的程序的具体意思还是不太明白，经过熟悉以及向教师及同学询问后，才渐渐的明白每一个语句的具体含义。4、设计总结在本次课程设计中，由于对matlab不是很熟悉，用着很生疏，只能通过查询相关资料自己学习，最后还是抑制了这个难题，能过正确运用matlab的根本操作，通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能提高自己的实际动手能力和独立思考的能力。同时在设计的过程中发现了自己的缺乏之处，对以前所学过的知识理解得不够深刻，掌握得不够结实。参考文献陈后金，?DSP技术及应用?，机械工业，2004年?基于matlab7.*的系统分析与设计?.楼顺天.**电子科技大学附录：〔源程序代码〕1、原始信号的时域及频域分析的源代码如下：[y,fs,nbits]=wavread('F:\a');sound(y,fs,nbits);%回放语音信号N=length(y);%求出语音信号的长度Y=fft(y,N);%傅里叶变换subplot(2,1,1);plot(y);title('原始信号波形');subplot(2,1,2);plot(abs(Y));title('原始信号频谱')2、声音加噪后的时域及频域的分析代码如下：[y,Fs,bits]=wavread('F:\a');%读出信号，采样率和采样位数。y=y(:,1);%取单声道作分析yl=length(y)%求语音信号长度yy=fft(y,yl);%傅里叶变换t=[0:1/8000:4zeros(1,yl-32001)]';m=0.07*sin(10000*pi*t);%产生噪声n=y+m;%参加噪声nl=length(n)%求语音信号长度nn=fft(n,nl);%傅里叶变换figure(1);subplot(2,1,1);plot(n);title('噪声信号波形')subplot(2,1,2);plot(y);title('原信号波形')figure(2);subplot(2,1,1);plot(abs(nn));title('噪声信号频谱');subplot(2,1,2);plot(abs(yy));title('原信号频谱');sound(n,fs)3、滤波器的幅频特性分析代码如下：Wp=3000*2*pi;Ws=3500*2*pi;Ap=0.3;As=50;fs=22050;wp=Wp/fs;ws=Ws/fs;%归一化通带截频，阻带截频wdelta=ws-wp;%过渡带宽%由阻带最小衰减指标确定选用汉明窗，并确定其阶数N=ceil(6.6*pi/wdelta)Nw=N;n=0:N-1;alpha=(N-1)/2;m=n-alpha+eps;%理想滤波器脉冲响应wc=(wp+ws)/2;%理想截频hd=sin(wc*m)./(pi*m);%采用汉明窗设计win=hamming(Nw);h=hd.*win';b=h;[H,f]=freqz(b,1,512,fs);%绘制滤波器的幅频响应图plot(f,20*log10(abs(H)))*label('频率');ylabel('幅度');title(’滤波器幅频特性’);gridon;4：滤波器的设计代码如下：[y,Fs,bits]=wavread('F:\a');%读出信号，采样率和采样位数。y=y(:,1);%取单声道作分析yl=length(y)；%求语音信号长度yy=fft(y,yl);%傅里叶变换t=[0:1/8000:4zeros(1,yl-32001)];m=0.07*sin(10000*pi*t);%产生噪声k=y+m;%参加噪声kl=length(k)；%求语音信号长度K=fft(k,kl);%傅里叶变换fb=800;fc=1000;As=100;Ap=1;fs=44000;wc=2*fc/fs;wb=2*fb/fs;[n,wn]=ellipord(wc,wb,Ap,As);[b,a]=ellip(n,Ap,As,wn);*=filter(b,a,k);s1=length(*);%求语音信号长度*=fft(*,s1);Figure(1);freqz(b,a,1000,8000);title('滤波器幅频相频图');figure(2);subplot(3,1,1);