circevod.m

来自「《数字信号处理及其MATLAB实现》程序集(2004)」· M 代码 · 共 12 行

function [xec, xoc] = circevod(x)
% 将信号分解为圆周偶和圆周奇两部分
% --------------------------------------------------------------
% [xec, xoc] = circecod(x)
%
if any(imag(x) ~= 0)
      error('x 非实数序列')
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));