From da87afc8d416d26cbfbd8310e0dd413198f5c700 Mon Sep 17 00:00:00 2001 From: Geoffrey Merck Date: Wed, 29 Apr 2020 20:41:41 +0200 Subject: [PATCH] Move helper script to separate folder --- .../FMGenerateFilterCoefficients.py | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) rename FMGenerateFilterCoefficients.py => Tools/FMGenerateFilterCoefficients.py (66%) diff --git a/FMGenerateFilterCoefficients.py b/Tools/FMGenerateFilterCoefficients.py similarity index 66% rename from FMGenerateFilterCoefficients.py rename to Tools/FMGenerateFilterCoefficients.py index 84ea62b..0580cdf 100644 --- a/FMGenerateFilterCoefficients.py +++ b/Tools/FMGenerateFilterCoefficients.py @@ -17,21 +17,22 @@ f2 = 2700 rp = 0.2 # scaling factor in bits, do not change ! -q = 15 +q = 0 # scaling factor as facor... scaling_factor = 2**q # let's generate a sequence of 2nd order IIR filters -#sos = signal.butter(2,[f1/fs*2,f2/fs*2],'pass',output='sos') -sos = signal.cheby1(3,rp,[f1/fs*2,f2/fs*2],'bandpass', output='sos') +#sos = signal.cheby1(3,rp,[f1, f2],'bandpass', output='sos', fs=fs) +#sos = signal.cheby1(1, rp, 2122, 'lowpass', output='sos', fs=fs) #deemphasis filter +sos = signal.cheby1(1, rp, 2122, 'highpass', output='sos', fs=fs) #deemphasis filter -sos = np.round((sos) * scaling_factor) +#sos = np.round((sos) * scaling_factor) # print coefficients for biquad in sos: for coeff in biquad: - print(int(coeff),",",sep="",end="") - #print((coeff),",",sep="",end="") + #print(int(coeff),",",sep="",end="") + print((coeff),",",sep="",end="") print("") # plot the frequency response