hc232.m

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

%《数字信号处理教程——MATLAB释义与实现》第二章例2.3.2程序hc232
% 实信号序列分解为奇部与偶部
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
n = [-2:10]; x = stepseq(-2,-2,10)-stepseq(7,-2,10);    % 给定的实序列
[x1,n1] = seqfold(x,n);		                            % 求出序列x(-n)
[xe,nxe] = seqadd(0.5*x,n,0.5*x1,n1);                   % 求分解出的偶序列
[xo,nxo] = seqadd(0.5*x,n,-0.5*x1,n1);                  % 求分解出的奇序列
subplot(2,2,1),stem(n,x),axis([-10,10,0,1.2])           % 画出实序列样本图
xlabel('n'),ylabel('x(n)'),title('矩形序列')
subplot(2,2,2),stem(nxe,xe),axis([-10,10,0,1.2])        % 画出偶序列样本图
xlabel('nxe'),ylabel('xe(n)'),title('偶部')
subplot(2,2,3),stem(n1,x1),axis([-10,10,0,1.2])         % 画出折叠后实序列样本图
xlabel('n1'),ylabel('x1(n)'),title('折叠后的矩形序列')
subplot(2,2,4),stem(nxo,xo),axis([-10,10,-0.5,0.5])
xlabel('nxe'),ylabel('xo(n)'),title('奇部')             % 画出偶序列样本图
set(gcf,'color','w')                                    % 置图形背景色为白