circevod.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 8 行

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