matlab代码傅里叶变换（matlab代码傅里叶变换）.docVIP

matlab代码傅里叶变换（matlab代码傅里叶变换） Fu Liye transform Matlab code: Img=imread (e:/4.jpg); Subplot (2,2,1), imshow (IMG); F=rgb2gray (IMG); Subplot (2,2,2), imshow (f); % two-dimensional Fu Liye transformation F=fft2 (f); FS=fftshift (F); % spectrum S=log (1+abs (FS)); Subplot (2,2,3); imshow (S, []); % of two-dimensional Fourier transform Fr=real (ifft2 (ifftshift (FS))); Ret=im2uint8 (mat2gray (FR)); Subplot (2,2,4), imshow (RET); Img=imread (e:/4.jpg); Subplot (2,2,1), imshow (IMG); F=rgb2gray (IMG); Subplot (2,2,2), imshow (f); The reference procedure: [i, lcmp]=imread (f:/tank.bmp);%======= read image display Imshow (I, lcmp); Ii=im2double (I);%===== image matrix type is converted to double (the image computation is not many integer) I1 = fft2 (II)%======; Fourier transform I2 =fftshift (I1);%====== will transform the image frequency four angle moving to the center (the original good processing portion in the four corner now mobile center, convenient behind) I3=log (ABS (i2));%===== display center frequency, and the log is to better display A program Figure, imshow (I3, []); Colormap (lcmp); Imwrite (I3, lcmp,f:/ffttank.bmp);%=== will be above I3 input to the ffttank file I5 = real (ifft2 (ifftshift (i2)));%=== frequency domain graph inverse transform and real part of airspace I6 = im2uint8 (mat2gray (i5));%=== the gray map Imwrite (I6, lcmp,f:/tank2.bmp,bmp);%=== with gray image and the original image to get the color reduction template Figure, imshow (I6); Title (inverse transform image); Program two RGB = imread (autumn.tif)%=== read image; I = rgb2gray (RGB);%== will turn it into a gray Imshow (I); J = DCT2 (I);%== two-dimensional DCT transform, the discrete cosine transform, Figure, imshow (log (ABS (J)), [])%== image shows after transformation, the energy is concentrated in the upper left corner Colormap (jet (64)), colorbar%==== established color template Figure; J (ABS (J) 10) = 0;%== I understand: because DCT transform as long as the preservation of the upper