hc462.m

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

%《数字信号处理教程——MATLAB释义与实现》第四章例4.6.2程序hc462
%  连续信号频谱计算中含直流和多个正余弦分量的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% Program 4.6.2 
wm=62.83;T=pi/wm;       % 选择采样周期T
N0=[8,20];          % 设定两种N值      
for r=1:2                       % 作两次循环计算
    N=N0(r);D=2*pi/(N*T);    % 取N,并求出模拟频率分辨率
    k=[0:N-1]+eps;            % 求频率下标，偏移微量
    X=sin(0.275*k*D)./sin(0.025*k*D); % 算出频谱Xd
    x=ifft(X);                 % 求Xd的IDFT，得出x(n)
    subplot(1,2,r),stem([0:N-1]*T,x,'.')    % 绘图
    set(gcf,'color','w');              % 置图形背景色为白
end