p_107_14.m

eda实例DSP上机程序

clc
clear
x1n=ones(1,8);
x2n=ones(1,20);

figure;
n1=0:length(x1n)-1;
subplot(2,2,1),stem(n1,x1n,'o');
n2=0:length(x2n)-1;
subplot(2,2,2),stem(n2,x2n,'o');

%%%%线性卷积
yn3=conv(x1n,x2n);
n3=0:length(yn3)-1;
subplot(2,2,3),stem(n3,yn3,'o');title('线性卷积');

%%%%%%%%%通过DFT计算卷积
N=20;
X1k1=fft(x1n,N);
X2k1=fft(x2n,N);
Yck1=X1k1.*X2k1;
ycn1=ifft(Yck1);
n1=0:length(ycn1)-1;
subplot(2,2,4),stem(n1,ycn1,'o');title('循环卷积');