audio-ml 1.0.0

package/dist/analysis/MFCCAnalyzer.js

@@ -0,0 +1,101 @@
+/**
+ * MFCCAnalyzer extracts Mel-Frequency Cepstral Coefficients (MFCCs) from PCM audio frames.
+ * MFCCs are widely used features in speech and audio processing, representing the short-term power spectrum of sound on a mel scale.
+ * This class uses FFT for spectral analysis, applies a mel filter bank, takes the log, and then applies a Discrete Cosine Transform (DCT).
+ * Usage: const mfcc = new MFCCAnalyzer({ sampleRate: 16000 }); mfcc.analyzeFrame(pcm)
+ */
+import FFT from 'fft.js';
+export class MFCCAnalyzer {
+    fft;
+    fftSize;
+    melBands;
+    numCoeffs;
+    sampleRate;
+    melFilters;
+    /**
+     * @param config - MFCC extraction configuration
+     */
+    constructor(config) {
+        this.fftSize = config.fftSize || 1024;
+        this.melBands = config.melBands || 26;
+        this.numCoeffs = config.numCoeffs || 13;
+        this.sampleRate = config.sampleRate;
+        this.fft = new FFT(this.fftSize);
+        this.melFilters = this.createMelFilterBank();
+    }
+    /**
+     * Extract MFCCs from a PCM frame.
+     * @param pcm - Input PCM audio frame (Float32Array, length = fftSize)
+     * @returns Array of MFCC coefficients (length = numCoeffs)
+     */
+    analyzeFrame(pcm) {
+        if (pcm.length !== this.fftSize)
+            throw new Error('PCM frame must be fftSize');
+        const spectrum = this.fft.createComplexArray();
+        this.fft.realTransform(spectrum, pcm);
+        this.fft.completeSpectrum(spectrum);
+        const mags = [];
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            const re = spectrum[2 * i];
+            const im = spectrum[2 * i + 1];
+            mags.push(Math.sqrt(re * re + im * im));
+        }
+        const melEnergies = this.melFilters.map(filter => filter.reduce((sum, w, k) => sum + w * mags[k], 0));
+        const logMel = melEnergies.map(e => Math.log(e + 1e-6));
+        return this.dct(logMel).slice(0, this.numCoeffs);
+    }
+    /**
+     * Create a mel filter bank matrix.
+     * @returns Array of mel filters (each filter is an array of weights for FFT bins)
+     */
+    createMelFilterBank() {
+        const filters = [];
+        const lowFreq = 0;
+        const highFreq = this.sampleRate / 2;
+        const melLow = this.hzToMel(lowFreq);
+        const melHigh = this.hzToMel(highFreq);
+        const melPoints = [];
+        for (let i = 0; i <= this.melBands + 2; i++) {
+            melPoints.push(melLow + (melHigh - melLow) * i / (this.melBands + 2));
+        }
+        const hzPoints = melPoints.map(mel => this.melToHz(mel));
+        const binPoints = hzPoints.map(hz => Math.floor((this.fftSize + 1) * hz / this.sampleRate));
+        for (let i = 0; i < this.melBands; i++) {
+            const filter = new Array(this.fftSize / 2).fill(0);
+            for (let j = binPoints[i]; j < binPoints[i + 1]; j++) {
+                filter[j] = (j - binPoints[i]) / (binPoints[i + 1] - binPoints[i]);
+            }
+            for (let j = binPoints[i + 1]; j < binPoints[i + 2]; j++) {
+                filter[j] = (binPoints[i + 2] - j) / (binPoints[i + 2] - binPoints[i + 1]);
+            }
+            filters.push(filter);
+        }
+        return filters;
+    }
+    /** Convert Hz to mel scale */
+    hzToMel(hz) {
+        return 2595 * Math.log10(1 + hz / 700);
+    }
+    /** Convert mel scale to Hz */
+    melToHz(mel) {
+        return 700 * (Math.pow(10, mel / 2595) - 1);
+    }
+    /**
+     * Discrete Cosine Transform (DCT-II) for MFCCs.
+     * @param vector - Input array (log-mel energies)
+     * @returns DCT coefficients
+     */
+    dct(vector) {
+        const N = vector.length;
+        const result = new Array(N).fill(0);
+        for (let k = 0; k < N; k++) {
+            let sum = 0;
+            for (let n = 0; n < N; n++) {
+                sum += vector[n] * Math.cos(Math.PI * k * (2 * n + 1) / (2 * N));
+            }
+            result[k] = sum;
+        }
+        return result;
+    }
+}
+//# sourceMappingURL=MFCCAnalyzer.js.map

