基于FFT实现自相关计算的C#代码单元测试失败排查

针对你遇到的FFT重构自相关后的测试失败问题，核心问题大概率出在归一化逻辑、Vector3分量处理、FFT卷积模式这三个方向，以下是具体排查点：

1. 错误的归一化导致AutoCorr(0)为1 #

你的测试显示归一化因子为1，说明代码中可能强制对自相关结果做了完全归一化（将所有结果除以lag=0的自相关值），但原TimeSeriesAnalysis类返回的是未归一化的原始自相关值（lag=0的结果等于序列的总能量，即所有Vector3元素的X/Y/Z分量平方和之和）。

常见错误代码示例： #

// 错误：强制归一化到[-1,1]范围
var autoCorr = InverseFFT(powerSpectrum);
autoCorr = autoCorr / autoCorr[0]; // 这一步会让lag=0的结果变成1

修复： #

如果原实现不需要归一化，直接移除归一化步骤；如果需要保留可选归一化，需增加开关参数，默认使用与原实现一致的未归一化模式：

var autoCorr = InverseFFT(powerSpectrum);
// 仅当需要归一化时才执行以下步骤，否则直接返回原始值
if (normalize)
{
    autoCorr = autoCorr / autoCorr[0];
}

2. 未正确拆分Vector3的分量进行独立计算 #

Vector3序列的自相关需要对X、Y、Z三个分量分别计算自相关，再将三个分量的结果组合成新的Vector3序列。如果你的代码将整个Vector3当作单一数值处理，或者没有拆分分量单独执行FFT流程，会导致结果完全不符合预期。

常见错误代码示例： #

// 错误：直接将Vector3序列传入FFT，未拆分分量
var fft = FFT.Forward(vector3Sequence.Select(v => v.X + v.Y + v.Z).ToArray());

修复： #

拆分三个分量分别计算自相关，再合并结果：

// 拆分X/Y/Z分量
var xVals = vector3Sequence.Select(v => v.X).ToArray();
var yVals = vector3Sequence.Select(v => v.Y).ToArray();
var zVals = vector3Sequence.Select(v => v.Z).ToArray();

// 分别对每个分量执行FFT自相关流程
var xCorr = ComputeAutoCorrViaFFT(xVals);
var yCorr = ComputeAutoCorrViaFFT(yVals);
var zCorr = ComputeAutoCorrViaFFT(zVals);

// 合并为Vector3序列
var result = xCorr.Zip(yCorr, zCorr, (x, y, z) => new Vector3(x, y, z)).ToList();

3. 未处理FFT的循环卷积混叠问题 #

用FFT计算自相关时，默认是循环卷积，而自相关需要的是线性卷积。如果不对序列进行补零扩展到2N-1长度（N为原序列长度），会导致循环卷积的混叠，使得lag较大的结果错误，甚至出现无效lag提示。

常见错误代码示例： #

// 错误：直接用原序列长度做FFT，未补零
var fft = FFT.Forward(sequence);
var powerSpectrum = fft.Select(v => v * Complex.Conjugate(v)).ToArray();
var autoCorr = FFT.Inverse(powerSpectrum);

修复： #

补零到2N-1长度后再执行FFT：

int n = sequence.Length;
int l = 2 * n - 1; // 线性卷积所需的最小长度

// 补零扩展序列
var padded = new Complex[l];
for (int i = 0; i < n; i++)
{
    padded[i] = sequence[i];
}

// 执行FFT流程
var fft = FFT.Forward(padded);
var powerSpectrum = fft.Select(v => v * Complex.Conjugate(v)).ToArray();
var autoCorr = FFT.Inverse(powerSpectrum);

// 取前n个值（对应lag 0到n-1）
var result = autoCorr.Take(n).Select(v => v.Real).ToArray();

4. Lag范围校验逻辑错误 #

无效lag提示说明代码对lag的输入范围判断有误，比如允许lag超过n-1（n为序列长度），或者计算后的自相关序列长度不足导致访问越界。需确保校验逻辑与原实现一致：

修复： #

public Vector3 AutoCorr(int lag)
{
    int n = _sequence.Count;
    if (lag < 0 || lag >= n)
    {
        throw new ArgumentOutOfRangeException(nameof(lag), "无效的lag值，范围应为0到序列长度-1");
    }
    // 后续逻辑
}

验证方向 #

• 先单独测试单个分量的自相关（比如只计算X分量），看lag=0的结果是否等于该分量的平方和，确认FFT流程的正确性。
• 对比原实现与FFT实现的中间结果（比如前几个lag的自相关值），定位差异点。

内容的提问来源于stack exchange，提问作者user366312