したがって、私の質問は Infinity を返すことから NaN を返すことに変わりました。FFT が常に Infinity を返す場合、これが役立つ場合があります ( http://gerrybeauregard.wordpress.com/2011/04/01/an-fft-in-c/#comment-196)。したがって、C#は負の数の平方根を取得しようとしているため、NaNを返していると思いますが、以下のコードに基づいて数値を負にするべきではありません.ポジティブ)。ただし、返される数値は負です。複数の変数を使用して re * re と im * im を取得し、2 つの結果を加算するという非効率的な方法を試してみましたが、結果も負です。Math.Absも駄目でした。FFT クラスの作成者 (上記のリンクを参照) に連絡し、次の返信を待っています。以下のコードの一部は、以前に行った AS3 バージョンから取得しました。ここで回答を得る前にクラスの作成者から回答を得た場合は、それを投稿します。どんな洞察も最も役に立ちます。これまで私を助けてくれたすべての人に感謝します. 私は C# を使用する AS3 プログラマーです (C# の方がはるかに優れているため)。私はユニティを使用しています。
private const uint LOGN = 11; // Log2 FFT Length
private const uint N = 1 << (int)LOGN; // FFT Length
private const uint BUF_LEN = N; // Audio buffer length
public FFT2 fft; // FFT Object
private double[] tempIm = new double[N]; // Temporary Imaginary Number array
private double[] m_mag = new double[N/2]; // Magnitude array
private double[] m_win = new double[N]; // Hanning Window
private int fftCount = 0; // How many times the FFT has been performed
private double SCALE = (double)20/System.Math.Log(10); // used to convert magnitude from FFT to usable dB
private double MIN_VALUE = (double)System.Double.MinValue;
...
// Hanning analysis window
for (int i = 0; i < N; i++) // for i < 2048
m_win[i] = (4.0/N) * 0.5*(1-Mathf.Cos(2*Mathf.PI*i/N)); // Hanning Vector [1] = 1 / 4595889085.750801
…</p>
// Perform FFT
fft.run(tempRe, tempIm);
fftCount++;
// Convert from Decibel to Magnitude
for (int i = 0; i < N/2; i++) {
double re = tempRe[i]; // get the Real FFT Number at position i
double im = tempIm[i]; // get the Imaginary FFT Number at position i
m_mag[i] = Math.Sqrt(re * re + im * im); // Convert magnitude to decibels
m_mag[i] = SCALE * Math.Log(m_mag[i] + MIN_VALUE);
if (fftCount == 50 && i == 400) print ("dB @ 399: " + m_mag[399]);
}
前のコードは次を出力します。
dB @ 400: NaN
-5.56725062513722E+33
ありがとうございました!!