ex3p4.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 24 行

n=-1:3;x=1:5;
k=0:500;
w=(pi/500)*k;
X=x*(exp(-j*pi/500)).^(n'*k);
magX=abs(X);
realX=real(X);
angX=angle(X);
imagX=imag(X);
subplot(2,2,1);plot(w/pi,magX);
title('magnitude part');xlabel('n')

subplot(2,2,2);plot(w/pi,realX);
title('real part');xlabel('n')

subplot(2,2,3);plot(w/pi,angX);
title('phase part');xlabel('n')

subplot(2,2,4);plot(w/pi,imagX);
title('imaginary part');xlabel('n')