hc953a.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第九章例9.5.3程序hc953
% 窗函数对频谱分析的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=128;
fr=[0.25,0.28];
R0=[128,256];
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')                                % 置图形背景色为白
subplot(3,2,1),axis([29,39,0,70])                   % 设定比例尺
subplot(3,2,2),axis([29,39,0,40])
subplot(3,2,3),axis([58,78,0,70])
subplot(3,2,4),axis([58,78,0,40])
subplot(3,2,5),axis([0.24,0.3,0,60])
subplot(3,2,6),axis([0.24,0.3,0,40])