package/dist/analysis/MFCCAnalyzer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"MFCCAnalyzer.js","sourceRoot":"","sources":["../../src/analysis/MFCCAnalyzer.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,GAAG,MAAM,QAAQ,CAAC;AAIzB,MAAM,OAAO,YAAY;IACf,GAAG,CAAM;IACT,OAAO,CAAS;IAChB,QAAQ,CAAS;IACjB,SAAS,CAAS;IAClB,UAAU,CAAS;IACnB,UAAU,CAAa;IAE/B;;OAEG;IACH,YAAY,MAAkB;QAC5B,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,IAAI,CAAC;QACtC,IAAI,CAAC,QAAQ,GAAG,MAAM,CAAC,QAAQ,IAAI,EAAE,CAAC;QACtC,IAAI,CAAC,SAAS,GAAG,MAAM,CAAC,SAAS,IAAI,EAAE,CAAC;QACxC,IAAI,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;QACpC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,mBAAmB,EAAE,CAAC;IAC/C,CAAC;IAED;;;;OAIG;IACH,YAAY,CAAC,GAAiB;QAC5B,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI,CAAC,OAAO;YAAE,MAAM,IAAI,KAAK,CAAC,2BAA2B,CAAC,CAAC;QAC9E,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC/C,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QACpC,MAAM,IAAI,GAAa,EAAE,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC1C,CAAC;QACD,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACtG,MAAM,MAAM,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QACxD,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACnD,CAAC;IAED;;;OAGG;IACK,mBAAmB;QACzB,MAAM,OAAO,GAAe,EAAE,CAAC;QAC/B,MAAM,OAAO,GAAG,CAAC,CAAC;QAClB,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC;QACrC,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QACrC,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,EAAE,CAAC;QACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,QAAQ,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,SAAS,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC;QACxE,CAAC;QACD,MAAM,QAAQ,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;QAC5F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnD,KAAK,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7E,CAAC;YACD,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACvB,CAAC;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,8BAA8B;IACtB,OAAO,CAAC,EAAU;QACxB,OAAO,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;IACzC,CAAC;IAED,8BAA8B;IACtB,OAAO,CAAC,GAAW;QACzB,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9C,CAAC;IAED;;;;OAIG;IACK,GAAG,CAAC,MAAgB;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACnE,CAAC;YACD,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAClB,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/dist/analysis/MelSpectrogramAnalyzer.d.ts

@@ -0,0 +1,33 @@
+export interface MelSpectrogramConfig {
+    sampleRate: number;
+    fftSize?: number;
+    melBands?: number;
+}
+export declare class MelSpectrogramAnalyzer {
+    private fft;
+    private fftSize;
+    private melBands;
+    private sampleRate;
+    private melFilters;
+    constructor(config: MelSpectrogramConfig);
+    /**
+     * Analyze a single frame of PCM samples and return mel spectrogram features.
+     * @param pcm - Input PCM audio frame (must be fftSize samples)
+     * @returns Array of log-mel energies
+     */
+    analyzeFrame(pcm: Float32Array): number[];
+    /**
+     * Create mel filter bank.
+     * @returns Array of mel filters (each filter is an array of weights for FFT bins)
+     */
+    private createMelFilterBank;
+    /**
+     * Convert Hz to mel scale.
+     */
+    private hzToMel;
+    /**
+     * Convert mel scale to Hz.
+     */
+    private melToHz;
+}
+//# sourceMappingURL=MelSpectrogramAnalyzer.d.ts.map

package/dist/analysis/MelSpectrogramAnalyzer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MelSpectrogramAnalyzer.d.ts","sourceRoot":"","sources":["../../src/analysis/MelSpectrogramAnalyzer.ts"],"names":[],"mappings":"AAQA,MAAM,WAAW,oBAAoB;IACnC,UAAU,EAAE,MAAM,CAAC;IACnB,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB,QAAQ,CAAC,EAAE,MAAM,CAAC;CACnB;AAED,qBAAa,sBAAsB;IACjC,OAAO,CAAC,GAAG,CAAM;IACjB,OAAO,CAAC,OAAO,CAAS;IACxB,OAAO,CAAC,QAAQ,CAAS;IACzB,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,UAAU,CAAa;gBAEnB,MAAM,EAAE,oBAAoB;IAQxC;;;;OAIG;IACH,YAAY,CAAC,GAAG,EAAE,YAAY,GAAG,MAAM,EAAE;IAyBzC;;;OAGG;IACH,OAAO,CAAC,mBAAmB;IA0B3B;;OAEG;IACH,OAAO,CAAC,OAAO;IAIf;;OAEG;IACH,OAAO,CAAC,OAAO;CAGhB"}

package/dist/analysis/MelSpectrogramAnalyzer.js

@@ -0,0 +1,86 @@
+/**
+ * MelSpectrogramAnalyzer computes mel spectrogram features from PCM audio frames.
+ * Mel spectrogram applies a mel filter bank to the FFT spectrum and returns log-mel energies.
+ * This is similar to MFCC but without the DCT step.
+ * Usage: const mel = new MelSpectrogramAnalyzer({ sampleRate: 16000 }); mel.analyzeFrame(pcm)
+ */
+import FFT from 'fft.js';
+export class MelSpectrogramAnalyzer {
+    fft;
+    fftSize;
+    melBands;
+    sampleRate;
+    melFilters;
+    constructor(config) {
+        this.fftSize = config.fftSize || 1024;
+        this.melBands = config.melBands || 40;
+        this.sampleRate = config.sampleRate;
+        this.fft = new FFT(this.fftSize);
+        this.melFilters = this.createMelFilterBank();
+    }
+    /**
+     * Analyze a single frame of PCM samples and return mel spectrogram features.
+     * @param pcm - Input PCM audio frame (must be fftSize samples)
+     * @returns Array of log-mel energies
+     */
+    analyzeFrame(pcm) {
+        if (pcm.length !== this.fftSize) {
+            throw new Error(`PCM frame must be fftSize (${this.fftSize}) samples`);
+        }
+        const spectrum = this.fft.createComplexArray();
+        this.fft.realTransform(spectrum, pcm);
+        this.fft.completeSpectrum(spectrum);
+        // Get magnitude spectrum
+        const mags = [];
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            const re = spectrum[2 * i];
+            const im = spectrum[2 * i + 1];
+            mags.push(Math.sqrt(re * re + im * im));
+        }
+        // Apply mel filter bank
+        const melEnergies = this.melFilters.map(filter => filter.reduce((sum, w, k) => sum + w * mags[k], 0));
+        // Log scaling
+        return melEnergies.map(e => Math.log(e + 1e-6));
+    }
+    /**
+     * Create mel filter bank.
+     * @returns Array of mel filters (each filter is an array of weights for FFT bins)
+     */
+    createMelFilterBank() {
+        const filters = [];
+        const lowFreq = 0;
+        const highFreq = this.sampleRate / 2;
+        const melLow = this.hzToMel(lowFreq);
+        const melHigh = this.hzToMel(highFreq);
+        const melPoints = [];
+        for (let i = 0; i <= this.melBands + 2; i++) {
+            melPoints.push(melLow + (melHigh - melLow) * i / (this.melBands + 2));
+        }
+        const hzPoints = melPoints.map(mel => this.melToHz(mel));
+        const binPoints = hzPoints.map(hz => Math.floor((this.fftSize + 1) * hz / this.sampleRate));
+        for (let i = 0; i < this.melBands; i++) {
+            const filter = new Array(this.fftSize / 2).fill(0);
+            for (let j = binPoints[i]; j < binPoints[i + 1]; j++) {
+                filter[j] = (j - binPoints[i]) / (binPoints[i + 1] - binPoints[i]);
+            }
+            for (let j = binPoints[i + 1]; j < binPoints[i + 2]; j++) {
+                filter[j] = (binPoints[i + 2] - j) / (binPoints[i + 2] - binPoints[i + 1]);
+            }
+            filters.push(filter);
+        }
+        return filters;
+    }
+    /**
+     * Convert Hz to mel scale.
+     */
+    hzToMel(hz) {
+        return 2595 * Math.log10(1 + hz / 700);
+    }
+    /**
+     * Convert mel scale to Hz.
+     */
+    melToHz(mel) {
+        return 700 * (Math.pow(10, mel / 2595) - 1);
+    }
+}
+//# sourceMappingURL=MelSpectrogramAnalyzer.js.map

