Discrepancy between convolution and filtering
显示 更早的评论
I've been trying to optimize the convolution process for my neural network to make it work on a big matrix. I figured out the logic, but the implementation is giving me different results from the conv2 function.
I've written a code to exemplify this:
% The image is a 28x28 picture, and the filter 5x5. The code i've written
% here is not the one i've optimized for matrix multiplication, but it
% shows the same logic where you sweep the image, multiply it by the
% filter and sum it.
image = data(:,:,1);
filter = rand(5, 5);
a = conv2(image, filter, 'valid');
b = zeros(24, 24);
for j=1:24
for i=1:24
b(i, j) = sum(filter .* image(i:i+4, j:j+4), 'all');
end
end
Although a and b are very similar:
They're usually off by a non-trivial amount:
Is there something i'm doing wrong here? Or maybe the conv2 algorithm does something to speed up the process that changes the result in a significant way?
Also, when the filter is all ones (like ones(5,5) instead of rand), a and b are identical.
采纳的回答
Image Analyst
2021-4-5
1 个投票
Your manual way is not flipping the kernel like conv2() does. So when the kernel is symmetric, the results will be the same and when the kernel is not symmetric, the results will be different. Try flipping the kernel in your "manual" way and you'll find they are the same.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Correlation and Convolution 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
