fgp341.m

digital signal processing-MATLAB

%《数字信号处理教程——MATKAB释义与实现》第三章演示程序fgp341
% 四点序列周期延拓的离散时间傅立叶变换
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
disp('四点序列周期延拓的离散时间傅立叶变换')
disp('可选四个不同的周期延拓数进行计算')
x0=[2,-1,1,1];Nx=length(x0);			% x0是4点行向量
Nw=1000;dw=2*pi/Nw;                     % 把2π分为Nw份，求频率分辨率dw
k=floor((-Nw/2+0.5):(Nw/2-0.5));        % k*dw是正负对称的Nw点频率向量
for r=1:4
    K=input('延拓周期数K=（建议依次取1,10,100,500四种） ')
    nx=0:(K*Nx-1); x=x0(mod(nx,Nx)+1);  % 延拓后的时域信号x
    X=x*exp(j*dw*nx'*k);                % 用公式(3.2.30)求x的DTFT
    % X=freqz(x,1,k*w);                 % 用freqz函数求x的DTFT
    subplot(4,1,r),plot(k*dw,abs(X)),grid,shg		% 画图
    ylabel('abs(X(\omega))'),xlabel('\omega')       % 标注，注意希腊字母如何标注
end
set(gcf,'color','w')                    % 置图形背景色为白