%《數字信號處理教程——MATLAB釋義與實現》第三章例3.6.4程序hc364
% 實序列的偶、奇分解和它的DFT的對稱性
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
n=0:9;x=10*(0.7).^n;N=length(x);    	% 給出x序列
[xe, xo] = circevod(x);                 % 調用循環奇偶分解函數求奇偶序列
Xe=fft(xe,N);Xo=fft(xo,N);              % 分別求奇偶分量的傅立葉變換
subplot(2,3,1),stem(n,xe,'.'),          % 以下畫圖
title('xe(n)')
subplot(2,3,2),stem(n,real(Xe),'.'),    % 繪循環偶序列頻譜幅特性
title('real(Xe(k))')
subplot(2,3,3),stem(n,imag(Xe),'.')     % 繪循環偶序列頻譜相特性
title('imag(Xe(k))')
subplot(2,3,4),stem(n,xo,'.'),
title('xo(n)')
subplot(2,3,5),stem(n,real(Xo),'.'),    % 繪循環奇序列頻譜幅特性
title('real(Xo(k))')
subplot(2,3,6),stem(n,imag(Xo),'.')     % 繪循環奇序列頻譜相特性
title('imag(Xo(k))')
set(gcf,'color','w')                    % 置圖形背景色為白