总结 这里应该有
假设a b的尺寸分别为 ha *wa hb*wb
imfilt(a,b,'conv','same') = conv2(a,b,'same')=
现将a、b扩充成 a(ha+hb-1,wa+wb-1)=0; b(ha+hb-1,wa+wb-1)=0;
然后执行 ifft2(fft2(a).*fft2(b)) 之后在其中找到与a维数匹配的矩阵 应该是从floor(hb/2)+1:floor(hb/2)+ha * floor(wb/2)+1:floor(wb/2)+wa
下面是仿真程序
%% 先定义两个矩阵
ha = 9;% a ,b 矩阵 行列数定义wb = 9;
hb = 3;
wb = 3;
a = rand(ha,wa);% 创建两个随机矩阵
b = rand(hb,wb);
c = a;% 把ab放到其他矩阵里面保留 因为计算ifft2的时候需要给ab矩阵扩充维数
d = b;
e = conv2(a,b,'same'); % 计算 a与b的卷积 结果返回的值域a矩阵一致
f = imfilter(a,b,'conv','same'); %用imfilter计算ab卷积,返回结果与a矩阵维数一致
c(ha+hb-1,wa+wb-1)=0; %扩充维数
d(ha+hb-1,wa+wb-1)=0;
g = ifft2(fft2(c).*fft2(d)); %用ifft2计算a b卷积 不过这里返回的维数是 ha+hb-1,wa+wb-1
h = g(floor(hb/2)+1:floor(hb/2)+ha,floor(wb/2)+1:floor(wb/2)+wa); %找出与a匹配的矩阵
disp(e);
disp(f);
disp(h); %可以看到结果是一样的
结果如下:
0.2888 0.8512 0.5708 0.1797
1.0546 1.2345 1.1799 0.3251
0.4022 0.9156 1.1181 0.3801
0.7640 0.6469 1.6873 0.6974
0.7060 1.6673 1.5857 0.8648
0.4711 1.6789 1.4895 0.7817
0.7303 1.7992 1.1466 0.9305
0.4810 1.7491 0.9976 0.9398
0.2471 1.0858 0.9932 0.7670
0.2888 0.8512 0.5708 0.1797
1.0546 1.2345 1.1799 0.3251
0.4022 0.9156 1.1181 0.3801
0.7640 0.6469 1.6873 0.6974
0.7060 1.6673 1.5857 0.8648
0.4711 1.6789 1.4895 0.7817
0.7303 1.7992 1.1466 0.9305
0.4810 1.7491 0.9976 0.9398
0.2471 1.0858 0.9932 0.7670
0.2888 0.8512 0.5708 0.1797
1.0546 1.2345 1.1799 0.3251
0.4022 0.9156 1.1181 0.3801
0.7640 0.6469 1.6873 0.6974
0.7060 1.6673 1.5857 0.8648
0.4711 1.6789 1.4895 0.7817
0.7303 1.7992 1.1466 0.9305
0.4810 1.7491 0.9976 0.9398
0.2471 1.0858 0.9932 0.7670