c105.m

来自「matlab例程MATLAB 6.x版的内建数据类型（Built-in data」· M 代码 · 共 17 行

R=32;
m=0:1:R-1;
f=(0.8).^m; 
h=(0.5).^m;
L=length(f)+length(h)-1;
k=0:1:L-1;
% Linear convolution: Direct approach
y1=conv(f,h);
subplot(2,1,2)
stem(k,y1)
% Linear Convolution:via DFT
FE=fft(f,L);
HE=fft(h,L);
y2=ifft(FE.*HE);
k=0:1:L-1;
subplot(2,1,1);stem(k,y2)
stem(k,y2)