example2_8.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 10 行

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');ylabel('amplitude');
subplot(2,1,2);plot(w/pi,angX);grid;
xlabel('Frequency');ylabel('angle');