hc461.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例4.6.1程序hc461
%  连续信号频谱计算中含直流和多个正余弦分量的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N0=[3,10];
for r=1:2
    T=1;N=N0(r);D=2*pi/N; 		% 频谱序列长度及频率分辨率
    kn=floor(-(N-1)/2:-1/2); 	% 负频率下标向量
    kp=floor(0:(N-1)/2); 		% 正频率下标向量
    w=[kp,kn]*D;     % 将负频率移到正频率的右方，形成新的频率排序
    Xd=2-exp(-j*w)+exp(-j*2*w)+exp(-j*3*w); % 按新的频率排序输入给定的频谱
    x=ifft(Xd);               	% 对循环对称的Xd求IFFT,得出循环对称的x 
    subplot(1,2,r),stem(T*[0:N-1],x)           % 绘图
    if r==1 title('N<序列长度，有泄漏')        % 用if语句写子图标题
    elseif r==2 title('N>序列长度，无泄漏'),end
    set(gcf,'color','w');              % 置图形背景色为白
    line([0,N]*T,[0,0])
end