example3_5.m

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

%例3-5程序，example3_5
n=0:7;xn=(0.9).^n;N=8;k=[-N/2:N/2];
subplot(1,3,1);stem(n,xn)
title('xn')
Xk=dft(xn,N)
magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1)]);       %DFT幅度
subplot(1,3,2);stem(k,magXk)
title('Xk值')
phaXk=angle(Xk)*180/pi;np=[0:7];
subplot(1,3,3);plot(np,phaXk)
title('相位')