hc455.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.5程序hc455
%  连续信号频谱计算中含直流和多个正余弦分量的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
for i= 1:2
    N=input('N= ');
    T=0.2;n=1:N;                    % 原始数据
    D=2*pi/(N*T);                   % 频率分辨率
    x=1+2*sin(0.2*pi*n*T)-3*cos(pi*n*T);        %  生成有限长的正弦序列
    X=T*fftshift(fft(x));X(1);      % 求x的FFT，移到对称位置
    k=floor(-(N-1)/2:(N-1)/2);      % 频率下标向量
    subplot(1,2,i),plot(k*D,abs(X))         % 绘图
    title('N= ,L=')
    xlabel('\Omega'),ylabel('|X(\Omega)|')
    v=axis;axis([-10,10,v(3),v(4)])  
end
set(gcf,'color','w');               % 置图形背景色为白