Apply a matrix of second order section digital filters to a timeseries.

applysos.m

function y = applysos(soscoef, x)
% Apply a matrix of SOS cofficients to some time series.
%
% Tobin Fricke
% 2012-02-07

if size(soscoef,2) ~= 6
  error('soscoef is not the right size');
end

y = x;

for ii=1:size(soscoef,1)
  y = filter(soscoef(ii,1:3), soscoef(ii,4:6), y);
end