% test_fft.m also see test_fft2.m clear format compact diary test_fft_m.out disp('test_fft.m running') a=[0.0,0.7071,1.0,0.7071,0.0,-0.7071,-1.0,-0.7071]; disp('a=') disp(a) b=fft(a); disp('b=fft(a)=') disp(b) disp('abs(b)=') disp(abs(b)) disp('fftshift(abs(b))=') disp(fftshift(abs(b))) c=ifft(b); disp('c=ifft(b)=') disp(c) absc=abs(c); disp('abs(c)=') disp(absc) diary off % just to get output from disp to a file