package/dist/analysis/MelSpectrogramAnalyzer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"MelSpectrogramAnalyzer.js","sourceRoot":"","sources":["../../src/analysis/MelSpectrogramAnalyzer.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,GAAG,MAAM,QAAQ,CAAC;AAQzB,MAAM,OAAO,sBAAsB;IACzB,GAAG,CAAM;IACT,OAAO,CAAS;IAChB,QAAQ,CAAS;IACjB,UAAU,CAAS;IACnB,UAAU,CAAa;IAE/B,YAAY,MAA4B;QACtC,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,IAAI,CAAC;QACtC,IAAI,CAAC,QAAQ,GAAG,MAAM,CAAC,QAAQ,IAAI,EAAE,CAAC;QACtC,IAAI,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;QACpC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,mBAAmB,EAAE,CAAC;IAC/C,CAAC;IAED;;;;OAIG;IACH,YAAY,CAAC,GAAiB;QAC5B,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI,CAAC,OAAO,EAAE,CAAC;YAChC,MAAM,IAAI,KAAK,CAAC,8BAA8B,IAAI,CAAC,OAAO,WAAW,CAAC,CAAC;QACzE,CAAC;QACD,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC/C,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAEpC,yBAAyB;QACzB,MAAM,IAAI,GAAa,EAAE,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC1C,CAAC;QAED,wBAAwB;QACxB,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAC/C,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CACnD,CAAC;QAEF,cAAc;QACd,OAAO,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IAClD,CAAC;IAED;;;OAGG;IACK,mBAAmB;QACzB,MAAM,OAAO,GAAe,EAAE,CAAC;QAC/B,MAAM,OAAO,GAAG,CAAC,CAAC;QAClB,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC;QACrC,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QACrC,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,EAAE,CAAC;QACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,QAAQ,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,SAAS,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC;QACxE,CAAC;QACD,MAAM,QAAQ,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;QAE5F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnD,KAAK,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7E,CAAC;YACD,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACvB,CAAC;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IAED;;OAEG;IACK,OAAO,CAAC,EAAU;QACxB,OAAO,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;IACzC,CAAC;IAED;;OAEG;IACK,OAAO,CAAC,GAAW;QACzB,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9C,CAAC;CACF"}

