为什么我的 C++ conv 函数的输出与 conv Matlab 调用中的输出不同？

Question

因此，我实现了自己的卷积函数并将其输出与 Matlab conv 函数之一进行了比较。

具体来说，我希望 conv( [0.1, 0.23, 0.25, 0.18, 0.09], [0, 0, 1, 2, 3, 4, 5, 0, 0], 'same') 调用的输出与调用 conv({0, 0, 1, 2, 3, 4, 5, 0, 0}, 2, 6, {0.1, 0.23, 0.25, 0.18, 0.09}, 5, 5, output); 后写入 output 的输出相同。

这是我的代码（它假设信号已经被填充，这就是为什么我有 _start 和 _stop 东西）

void conv(double* signal, int conv_start, int conv_stop, double* kernel, int kernel_len, int output_len, double* output){   
    int halfKernel = floor(kernel_len/2.0);
    for (int i = 0; i<output_len; i++) output[i] = 0;
    for (int c = conv_start; c<=conv_stop; c++){
        for (int k = -halfKernel; k <=halfKernel; k++){ 
            output[c-conv_start] += kernel[k+halfKernel]*signal[c+k];
        }
    }
}

Matlab 函数的输出是： 1.0100 1.7700 2.6200 2.8700 2.2400

而我的是： 0.880000 1.630000 2.480000 2.790000 2.470000.

我也用上面的输入手工完成了计算，然后我得到了与我自己的实现相同的结果。那么，这是一个概念错误还是 Matlab 函数没有按我认为的方式执行？

Answer 1

首先，既然你写了这个函数，你的手算结果是一样的，这是有道理的。此外，它暗示这确实是一个概念错误。

在卷积中，内核应该围绕中间反射（或者如果你愿意的话翻转）。所以你可以这样做：

output[c-conv_start] += kernel[(kernel_len -1) - (k+halfKernel)]*signal[c+k];

为什么我的 C++ conv 函数的输出与 conv Matlab 调用中的输出不同？

Why is the output of my C++ conv-function not the same as in the conv Matlab call?

c++

algorithm

matlab

convolution