p_106_3.m

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

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;
subplot(2,2,3),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;
subplot(2,2,4),stem(n1,ycn,'o');title('循环卷积');