各類濾波器的MATLAB程序清單_第1頁
各類濾波器的MATLAB程序清單_第2頁
各類濾波器的MATLAB程序清單_第3頁
各類濾波器的MATLAB程序清單_第4頁
各類濾波器的MATLAB程序清單_第5頁
已閱讀5頁，還剩8頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

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 低通、同態(tài)濾波程序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、 同態(tài)濾波I=rgb2gray(imread('fabric00.b

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論