高通滤波matlab程序：clear all;I1=imread(E:matlab图片硬币.jpg);subplot(1,2,1);imshow(I1);title(原始图像);m=double(I1);f=fft2(m);f=fftshift(f);N1,N2=size(f); %返回矩阵的行和列n1=round(N1/2);n2=round(N2/2);n=2;d0=50; %滤波器截止频率,滤波半径for i=1:N1 for j=1:N2 d=sqrt(i-n1)2+(j-n2)2); %计算高通滤波转换函数 if d=d0 h=0+0.5; else h=1+0.5; end y(i,j)=h*f(i,j); endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A);subplot(1,2,2);imshow(B);title(高通滤波);低通滤波matlab程序：clear all;I1=imread(E:matlab图片硬币.jpg);subplot(1,2,1);imshow(I1);title(原始图像);m=double(I1);f=fft2(m);f=fftshift(f);N1,N2=size(f); %返回矩阵的行和列n1=round(N1/2);n2=round(N2/2);n=2;d0=50; %滤波器截止频率,滤波半径for i=1:N1 for j=1:N2 d=sqrt(i-n1)2+(j-n2)2); %计算低通滤波转换函数 if d=d0 h=1; else h=0; end y(i,j)=h*f(i,j); endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A);subplot(1,2,2);imshow(B);title(低通滤波);