package/dist/analysis/PLPAnalyzer.d.ts

@@ -0,0 +1,47 @@
+export interface PLPConfig {
+    sampleRate: number;
+    fftSize?: number;
+    order?: number;
+    numBands?: number;
+}
+export declare class PLPAnalyzer {
+    private fft;
+    private fftSize;
+    private order;
+    private numBands;
+    private sampleRate;
+    private barkFilters;
+    /**
+     * @param config - Configuration with sample rate, optional FFT size, LPC order, and number of critical bands
+     */
+    constructor(config: PLPConfig);
+    /**
+     * Extract PLP coefficients from a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns Array of PLP coefficients
+     */
+    analyzeFrame(pcm: Float32Array): number[];
+    /**
+     * Create a Bark scale filter bank.
+     * @returns Array of Bark filters
+     */
+    private createBarkFilterBank;
+    /**
+     * Convert Hz to Bark scale.
+     */
+    private hzToBark;
+    /**
+     * Equal loudness pre-emphasis weight.
+     */
+    private equalLoudnessWeight;
+    /**
+     * Convert Bark domain loudness to autocorrelation via inverse FFT.
+     * Interpolates Bark domain loudness back to frequency domain, then computes autocorrelation.
+     */
+    private barkToAutocorr;
+    /**
+     * Levinson-Durbin recursion for LPC coefficients.
+     */
+    private levinsonDurbin;
+}
+//# sourceMappingURL=PLPAnalyzer.d.ts.map

package/dist/analysis/PLPAnalyzer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PLPAnalyzer.d.ts","sourceRoot":"","sources":["../../src/analysis/PLPAnalyzer.ts"],"names":[],"mappings":"AAQA,MAAM,WAAW,SAAS;IACxB,UAAU,EAAE,MAAM,CAAC;IACnB,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB,KAAK,CAAC,EAAE,MAAM,CAAC;IACf,QAAQ,CAAC,EAAE,MAAM,CAAC;CACnB;AAED,qBAAa,WAAW;IACtB,OAAO,CAAC,GAAG,CAAM;IACjB,OAAO,CAAC,OAAO,CAAS;IACxB,OAAO,CAAC,KAAK,CAAS;IACtB,OAAO,CAAC,QAAQ,CAAS;IACzB,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,WAAW,CAAa;IAEhC;;OAEG;gBACS,MAAM,EAAE,SAAS;IAS7B;;;;OAIG;IACH,YAAY,CAAC,GAAG,EAAE,YAAY,GAAG,MAAM,EAAE;IAuCzC;;;OAGG;IACH,OAAO,CAAC,oBAAoB;IAkC5B;;OAEG;IACH,OAAO,CAAC,QAAQ;IAIhB;;OAEG;IACH,OAAO,CAAC,mBAAmB;IAM3B;;;OAGG;IACH,OAAO,CAAC,cAAc;IAkDtB;;OAEG;IACH,OAAO,CAAC,cAAc;CAqBvB"}

package/dist/analysis/PLPAnalyzer.js

