The inverse discrete Fourier transform Fn−1:Rn→Rn is given by the formula
x=Fn−1y, where xl=j=0∑n−1ωnjlyj,l=0,…,n−1
Here, ωn=exp(2πi/n) is the primitive root of unity of degree n and n=2p,p=1,2,…
(i) Show how to assemble x=F2m−1y in a small number of operations if the Fourier transforms of the even and odd parts of y,
x(E)=Fm−1y(E),x(O)=Fm−1y(O)
are already known.
(ii) Describe the Fast Fourier Transform (FFT) method for evaluating x, and draw a relevant diagram for n=8.
(iii) Find the costs of the FFT method for n=2p (only multiplications count).
(iv) For n=4 use the FFT method to find x=F4−1y when:
(a) y=(1,1,−1,−1),
(b) y=(1,−1,1,−1).