(整理)各类滤波器的MATLAB程序清单._第1页
(整理)各类滤波器的MATLAB程序清单._第2页
(整理)各类滤波器的MATLAB程序清单._第3页
(整理)各类滤波器的MATLAB程序清单._第4页
(整理)各类滤波器的MATLAB程序清单._第5页
已阅读5页，还剩7页未读， 继续免费阅读

下载本文档

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

文档简介

1、精品文档 各类滤波器的MATLAB程序 一、理想低通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); Hd=on es(size(IA); r=sqrt(f1.A2+f2.A2); Hd(r0.2)=0; Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure s

2、urf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); o o 二、理想高通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); Hd=on es(size(IA); r=sqrt(f1.A2+f2.A2); Hd(r0.2)=0; Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); la=real(ifft2(Ya); figure subplot(2,2,1),imshow(

3、ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); o o 三、Butterworth低通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); D=0.3; r=f1.A2+f2.A2; n=4; for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1/(tA n+1

4、); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 四、Butterworth高通滤波器 IA=imread(le na.bmp); f1,f2=freqspa

5、ce(size(IA),meshgrid); D=0.3; r=f1.A2+f2.A2; n=4; for i=1:size(IA,1) for j=1:size(IA,2) t=(D*D)/r(i,j); Hd(i,j)=1/(tA n+1); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf

6、(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 五、高斯低通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); D=100/size(IA,1); r=f1.A2+f2.A2; Hd=on es(size(IA); for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=exp(-t); end end Y=fft2(doub

7、le(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 六、高斯高通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(

8、IA),meshgrid); %D=100/size(IA,1); D=0.3; r=f1.A2+f2.A2; for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure

9、 surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 七、梯形低通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); %D=100/size(IA,1); D0=0.1; D1=0.4; r=sqrt(f1.A2+f2.A2); Hd=zeros(size(IA); Hd(r=D0 end end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd

10、; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 精品文档 0 0 八、梯形高通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); %D=100/size

11、(IA,1); D0=0.1; D1=0.4; r=sqrt(f1.A2+f2.A2); Hd=on es(size(IA); Hd(r=D0 end end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight

12、 in g,pho ng); 0 0 九、用其他方法编写的理想低通、理想高通、Butterworth低通、同态滤波程序 1、理想低通 i1=imread(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); d0=50; uO=floor(N1/2)+1; vO=floor(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); if d=d0 h=1; else h=0; end y(i,j)=g(i,j)*h；

13、end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); 2、理想高通 i1=imread(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); n=2; d0=10; uO=floor(N1/2)+1; vO=flo

14、or(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); if d=d0 h=0; else h=1; end y(i,j)=g(i,j)*h； end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); -.h Ji :“ .J1 3、Butterworth 低通 i1=imr

15、ead(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); n=2; d0=50; uO=floor(N1/2)+1; vO=floor(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); h=1/(1+(d/d0)A(2* n); y(i,j)=g(i,j)*h; end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); 4、同态滤波 l=rgb2gray(imread(fabric00.bmp); M,N=size(l); T=double(l); L=l

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论