Tuesday, 26 May 2015

Calculate IFFT using only FFT?



Say you start out with $A=\left\{3,−1,0,0\right\}$, run fft() and end up with $B=\left\{2,3+j,4,3-j\right\}$.


What specific steps would need to be done to obtain $A=\left\{3,−1,0,0\right\}$ again using some pre-calculation on $B$, then fft(), then maybe post-calculation steps before obtaining $A$? I have heard that this is possible but have not been able to work out how.



Answer



We want to do inverse DFT: $$ x(k) = \frac{1}{\sqrt{N}}\sum_n y(n) e^{j 2\pi kn/N}, $$ but can only do forward DFT. So let's try to make it "look" like a forward DFT. Note that we want a minus sign in the complex exponential. We can introduce it via complex conjugation: $$ \overline{x(k)} = \frac{1}{\sqrt{N}}\sum_n \overline{y(n)} e^{-j 2\pi kn/N}, $$ That is starting to look like a forward DFT now! So the recipe is:



  1. Complex conjugate the given sequence that we want to inverse DFT

  2. Calculate its forward DFT

  3. Calculate complex conjugate of the result.



That gives you the inverse DFT of the original sequence.


(Note that I cheated a little bit on the scaling factor to make it easier!)


No comments:

Post a Comment

readings - Appending 内 to a company name is read ない or うち?

For example, if I say マイクロソフト内のパートナーシップは強いです, is the 内 here read as うち or ない? Answer 「内」 in the form: 「Proper Noun + 内」 is always read 「ない...