fft_2dbyfft_1d.m

通过计算晶体结构来确定其光谱分布

%Test the feasibility of 2D_FFT by 1D_FFT
clear;
clear global;
close all;
N=64;
data=rand(N)*1e2;
A=fft2(data);
A=fftshift(A);
hN=0;

hN=hN+1;
figure(hN);
surfl(real(A));
%colormap('gray');
title('Real A');

hN=hN+1;
figure(hN);
surfl(imag(A));
%colormap('gray');
title('Imagery A');

B=zeros(N);
C=zeros(N);
for m=1:N
    B(m,:)=fft(data(m,:));
end
for m=1:N
    C(:,m)=fft(B(:,m));
end

B=fftshift(C);
hN=hN+1;
figure(hN);
surfl(real(B));
%colormap('gray');
title('Real B');

hN=hN+1;
figure(hN);
surfl(imag(B));
title('Imagery B');

F=fft(data);
G=F.';
H=fft(G);
D=H.';
D=fftshift(D);

hN=hN+1;
figure(hN);
surfl(real(D));
%colormap('gray');
title('Real D');

hN=hN+1;
figure(hN);
surfl(imag(D));
%colormap('gray');
title('Imagery D');