hc341.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第三章例3.4.1程序hc341
% 周期序列的离散傅立叶级数
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
n=-15:16;                       % 时域自变量序列
x=[ones(1,4),zeros(1,4)];       % 输入序列
xtide=x(mod(n,8)+1);            % 延拓为周期序列
subplot(2,1,1),stem(n,xtide,'.')    % 画出时域序列
xlabel('n'),ylabel('xtide(n)')
k=[-8:15]+eps;                  % 频域自变量序列
X=exp(-j*3/8*pi*k).*sin(pi/2*k)./sin(pi/8*k); % 周期序列频谱的解析解
subplot(2,1,2),stem(k,abs(X),'.')  % 画出傅立叶级数系数
xlabel('k'),ylabel('abs(X(k))')
set(gcf,'color','w')              % 置图形背景色为白