example2_1.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 14 行

%例2-1 程序:example2_1
subplot(1,1,1)
n=0:100;x=(0.5).^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);
axis([-2,2,0,2])
xlabel('频率（以pi为单位）');ylabel('X的绝对值')
title('幅度')
subplot(2,1,2);plot(w/pi,angX);
axis([-2,2,-1,1])
xlabel('频率（以pi为单位）');ylabel('弧度/pi')
title('相角')