hc364.m

来自「这个是西安电子科技大学的陈怀琛教授编写的《数字信号处理教程——MATLAB释义与」· M 代码 · 共 21 行

%《数字信号处理教程——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')                    % 置图形背景色为白