polyp1.m

来自「里面囊括了基于matlab滤波器设计的各种.m文件」· M 代码 · 共 24 行

 % A = polyp1(x,M)
% Find the type-I polyphase components of signal x.
% x is a column vector. The decimation factor is M.
% The result will be put in the columns of matrix A. 
% Also see POLYP2

function A = polyp1(x,M)

x = x(:);
l = dsample(x,M);
l = l(:);
A(:,1) = l;
N = size(l,1);

for I = 1:M-1
  l = [zeros(M-I,1);x];
  l = dsample(l,M);
  l = l(2:size(l,1));
  if size(l,2)-N < 0
    A(:,I+1) = [l;zeros(N-size(l,1),1)];
  else
    A(:,I+1) = l;
  end 
end