卷积的Matlab代码实现

《数字图像处理》实验课上老师让尽量给出自己的源码，查阅资料初步了解了卷积的基本原理，所以针对实验课上的例子给出了自己的卷积变换代码。总之，卷积变换代码并不复杂，前提是了解了卷积的基本原理，这里有一个参考链接感觉讲的很棒，请查阅过后再阅读我的代码。

博客链接：https://blog.csdn.net/jinv5/article/details/52874880

%%拉普拉斯算子增强（卷积在下面）
clear all
load imdemos circuit%在图像处理工具箱中加载'circuit'图像
I = circuit;
I = double(I);
imshow(I, []);
h = [0 1 0;1 -4 1; 0 1 0];
J = conv2(I, h, 'same');%以h为卷积核对I进行卷积，返回与I相同大小的中心部分
K = I - J;
figure, imshow(J, []);
figure, imshow(K, []);
%%自己编码实现卷积运算
%先补0
myF = zeros(2, 4 + size(I, 2));
myF(3:(size(I, 1) + 2), 1:2) = 0;
myF(3:(size(I, 1) + 2),3:(size(I, 1) + 2)) = I;
myF(3:(size(I, 1) + 2),size(I, 1) + 3:size(I, 1) + 4) = 0;
myF(size(I, 1) + 3:size(I, 1) + 4, 1:size(I, 1) + 4) = 0;
%进行卷积 卷积顺序改变不会影响卷积的结果
P = ones(4 + size(I, 1), 4 + size(I, 2));
for i = 3:(2 + size(I, 1))
    for j = 3:(2 + size(I, 2))
        for m = 1:3
            for n = 1:3
                P(i, j) = P(i, j) + myF(i + m-2, j + n-2)*h(m, n);
            end
        end
    end
end
P = P(3:size(I, 1) + 2, 3:size(I, 2));

figure, imshow(P, []);

效果展示：

原图I：                                   利用matlab卷积函数conv2(I, h, 'same')求得的结果：      自己编码求得的卷积变换结果：