convt_fft.m

来自「卷积实验代码 实现多种功能 包括线性卷积和圆周卷积的动态演示」· M 代码 · 共 11 行

function y=circonvt_fft(x1,x2)
M=length(x1);
N=length(x2);
L=M+N-1;
n=0:L-1;
x1=[x1 zeros(1,L-M)];
x2=[x2 zeros(1,L-N)];
y1=fft(x1);
y2=fft(x2);
[yk,n]=sigmult(y1,n,y2,n);
y=ifft(yk);