evenodd.m

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

function [xe,x0,m]=evenodd(x,n)
%if any(imag(x)=0)
%   error('x is not a real sequency')
%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
      xe=0.5*(x+fliplr(x));
      x0=0.5*(x-fliplr(x));