%《數(shù)字信號處理教程——MATLAB釋義與實現(xiàn)》第三章例3.6.3程序hc363
% 實序列的偶、奇分解和它的DFT的對稱性
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%

N=8; x=[ones(1,6),zeros(1,2)];nx=0:N-1; % 給定實序列x
xr =x(mod(-nx,N)+1);                    % 求x的循環(huán)折疊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)              % 繪循環(huán)折疊序列圖
title('循環(huán)折疊序列xr')
xlabel('n'),ylabel('xr(n)')
subplot(2,2,2),stem(nx,xe)              % 繪循環(huán)偶序列圖
title('偶序列xe')
xlabel('n'),ylabel('xe(n)')
subplot(2,2,4),stem(nx,xo)              % 繪循環(huán)奇疊序列圖
title('奇序列xo')
xlabel('n'),ylabel('xo(n)')
set(gcf,'color','w')                    % 置圖形背景色為白