なぜ私のC++ conv関数の出力はConv Matlab呼び出しと同じではありませんか？

Question

私は独自の畳み込み関数を実装し、その出力をMatlab conv関数のものと比較しました。

具体的には、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);を呼び出した後、conv([0.1, 0.23, 0.25, 0.18, 0.09], [0, 0, 1, 2, 3, 4, 5, 0, 0], 'same')コールの出力がoutputに書き込まれたものと同じになるようにします。ここで

は私のコード（それは私が_startと_stopブツを持っている理由である、信号がすでに埋められていることasssumes）

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];