高通滤波matlab程序:clearall;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;%滤波器截止频率,滤波半径fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%计算高通滤波转换函数ifd=d0h=0+0.5;elseh=1+0.5;endy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('高通滤波');低通滤波matlab程序:clearall;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;%滤波器截止频率,滤波半径fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%计算低通滤波转换函数ifd=d0h=1;elseh=0;endy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('低通滤波');