@@ -0,0 +1,176 @@
+/**
+ * PLPAnalyzer extracts Perceptual Linear Prediction (PLP) coefficients from PCM audio frames.
+ * PLP is used in speech processing to model the auditory spectrum using perceptual principles.
+ * This implementation includes critical band analysis, equal loudness pre-emphasis, and LPC analysis.
+ * Usage: const plp = new PLPAnalyzer({ sampleRate: 16000 }); plp.analyzeFrame(pcm)
+ */
+import FFT from 'fft.js';
+export class PLPAnalyzer {
+    fft;
+    fftSize;
+    order;
+    numBands;
+    sampleRate;
+    barkFilters;
+    /**
+     * @param config - Configuration with sample rate, optional FFT size, LPC order, and number of critical bands
+     */
+    constructor(config) {
+        this.fftSize = config.fftSize || 512;
+        this.order = config.order || 12;
+        this.numBands = config.numBands || 18;
+        this.sampleRate = config.sampleRate;
+        this.fft = new FFT(this.fftSize);
+        this.barkFilters = this.createBarkFilterBank();
+    }
+    /**
+     * Extract PLP coefficients from a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns Array of PLP coefficients
+     */
+    analyzeFrame(pcm) {
+        // Ensure correct frame size
+        if (pcm.length !== this.fftSize) {
+            throw new Error(`PCM frame must be fftSize (${this.fftSize}) samples, got ${pcm.length}`);
+        }
+        // Step 1: Compute power spectrum
+        const spectrum = this.fft.createComplexArray();
+        this.fft.realTransform(spectrum, pcm);
+        this.fft.completeSpectrum(spectrum);
+        const powerSpectrum = [];
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            const re = spectrum[2 * i];
+            const im = spectrum[2 * i + 1];
+            powerSpectrum.push(re * re + im * im);
+        }
+        // Step 2: Critical band analysis (Bark scale)
+        const barkEnergies = this.barkFilters.map(filter => filter.reduce((sum, w, k) => sum + w * powerSpectrum[k], 0));
+        // Step 3: Equal loudness pre-emphasis
+        const preemphasized = barkEnergies.map((energy, i) => {
+            const bark = i * 1.0; // Approximate bark value
+            const equalLoudness = this.equalLoudnessWeight(bark);
+            return energy * equalLoudness;
+        });
+        // Step 4: Intensity-loudness conversion (cube root)
+        const loudness = preemphasized.map(e => Math.pow(e + 1e-12, 1.0 / 3.0));
+        // Step 5: Inverse FFT to get autocorrelation
+        const autocorr = this.barkToAutocorr(loudness);
+        // Step 6: LPC analysis using Levinson-Durbin
+        return this.levinsonDurbin(autocorr, this.order);
+    }
+    /**
+     * Create a Bark scale filter bank.
+     * @returns Array of Bark filters
+     */
+    createBarkFilterBank() {
+        const filters = [];
+        const nyquist = this.sampleRate / 2;
+        const maxBark = this.hzToBark(nyquist);
+        for (let band = 0; band < this.numBands; band++) {
+            const filter = new Array(this.fftSize / 2).fill(0);
+            const centerBark = (band + 0.5) * maxBark / this.numBands;
+            for (let i = 0; i < this.fftSize / 2; i++) {
+                const freq = (i * this.sampleRate) / this.fftSize;
+                const bark = this.hzToBark(freq);
+                const distance = Math.abs(bark - centerBark);
+                // Triangular filter shape
+                if (distance < 1.0) {
+                    filter[i] = 1.0 - distance;
+                }
+            }
+            // Normalize
+            const sum = filter.reduce((a, b) => a + b, 0);
+            if (sum > 0) {
+                for (let i = 0; i < filter.length; i++) {
+                    filter[i] /= sum;
+                }
+            }
+            filters.push(filter);
+        }
+        return filters;
+    }
+    /**
+     * Convert Hz to Bark scale.
+     */
+    hzToBark(hz) {
+        return 13 * Math.atan(0.00076 * hz) + 3.5 * Math.atan(Math.pow(hz / 7500, 2));
+    }
+    /**
+     * Equal loudness pre-emphasis weight.
+     */
+    equalLoudnessWeight(bark) {
+        // Simplified equal loudness curve
+        const w = 1.0 / (1.0 + Math.pow(bark / 15.0, 2));
+        return w;
+    }
+    /**
+     * Convert Bark domain loudness to autocorrelation via inverse FFT.
+     * Interpolates Bark domain loudness back to frequency domain, then computes autocorrelation.
+     */
+    barkToAutocorr(loudness) {
+        // Interpolate Bark domain loudness back to frequency domain
+        const nyquist = this.sampleRate / 2;
+        const maxBark = this.hzToBark(nyquist);
+        const freqSpectrum = new Array(this.fftSize / 2).fill(0);
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            const freq = (i * this.sampleRate) / this.fftSize;
+            const bark = this.hzToBark(freq);
+            const bandIndex = (bark / maxBark) * this.numBands;
+            // Linear interpolation
+            const lowerBand = Math.floor(bandIndex);
+            const upperBand = Math.min(Math.ceil(bandIndex), this.numBands - 1);
+            const t = bandIndex - lowerBand;
+            if (lowerBand >= 0 && lowerBand < loudness.length && upperBand >= 0 && upperBand < loudness.length) {
+                freqSpectrum[i] = loudness[lowerBand] * (1 - t) + loudness[upperBand] * t;
+            }
+            else if (lowerBand >= 0 && lowerBand < loudness.length) {
+                freqSpectrum[i] = loudness[lowerBand];
+            }
+        }
+        // Pad to fftSize for IFFT (mirror for real signal)
+        const padded = new Array(this.fftSize).fill(0);
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            padded[i] = freqSpectrum[i];
+            if (i > 0) {
+                padded[this.fftSize - i] = freqSpectrum[i]; // Mirror for real IFFT
+            }
+        }
+        // Inverse FFT to get autocorrelation
+        const input = this.fft.createComplexArray();
+        const output = this.fft.createComplexArray();
+        for (let i = 0; i < padded.length; i++) {
+            input[2 * i] = padded[i];
+            input[2 * i + 1] = 0;
+        }
+        this.fft.inverseTransform(output, input);
+        // Extract autocorrelation (first order+1 values)
+        const autocorr = new Array(this.order + 1).fill(0);
+        for (let i = 0; i <= this.order; i++) {
+            autocorr[i] = output[2 * i] / this.fftSize;
+        }
+        return autocorr;
+    }
+    /**
+     * Levinson-Durbin recursion for LPC coefficients.
+     */
+    levinsonDurbin(R, order) {
+        const a = new Array(order + 1).fill(0);
+        const e = new Array(order + 1).fill(0);
+        a[0] = 1;
+        e[0] = R[0];
+        for (let i = 1; i <= order; i++) {
+            let acc = 0;
+            for (let j = 1; j < i; j++) {
+                acc += a[j] * R[i - j];
+            }
+            const k = (R[i] - acc) / (e[i - 1] + 1e-12);
+            a[i] = k;
+            for (let j = 1; j < i; j++) {
+                a[j] = a[j] - k * a[i - j];
+            }
+            e[i] = (1 - k * k) * e[i - 1];
+        }
+        return a.slice(1);
+    }
+}
+//# sourceMappingURL=PLPAnalyzer.js.map

