我想写一个这个公式的矢量化版本:
总和从1
到j
,从1
到k
。
我用循环编写了这个版本:
function [ F] = calculation( A,B)
C=zeros(size(A)+2);
for j=1:size(C,1)
for k=1:size(C,2)
for p=1:j
for q=1:k
C(j,k)=C(j,k)+A(p,q)*B(j-p+1,k-q+1);
end
end
end
end
end
你能帮我把它矢量化吗?
正如您将注意到的,它是conv2
的公式。此函数将为语法F=conv2(A,B)
显示一个大小为F
的输出size(F)=size(A)+size(B)-1;
。
我的第二个问题是,是否有办法计算卷积但仅针对区域F(1:size(A,1)+2,1:size(A,2)+2)
并且不计算其他分量(以加快计算速度)?