ButterParam - Fast parameters for a digital Butterworth filter
The calculation of filter parameters with BUTTER is expensive. Therefore ButterParam stores calculated parameters in a lookup-table to accelerate repeated calls with the same input.
x = rand(10000, 1); % test data
[B, A] = butter(3, 0.6);
y = filter(B, A, x);
The relation of the processing times is: BUTTER=68% / FILTER=32%!
With [B, A] = ButterParam(3, 0.6) we get: ButterParam=14% / FILTER = 86%, and the both calls together use 35% of the version with BUTTER.
Input and output equal the BUTTER command of the Signal-Processing-Toolbox
when used with 3 inputs and 2 outputs:
[B, A] = ButterParam(N, Wn, Pass)
N: Order of the filter.
Wn: Cutoff frequency, scalar or [1x2].
Pass: Optional, 'low', 'high', 'bandpass', 'stop'.
B, A: Coefficients for FILTER.
New parameters are calculated by BUTTER (Signal-Processing-Toolbox).
The list of calculated parameters can be saved to a MAT file for future access.
If all used parameters are found in this MAT file, the SPT is not needed to
get the Butterworth filter parameters.
Run the included unittest function uTest_ButterParam to check validity and speed.