p_109_24.m

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

clc
clear
x1n=[2 1 1 2];
x2n=[1 -1 -1 1];
%%%%线性卷积
yn=conv(x1n,x2n);
n=0:length(yn)-1
figure;
subplot(2,2,1),stem(n,yn,'o');
%%%%%%%%%通过DFT计算卷积
M1=length(x1n);
M2=length(x2n);
N=M1+M2-1;
X1k=fft(x1n,N);
X2k=fft(x2n,N);
Yck=X1k.*X2k;
ycn=ifft(Yck);
n1=0:length(ycn)-1
subplot(2,2,2),stem(n1,ycn,'o');