hc363.m

某本matlab书籍中的代码

%《数字信号处理教程——MATLAB释义与实现》第三章例3.6.3程序hc363
% 实序列的偶、奇分解和它的DFT的对称性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%

N=8; x=[ones(1,6),zeros(1,2)];nx=0:N-1; % 给定实序列x
xr =x(mod(-nx,N)+1);                    % 求x的循环折叠xr
xe=0.5*(x+xr); xo=0.5*(x-xr);           % 分解为偶、奇序列xe,xo
X=fft(x,8), Xe=fft(xe,8), Xo=fft(xo,8)  % 求x,xe,xo的DFT
Error1=(max(abs(real(X)-Xe)))           % 岩盐X的实部与DFT[xe]的误差
Error2=(max(abs(j*imag(X)-Xo)))         % 岩盐X的虚部与DFT[xo]的误差
subplot(2,2,1),stem(nx,x)               % 绘初始序列图
xlabel('n'),ylabel('x(n)')
title('初始序列x')
subplot(2,2,3),stem(nx,xr)              % 绘循环折叠序列图
title('循环折叠序列xr')
xlabel('n'),ylabel('xr(n)')
subplot(2,2,2),stem(nx,xe)              % 绘循环偶序列图
title('偶序列xe')
xlabel('n'),ylabel('xe(n)')
subplot(2,2,4),stem(nx,xo)              % 绘循环奇叠序列图
title('奇序列xo')
xlabel('n'),ylabel('xo(n)')
set(gcf,'color','w')                    % 置图形背景色为白