freq2time: use czt to avoid taking conjugate x2

garrettj403 · Mar 22, 2021 · 7c34843 · 7c34843
1 parent 2406cd6
commit 7c34843
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/czt.py b/czt.py
@@ -284,16 +284,16 @@ def freq2time(f, X, t=None):
  dt = tspan / (nt - 1) # more accurate than t[1] - t[0]
 
  # Starting point
- A = np.exp(2j * np.pi * t[0] * df)
+ A = np.exp(-2j * np.pi * t[0] * df)
 
  # Step
- W = np.exp(-2j * np.pi * df * dt)
+ W = np.exp(2j * np.pi * df * dt)
 
  # Phase correction
  phase = np.exp(2j * np.pi * f[0] * t)
 
  # Time-domain transform
- time_data = iczt(X, N=nt, W=W, A=A) * phase
+ time_data = czt(X, nt, W=W, A=A) * phase / nf
 
  return t, time_data