package/dist/analysis/PLPAnalyzer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"PLPAnalyzer.js","sourceRoot":"","sources":["../../src/analysis/PLPAnalyzer.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,GAAG,MAAM,QAAQ,CAAC;AASzB,MAAM,OAAO,WAAW;IACd,GAAG,CAAM;IACT,OAAO,CAAS;IAChB,KAAK,CAAS;IACd,QAAQ,CAAS;IACjB,UAAU,CAAS;IACnB,WAAW,CAAa;IAEhC;;OAEG;IACH,YAAY,MAAiB;QAC3B,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,GAAG,CAAC;QACrC,IAAI,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,IAAI,EAAE,CAAC;QAChC,IAAI,CAAC,QAAQ,GAAG,MAAM,CAAC,QAAQ,IAAI,EAAE,CAAC;QACtC,IAAI,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;QACpC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,oBAAoB,EAAE,CAAC;IACjD,CAAC;IAED;;;;OAIG;IACH,YAAY,CAAC,GAAiB;QAC5B,4BAA4B;QAC5B,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI,CAAC,OAAO,EAAE,CAAC;YAChC,MAAM,IAAI,KAAK,CAAC,8BAA8B,IAAI,CAAC,OAAO,kBAAkB,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC;QAC5F,CAAC;QACD,iCAAiC;QACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC/C,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAEpC,MAAM,aAAa,GAAa,EAAE,CAAC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/B,aAAa,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACxC,CAAC;QAED,8CAA8C;QAC9C,MAAM,YAAY,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CACjD,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAC5D,CAAC;QAEF,sCAAsC;QACtC,MAAM,aAAa,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACnD,MAAM,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,yBAAyB;YAC/C,MAAM,aAAa,GAAG,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC;YACrD,OAAO,MAAM,GAAG,aAAa,CAAC;QAChC,CAAC,CAAC,CAAC;QAEH,oDAAoD;QACpD,MAAM,QAAQ,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QAExE,6CAA6C;QAC7C,MAAM,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,CAAC,CAAC;QAE/C,6CAA6C;QAC7C,OAAO,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;IACnD,CAAC;IAED;;;OAGG;IACK,oBAAoB;QAC1B,MAAM,OAAO,GAAe,EAAE,CAAC;QAC/B,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC;QACpC,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC;QAEvC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,IAAI,CAAC,QAAQ,EAAE,IAAI,EAAE,EAAE,CAAC;YAChD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnD,MAAM,UAAU,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC;YAE1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;gBAClD,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;gBACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,UAAU,CAAC,CAAC;gBAE7C,0BAA0B;gBAC1B,IAAI,QAAQ,GAAG,GAAG,EAAE,CAAC;oBACnB,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,QAAQ,CAAC;gBAC7B,CAAC;YACH,CAAC;YAED,YAAY;YACZ,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9C,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;gBACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBACvC,MAAM,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;gBACnB,CAAC;YACH,CAAC;YAED,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACvB,CAAC;QAED,OAAO,OAAO,CAAC;IACjB,CAAC;IAED;;OAEG;IACK,QAAQ,CAAC,EAAU;QACzB,OAAO,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC;IAED;;OAEG;IACK,mBAAmB,CAAC,IAAY;QACtC,kCAAkC;QAClC,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;QACjD,OAAO,CAAC,CAAC;IACX,CAAC;IAED;;;OAGG;IACK,cAAc,CAAC,QAAkB;QACvC,4DAA4D;QAC5D,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC;QACpC,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC;QACvC,MAAM,YAAY,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAEzD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;YAClD,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YACjC,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;YAEnD,uBAAuB;YACvB,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;YACxC,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,EAAE,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;YACpE,MAAM,CAAC,GAAG,SAAS,GAAG,SAAS,CAAC;YAEhC,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,GAAG,QAAQ,CAAC,MAAM,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;gBACnG,YAAY,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;YAC5E,CAAC;iBAAM,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;gBACzD,YAAY,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,SAAS,CAAC,CAAC;YACxC,CAAC;QACH,CAAC;QAED,mDAAmD;QACnD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;YAC5B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACV,MAAM,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,uBAAuB;YACrE,CAAC;QACH,CAAC;QAED,qCAAqC;QACrC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC5C,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACzB,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACvB,CAAC;QACD,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAEzC,iDAAiD;QACjD,MAAM,QAAQ,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QAC7C,CAAC;QAED,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED;;OAEG;IACK,cAAc,CAAC,CAAW,EAAE,KAAa;QAC/C,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACvC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACzB,CAAC;YACD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;YAC5C,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC7B,CAAC;YACD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAChC,CAAC;QAED,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;CACF"}

