Cascade two filters together to get a new filter. usage: [a,b] = abcascade(a1,b1,a2,b2) M Hewitson 03-04-07 $Id: abcascade.html,v 1.14 2008/03/31 10:27:32 hewitson Exp $
0001 function [ao,bo] = abcascade(a1,b1,a2,b2) 0002 0003 % Cascade two filters together to get a new filter. 0004 % 0005 % usage: [a,b] = abcascade(a1,b1,a2,b2) 0006 % 0007 % M Hewitson 03-04-07 0008 % 0009 % $Id: abcascade.html,v 1.14 2008/03/31 10:27:32 hewitson Exp $ 0010 % 0011 0012 n = length(a1); 0013 m = length(a2); 0014 N = m+n-1; 0015 0016 alpha = zeros(1, N); 0017 beta = zeros(1, N); 0018 0019 for i=0:n-1 0020 for j=0:m-1 0021 alpha(i+j+1) = alpha(i+j+1) + a1(i+1)*a2(j+1); 0022 beta(i+j+1) = beta(i+j+1) + b1(i+1)*b2(j+1); 0023 end 0024 end 0025 0026 ao = alpha; 0027 bo = beta;