analytic wafeforms

analytic_wfm.py

from numpy import linspace
from math import pi, sqrt, sin

#Heavyside step function
H = lambda t: 1 if t > 0 else 0

# pure sine
ACV_A1 = lambda T, Hz = 50: 1000 * sqrt(2) * np.sin(2*pi*Hz * T)
# 
ACV_A2 = lambda T, Hz = 50: 1000 * sqrt(2) * np.sin(2*pi*Hz * T) + 500
# 
ACV_A3 = lambda T, Hz = 50: 1000 * sqrt(2) * (np.sin(2*pi*Hz * T) + 
                            0.05 * np.sin(2*pi*Hz * T * 4 + pi * 2 / 3))
#     
ACV_A4 = lambda T, Hz = 50:( 1000 * sqrt(2) * (np.sin(2*pi*Hz * T) + 
                            0.07 * np.sin(2*pi*Hz * T * 5 + pi * 22 / 18)))

# Realistic triangle
ACV_A5 = lambda T, Hz = 50:( 1000 * sqrt(2) * (np.sin(2*pi*Hz * T) + 
                            0.05 * np.sin(2*pi*Hz * T * 3 - pi) + 
                            0.05 * np.sin(2*pi*Hz * T * 5) + 
                            0.02 * np.sin(2*pi*Hz * T * 7 - pi) +
                            0.01 * np.sin(2*pi*Hz * T * 9)))
# 
ACV_A6 = lambda T, Hz = 50:( 1000 * sqrt(2) * (np.sin(2*pi*Hz * T) + 
                            0.02 * np.sin(2*pi*Hz * T * 3 - pi) + 
                            0.02 * np.sin(2*pi*Hz * T * 5) + 
                            0.0015 * np.sin(2*pi*Hz * T * 7 - pi) + 
                            0.009 * np.sin(2*pi*Hz * T * 9)))

#A7 & A8 convert so that a input of 16*pi corresponds to a input 0.25 in the current verion
ACV_A7 = lambda T: [1000*sqrt(2)*sin(100*pi*t)*(0.9*t/5*H(5-t) + H(t-5)*H(10-t)*(0.9+0.1*(t-5)/5)) for t in T]
ACV_A8 = lambda T: [1000*sqrt(2)*sin(t)*t/(10*pi)*H(10-t/(100*pi)) for t in T]



if __name__ == "__main__":
    #create 1 period triangle
    x = linspace(0,0.2*pi,8000)
    y = ACV_A5(x)