LPC: To Spectrum (slice)...

You can choose this command after selecting 1 or more LPC objects.


Time (s)
the time at which the Spectrum should be calculated.
Minimum frequency resolution (Hz)
successive frequencies in the Spectrum will be maximally this distance apart.
Bandwidth reduction (Hz)
formants with small bandwidths show up very well as peaks in the spectrum because the poles lie close to the contour along which the spectrum is computed (the unit circle in the z-plane). Peak enhancement can be realized by computing the spectrum in the z-plane along a contour of radius r = exp (– π · bandwidthReduction / samplingFrequency). This technique is also called off-axis spectrum computation. Negative values evaluate the spectrum on a contour outside the unit circle and therefore result in a flattened spectrum.
De-emphasis frequency (Hz)
Performs de-emphasis when frequency is in the interval (0, Nyquist frequency)


The Spectrum at time t will be calculated from the nearest LPC_Frame according to:

Spectrum (f) = √(gain·T/df) / (1 + ∑k=1..numberOfCoefficients akzk),

where T is the sampling period and z = exp (–2 π i f T) and df is the distance in Hz between two successive components in the Spectrum.

1. Allocate a large enough buffer[1..nfft] to perform an FFT analysis.
2. Make the first value of the buffer 1 and copy the prediction coefficients a into the buffer. This results in buffer values: (1, a1, ..., anumberOfCoefficients, 0, ..., 0).
3. If De-emphasis frequency is in the interval (0, nyquistFrequency) then "multiply" the buffer with (1 - b z–1), where b = exp (– π deEmphasisFrequency T). This results in buffer values: (1, a1b, a2b·a1, ..., anumberOfCoefficientsb·anumberOfCoefficients–1, –b·anumberOfCoefficients, 0, ..., 0). Note that the number of values in the buffer that differ from 0 has increased by one.
4. If Bandwidth reduction is greater than 0 then multiply corresponding values in the buffer by gi–1 where g = exp (2π bandwidthReduction T / nfft), and i is the position index in the buffer. i runs from 1 to numberOfCoefficients+1+t, where t equals 1 when de-emphasis was performed, else 0.
5. Calculate the FFT spectrum of the buffer with the coefficients. This results in complex amplitudes (aj,bj), j=1..nfft/2+1.
6. Calculate the LPC Spectrum by taking the inverse of the FFT spectrum, i.e., each complex amplitude becomes (aj,bj)–1 = (aj,–bj) / (aj2 + bj2)
7. Multiply all values with the scale factor √(gain·T/df).

Links to this page

© djmw, November 20, 2007