2d:null 2e:["$","title",null,{"children":"audioprocessor | Korakot C. | Observable"}] 33:Ta1f,AudioProcessor = { var AudioProcessor = {} AudioProcessor.F1_MIN = 200 AudioProcessor.F1_MAX = 1000 AudioProcessor.F2_MIN = 700 AudioProcessor.F2_MAX = 3000 AudioProcessor.F3_MIN = 1500 AudioProcessor.F3_MAX = 5000 AudioProcessor.DEFAULT_MAX_FORMANT_MALE = 5000 AudioProcessor.DEFAULT_MAX_FORMANT_FEMALE = 5500 AudioProcessor.DEFAULT_MAX_FORMANT_CHILD = 8000 AudioProcessor.HANNING_SHIFT = 32 AudioProcessor._HANNING_WINDOW = { 61: new Float32Array([ 0. , 0.00273905, 0.0109262 , 0.02447174, 0.04322727, 0.0669873 , 0.0954915 , 0.12842759, 0.1654347 , 0.20610737, 0.25 , 0.29663168, 0.3454915 , 0.39604415, 0.44773577, 0.5 , 0.55226423, 0.60395585, 0.6545085 , 0.70336832, 0.75 , 0.79389263, 0.8345653 , 0.87157241, 0.9045085 , 0.9330127 , 0.95677273, 0.97552826, 0.9890738 , 0.99726095, 1. , 0.99726095, 0.9890738 , 0.97552826, 0.95677273, 0.9330127 , 0.9045085 , 0.87157241, 0.8345653 , 0.79389263, 0.75 , 0.70336832, 0.6545085 , 0.60395585, 0.55226423, 0.5 , 0.44773577, 0.39604415, 0.3454915 , 0.29663168, 0.25 , 0.20610737, 0.1654347 , 0.12842759, 0.0954915 , 0.0669873 , 0.04322727, 0.02447174, 0.0109262 , 0.00273905, 0. ]), 75: new Float32Array([ 0. , 0.00180126, 0.00719204, 0.01613353, 0.02856128, 0.04438575, 0.06349294, 0.08574518, 0.11098212, 0.13902195, 0.16966264, 0.20268341, 0.23784636, 0.27489813, 0.31357176, 0.35358861, 0.39466037, 0.43649109, 0.4787794 , 0.5212206 , 0.56350891, 0.60533963, 0.64641139, 0.68642824, 0.72510187, 0.76215364, 0.79731659, 0.83033736, 0.86097805, 0.88901788, 0.91425482, 0.93650706, 0.95561425, 0.97143872, 0.98386647, 0.99280796, 0.99819874, 1. , 0.99819874, 0.99280796, 0.98386647, 0.97143872, 0.95561425, 0.93650706, 0.91425482, 0.88901788, 0.86097805, 0.83033736, 0.79731659, 0.76215364, 0.72510187, 0.68642824, 0.64641139, 0.60533963, 0.56350891, 0.5212206 , 0.4787794 , 0.43649109, 0.39466037, 0.35358861, 0.31357176, 0.27489813, 0.23784636, 0.20268341, 0.16966264, 0.13902195, 0.11098212, 0.08574518, 0.06349294, 0.04438575, 0.02856128, 0.01613353, 0.00719204, 0.00180126, 0. ]) }; return AudioProcessor }34:T4c3,initFilterBanks = function (Nspec, noFilterBanks, minFreq, maxFreq, sampleRate) { var filterBanks = []; var NFFT = Nspec * 2; var totalFilters = noFilterBanks; var minMel = freqToMel(minFreq); var maxMel = freqToMel(maxFreq); var dMel = (maxMel - minMel) / (noFilterBanks + 1); var melSpacing = []; var fftFreqs2Mel = []; var lower = []; var center = []; // Init melSpacing for (var i = 0; i < noFilterBanks + 2; i++) { var mel = minMel + i * dMel; melSpacing.push(mel); } // Init fftFreqs2Mel for (var i = 0; i < Nspec; i++) { var fftFreq = i * sampleRate / NFFT; var fftFreq2Mel = freqToMel(fftFreq); fftFreqs2Mel.push(fftFreq2Mel); } // Init lower for (var i = 0; i < noFilterBanks; i++) { lower.push(melSpacing[i]); } // Init center for (var i = 1; i < noFilterBanks + 1; i++) { center.push(melSpacing[i]); } // Prepare the mel scale filterbank for (var i = 0; i < totalFilters; i++) { var fBank = []; for (var j = 0; j < Nspec; j++) { var val = Math.max(0.0, (1 - Math.abs(fftFreqs2Mel[j] - center[i]) / (center[i] - lower[i]))); fBank.push(val); } filterBanks.push(fBank); } return filterBanks; };35:T4b7,// AudioProcessor.getMFCCs getMFCCs = function (options) { var toLM = options.toLinearMagnitude === undefined ? true : !!options.toLinearMagnitude; var fft = options.fft; if (toLM) { for (var j = 0; j < fft.length; j++) { fft[j] = AudioProcessor.decibelsToLinear(fft[j]); } } else { // We need to ensure that these are all positive values var tmpFFT = []; for (var i = 0; i < fft.length; i++) { tmpFFT[i] = Math.abs(fft[i]); } fft = tmpFFT; } var filterBanks = [], numFilterBanks = options.filterBanks || DEFAULT_FILTER_BANKS, NFFT = options.fftSize, minFreq = options.minFreq || DEFAULT_MIN_HZ, maxFreq = options.maxFreq || DEFAULT_MAX_HZ, sampleRate = options.sampleRate; // Build the mel scale filter banks filterBanks = initFilterBanks(NFFT / 2, numFilterBanks, minFreq, maxFreq, sampleRate); // 1. Map the spectrum to the mel scale (apply triangular filters) // 2. Take the log of it // Performance is better when both of these are done together var logMelSpectrum = logMel(fft, filterBanks); // Perform the Discrete Cosine Transformation var postDCT = dct(logMelSpectrum); return postDCT; };36:T43a,getPeaks = function (smoothedArray, sampleRate, fftSize) { var peaks = []; var previousNum; var currentNum; var nextNum; for (var i = 0; i < smoothedArray.length; i++) { var hz = toFrequency(i, sampleRate, fftSize); var formant = peaks.length + 1; switch (formant) { case 1: if (hz < AudioProcessor.F1_MIN) { continue; } break; case 2: if (hz < AudioProcessor.F2_MIN || hz - peaks[0] < MIN_DIFF_F1_F2) { continue; } break; case 3: if (hz < AudioProcessor.F3_MIN || hz - peaks[1] < MIN_DIFF_F2_F3) { continue; } break; default: return null; } previousNum = smoothedArray[i - 1] || 0; currentNum = smoothedArray[i] || 0; nextNum = smoothedArray[i + 1] || 0; if (currentNum > previousNum && currentNum > nextNum) { if (peakHeight(i, smoothedArray) >= MIN_PEAK_HEIGHT) { peaks.push(hz); if (formant === 3) { return peaks; } } } } return peaks; };37:T910,MathUtils = { var MathUtils = {}; MathUtils.nextPow2 = function(v) { v -= 1; var p = 2; while (v >>= 1) {p <<= 1;} return p; } MathUtils.predict = function(features, weights) { if (features.length !== weights.length) { throw new Error("features and weights must be equal in " + "length when using the regression prediction method. " + "Features length: " + features.length + ". Weights length: " + weights.length); } var sum = 0; for (var i = 0; i < features.length; i++) { sum += features[i] * weights[i]; } return sum; } MathUtils.mapToScale = function(value, sourceMin, sourceMax, destMin, destMax) { return (value - sourceMin) / (sourceMax - sourceMin) * (destMax - destMin) + destMin; } MathUtils.arrayAbs = function(y) { for (var i = 0; i < y.length; i++) { y[i] = Math.abs(y[i]); } return y; }; MathUtils.negArrayAddValue = function(y, value) { for (var i = 0; i < y.length; i++) { y[i] = -y[i] + value; } return y; }; MathUtils.addToArray = function(y, value) { for (var i = 0; i < y.length; i++) { y[i] = y[i] + value; } return y; }; MathUtils.concatenate = function(args) { var p = new Array(); for (var i = 0; i < arguments.length; i++) { for (var j = 0; j < arguments[i].length; j++) { p.push(arguments[i][j]); } } return p; }; MathUtils.flipArray = function(y) { var p = new Array(); for (var i = y.length - 1; i > -1; i--) { p.push(y[i]); } return p; }; MathUtils.pinv = function(A) { var z = numeric.svd(A), foo = z.S[0]; var U = z.U, S = z.S, V = z.V; var m = A.length, n = A[0].length, tol = Math.max(m, n) * numeric.epsilon * foo, M = S.length; var i, Sinv = new Array(M); for (i = M - 1; i !== -1; i--) { if (S[i] > tol) Sinv[i] = 1 / S[i]; else Sinv[i] = 0; } return numeric.dot(numeric.dot(V, numeric.diag(Sinv)), numeric.transpose(U)) }; return MathUtils }38:T66f,// AudioProcessor.savitzkyGolay savitzkyGolay = function(y, window_size, order) { var windowSize = Math.abs(parseInt(window_size, 10)); order = Math.abs(parseInt(order, 10)); var order_range = order + 1; var half_window = (windowSize - 1) / 2; var b = new Array(); for (var k = -half_window; k < half_window + 1; k++) { var row = new Array(); for (var i = 0; i < order_range; i++) { row.push(Math.pow(k, i)); } b.push(row); } //This line needs to be changed if you use something other than 0 for derivative var temp = MathUtils.pinv(b); var m = temp[0]; //if you take a look at firstvals in the python code, and then at this code you'll see that I've only broken firstvals down into different parts such as first taking a sub array, flipping it, and so on var yTemp = new Array(); yTemp = y.subarray ? y.subarray(1, half_window + 1) : y.slice(1, half_window + 1); yTemp = MathUtils.flipArray(yTemp); yTemp = MathUtils.addToArray(yTemp, -y[0]); yTemp = MathUtils.arrayAbs(yTemp); yTemp = MathUtils.negArrayAddValue(yTemp, y[0]); var firstvals = yTemp; //Same thing was done for lastvals var yTemp2 = new Array(); yTemp2 = y.subarray ? y.subarray(-half_window - 1, -1) : y.slice(-half_window - 1, -1); yTemp2 = MathUtils.flipArray(yTemp2); yTemp2 = MathUtils.addToArray(yTemp2, -y[y.length - 1]); yTemp2 = MathUtils.arrayAbs(yTemp2); yTemp2 = MathUtils.addToArray(yTemp2, y[y.length - 1]); var lastvals = yTemp2; y = MathUtils.concatenate(firstvals, y, lastvals); m = MathUtils.flipArray(m); var result = new Array(); result = convolve(m, y); return result; };39:[] 31:{"id":"368ba7012dca3d55","slug":"audioprocessor","title":"audioprocessor","thumbnail":"f5d44d86ace597c9e006495b5a738f03f134c63b0ebe9b3d96916c00e8a1a92c","access_level":"public","create_time":"$D2019-05-24T10:47:16.753Z","update_time":"$D2019-05-24T12:22:04.266Z","view_time":null,"first_public_time":"$D2019-05-24T12:05:43.553Z","has_unpublished_edits":null,"fork":false,"stars":0,"starred":false,"owner":{"type":"team","id":"a8da4966abe85e60","login":"kora","name":"Korakot C.","avatar_url":"https://avatars-next.observableusercontent.com/c22a27a09015038f899418e3a73866ded057c2a02c93b43488c6a3e985d25824","tier":"starter_2024"},"version":180,"publish_version":180,"role":null,"authored":false,"authors":[{"type":"individual","id":"5a8371b5bb129d9a","login":"kora","name":"Korakot C.","avatar_url":"https://avatars-next.observableusercontent.com/c22a27a09015038f899418e3a73866ded057c2a02c93b43488c6a3e985d25824","tier":null}],"members":[{"type":"individual","id":"5a8371b5bb129d9a","login":"kora","name":"Korakot C.","avatar_url":"https://avatars-next.observableusercontent.com/c22a27a09015038f899418e3a73866ded057c2a02c93b43488c6a3e985d25824","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":2,"collection":null,"type":"notebook","body":{"title":"audioprocessor","cells":[{"id":0,"value":"md`# audioprocessor\n\nพยายามแกะจาก [audioprocessor.js](https://github.com/BYU-ODH/apeworm/blob/master/src/modules/utils/audioprocessor.js)\nเพื่อหา formant`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":24,"value":"DEFAULT_MIN_HZ = 0","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":27,"value":"DEFAULT_MAX_HZ = 8000","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":30,"value":"DEFAULT_FILTER_BANKS = 40","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":33,"value":"WINDOW_SIZES = [\n 75,\n 61\n ]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":18,"value":"MIN_DIFF_F1_F2 = 150","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":20,"value":"MIN_DIFF_F2_F3 = 500","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":36,"value":"MIN_PEAK_HEIGHT = 0.1","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":12,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":50,"value":"// AudioProcessor.dct\ndct = function(values) {\n var result = [];\n for (var i = 0; i < values.length; i++) {\n var val = 0;\n for (var j = 0; j < values.length; j++) {\n val += values[j] * Math.cos(i * (j + 0.5) * Math.PI / values.length);\n }\n\n // Perform scaling used by matlab implementation of dct\n if (i == 0) {\n val /= Math.sqrt(2.0);\n }\n val *= Math.sqrt(2.0 / values.length);\n\n result.push(val);\n }\n\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":58,"value":"// AudioProcessor.inverseDct \ninverseDct = function(values) {\n var result = [];\n for (var i = 0; i < values.length; i++) {\n // Add the dc term\n var val = values[0] / Math.sqrt(2);\n\n // Add ac terms\n for (var j = 1; j < values.length; j++) {\n val += values[j] * Math.cos(j * (i + 0.5) * Math.PI / values.length);\n }\n\n // Perform scaling used by matlab implementation of dct\n val *= Math.sqrt(2.0 / values.length);\n\n result.push(val);\n }\n\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":61,"value":"// AudioProcessor.computeDft \ncomputeDft = function (inreal, inimag) {\n var n = inreal.length;\n var outreal = new Array(n);\n var outimag = new Array(n);\n for (var k = 0; k < n; k++) { // For each output element\n var sumreal = 0;\n var sumimag = 0;\n for (var t = 0; t < n; t++) { // For each input element\n var angle = 2 * Math.PI * t * k / n;\n sumreal += inreal[t] * Math.cos(angle) + inimag[t] * Math.sin(angle);\n sumimag += -inreal[t] * Math.sin(angle) + inimag[t] * Math.cos(angle);\n }\n outreal[k] = sumreal;\n outimag[k] = sumimag;\n }\n return [outreal, outimag];\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":64,"value":"// AudioProcessor.inverseDft \ninverseDft = function (inreal, inimag) {\n var n = inreal.length;\n var outreal = new Array(n);\n var outimag = new Array(n);\n for (var k = 0; k < n; k++) { // For each output element\n var sumreal = 0;\n var sumimag = 0;\n for (var t = 0; t < n; t++) { // For each input element\n var angle = 2 * Math.PI * t * k / n;\n sumreal += inreal[t] * Math.cos(angle) + inimag[t] * Math.sin(angle);\n sumimag += -inreal[t] * Math.sin(angle) + inimag[t] * Math.cos(angle);\n }\n outreal[k] = sumreal;\n outimag[k] = sumimag;\n }\n return [outreal, outimag];\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":75,"value":"freqToMel = function(frequency) {\n return 1127.01048 * Math.log(1.0 + frequency / 700);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":81,"value":"melToFreq = function(mel) {\n return 700 * (Math.exp(mel/1127.01048) - 1);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":73,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":69,"value":"logMel = function(spectrum, filterBanks) {\n var result = [];\n for (var i = 0; i < filterBanks.length; i++) {\n var val = 0.0;\n for (var j = 0; j < filterBanks[i].length; j++) {\n val += (filterBanks[i][j]) * spectrum[j];\n }\n result.push(Math.log10(val));\n }\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":67,"value":"$35","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":89,"value":"// AudioProcessor.toFrequency\ntoFrequency = function(position, sampleRate, fftSize) {\n\n var nyquist = sampleRate / 2;\n\n var totalBins = fftSize / 2;\n\n return position * (nyquist / totalBins);\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":93,"value":"peakHeight = function (index, values) {\n var peak = values[index],\n lheight = null,\n rheight = null;\n\n var prev = null;\n // check the left\n for (var i = index - 1; i >= 0; i--) {\n if (prev !== null && values[i] > prev) {\n break;\n }\n prev = values[i];\n lheight = peak - prev;\n }\n\n prev = null;\n // check the right\n for (var i = index + 1; i < values.length; i++) {\n if (prev !== null && values[i] > prev) {\n break;\n }\n prev = values[i];\n rheight = peak - prev;\n }\n\n var result;\n if (lheight === null) {\n result = +rheight;\n }\n else if (rheight === null) {\n result = +lheight;\n }\n else {\n result = lheight < rheight ? lheight : rheight;\n }\n\n if (result < 0) {\n return 0;\n }\n return result;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":11,"value":"$36","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":85,"value":"// AudioProcessor.getFormants\ngetFormants = function (data, sampleRate) {\n var that = this;\n\n if (arguments.length !== 2) {\n throw new Error(\"Invalid arguments. Function must be called as \" +\n \" getFormants(data, sampleRate)\");\n }\n\n var fftSize = data.length * 2;\n\n for (var i = 0; i < WINDOW_SIZES.length; i++) {\n var smooth = hann(data, WINDOW_SIZES[i]) // 75 and 61\n .slice(AudioProcessor.HANNING_SHIFT); // 32\n var formants = getPeaks(smooth, sampleRate, fftSize);\n\n if (formants[0] < AudioProcessor.F1_MIN || formants[0] > AudioProcessor.F1_MAX || formants[0] >= formants[1] ||\n formants[1] < AudioProcessor.F2_MIN || formants[1] > AudioProcessor.F2_MAX || formants[1] >= formants[2] ||\n formants[2] < AudioProcessor.F3_MIN || formants[2] > AudioProcessor.F3_MAX) {\n continue;\n }\n else {\n return formants;\n }\n }\n return []; // no good formants found\n};","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":102,"value":"MathUtils_findPeaks = function (values, numPeaks) {\n var peaks = [];\n var previousNum;\n var currentNum;\n var nextNum;\n\n for (var i = 1; i < values.length - 1 && peaks.length < numPeaks; i++) {\n\n previousNum = values[i - 1] || -Infinity;\n currentNum = values[i];\n nextNum = values[i + 1] || -Infinity;\n\n if (currentNum > previousNum && currentNum > nextNum) {\n peaks.push(i);\n }\n }\n return peaks;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":99,"value":"//AudioProcessor.getFormantsFromCepstrum\ngetFormantsFromCepstrum = function(cepstrum, options) {\n\n var numFormants = options.numFormants || 3;\n var sampleRate = options.sampleRate || 44100;\n var fftSize = options.fftSize || 1024;\n var cutoff = options.cutoff || Math.round(fftSize / 8);\n\n var formants = [];\n\n // Lifter the cepstrum\n for (var i = cutoff; i < cepstrum.length; i++) {\n cepstrum[i] = 0;\n }\n\n // Inverse dct to get the envelope\n var envelope = inverseDct(cepstrum);\n\n // Find peaks in envelope\n var peakPositions = MathUtils_findPeaks(envelope, numFormants);\n\n for (var i = 0; i < numFormants; i++) {\n var frequency = toFrequency(peakPositions[i], sampleRate, fftSize)\n formants.push(frequency);\n }\n\n return formants;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":115,"value":"linearToDecibel = function(magnitude) {\n return Math.log10(magnitude) / 0.05;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":110,"value":"decibelsToLinear = function(dB) {\n return Math.pow(10, 0.05 * dB);\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":108,"value":"// AudioProcessor.getCepstrum\ngetCepstrum = function(fft, options) {\n\n var toLM = options.toLinearMagnitude === undefined ? true : !!options.toLinearMagnitude;\n\n if (toLM) {\n for (var j = 0; j < fft.length; j++) {\n fft[j] = decibelsToLinear(fft[j]);\n }\n }\n else {\n // we need to ensure that these are all positive values\n var tmpFFT = [];\n for (var i = 0; i < fft.length; i++) {\n tmpFFT[i] = Math.abs(fft[i]);\n }\n fft = tmpFFT;\n }\n\n // Compute the log of the spectrum\n var logSpectrum = [];\n for (var i = 0; i < fft.length; i++) {\n logSpectrum.push(Math.log(fft[i]));\n }\n\n // Perform the Discrete Cosine Transformation\n var postDCT = dct(logSpectrum);\n\n return postDCT;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":121,"value":"//AudioProcessor.hann\nhann = function(vals, window_size) {\n if (typeof AudioProcessor._HANNING_WINDOW[window_size] === 'undefined') {\n throw new Error('No precomputed Hanning Window values found for ' +\n window_size);\n }\n\n var s = [];\n\n for (var i = window_size - 1; i > 0; i--) {\n s.push(vals[i]);\n }\n for (var i = 0; i < vals.length; i++) {\n s.push(vals[i]);\n }\n for (var i = vals.length - 1; i > vals.length - window_size; i--) {\n s.push(vals[i]);\n }\n\n var w = AudioProcessor._HANNING_WINDOW[window_size];\n\n var sum = 0;\n var wMorph = [];\n for (var i = 0; i < w.length; i++) {\n sum += w[i];\n }\n for (var i = 0; i < w.length; i++) {\n wMorph[i] = w[i] / sum;\n }\n return convolve(wMorph, s);\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":124,"value":"// AudioProcessor.convolve\nconvolve = function (m, y) {\n var result = new Array(),\n first = null,\n second = null;\n\n if (m.length > y.length) {\n first = y;\n second = m;\n }\n else {\n first = m;\n second = y;\n }\n var size = second.length - first.length + 1;\n for (var i = 0; i < size; i++) {\n var newNum = 0,\n len = first.length;\n\n for (var j = 0; j < first.length; j++) {\n newNum = newNum + first[len - 1 - j] * second[j + i];\n }\n result.push(newNum);\n }\n return result;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":137,"value":"numeric = require('https://bundle.run/numeric@1.2.6')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":136,"value":"$37","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":156,"value":"$38","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":159,"value":"md`เรียกใช้ (import) จากหน้าอื่นได้เลยนะ\n\n import {getFormant} from \"@korakot/audioprocessor\"\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":170,"value":"data = {\n var data = new Float32Array(512)\n return data\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":174,"value":"// ทดลองดู\ngetFormant(data, 22050) // data ต้องเป็น frequency ตะหาก","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W39","stdlib":"1","titling":"auto"}} 32:"$31:body"