Skip to content

Instantly share code, notes, and snippets.

View npyoung's full-sized avatar

Noah Young npyoung

  • IndustrialNext
  • San Francisco
View GitHub Profile
@npyoung
npyoung / joy.py
Created February 24, 2022 23:45
Apply Joy Division to your timeseries data
import numpy as np
from scipy.signal.windows import hamming
class Joy:
def __truediv__(cls, arr):
T, N = arr.shape
spacer = 2 * np.arange(N)
t = np.arange(T)
window = hamming(T)
normed = (arr - arr.mean(0)) / (1e-3 + arr.std(0))