firofb.m

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

 % FIROFB --- Orthogonal Maximally decimation FIR Filter Bank
%             computation(1D)
% y = firofb(x,A,M)
% x, y are input and output respectively. A is N-by-M matrix,
% whose columns are synthesis filte bank, their time-reverse 
% are analysis 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 = firofb(x,A,M)

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

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