hc951.m

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

%《数字信号处理教程——MATLAB释义与实现》第九章例9.5.1程序hc951
% 正余弦信号的频谱分析
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=32; f=10;Fs=64;                       % 给定原始参数
n=0:N-1;g1=cos(2*pi*f*n/Fs);k=n;        % 生成有限长正弦信号序列
G1=fft(g1,N);                           % 求频谱
subplot(2,1,1),
stem(k,abs(G1),'.'),hold on             % 画出离散频谱样本
plot(N*[0:1023]/1024,abs(fft(g1,1024)),':')   % 画出频谱曲线
title('(a)N=32; f=10;Fs=64')
xlabel('k'),ylabel('|X(k)|')
N=32; f=11;Fs=64;                       % 改变原始参数
n=0:N-1;g1=cos(2*pi*f*n/Fs);k=n;        % 生成另一个有限长正弦信号序列
G1=fft(g1,N);                           % 求频谱
subplot(2,1,2),
stem(k,abs(G1),'.'),hold on             % 画出离散频谱样本
plot(N*[0:1023]/1024,abs(fft(g1,1024)),':')  % 画出频谱曲线
title('(b)N=32; f=11;Fs=64')
xlabel('k'),ylabel('|X(k)|')
set(gcf,'color','w')                 % 置图形背景色为白