modusa 0.4.27__tar.gz → 0.4.29__tar.gz

{modusa-0.4.27 → modusa-0.4.29}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: modusa
-Version: 0.4.27
+Version: 0.4.29
 Summary: A modular signal analysis python library.
 Author-Email: Ankit Anand <ankit0.anand0@gmail.com>
 License: MIT

{modusa-0.4.27 → modusa-0.4.29}/pyproject.toml

@@ -16,7 +16,7 @@ dependencies = [
 ]
 requires-python = ">=3.11"
 readme = "README.md"
-version = "0.4.27"
+version = "0.4.29"
 
 [project.license]
 text = "MIT"

{modusa-0.4.27 → modusa-0.4.29}/src/modusa/__init__.py

@@ -8,4 +8,6 @@ from modusa.tools import play, convert, record, save
 from modusa.tools import download
 from modusa.tools import load, load_ann
 
-__version__ = "0.4.27"
+from modusa.tools import synth_f0
+
+__version__ = "0.4.29"

{modusa-0.4.27 → modusa-0.4.29}/src/modusa/tools/__init__.py

@@ -9,4 +9,6 @@ from .ann_loader import load_ann
 from .audio_recorder import record
 
 from .plotter import Fig as fig
-from .plotter import dist_plot, hill_plot, plot
+from .plotter import dist_plot, hill_plot, plot
+
+from .synth import synth_f0

modusa-0.4.29/src/modusa/tools/synth.py

@@ -0,0 +1,55 @@
+#!/usr/bin/env python3
+
+#---------------------------------
+# Author: Ankit Anand
+# Date: 22/10/25
+# Email: ankit0.anand0@gmail.com
+#---------------------------------
+
+import numpy as np
+
+def synth_f0(f0, f0t, sr, nharm=0):
+	"""
+	Synthesize f0 contour so that you can
+	hear it back.
+
+	Parameters
+	----------
+	f0: ndarray
+		- Fundamental frequency (f0) contour in Hz.
+	f0t: ndarray
+		- Timestamps in seconds
+	sr: int
+		- Sampling rate in Hz for the synthesized audio.
+	nharm: int
+		- Number of harmonics
+		- Default: 0 => Only fundamental frequency (No harmonics)
+
+	Returns
+	-------
+	ndarray
+		- Syntesized audio.
+	sr
+		- Sampling rate of the synthesized audio
+	"""
+	
+	# Create new time axis
+	t = np.arange(0, f0t[-1], 1 / sr)
+	
+	# Interpolate the f0 to match the sampling time points
+	f0_interp = np.interp(t, f0t, f0)
+	
+	# Compute phase by integrating frequency over time
+	phase = 2 * np.pi * np.cumsum(f0_interp) / sr
+	
+	# Start with fundamental
+	y = np.sin(phase)
+	
+	# Add harmonics if requested
+	for n in range(2, nharm + 2):  # from 2nd to (nharm+1)th harmonic
+		y += np.sin(n * phase) / n**2  # dividing by n to reduce harmonic amplitude
+	
+	# Normalize output to avoid clipping
+	y /= np.max(np.abs(y))
+	
+	return y, sr