基于matlab 的巴特沃斯低通滤波器的实现_第1页
基于matlab 的巴特沃斯低通滤波器的实现_第2页
基于matlab 的巴特沃斯低通滤波器的实现_第3页
基于matlab 的巴特沃斯低通滤波器的实现_第4页
基于matlab 的巴特沃斯低通滤波器的实现_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、基于matlab 的巴特沃斯低通滤波器的实现一 、课程设计的目的运用MATLAB实现巴特沃斯低通滤波器的设计以及相应结果的显示，另外还对多种低通滤波窗口进行了比较。二、课程设计的基本要求1）熟悉和掌握MATLAB 的基本应用技巧。2）学习和熟悉MATLAB相关函数的调用和应用。3）学会运用MATLAB实现低通滤波器的设计并进行结果显示。 三、双线性变换实现巴特沃斯低通滤波器的技术指标：1.采样频率10Hz。2.通带截止频率fp=0.2*pi Hz。3.阻带截止频率fs=0.3*pi Hz。4.通带衰减小于1dB，阻带衰减大于20dB四、使用双线性变换法由模拟滤波器原型设计数字滤波器程序代码：T

2、=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(Rp/10)-1)(1/(2*N);z,p,k = buttap(N); %获取零极点参数

3、p = p * OmegaC ;k = k*OmegaCN;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);stem(x,'.'); title ('输入信号');y=filter(b,a,x);

4、subplot(2,2,2);stem(y,'.') ;title('滤波之后的信号');figure(3) ;stem(y,'.') title（'输出的信号')）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),'-')legend('矩形窗','汉宁窗','布莱克曼窗','汉明窗');结果显示;200-20-40-60-80-100-120-1

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论