此傅立叶变换实现有什么问题
|
我正在尝试实现离散的傅立叶变换，但是没有用。我可能在某个地方写了一个错误，但是我还没有找到它。
根据以下公式：
该函数执行第一个循环，循环到X0-Xn-1 ...
public Complex[] Transform(Complex[] data, bool reverse)
{
var transformed = new Complex[data.Length];
for(var i = 0; i < data.Length; i++)
{
//I create a method to calculate a single value
transformed[i] = TransformSingle(i, data, reverse);
}
return transformed;
}
private Complex TransformSingle(int k, Complex[] data, bool reverse)
{
var sign = reverse ? 1.0: -1.0;
var transformed = Complex.Zero;
var argument = sign*2.0*Math.PI*k/data.Length;
for(var i = 0; i < data.Length; i++)
{
transformed += data[i]*Complex.FromPolarCoordinates(1, argument*i);
}
return transformed;
}
var sign = reverse ? 1.0: -1.0;
-1
-1
var argument = sign*2.0*Math.PI*k/data.Length;
transformed += data[i]*Complex.FromPolarCoordinates(1, argument*i);
没有找到相关结果
已邀请:
2 个回复
穗揣
什么时候应该更像是：除非您将其包裹在方法ѭ9中 更新： 我在AForge.NET Framework库中找到了以下代码，它显示了您的代码中未处理的其他Cos / Sin操作正在执行。可以在Sources \\ Math \\ FourierTransform.cs：DFT方法的完整上下文中找到此代码。for ( int i = 0; i < n; i++ ) { dst[i] = Complex.Zero; arg = - (int) direction * 2.0 * System.Math.PI * (double) i / (double) n; // sum source elements for ( int j = 0; j < n; j++ ) { cos = System.Math.Cos( j * arg ); sin = System.Math.Sin( j * arg ); dst[i].Re += ( data[j].Re * cos - data[j].Im * sin ); dst[i].Im += ( data[j].Re * sin + data[j].Im * cos ); } }祁澈盘哼颗
这只是归一化的约定，只有幅度会发生变化，因此您的结果也不会那么糟（如果您没有在逆变换中忘记1 / N）。 干杯