hc952.m

数字信号处理,对与初雪着很有用,注释非常详细,调试都通过了!

%《数字信号处理教程——MATLAB释义与实现》第九章例9.5.1程序hc951
% 信号长度和DFT长度对频谱分析的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=input('信号长度N=(例9.5.2中取16)');           % 给定信号长度
R=input('DFT 长度R=(例9.5.2中取16,32,64,128)');           % 给定DFT长度
fr=input('[f1,f2]= (例9.5.2中取[0.22,0.34])');           % 输入两个正弦频率
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;                     % 求其DFT
stem(k,abs(X),'.');grid on              % 画图
title(['N=',num2str(N),',R=',num2str(R)])
set(gcf,'color','w')                    % 置图形背景色为白