circevod.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 9 行

function [xev,xod]=circevod(x)
if any(imag(x)~=0)
   error('The sequence is not real.')
end
N=length(x);
n=0:(N-1);
xev=0.5*(x+x(sigmod(-n,N)+1));
xod=0.5*(x-x(sigmod(-n,N)+1));