fig312.m

高西全2009年最新版<数字信号处理>源码程序

%图3.2.1绘图程序：fig321.m
M=24;
xn=[1 2 3 4 5 6];
n=-M/3:2*M/3-1;
x=zeros(1,M);
for l=1:6,
    x(8+l)=xn(l);
end
subplot(5,2,1)
stem(n,x,'.');axis([-8,16,0,7])
title('(a) x(n)');xlabel('n');ylabel('x(n)')
b=[xn 0 0];
ny=0:M-1;
y8=b(mod(ny,8)+1);
subplot(5,2,5)
stem(n,y8,'.');axis([-8,16,0,7])
title('(b) x(n)的周期延拓序列,N=8');xlabel('n');ylabel('x8(n)')
xn(1)=xn(1)+xn(5);xn(2)=xn(2)+xn(6);
y5=xn(mod(ny,4)+1);
subplot(5,2,9)
stem(n,y5,'.');axis([-8,16,0,10])
title('(c) x(n)的周期延拓序列,N=4');xlabel('n');ylabel('x4(n)')