ABCASCADE Cascade two filters together to get a new filter. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DESCRIPTION: ABCASCADE Cascade two filters together to get a new filter. CALL: [a,b] = abcascade(a1,b1,a2,b2) VERSION: $Id: abcascade.m,v 1.2 2008/07/22 10:24:00 ingo Exp $ HISTORY: 03-04-2007 M Hewitson Creation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0001 % ABCASCADE Cascade two filters together to get a new filter. 0002 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0003 % 0004 % DESCRIPTION: ABCASCADE Cascade two filters together to get a new filter. 0005 % 0006 % CALL: [a,b] = abcascade(a1,b1,a2,b2) 0007 % 0008 % VERSION: $Id: abcascade.m,v 1.2 2008/07/22 10:24:00 ingo Exp $ 0009 % 0010 % HISTORY: 03-04-2007 M Hewitson 0011 % Creation 0012 % 0013 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0014 0015 function [ao,bo] = abcascade(a1,b1,a2,b2) 0016 0017 n = length(a1); 0018 m = length(a2); 0019 N = m+n-1; 0020 0021 alpha = zeros(1, N); 0022 beta = zeros(1, N); 0023 0024 for i=0:n-1 0025 for j=0:m-1 0026 alpha(i+j+1) = alpha(i+j+1) + a1(i+1)*a2(j+1); 0027 beta(i+j+1) = beta(i+j+1) + b1(i+1)*b2(j+1); 0028 end 0029 end 0030 0031 ao = alpha; 0032 bo = beta; 0033 end 0034