40/CH8/EX8.18/Exa_8_18.sce


1
2
3
4
5
6
7
8
9

//The DFT from the matrix formulation
xn=[1;2;1;0];
w=exp(-%i*%pi/2);
for i=1:4
    for j=1:4
        WN(i,j)=w^((i-1)*(j-1));
    end
end
XDFT=WN*xn