矩阵卷积(matlab-函数).docxVIP

function[hp]=juanji(f,g)

%此函数只为计算两任意二维矩阵的卷积

%命令格式：C=juanji(A,B)C为AB的卷积;

%若A为m*n矩阵，B为p*q矩阵，则C为(m+p-1)*(n+q-1)矩阵;

[m1n1]=size(f);

[m2n2]=size(g);

M=m1+m2-1;

N=n1+n2-1;

fp=zeros(M,N);

gp=zeros(M,N);

fp(1:m1,1:n1)=f;

gp(1:m2,1:n2)=g;

Fp=reshape(fp,[M*N1]);%fp必须转置reshape把原矩阵的各列连接为列矩阵

fori=1:M

G(:,:,i)=gallery(circul,gp(i,:));%生成循环矩阵

end

temp=gallery(circul,1:M);

fori=1:M

forj=1:M

a1=(i-1)*N+1;

a2=i*N;

b1=(j-1)*N+1;

b2=j*N;

c=temp(i,j);

Gb(a1:a2,b1:b2)=G(:,:,temp(i,j));

end

end

Hp=Gb*Fp;

hp=reshape(Hp,[MN]);

end