circevod.m

来自「很多matlab的源代码」· M 代码 · 共 12 行

function [xec,xoc] = circevod(x)
% signal decomposition into circular-even and circular-odd parts
% ---------------------------------------------------------------
% [xec,xoc] = circevod(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));