matlab卷积代码,卷积的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')求得的结果：      自己编码求得的卷积变换结果：