ep312.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理(第三版)》第三章例3.1.2程序ep312.m
% 西安电子科技大学出版社出版  高西全  丁玉美合著 2008年5月
% DFT的MATLB计算
xn=[1 1 1 1];           %输入时域序列向量xn=R8(n)
Xk16=fft(xn,16);        %计算xn的16点DFT
Xk32=fft(xn,32);        %计算xn的32点DFT
%以下为绘图部分
k=0:15;wk=2*k/16;       %产生16点DFT对应的采样点频率(关于π归一化值)
subplot(3,2,1);stem(wk,abs(Xk16),'.'); %绘制16点DFT的幅频特性图
title('(a)16点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度')
subplot(3,2,5);stem(wk,angle(Xk16),'.'); %绘制16点DFT的相频特性图
line([0,2],[0,0]);title('(b)16点DFT的相频特性图')
xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])
k=0:31;wk=2*k/32;       %产生32点DFT对应的采样点频率(关于π归一化值)
subplot(3,2,2);stem(wk,abs(Xk32),'.'); %绘制32点DFT的幅频特性图
title('(c)32点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度')
subplot(3,2,6);stem(wk,angle(Xk32),'.'); %绘制32点DFT的相频特性图
line([0,2],[0,0]);title('(d)32点DFT的相频特性图');
xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])