ex3p3.m

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

w=[0:1:500]*pi/500;
X=exp(j*w)./(exp(j*w)-0.5*ones(1,501));
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')