Second-order (biquadratic) IIR digital filtering
Note: If either input to sosfilt is single precision, filtering is implemented using single-precision arithmetic. The output, y, is single precision.
sos represents the second-order section digital filter H(z)
by an L-by-6 matrix containing the coefficients of each second-order section in its rows.
If x is a matrix, sosfilt applies the filter to each column of x independently. The output y is a matrix of the same size, containing the filtered data corresponding to each column of x.
If x is a multidimensional array, sosfilt filters along the first nonsingleton dimension. The output y is a multidimensional array of the same size as x, containing the filtered data corresponding to each row and column of x.
The second order sections matrix,sos, the input signal,x, or both can be double or single precision. If at least one input is single precision, filtering is done with single precision arithmetic.
y = sosfilt(sos,x,dim) operates along the dimension dim.
 Orfanidis, S.J., Introduction to Signal Processing, Prentice-Hall, Englewood Cliffs, NJ, 1996.