sigevenodd.m

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

function [xeven,xodd,m]=sigevenodd(x,n)
%Real signal decomposition into even and odd parts
if (imag(x)~=0)
   error('x is not a real sequence')
end

m=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;
xeven=0.5*(x+fliplr(x));
xodd =0.5*(x-fliplr(x));