2017-11-12 1 views
0

私はPythonの高速フーリエ変換に問題があります。私は信号の低周波と高周波のパワーを抽出したいと思います。私は、次の例を見てきました:因子nによるフーリエ変換の除算

Y = np.fft.fft(signal)/n #Calculate FFT 
Y = Y[range(n/2)] #Return one side of the FFT 
lf = np.trapz(abs(Y[(frq>=0.04) & (frq<=0.15)])) 
hf = np.trapz(abs(Y[(frq>=0.16) & (frq<=0.5)])) #Do the same for 0.16-0.5Hz (HF) 

signalは、値の配列だけ(のは10000個の値を言ってみましょう)です。

ここでの問題は、私がnがすべきではないことです。 nはちょうどsignalの長さですか?

答えて

2

いいえ、nは必ずしも信号の長さである必要はありません。 nはFFTの長さであり、必要とされる周波数分解能および時間位置特定によって決定され、これはまた、定常性および/または信号の他の特性、ならびにウィンドウ処理またはゼロパディングが使用されるかどうか。

関連する問題