function [xe, xo, m] = evenodd(x,n)% Real signal decomposition into even and odd parts% -------------------------------------------------% [xe, xo, m] = evenodd(x,n)%if any(imag(x) ~= 0)      error('x is not a real sequence')endm = -fliplr(n);m1 = min([m,n]); m2 = max([m,n]); m = m1:m2;nm = n(1)-m(1); n1 = 1:length(n);x1 = zeros(1,length(m));x1(n1+nm) = x; x = x1;xe = 0.5*(x + fliplr(x));xo = 0.5*(x - fliplr(x));