circevod.m

来自「matlab例子」· M 代码 · 共 12 行

function [xec,xoc]=circevod(x)
%将信号分解为循环偶和循环奇两部分
%----------------------------
% [xec,xoc]=circevod(x)
%
if any(imag(x)~=0)
   error('x非实数序列')
end
N=length(x);n=0:(N-1);
xec=0.5*(x+x(mod(-n,N)+1));
xoc=0.5*(x-x(mod(-n,N)+1));