ASCEE/lasp
ASCEE/lasp
5
1
Fork 0
Code Issues 5 Pull Requests Actions Releases 14 Activity

Reversed poles of weighting filters to make filter stable

Browse Source
This commit is contained in:
Anne de Jong 2020-01-20 20:13:24 +01:00
parent b3197a3dee
commit b4ed3138d4
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lasp/filter/soundpressureweighting.py
View File

@ -138,7 +138,7 @@ class SPLFilterDesigner:
p1 = 2*np.pi*self.f1 p1 = 2*np.pi*self.f1
p4 = 2*np.pi*self.f4 p4 = 2*np.pi*self.f4
zeros_analog = [0,0] zeros_analog = [0,0]
poles_analog = [p1, p1, p4, p4] poles_analog = [-p1, -p1, -p4, -p4]
k_analog = p4**2/self._C_uncor(self.fr) k_analog = p4**2/self._C_uncor(self.fr)
z, p, k = bilinear_zpk(zeros_analog, poles_analog, k_analog, fs) z, p, k = bilinear_zpk(zeros_analog, poles_analog, k_analog, fs)
@ -159,7 +159,7 @@ class SPLFilterDesigner:
p3 = 2*np.pi*self.f3 p3 = 2*np.pi*self.f3
p4 = 2*np.pi*self.f4 p4 = 2*np.pi*self.f4
zeros_analog = [0,0,0,0] zeros_analog = [0,0,0,0]
poles_analog = [p1,p1, p2,p3,p4, p4] poles_analog = [-p1,-p1,-p2,-p3,-p4,-p4]
k_analog = p4**2/self._A_uncor(self.fr) k_analog = p4**2/self._A_uncor(self.fr)
z, p, k = bilinear_zpk(zeros_analog, poles_analog, k_analog, fs) z, p, k = bilinear_zpk(zeros_analog, poles_analog, k_analog, fs)