hc953.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第九章例9.5.3程序hc953
% 窗函数对频谱分析的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=128;                      % 给定信号长度N
fr=[0.25,0.28];             % 给定信号频率fr
R0=[128,256];               % 给定DFT长度R
for i=1:2                   % 计算两组参数下的结果
    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,2*i-1),stem(k,abs(X),'.');grid on  % 画频谱图
    title(['(a) N=',num2str(N),',R=',num2str(R),',矩形窗DFT'])
    ylabel('|X(k)|')
    x1=x.*hamming(N)';                              % 对信号序列加哈明窗
    X1=fft(x1,R);k=0:R-1;                           % 求离散频谱
    subplot(3,2,2*i),stem(k,abs(X1),'.');grid on    % 画频谱图
    title(['(b) N=',num2str(N),',R=',num2str(R),',哈明窗DFT'])
    ylabel('|X(k)|')
end
subplot(3,2,5),plot(k/R,abs(X)),grid on             % 画X的连续频谱
title(['(e) N=',num2str(N),',R=',num2str(R),',矩形窗DTFT'])
xlabel('f=k/R'),ylabel('|X(f)|')
subplot(3,2,6),plot(k/R,abs(X1)),grid on            % 画X1的连续频谱
title(['(f) N=',num2str(N),',R=',num2str(R),',哈明窗DTFT'])
xlabel('f=k/R'),ylabel('|X(f)|')
set(gcf,'color','w')                                % 置图形背景色为白