firbfb.m

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

 % FIRBFB --- Biorthogonal Maximally decimation FIR Filter Bank
%              Computation(1D)
% y = firbfb(x,A,T,M)
% x, y are input and output respectively. A is N1-by-M matrix;
% whose columns are analysis filter bank. T is N2-by-M matrix,
% whose columns are synthesis filter bank. M is the channel
% numbers decimation and interpolation factor. Compare output
% with input, we may know if the system is PR.
%

function y = firbfb(x,A,T,M)

if M ~= size(A,2)
  error('M must be equal to the number of columns of A');
end

if M ~= size(T,2)
  error('M must be equal to the number of columns of T');
end

%A = antieye(size(A,1))*A;

for I = 1:M
 l = filter(A(:,I),1,x);
% l = conv(A(:,I),x);
 l = dsample(l,M);
 l = usample(l,M);
% l = conv(T(:,I),l);
 l = filter(T(:,I),1,l);
 l = l(:)';
 f(I,:) = l;
end 
y = sum(f);
dplot(y,'g');hold on;dplot(x,'r');