ms11_5.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 11 行

xn1=[0,1,2,3,0,0];         %建立xn1序列
xn2=[1,1,1,0,0,0];         %建立xn2序列
N=length(xn1);      
n=0:N-1;k=0:N-1;
Xk1=xn1*(exp(-j*2*pi/N)).^(n'*k);
Xk2=xn2*(exp(-j*2*pi/N)).^(n'*k);
Yk=Xk1.*Xk2;
yn=Yk*(exp(j*2*pi/N)).^(n'*k)/N;
stem(n,yn,'k');