fftsymmetrytest.m

来自「通过计算晶体结构来确定其光谱分布」· M 代码 · 共 17 行

clear;
close all;
N=100;
A=rand(1,N);
B=fft(A);
C=fftshift(B);
t=mod(N,2);
if t==0     %if N is even
    x=(-N/2:N/2-1);
else        %if N is odd
    t1=floor(N/2);
    x=(-t1:t1);
end
figure(1);
plot(x,real(C),'b',x,imag(C),'r');
grid on;
legend('Real part','Imagery part');