hc452.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.2程序hc452
% 连续信号频谱计算中有相邻谐振峰的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
T0=[0.6,0.15,0.15,0.15];   % 各次计算拟采用的T,编成向量T0
N0=[256,256,256,2048];     % 各次计算拟采用的N,编成向量N0
for r=1:4
    T=T0(r);N=N0(r); n=0:N-1;      % 根据计算顺序选用T及N
    D=2*pi/(N*T);                % 频率分辨率
    xa=exp(-0.01*n*T).*cos(2*n*T)+2*exp(-0.015*n*T).*cos(2.1*n*T); %样本序列xa
    Xa=T*fftshift(fft(xa));Xa(1)       % 求xa的FFT，移到对称位置
    k=floor(-(N-1)/2:(N-1)/2);      % 乃奎斯特频率下标D=2*pi/(N*T);
    set(gcf,'color','w')              % 置图形背景色为白
    subplot(2,2,r),plot(k*D,abs(Xa))    % 绘制幅频特性
    xlabel('\Omega'),ylabel('|Xa(\Omega)|')
    switch r                                % 用switch语句写子图标题
    case 1, title('(a)T=0.6,N=256')
    case 2, title('(b)T=0.15,N=256')
    case 3, title('(c)T=0.15,N=256,频率轴放大')
    case 4, title('(d)T=0.15,N=2048,频率轴放大')
    otherwise 
    end    
    if r>2 axis([1,3,0,60]), end            % 不同r定不同比例尺
end