40/CH8/EX8.21/Exa_8_21.sce


1
2
3
4
5
6
7
8
9
10
11

//A four point decimation-in-time FFT algorithm
x=[1 2 1 0];
w=-%i;
xdft=[0 0 0 0];
for i=1:4
    for j=1:4
        xdft(i)=xdft(i)+x(j)*w^((i-1)*(j-1));
    end
end
XDFT=dft(x,-1);
xdft,XDFT