hc324.m

来自「数字信号处理教程：MATLAB释义与实现 书中实例码源」· M 代码 · 共 14 行

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.4程序hc324
% 信号序列的DTFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=[2,-1,1,1];nx=0:3;                % 给定x序列
K=64;dw=2*pi/K;                     % 确定频率分辨率dw
k=floor((-K/2+0.5):(K/2-0.5));      % 设定对称的频率向量下标
X=x*exp(j*dw*nx'*k);                % 求序列x的DTFT
subplot(2,1,1),plot(k*dw,abs(X))    % 绘图
xlabel('频率'),ylabel('幅频特性')
subplot(2,1,2),plot(k*dw,angle(X))
xlabel('频率'),ylabel('相频特性')
set(gcf,'color','w')                % 置图形背景色为白