@audio/spectral-ltas 1.0.0

package/README.md

@@ -0,0 +1,5 @@
+# @audio/spectral-ltas
+
+> LTAS — long-term average spectrum (Welch), critical-band smoothed; adaptive-EQ + match-EQ substrate
+
+Planned — not implemented yet. See the umbrella README for status and sources.

package/ltas.js

@@ -0,0 +1,25 @@
+// LTAS — long-term average spectrum via Welch's method: hann-windowed frames,
+// power-averaged, returned as RMS magnitude per bin. The adaptive-EQ / match-EQ substrate.
+
+import { fft } from 'fourier-transform'
+
+export default function ltas (data, { frameSize = 4096, hop = frameSize / 2 } = {}) {
+	let half = frameSize / 2
+	let win = new Float64Array(frameSize)
+	for (let i = 0; i < frameSize; i++) win[i] = 0.5 - 0.5 * Math.cos(2 * Math.PI * i / (frameSize - 1))
+
+	let acc = new Float64Array(half + 1)
+	let frames = 0
+	let buf = new Float64Array(frameSize)
+	for (let pos = 0; pos + frameSize <= data.length; pos += hop) {
+		for (let i = 0; i < frameSize; i++) buf[i] = data[pos + i] * win[i]
+		let [re, im] = fft(buf)
+		for (let k = 0; k <= half; k++) acc[k] += re[k] * re[k] + im[k] * im[k]
+		frames++
+	}
+	if (!frames) return new Float32Array(half + 1)
+
+	let out = new Float32Array(half + 1)
+	for (let k = 0; k <= half; k++) out[k] = Math.sqrt(acc[k] / frames)
+	return out
+}

package/package.json

@@ -0,0 +1,33 @@
+{
+  "name": "@audio/spectral-ltas",
+  "version": "1.0.0",
+  "description": "LTAS — long-term average spectrum via Welch's method (adaptive-EQ / match-EQ substrate)",
+  "type": "module",
+  "sideEffects": false,
+  "main": "ltas.js",
+  "exports": {
+    ".": "./ltas.js",
+    "./package.json": "./package.json"
+  },
+  "files": [
+    "ltas.js"
+  ],
+  "dependencies": {
+    "fourier-transform": "^2.2.0"
+  },
+  "keywords": [
+    "audio",
+    "dsp",
+    "spectral",
+    "ltas",
+    "welch"
+  ],
+  "license": "MIT",
+  "author": "Dmitry Iv. <dfcreative@gmail.com>",
+  "engines": {
+    "node": ">=18"
+  },
+  "publishConfig": {
+    "access": "public"
+  }
+}