96-04-01.m

来自「Matlab數位影音處理 影像方面的處理」· M 代码 · 共 18 行

clear all;
a=zeros(64);
b=ones(64);
f1=[a,b;b,a];
f=[f1,f1,f1,f1;f1,f1,f1,f1;f1,f1,f1,f1;f1,f1,f1,f1];
Fi=fft2(f1,512,512);
F=abs(Fi);
Fc=fftshift(F);
Fl=log(1+abs(Fc));
figure
subplot(2,2,1);
imshow(f,[])
subplot(2,2,2);
imshow(Fc,[])
subplot(2,2,3);
imshow(F,[])
subplot(2,2,4);
imshow(Fl,[])