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

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

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)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論