package/dist/analysis/RMSEAnalyzer.d.ts

@@ -0,0 +1,21 @@
+/**
+ * RMSEAnalyzer computes the Root Mean Square Energy (RMSE) of a PCM audio frame.
+ * RMSE is a simple measure of signal energy, useful for detecting silence, loudness, or dynamic range.
+ * Usage: const rmse = new RMSEAnalyzer({ sampleRate: 16000 }); rmse.analyzeFrame(pcm)
+ */
+export interface RMSEConfig {
+    sampleRate: number;
+}
+export declare class RMSEAnalyzer {
+    /**
+     * @param _config - Configuration with sample rate
+     */
+    constructor(_config: RMSEConfig);
+    /**
+     * Compute the root mean square energy for a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns RMSE value
+     */
+    analyzeFrame(pcm: Float32Array): number;
+}
+//# sourceMappingURL=RMSEAnalyzer.d.ts.map

package/dist/analysis/RMSEAnalyzer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"RMSEAnalyzer.d.ts","sourceRoot":"","sources":["../../src/analysis/RMSEAnalyzer.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AACH,MAAM,WAAW,UAAU;IAAG,UAAU,EAAE,MAAM,CAAA;CAAE;AAElD,qBAAa,YAAY;IACvB;;OAEG;gBACS,OAAO,EAAE,UAAU;IAE/B;;;;OAIG;IACH,YAAY,CAAC,GAAG,EAAE,YAAY,GAAG,MAAM;CAOxC"}

package/dist/analysis/RMSEAnalyzer.js

@@ -0,0 +1,20 @@
+export class RMSEAnalyzer {
+    /**
+     * @param _config - Configuration with sample rate
+     */
+    constructor(_config) {
+    }
+    /**
+     * Compute the root mean square energy for a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns RMSE value
+     */
+    analyzeFrame(pcm) {
+        let sum = 0;
+        for (let i = 0; i < pcm.length; i++) {
+            sum += pcm[i] * pcm[i];
+        }
+        return Math.sqrt(sum / pcm.length);
+    }
+}
+//# sourceMappingURL=RMSEAnalyzer.js.map

package/dist/analysis/RMSEAnalyzer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RMSEAnalyzer.js","sourceRoot":"","sources":["../../src/analysis/RMSEAnalyzer.ts"],"names":[],"mappings":"AAOA,MAAM,OAAO,YAAY;IACvB;;OAEG;IACH,YAAY,OAAmB;IAC/B,CAAC;IACD;;;;OAIG;IACH,YAAY,CAAC,GAAiB;QAC5B,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACzB,CAAC;QACD,OAAO,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;IACrC,CAAC;CACF"}

package/dist/analysis/SpectralBandwidthAnalyzer.d.ts

@@ -0,0 +1,22 @@
+export interface SpectralBandwidthConfig {
+    /** Sample rate of the audio (Hz) */
+    sampleRate: number;
+    /** FFT size (default: 1024) */
+    fftSize?: number;
+}
+export declare class SpectralBandwidthAnalyzer {
+    private fft;
+    private fftSize;
+    private sampleRate;
+    /**
+     * @param config - Configuration with sample rate and optional FFT size
+     */
+    constructor(config: SpectralBandwidthConfig);
+    /**
+     * Compute the spectral bandwidth for a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns Spectral bandwidth in Hz
+     */
+    analyzeFrame(pcm: Float32Array): number;
+}
+//# sourceMappingURL=SpectralBandwidthAnalyzer.d.ts.map

package/dist/analysis/SpectralBandwidthAnalyzer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpectralBandwidthAnalyzer.d.ts","sourceRoot":"","sources":["../../src/analysis/SpectralBandwidthAnalyzer.ts"],"names":[],"mappings":"AAOA,MAAM,WAAW,uBAAuB;IACtC,oCAAoC;IACpC,UAAU,EAAE,MAAM,CAAC;IACnB,+BAA+B;IAC/B,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED,qBAAa,yBAAyB;IACpC,OAAO,CAAC,GAAG,CAAM;IACjB,OAAO,CAAC,OAAO,CAAS;IACxB,OAAO,CAAC,UAAU,CAAS;IAC3B;;OAEG;gBACS,MAAM,EAAE,uBAAuB;IAK3C;;;;OAIG;IACH,YAAY,CAAC,GAAG,EAAE,YAAY,GAAG,MAAM;CAyBxC"}

package/dist/analysis/SpectralBandwidthAnalyzer.js

