From bbeb6795abeea51514b478b28b6379c4cd69dd91 Mon Sep 17 00:00:00 2001
From: millen1m <maxim.millen@gmail.com>
Date: Fri, 9 Nov 2018 18:56:04 +0000
Subject: [PATCH] added bandwidth calculator.

---
 eqsig/__about__.py |  2 +-
 eqsig/measures.py  | 18 ++++++++++++++++++
 2 files changed, 19 insertions(+), 1 deletion(-)

diff --git a/eqsig/__about__.py b/eqsig/__about__.py
index b2a26e1..466f5a7 100644
--- a/eqsig/__about__.py
+++ b/eqsig/__about__.py
@@ -1 +1 @@
-__version__ = "0.5.28"
+__version__ = "0.5.29"
diff --git a/eqsig/measures.py b/eqsig/measures.py
index 4baccaf..5c52832 100644
--- a/eqsig/measures.py
+++ b/eqsig/measures.py
@@ -181,3 +181,21 @@ def max_fa_period(asig):
     max_index = np.argmax(asig.fa_spectrum)
     max_period = 1. / asig.fa_frequencies[max_index]
     return max_period
+
+
+def calc_bandwidth_f_min(asig, ratio=0.707):
+    fas1_smooth = asig.smooth_fa_spectrum
+    max_fas1 = max(fas1_smooth)
+    lim_fas = max_fas1 * ratio
+    ind2 = np.where(fas1_smooth > lim_fas)
+    min_freq = asig.smooth_fa_frequencies[ind2[0][0]]
+    return min_freq
+
+
+def calc_bandwidth_f_max(asig, ratio=0.707):
+    fas1_smooth = asig.smooth_fa_spectrum
+    max_fas1 = max(fas1_smooth)
+    lim_fas = max_fas1 * ratio
+    ind2 = np.where(fas1_smooth > lim_fas)
+    max_freq = asig.smooth_fa_frequencies[ind2[0][-1]]
+    return max_freq