ex3p5.m

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

n=0:10;x=(0.9*exp(j*pi/3)).^n;
k=-200:200;w=(pi/100)*k;
X=x*(exp(-j*pi/100)).^(n'*k);
magX=abs(X);angX=angle(X);
subplot(2,1,1);plot(w/pi,magX);grid
xlabel('frequency in units of pi');ylabel('|X|')
title('Magnitude part')
subplot(2,1,2);plot(w/pi,angX/pi);grid
xlabel('frequency in units of pi');ylabel('radians/pi')
title('angel Part')