built-in 2D convolution vs for-loop convolution differences?

Question

Matthew Lang 2013-5-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75073-built-in-2d-convolution-vs-for-loop-convolution-differences

评论： Burak Karadeniz 2021-3-20

Can someone tell me what's up with the following code? Why do I get different results implementing my own convolution w/ for loops vs the conv2 function?

x = 0:.1:pi;
A = rand(100,100);
M = sin(x).^2'*sin(x).^2; 
figure(1);
subplot(3,2,1);imagesc(x,x,M);colorbar('vert');title('Convolution kernel');
subplot(3,2,2);imagesc(A);colorbar('vert');title('Data matrix');
res = conv2(A,M,'valid');
subplot(3,2,3);imagesc(res);colorbar('vert');title('conv2(A,M)');
res2 = zeros(size(res));
for ii = 1:size(A,1)-size(M,1)+1
    for jj = 1:size(A,2)-size(M,2)+1
        res2(ii,jj) = sum(sum(M.*A(ii:ii+size(M,1)-1,jj:jj+size(M,2)-1)));
    end
end
subplot(3,2,4);imagesc(res2);colorbar('vert');title('for loop conv2(A,M)');
subplot(3,2,5);imagesc(res-res2);colorbar('vert');title('difference');
subplot(3,2,6);imagesc(100*(res-res2)./res);colorbar('vert');title('% error');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Iman Ansari 2013-5-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75073-built-in-2d-convolution-vs-for-loop-convolution-differences#answer_84790

在 MATLAB Online 中打开

Hi. You need to rotate your kernel 180 degree:

x = 0:.1:pi;
A = rand(100,100);
M = sin(x).^2'*sin(x).^2; 
figure(1);
subplot(3,2,1);imagesc(x,x,M);colorbar('vert');title('Convolution kernel');
subplot(3,2,2);imagesc(A);colorbar('vert');title('Data matrix');
res = conv2(A,M,'valid');
subplot(3,2,3);imagesc(res);colorbar('vert');title('conv2(A,M)');
res2 = zeros(size(res));
M=rot90(M,2);
for ii = 1:size(A,1)-size(M,1)+1
    for jj = 1:size(A,2)-size(M,2)+1
        res2(ii,jj) = sum(sum(M.*A(ii:ii+size(M,1)-1,jj:jj+size(M,2)-1)));
    end
end
subplot(3,2,4);imagesc(res2);colorbar('vert');title('for loop conv2(A,M)');
subplot(3,2,5);imagesc(res-res2);colorbar('vert');title('difference');
subplot(3,2,6);imagesc(100*(res-res2)./res);colorbar('vert');title('% error');