From 45d22d3a13dd93a6fbf3a9f5ac5e5d2311b6953b Mon Sep 17 00:00:00 2001
From: "Michele.Maris" <michele.maris@inaf.it>
Date: Tue, 1 Aug 2023 11:40:00 +0200
Subject: [PATCH] a colored_noise.py
---
src/yapsut/colored_noise.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/yapsut/colored_noise.py b/src/yapsut/colored_noise.py
index 907c72c..055d7f4 100644
--- a/src/yapsut/colored_noise.py
+++ b/src/yapsut/colored_noise.py
@@ -92,10 +92,10 @@ Keywords:
#
self._freq=np.fft.rfftfreq(self.N)
#
- pl=self.alpha/2
+ pl=self.alpha
with np.errstate(divide='ignore'):
x=self.fknee/np.where(self._freq==0, np.inf , self._freq)
- self._S=(1+x)**pl
+ self._S=(1+x**pl)**0.5
#
if self._zero_policy=='0' :
self._S[0]=0.
@@ -104,7 +104,7 @@ Keywords:
elif self._zero_policy=='M' :
self._S[0]=self.wn_mean
elif self._zero_policy=='100x1' :
- self._S[0]=1 if fknee>ff[1] else (1+100*self.fknee/self._freq[1])**pl
+ self._S[0]=1 if fknee>ff[1] else (1+100*(self.fknee/self._freq[1])**pl)**0.5
else : # self._zero_policy=='I' :
# left things as they are
pass
--
GitLab