资源描述：
《基于matlab-的巴特沃斯低通滤波器的实现.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、基于matlab的巴特沃斯低通滤波器的实现一、课程设计的目的运用MATLAB实现巴特沃斯低通滤波器的设计以及相应结果的显示，另外还对多种低通滤波窗口进行了比较。二、课程设计的基本要求1）熟悉和掌握MATLAB的基本应用技巧。2）学习和熟悉MATLAB相关函数的调用和应用。3）学会运用MATLAB实现低通滤波器的设计并进行结果显示。三、双线性变换实现巴特沃斯低通滤波器的技术指标：1.采样频率10Hz。2.通带截止频率fp=0.2*piHz。3.阻带截止频率fs=0.3*piHz。4.通带衰减小于1dB，阻带衰减大于20dB四、使用双
2、线性变换法由模拟滤波器原型设计数字滤波器程序代码：T=0.1;FS=1/T;fp=0.2*pi;fs=0.3*pi;wp=fp/FS*2*pi;ws=fs/FS*2*pi;Rp=1;%通带衰减As=15;%阻带衰减OmegaP=(2/T)*tan(wp/2);%频率预计OmegaS=(2/T)*tan(ws/2);%频率预计%设计巴特沃斯低通滤波器原型N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS)));OmegaC=OmegaP/((10^
3、(Rp/10)-1)^(1/(2*N)));[z,p,k]=buttap(N);%获取零极点参数p=p*OmegaC;k=k*OmegaC^N;B=real(poly(z));b0=k;cs=k*B;ds=real(poly(p));[b,a]=bilinear(cs,ds,FS);%双线性变换figure(1);%绘制结果freqz(b,a,512,FS);%进行滤波验证figure(2);%绘制结果f1=50;f2=250;n=0:63;x=sin(2*pi*f1*n)+sin(2*pi*f2*n);subplot(2,2,1
4、);stem(x,'.');title('输入信号');y=filter(b,a,x);subplot(2,2,2);stem(y,'.');title('滤波之后的信号');figure(3);stem(y,'.')title（'输出的信号')）显示结果：1.滤波器的幅频与相频响应2.输入与滤波后信号显示3.输出的信号波形五、各种滤波窗口的比较：代码：wc=0.2*pi;N=33;tao=(N-1)/2;n=[0:(N-1)];m=n-tao+0.000001;hd=sin(wc*m)./(pi*m);wd1=boxcar(N)
5、';b1=hd.*wd1;wd2=hanning(N)';b2=hd.*wd2;wd3=blackman(N)';b3=hd.*wd3;wd4=hamming(N)';b4=hd.*wd4;[h1,w]=freqz(b1,1);[h2,w]=freqz(b2,1);[h3,w]=freqz(b3,1);[h4,w]=freqz(b4,1);plot(w,20*log10(abs(h1)),w,20*log10(abs(h2)),':',w,20*log10(abs(h3)),'-.',w,20*log10(abs(h4)),'--
6、')legend('矩形窗','汉宁窗','布莱克曼窗','汉明窗');结果显示;六、总结及心得:基于matlab的巴特沃斯低通滤波器的实现姓名：学号：S20100专业：日期：2010年06月08日
