From 5bba29570edb05a58845236134ab5f917bb0a70d Mon Sep 17 00:00:00 2001
From: Anthony Wang
Date: Tue, 21 Mar 2023 22:56:34 -0400
Subject: Implement saw wave and honk waveform

---
 yue.py | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

(limited to 'yue.py')

diff --git a/yue.py b/yue.py
index 5ded276..7435d48 100644
--- a/yue.py
+++ b/yue.py
@@ -8,6 +8,13 @@ bitrate = 44100
 music = []
 
 
+def saw(x):
+    """
+    Sawtooth wave
+    """
+    return x / 2 / math.pi % 1
+
+
 def freq(octave, step):
     """
     Returns the frequency of a note
@@ -45,6 +52,16 @@ def seething(f, t):
     return Y
 
 
+def honk(f, t):
+    """
+    Returns the intensity of the "honk" waveform of frequency f sampled at time t
+    """
+    w = 2 * math.pi * f
+    Y = sum([(saw(i * w * t) / 3 + math.sin(i * w * t)) * math.exp(-0.002 * w * t) for i in range(1, 6)]) / 2
+    Y = math.copysign(abs(Y) ** (1 / 2), Y)
+    return Y
+
+
 def at(t):
     """
     Returns the total intensity of music sampled at time t
-- 
cgit v1.2.3-70-g09d2