6-33.m

来自「MATLAB科学计算与工程分析源代码源程序5」· M 代码 · 共 10 行

clear
A=magic(4);            %生成4*4的魔幻矩阵
B=ones(3);              %生成3*3的全1矩阵
A(6,6)=0;               %将A用零补全为（4+3-1）*（4+3-1）维
B(6,6)=0;               %将B用零补全为（4+3-1）*（4+3-1）维
C=ifft2(fft2(A).*fft2(B));   %对A、B进行二维快速Fourier变换，并将结果相乘，对%乘积进行二维逆快速Fourier变换，得到卷积
A=magic(4);
B=ones(3);
D=conv2(A,B