% MTXIIRRESP calculate iir resp by matrix product%% A contains numerators coefficients (a row for each filter)% A contains denominators coefficients (a row for each filter)%% NOTE: A and B should have the same size, zero pad if necessary%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function rsp = mtxiirresp2(A,B,freq,fs) T = 1/fs; % build Z matrix Z = ones(size(A,2),numel(freq)); zz = exp(-2.*pi.*1i.*T.*freq); for jj=2:size(Z,1) Z(jj,:) = zz.^(jj-1); end % get numerator and denominator num = A*Z; den = B*Z; % get response rspm = num./den; if size(A,1)>1 rsp = sum(rspm,1); else rsp = rspm; endend