fft에서 ifft를 계산하는 방법은 무엇입니까?

Question

저는 기본 주파수를 실시간으로 얻고 고역 통과 필터를 구현하기 위해 fft를했습니다.

이제 필터를 적용한 후 .wav 파일로 녹음 할 수있게하려고합니다.

먼저 fft를 반전시켜야하며 이것이 내 질문입니다. 이 작업을 수행하기위한 단계는 무엇입니까?

이 project에 정의 된 FFT를 사용합니다. 정확한 FFT 및 IFFT 정의에 따라 차이는 상수이다

using System; 
using System.Collections.Generic; 
using System.Text; 

namespace SoundLog 
{ 
    public class FourierTransform 
    { 
     static private int n, nu; 

     static private int BitReverse(int j) 
     { 
      int j2; 
      int j1 = j; 
      int k = 0; 
      for (int i = 1; i <= nu; i++) 
      { 
       j2 = j1/2; 
       k = 2 * k + j1 - 2 * j2; 
       j1 = j2; 
      } 
      return k; 
     } 

     static public double[] FFT(ref double[] x) 
     { 
      // Assume n is a power of 2 
      n = x.Length; 
      nu = (int)(Math.Log(n)/Math.Log(2)); 
      int n2 = n/2; 
      int nu1 = nu - 1; 
      double[] xre = new double[n]; 
      double[] xim = new double[n]; 
      double[] magnitude = new double[n2]; 
      double[] decibel = new double[n2]; 
      double tr, ti, p, arg, c, s; 
      for (int i = 0; i < n; i++) 
      { 
       xre[i] = x[i]; 
       xim[i] = 0.0f; 
      } 
      int k = 0; 
      for (int l = 1; l <= nu; l++) 
      { 
       while (k < n) 
       { 
        for (int i = 1; i <= n2; i++) 
        { 
         p = BitReverse(k >> nu1); 
         arg = 2 * (double)Math.PI * p/n; 
         c = (double)Math.Cos(arg); 
         s = (double)Math.Sin(arg); 
         tr = xre[k + n2] * c + xim[k + n2] * s; 
         ti = xim[k + n2] * c - xre[k + n2] * s; 
         xre[k + n2] = xre[k] - tr; 
         xim[k + n2] = xim[k] - ti; 
         xre[k] += tr; 
         xim[k] += ti; 
         k++; 
        } 
        k += n2; 
       } 
       k = 0; 
       nu1--; 
       n2 = n2/2; 
      } 
      k = 0; 
      int r; 
      while (k < n) 
      { 
       r = BitReverse(k); 
       if (r > k) 
       { 
        tr = xre[k]; 
        ti = xim[k]; 
        xre[k] = xre[r]; 
        xim[k] = xim[r]; 
        xre[r] = tr; 
        xim[r] = ti; 
       } 
       k++; 
      } 
      for (int i = 0; i < n/2; i++) 
       //magnitude[i] = (float)(Math.Sqrt((xre[i] * xre[i]) + (xim[i] * xim[i]))); 
       decibel[i] = 10.0 * Math.Log10((float)(Math.Sqrt((xre[i] * xre[i]) + (xim[i] * xim[i])))); 
      //return magnitude; 
      return decibel; 
     } 
    } 
} 

Answer 1

FFTW과 같이 주위에 정말 좋은 fft 구현이 많이 있는데이 중 하나를 사용하는 것이 좋습니다. 그들은 또한 ifft와 함께 온다. 구현 된대로 당신은 엄청나게 천천히 움직일 것입니다.

Answer 2

: 여기

는 코드입니다. 귀하의 경우 상수를 결정하는 가장 쉬운 방법은 아마도 시험입니다 & 오류입니다.