各类滤波器的MATLAB程序清单Word文档格式.docx

各类滤波器的MATLAB程序清单Word文档格式.docx

《各类滤波器的MATLAB程序清单Word文档格式.docx》由会员分享，可在线阅读，更多相关《各类滤波器的MATLAB程序清单Word文档格式.docx（15页珍藏版）》请在冰豆网上搜索。

Y=fft2（double（IA））;

Y=fftshift（Y）;

Ya=Y.*Hd;

Ya=ifftshift（Ya）;

Ia=ifft2（Ya）;

figure

subplot（2,2,1）,imshow（uint8（IA））;

subplot（2,2,2）,imshow（uint8（Ia））;

surf（Hd,'

Facecolor'

'

interp'

Edgecolor'

none'

Facelighting'

phong'

二、理想高通滤波器

Hd（r<

Ia=real（ifft2（Ya））;

三、Butterworth低通滤波器

D=0.3;

r=f1.^2+f2.^2;

n=4;

fori=1:

size（IA,1）

forj=1:

size（IA,2）

t=r（i,j）/（D*D）;

Hd（i,j）=1/（t^n+1）;

end

end

四、Butterworth高通滤波器

t=（D*D）/r（i,j）;

五、高斯低通滤波器

IB=imread（'

babarra.bmp'

D=100/size（IA,1）;

Hd（i,j）=exp（-t）;

六、高斯高通滤波器

%D=100/size（IA,1）;

Hd（i,j）=1-exp（-t）;

七、梯形低通滤波器

D0=0.1;

D1=0.4;

Hd=zeros（size（IA））;

D0）=1;

ifr（i,j）>

=D0&

r（i,j）<

=D1

Hd（i,j）=（D1-r（i,j））/（D1-D0）;

八、梯形高通滤波器

D1）=0;

Hd（i,j）=（D0-r（i,j））/（D0-D1）;

九、用其他方法编写的理想低通、理想高通、Butterworth低通、同态滤波程序

1、理想低通

i1=imread（'

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;

N1

N2

d=sqrt（（i-u0）^2+（j-v0）^2）;

ifd<

=d0

h=1;

else

h=0;

y（i,j）=g（i,j）*h;

y=ifftshift（y）;

E1=ifft2（y）;

E2=real（E1）;

subplot（2,2,1）,imshow（uint8（i1））;

subplot（2,2,2）,imshow（uint8（i2））;

subplot（2,2,3）,imshow（uint8（E2））;

2、理想高通

n=2;

d0=10;

elseh=1;

3、Butterworth低通

h=1/（1+（d/d0）^（2*n））;

4、同态滤波

I=rgb2gray（imread（'

fabric00.bmp'

））;

[M,N]=size（I）;

T=double（I）;

L=log（T）;

F=fft2（L）;

A=2;

B=0.3;

M

N

D（i,j）=（（i-M/2）^2+（j-N/2）^2）;

c=1.1;

%锐化参数

D0=max（M,N）;

H=（A-B）*（1-exp（c*（-D/（D0^2））））+B;

F=F.*H;

F=ifft2（F）;

Y=exp（F）;

subplot（1,2,1）,imshow（I）;

subplot（1,2,2）,imshow（uint8（real（Y）））;

一十、Gabor滤波器