example1_8b.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

n=[-2:10];
x=[1:7,6:-1:1];
[x21,n21]=sigfold(x,n); %x(-n)
[x21,n21]=sigshift(x21,n21,3);   %x(-(n-3))
[x22,n22]=sigshift(x,n,2);
[x22,n22]=sigmult(x,n,x22,n22);
[x2,n2]=sigadd(x21,n21,x22,n22);
stem(n2,x2);
title('eg1.8b');
xlabel('n');ylabel('x2(n)');