circevod.m

来自「用Matlab实现重叠保留法计算循环卷积的程序实例」· M 代码 · 共 13 行

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