hc432.m

来自「这是数字信号处理——matlab释义与实现一书的源代码」· M 代码 · 共 16 行

%《数字信号处理教程——MATLAB释义与实现》第四章例4.3.2程序hc432
% 用FFT计算离散时间傅立叶变换DTFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
C=[20,1024];T=0.5;  % 设定两种补零后长度
for i=[1,2]
    N=C(i);D=2*pi/(N*T);                    % 根据N值确定频率分辨率
    x=[1 1 1 1 1 1 zeros(1,N-11) 1 1 1 1 1];% 根据N值确定序列
    X=fftshift(fft(x,N));                   % 作DFT，移到对称位置
    k=floor(-N/2+0.5:N/2-0.5);              % 确定对称的频率位置向量
    subplot(1,2,i),plot(k*D,X)              % 绘制幅频特性
    xlabel('频率'),ylabel('幅度')
    line([-10,10],[0,0])
end
set(gcf,'color','w')              % 置图形背景色为白