function [V,D] = DFTPolyEntries2(K,N); % [V,D] = DFTPolyEntries2(K,N); % % Calculates necessary indices for polyphase implementation of an % oversampled DFT filter bank, directly applying eqn (44) in [1]. % [1]. % % Input parameter: % K number of channels % N decimation ratio % % Output parameter: % V sample block with polyphase indices % D sample block with delays % % References: % [1] Cvetkovic Z and Vetterli M: "Oversampled Filter Banks". % To appear in IEEE Trans SP. % % St.Weiss, University of Strathclyde, 6.7.1997 M = lcm(K,N); L = M/N; J = M/K; B = N/J; % # of blocks K1 = []; N1 = []; for j = 1:J, K1 = [K1 diag(ones(1,K))]; end; for l = 1:L, N1 = [N1; diag(ones(1,N))]; end; V1 = K1*diag(0:M-1)*N1; V = V1(1:B:K,1:B:N); D = 1;