p_109_25.m

来自「eda实例DSP上机程序」· M 代码 · 共 24 行

clc
clear
x1n=ones(1,6);
x2n0=ones(1,8);
n=0:length(x2n0)-1;
x2n=n.*x2n0;
X1k1=fft(x1n,8);
X2k1=fft(x2n,8);
Yck1=X1k1.*X2k1;
ycn1=ifft(Yck1);
n1=0:length(ycn1)-1;
subplot(2,2,1),stem(n1,ycn1,'o');

X1k2=fft(x1n,16);
X2k2=fft(x2n,16);
Yck2=X1k2.*X2k2;
ycn2=ifft(Yck2);
n2=0:length(ycn2)-1;
subplot(2,2,2),stem(n2,ycn2,'o');

yn3=conv(x1n,x2n);
n3=0:length(yn3)-1;
subplot(2,2,4),stem(n3,yn3,'o');