circevod.m

来自「DSP2 using matlab v4」· M 代码 · 共 23 行

function [xec, xoc] = circevod(x)

% signal decomposition into circular-even and circular-odd parts

% --------------------------------------------------------------

% [xec, xoc] = circecod(x)

%

if any(imag(x) ~= 0)

      error('x is not a real sequence')

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));