手動で逆FFTを計算する

Question

fft - 出力の逆FFTを手動で計算しようとしています。私は最初にfftを使ってデータセットのFFTを計算する次のスクリプトを使用しています。私は手動で逆FFTを見つけようとしますが、それはifftの結果と似ていません。

あなたは自分のエラーを発見できますか？私は単に

data = [ 
    -0.0005 
    -0.0004 
    -0.0003 
    -0.0002 
    -0.0001 
    -0.0000 
    0.0001 
    0.0001 
    0.0001 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0002 
    0.0003 
    0.0004 
    0.0005 
    0.0006 
    0.0007 
    0.0009 
    0.0010 
    0.0011 
    0.0011 
    0.0012 
    0.0011 
    0.0011 
    0.0011 
    0.0010 
    0.0011]; 


delta = 0.0125; 
fs = 1/delta; 
x  = (0:1:length(data)-1)/fs; 

X=fft(data); 

%find fft 
N=length(data); 
ws = 2*pi/N; 
wnorm = -pi:ws:pi; 
wnorm = wnorm(1:length(x)); 
w = wnorm*fs; 
figure(2) 
plot(w/(2*pi),abs(fftshift(X))) 


%find inverse fft manually 
for m=1:length(X) 
    for k=1:length(data) 
     X_real(m) = X(k)*exp(i*k*ws*(m-1)); 
    end 
end 

figure(3) 
plot(1:length(data), abs(X_real), 1:length(data), ifft(X)) 

Answer 1

https://en.wikipedia.org/wiki/Fast_Fourier_transform#Definition_and_speedは、以下のようなループのためにあなたを変更してください、FFTの標準逆式がここで提示し使用しています。

for m=1:length(X) 
    for k=1:length(data) 
     temp(k) = X(k)*exp(i*(m-1)*ws*(k-1)); 
    end 
    X_real(m)=(1/N)*sum(temp); 
end 

figure(3) 
plot(1:length(data), real(X_real)) 

ifftの方程式は、hereで見つけることができます。

あなたは2つのことを逃しました。

1つは正規化、もう1つは加算です。