IIR filters code & examples¶

Table of Contents¶

24th-order Bessel filter, lowpass 1khz
12th-order Bessel filter, lowpass 1khz
6th-order Bessel filter, lowpass 1khz
24th-order Butterworth filter, lowpass 1khz
12th-order Butterworth filter, lowpass 1khz
12th-order Butterworth filter, highpass 1khz
12th-order Butterworth filter, bandpass
12th-order Butterworth filter, bandstop
4th-order Bessel filter, bandpass
8th-order Chebyshev type I filter, lowpass
8th-order Chebyshev type II filter, lowpass

24th-order Bessel filter, lowpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(bessel<fbase>(24), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("bessel_lowpass24", output, options + ", title='24th-order Bessel filter, lowpass 1khz'");

Result

bessel_lowpass24

12th-order Bessel filter, lowpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(bessel<fbase>(12), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("bessel_lowpass12", output, options + ", title='12th-order Bessel filter, lowpass 1khz'");

Result

bessel_lowpass12

6th-order Bessel filter, lowpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(bessel<fbase>(6), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("bessel_lowpass6", output, options + ", title='6th-order Bessel filter, lowpass 1khz'");

Result

bessel_lowpass6

24th-order Butterworth filter, lowpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(butterworth<fbase>(24), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_lowpass24", output,
              options + ", title='24th-order Butterworth filter, lowpass 1khz'");

Result

butterworth_lowpass24

12th-order Butterworth filter, lowpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(butterworth<fbase>(12), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_lowpass12", output,
              options + ", title='12th-order Butterworth filter, lowpass 1khz'");

Result

butterworth_lowpass12

12th-order Butterworth filter, highpass 1khz¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_highpass(butterworth<fbase>(12), 1000, 48000);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_highpass12", output,
              options + ", title='12th-order Butterworth filter, highpass 1khz'");

Result

butterworth_highpass12

12th-order Butterworth filter, bandpass¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_bandpass(butterworth<fbase>(12), 0.1, 0.2);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_bandpass12", output,
              options + ", title='12th-order Butterworth filter, bandpass'");

Result

butterworth_bandpass12

12th-order Butterworth filter, bandstop¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_bandstop(butterworth<fbase>(12), 0.1, 0.2);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_bandstop12", output,
              options + ", title='12th-order Butterworth filter, bandstop'");

Result

butterworth_bandstop12

4th-order Butterworth filter, bandpass¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_bandpass(butterworth<fbase>(4), 0.005, 0.9);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("butterworth_bandpass4", output, options + ", title='4th-order Butterworth filter, bandpass'");

Result

butterworth_bandpass4

8th-order Chebyshev type I filter, lowpass¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(chebyshev1<fbase>(8, 2), 0.09);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("chebyshev1_lowpass8", output,
              options + ", title='8th-order Chebyshev type I filter, lowpass'");

Result

chebyshev1_lowpass8

8th-order Chebyshev type II filter, lowpass¶

Code

const std::string options = "phaseresp=True, log_freq=True, freq_dB_lim=(-160, 10), padwidth=8192";
constexpr size_t maxorder = 32;
univector<fbase, 1024> output;

zpk<fbase> filt                       = iir_lowpass(chebyshev2<fbase>(8, 80), 0.09);
std::vector<biquad_params<fbase>> bqs = to_sos(filt);
output                                = biquad<maxorder>(bqs, unitimpulse());
plot_save("chebyshev2_lowpass8", output,
              options + ", title='8th-order Chebyshev type II filter, lowpass'");

Result

chebyshev2_lowpass8