p_106_3_1.m

eda实例DSP上机程序

clc
clear
x1n=[ones(1,5),zeros(1,5)];
x2n=[ones(1,5),-ones(1,5)];


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');

%%%%线性卷积
yn=conv(x1n,x2n);
n=0:length(yn)-1;
figure;
subplot(1,1,1),stem(n,yn,'o');title('线性卷积');

%%%%%%%%%通过DFT计算卷积
N=10;
X1k=fft(x1n,N);
X2k=fft(x2n,N);
Yck=X1k.*X2k;
ycn=ifft(Yck);
n1=0:length(ycn)-1;
figure;
subplot(1,1,1),stem(n1,ycn,'o');title('循环卷积');