矩阵卷积、矩阵相乘以及频域相乘之间的关系

通过查阅了网上很多的资料,在此做一个总结,有错误之处,还请评论指出,谢谢!

设矩阵A与矩阵B,其中矩阵B为卷积模板,B1是卷积模板翻转180度,FA是矩阵A在频域下的矩阵,FB是矩阵B在频域下的矩阵。
结论:矩阵A与矩阵B1相乘=矩阵A与矩阵B的卷积=矩阵FA与矩阵FB相乘。
(PS:上面说的两个相乘是不同的,具体可以从下面的讲述中可以清楚的明白)

1)验证:矩阵A与矩阵B1相乘=矩阵A与矩阵B的卷积
假设矩阵A为5×5的矩阵,矩阵B为3×3的矩阵。
在这里插入图片描述
在这里插入图片描述
B1是B旋转180度的矩阵,B1=rot90(rot90(B))
在这里插入图片描述
矩阵A与矩阵B1的相乘过程如下图所示:
在这里插入图片描述
依次计算,最后的结果为:
在这里插入图片描述
通过matlab来证明conv2(A,B)与上述操作的结果一致:
(关于im2col函数的作用可以参考:https://blog.csdn.net/dwyane12138/article/details/78449898)
为了让im2col()函数更明显,可以输出A2。

A=[1 2 4 5 6;6 8 9 1 3;6 8 9 2 4;7 3 5 7 5;1 5 8 9 3];%5*5
B=[3 4 5;6 7 8;3 2 1];%卷积核 3*3
C1=conv2(A,B,'full') %直接用卷积函数conv2的卷积结果
A1=padarray(A,[2 2],'both');%用0扩充A
A2=im2col(A1,[3 3]);%块数据向量化
B1=B(:);%变成一个列向量
B2=B1(end:-1:1);%倒序
B2=B2';%转置
C2=reshape(B2*A2,[7 7])%卷积转换为乘积并且转化为结果数据

结果如下图所示:
在这里插入图片描述
在这里插入图片描述
运算结果成功证明:矩阵A与矩阵B1相乘=矩阵A与矩阵B的卷积。

2)证明:矩阵A与矩阵B的卷积=矩阵FA与矩阵FB相乘

A=[1 2 4 5 6;6 8 9 1 3;6 8 9 2 4;7 3 5 7 5;1 5 8 9 3];%5*5
B=[3 4 5;6 7 8;3 2 1];%卷积核 3*3
C1=conv2(A,B,'full') %直接用卷积函数conv2的卷积结果
FA=fft2(A,7,7);
FB=fft2(B,7,7);
%FC=FA.*FB;  频域乘积,结果就是下面的for循环的操作,从这里可以看出,上面提到的两个相乘的方法是不同的。
for i=1:5
    for j=1:5
        FC(i,j)=FA(i,j)*FB(i,j);
    end
end
C3=ifft2(FC); %频域转换

结果入下图所示:
在这里插入图片描述
在这里插入图片描述
运算结果成功证明:矩阵A与矩阵B的卷积=矩阵FA与矩阵FB相乘。

综上所述:矩阵A与矩阵B1相乘=矩阵A与矩阵B的卷积=矩阵FA与矩阵FB相乘。

  • 4
    点赞
  • 32
    收藏
    觉得还不错? 一键收藏
  • 3
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值