marcel/pylot
3
0
Fork 0
Code Issues 7 Pull Requests 1 Releases Wiki Activity

Bugfix in fitSourceModel: Be sure, to set trial window aroun initial corner frequency within available data.

Browse Source
This commit is contained in:
Ludger Küperkoch 2016-10-27 16:06:45 +02:00
parent 41eebfbcda
commit 9e5b3f6916
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
pylot/core/analysis/magnitude.py
View File

@ -628,12 +628,13 @@ def fitSourceModel(f, S, fc0, iplot, verbosity=False):
STD = [] STD = []
# get window around initial corner frequency for trials # get window around initial corner frequency for trials
# left side of initial corner frequency # left side of initial corner frequency
fcstopl = fc0 - max(1, fc0 / 2) fcstopl = max(f[0], fc0 - max(1, fc0 / 2))
il = np.where(f <= fcstopl) il = np.where(f <= fcstopl)
il = il[0][np.size(il) - 1] il = il[0][np.size(il) - 1]
# right side of initial corner frequency # right side of initial corner frequency
fcstopr = fc0 + (fc0 / 2) fcstopr = min(fc0 + (fc0 / 2), f[len(f) - 1])
ir = np.where(f >= fcstopr) ir = np.where(f >= fcstopr)
print fcstopl, fcstopr
# check, if fcstopr is available # check, if fcstopr is available
if np.size(ir) == 0: if np.size(ir) == 0:
fcstopr = fc0 fcstopr = fc0