hc952a.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第九章例9.5.1a程序hc951a
% 信号长度和DFT长度对频谱分析的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=16;                   % 给定信号长度
fr=[0.22,0.34];         % 输入两个正弦频率
R0=[16,32,64,128];      % 给定DFT长度
for i=1:4               % 计算四组参数下的结果
    R=R0(i);
    n=0:N-1;
    x=0.5*sin(2*pi*n*fr(1))+sin(2*pi*n*fr(2));      % 列出信号序列
    X=fft(x,R);k=0:R-1;                             % 求离散频谱
    subplot(3,2,i),stem(k,abs(X),'.');grid on       % 绘图
    title(['(a).N=',num2str(N),',R=',num2str(R)])
    xlabel(''),ylabel('DFT幅度')
end
subplot(3,1,3),plot(k/R,abs(X)),grid on             % 求连续频谱
title('x(n)的DTFT幅度')
xlabel('数字频率f=k/R'),ylabel('DTFT幅度')
set(gcf,'color','w')                    % 置图形背景色为白