各类滤波器的MATLAB程序清单_第1页
各类滤波器的MATLAB程序清单_第2页
各类滤波器的MATLAB程序清单_第3页
各类滤波器的MATLAB程序清单_第4页
各类滤波器的MATLAB程序清单_第5页
免费预览已结束，剩余8页可下载查看

下载本文档

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

文档简介

1、各类滤波器的MATLA醍序清单各类滤波器的 MATLAB 程序一、 理 想低通滤波器IA=imread('lena.bmp'); f1,f2=freqspace(size(IA),'meshgrid');Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r>0.2)=0;Y=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(

2、uint8(Ia);figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')二、理想高通滤波器IA=imread('lena.bmp'); f1,f2=freqspace(size(IA),'meshgrid');Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r<0.2)=0;Y=fft2(double(IA);Y=fftshift(Y

3、);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')三、 Butterworth 低通滤波器IA=imread('lena.bmp'); f1,f2=f

4、reqspace(size(IA),'meshgrid'); D=0.3;r=f1.2+f2.2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=r(i,j)/(D*D);Hd(i,j)=1/(tn+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(uint8(IA); subplot(2,2,2),imshow(uint8(Ia); figuresur

5、f(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')四、 Butterworth 高通滤波器IA=imread('lena.bmp');f1,f2=freqspace(size(IA),'meshgrid');D=0.3;r=f1.2+f2.2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=(D*D)/r(i,j);Hd(i,j)=1/(tn+1)

6、;endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figure五、高 斯低通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');D=100/size(IA,1);r=f1

7、.2+f2.2;Hd=ones(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);endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,'Facecolor','interp',&#

8、39;Edgecolor','none','Facelighting','phong')六、高 斯高通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D=0.3;r=f1.2+f2.2;for i=1:size(IA,1)for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t);en

9、dendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figure七、梯 形低通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1

10、;D1=0.4; r=sqrt(f1.2+f2.2); Hd=zeros(size(IA);Hd(r<D0)=1;for i=1:size(IA,1)for j=1:size(IA,2) if r(i,j)>=D0 & r(i,j)<=D1 Hd(i,j)=(D1-r(i,j)/(D1-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(uint8(IA); subpl

11、ot(2,2,2),imshow(uint8(Ia); figure surf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')八、梯 形高通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1;D1

12、=0.4;r=sqrt(f1.2+f2.2);Hd=ones(size(IA);Hd(r<D1)=0;for i=1:size(IA,1)for j=1:size(IA,2)if r(i,j)>=D0 & r(i,j)<=D1Hd(i,j)=(D0-r(i,j)/(D0-D1);endendendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(u

13、int8(Ia);figure九、用 其他方法编写的理想低通、理想高通、Butterworth 低通、同态滤波程序1、 理想低通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1);f=double(i2);k=fft2(f);g=fftshift(k);N1,N2=size(g);d0=50;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2 d=sqrt(i-u0)2+(j-v0)2);if d<=d0h=1;elseh=0;e

14、ndy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2,1),imshow(uint8(i1);subplot(2,2,2),imshow(uint8(i2);subplot(2,2,3),imshow(uint8(E2);2、 理想高通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f);g=fftshift(k);N1,N2=size(g)

15、;n=2;d0=10;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2d=sqrt(i-u0)2+(j-v0)2);if d<=d0h=0;else h=1;endy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figure subplot(2,2,1),imshow(uint8(i1); subplot(2,2,2),imshow(uint8(i2); subplot(2,2,3),imshow(uint8(E2);3、 Butterworth 低通i1

16、=imread('lena.bmp'); i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f); g=fftshift(k);N1,N2=size(g);n=2; d0=50; u0=floor(N1/2)+1; v0=floor(N2/2)+1;for i=1:N1for j=1:N2 d=sqrt(i-u0)2+(j-v0)2); h=1/(1+(d/d0)(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(uint8(i1); subplot(2,2,2),imshow(uint8(i2); subplot(2,2,3),imshow(uint8(E2);4、 同态滤波I=rgb2gray(imread('fabric00.b

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

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

最新文档

评论

0/150

提交评论