@@ -0,0 +1,51 @@
+/**
+ * SpectralBandwidthAnalyzer computes the spectral bandwidth of a PCM audio frame.
+ * The spectral bandwidth is the spread of the spectrum around the centroid, related to the perceived timbre of a sound.
+ * Usage: const sb = new SpectralBandwidthAnalyzer({ sampleRate: 16000 }); sb.analyzeFrame(pcm)
+ */
+import FFT from 'fft.js';
+export class SpectralBandwidthAnalyzer {
+    fft;
+    fftSize;
+    sampleRate;
+    /**
+     * @param config - Configuration with sample rate and optional FFT size
+     */
+    constructor(config) {
+        this.fftSize = config.fftSize || 1024;
+        this.fft = new FFT(this.fftSize);
+        this.sampleRate = config.sampleRate;
+    }
+    /**
+     * Compute the spectral bandwidth for a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns Spectral bandwidth in Hz
+     */
+    analyzeFrame(pcm) {
+        const spectrum = this.fft.createComplexArray();
+        this.fft.realTransform(spectrum, pcm);
+        this.fft.completeSpectrum(spectrum);
+        let mags = [];
+        for (let i = 0; i < this.fftSize / 2; i++) {
+            const re = spectrum[2 * i];
+            const im = spectrum[2 * i + 1];
+            mags.push(Math.sqrt(re * re + im * im));
+        }
+        // Spectral centroid
+        let num = 0, denom = 0;
+        for (let i = 0; i < mags.length; i++) {
+            num += i * mags[i];
+            denom += mags[i];
+        }
+        if (denom === 0)
+            return 0;
+        const centroid = (num / denom);
+        // Bandwidth
+        let spread = 0;
+        for (let i = 0; i < mags.length; i++) {
+            spread += mags[i] * Math.pow(i - centroid, 2);
+        }
+        return Math.sqrt(spread / denom) * (this.sampleRate / 2) / mags.length;
+    }
+}
+//# sourceMappingURL=SpectralBandwidthAnalyzer.js.map

package/dist/analysis/SpectralBandwidthAnalyzer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpectralBandwidthAnalyzer.js","sourceRoot":"","sources":["../../src/analysis/SpectralBandwidthAnalyzer.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AACH,OAAO,GAAG,MAAM,QAAQ,CAAC;AASzB,MAAM,OAAO,yBAAyB;IAC5B,GAAG,CAAM;IACT,OAAO,CAAS;IAChB,UAAU,CAAS;IAC3B;;OAEG;IACH,YAAY,MAA+B;QACzC,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,IAAI,CAAC;QACtC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;IACtC,CAAC;IACD;;;;OAIG;IACH,YAAY,CAAC,GAAiB;QAC5B,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;QAC/C,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QACpC,IAAI,IAAI,GAAa,EAAE,CAAC;QACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1C,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC1C,CAAC;QACD,oBAAoB;QACpB,IAAI,GAAG,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,CAAC;QACD,IAAI,KAAK,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC;QAC1B,MAAM,QAAQ,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC;QAC/B,YAAY;QACZ,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,EAAE,CAAC,CAAC,CAAC;QAChD,CAAC;QACD,OAAO,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;IACzE,CAAC;CACF"}

package/dist/analysis/SpectralCentroidAnalyzer.d.ts

@@ -0,0 +1,22 @@
+export interface SpectralCentroidConfig {
+    /** Sample rate of the audio (Hz) */
+    sampleRate: number;
+    /** FFT size (default: 1024) */
+    fftSize?: number;
+}
+export declare class SpectralCentroidAnalyzer {
+    private fft;
+    private fftSize;
+    private sampleRate;
+    /**
+     * @param config - Configuration with sample rate and optional FFT size
+     */
+    constructor(config: SpectralCentroidConfig);
+    /**
+     * Compute the spectral centroid for a PCM frame.
+     * @param pcm - Input PCM audio frame
+     * @returns Spectral centroid in Hz
+     */
+    analyzeFrame(pcm: Float32Array): number;
+}
+//# sourceMappingURL=SpectralCentroidAnalyzer.d.ts.map

package/dist/analysis/SpectralCentroidAnalyzer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpectralCentroidAnalyzer.d.ts","sourceRoot":"","sources":["../../src/analysis/SpectralCentroidAnalyzer.ts"],"names":[],"mappings":"AAOA,MAAM,WAAW,sBAAsB;IACrC,oCAAoC;IACpC,UAAU,EAAE,MAAM,CAAC;IACnB,+BAA+B;IAC/B,OAAO,CAAC,EAAE,MAAM,CAAC;CAClB;AAED,qBAAa,wBAAwB;IACnC,OAAO,CAAC,GAAG,CAAM;IACjB,OAAO,CAAC,OAAO,CAAS;IACxB,OAAO,CAAC,UAAU,CAAS;IAC3B;;OAEG;gBACS,MAAM,EAAE,sBAAsB;IAK1C;;;;OAIG;IACH,YAAY,CAAC,GAAG,EAAE,YAAY,GAAG,MAAM;CAkBxC"}