autoBarkSpectrogram Sound_to_BarkSpectrogram (Sound me, double analysisWidth, double dt, double f1_bark, double fmax_bark, double df_bark) {
try {
double nyquist = 0.5 / my dx, samplingFrequency = 2 * nyquist;
double windowDuration = 2 * analysisWidth; /* gaussian window */
double zmax = NUMhertzToBark2 (nyquist);
double fmin_bark = 0;
// Check defaults.
if (f1_bark <= 0) {
f1_bark = 1;
}
if (fmax_bark <= 0) {
fmax_bark = zmax;
}
if (df_bark <= 0) {
df_bark = 1;
}
fmax_bark = MIN (fmax_bark, zmax);
long numberOfFilters = lround ( (fmax_bark - f1_bark) / df_bark);
if (numberOfFilters <= 0) {
Melder_throw (U"The combination of filter parameters is not valid.");
}
long numberOfFrames; double t1;
Sampled_shortTermAnalysis (me, windowDuration, dt, & numberOfFrames, & t1);
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoBarkSpectrogram thee = BarkSpectrogram_create (my xmin, my xmax, numberOfFrames, dt, t1, fmin_bark, fmax_bark, numberOfFilters, df_bark, f1_bark);
autoMelderProgress progess (U"BarkSpectrogram analysis");
for (long iframe = 1; iframe <= numberOfFrames; iframe++) {
double t = Sampled_indexToX (thee.get(), iframe);
Sound_into_Sound (me, sframe.get(), t - windowDuration / 2.0);
Sounds_multiply (sframe.get(), window.get());
Sound_into_BarkSpectrogram_frame (sframe.get(), thee.get(), iframe);
if (iframe % 10 == 1) {
Melder_progress ( (double) iframe / numberOfFrames, U"BarkSpectrogram analysis: frame ",
iframe, U" from ", numberOfFrames, U".");
}
}
_Spectrogram_windowCorrection ((Spectrogram) thee.get(), window -> nx);
return thee;
} catch (MelderError) {
Melder_throw (me, U": no BarkSpectrogram created.");
}
}
BarkFilter Sound_to_BarkFilter (Sound me, double analysisWidth, double dt,
double f1_bark, double fmax_bark, double df_bark) {
try {
double t1, nyquist = 0.5 / my dx, samplingFrequency = 2 * nyquist;
double windowDuration = 2 * analysisWidth; /* gaussian window */
double zmax = NUMhertzToBark2 (nyquist);
double fmin_bark = 0;
long nt, frameErrorCount = 0;
// Check defaults.
if (f1_bark <= 0) {
f1_bark = 1;
}
if (fmax_bark <= 0) {
fmax_bark = zmax;
}
if (df_bark <= 0) {
df_bark = 1;
}
fmax_bark = MIN (fmax_bark, zmax);
long nf = floor ( (fmax_bark - f1_bark) / df_bark + 0.5);
if (nf <= 0) {
Melder_throw ("The combination of filter parameters is not valid.");
}
Sampled_shortTermAnalysis (me, windowDuration, dt, & nt, & t1);
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoBarkFilter thee = BarkFilter_create (my xmin, my xmax, nt, dt, t1,
fmin_bark, fmax_bark, nf, df_bark, f1_bark);
autoMelderProgress progess (L"BarkFilter analysis");
for (long i = 1; i <= nt; i++) {
double t = Sampled_indexToX (thee.peek(), i);
Sound_into_Sound (me, sframe.peek(), t - windowDuration / 2);
Sounds_multiply (sframe.peek(), window.peek());
if (! Sound_into_BarkFilter_frame (sframe.peek(), thee.peek(), i)) {
frameErrorCount++;
}
if ( (i % 10) == 1) {
Melder_progress ( (double) i / nt, L"BarkFilter analysis: frame ",
Melder_integer (i), L" from ", Melder_integer (nt), L"."); therror
}
}
if (frameErrorCount > 0) {
Melder_warning (L"Analysis results of ", Melder_integer (frameErrorCount), L" frame(s) out of ",
Melder_integer (nt), L" will be suspect.");
}
double ref = FilterBank_DBREF * gaussian_window_squared_correction (window -> nx);
NUMdmatrix_to_dBs (thy z, 1, thy ny, 1, thy nx, ref, FilterBank_DBFAC, FilterBank_DBFLOOR);
return thee.transfer();
} catch (MelderError) {
LPC LPC_and_Sound_to_LPC_robust (LPC thee, Sound me, double analysisWidth, double preEmphasisFrequency, double k,
int itermax, double tol, int wantlocation) {
struct huber_struct struct_huber = { 0 };
try {
double t1, samplingFrequency = 1.0 / my dx, tol_svd = 0.000001;
double location = 0, windowDuration = 2 * analysisWidth; /* Gaussian window */
long nFrames, frameErrorCount = 0, iter = 0;
long p = thy maxnCoefficients;
if (my xmin != thy xmin || my xmax != thy xmax) {
Melder_throw ("Time domains differ.");
}
if (my dx != thy samplingPeriod) {
Melder_throw ("Sampling intervals differ.");
}
if (floor (windowDuration / my dx) < p + 1) {
Melder_throw ("Analysis window too short.");
}
Sampled_shortTermAnalysis (me, windowDuration, thy dx, & nFrames, & t1);
if (nFrames != thy nx || t1 != thy x1) {
Melder_throw ("Incorrect retrieved analysis width");
}
autoSound sound = Data_copy (me);
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoLPC him = Data_copy (thee);
huber_struct_init (&struct_huber, windowDuration, p, samplingFrequency, location, wantlocation);
struct_huber.k = k;
struct_huber.tol = tol;
struct_huber.tol_svd = tol_svd;
struct_huber.itermax = itermax;
autoMelderProgress progess (L"LPC analysis");
Sound_preEmphasis (sound.peek(), preEmphasisFrequency);
for (long i = 1; i <= nFrames; i++) {
LPC_Frame lpc = (LPC_Frame) & thy d_frames[i];
LPC_Frame lpcto = (LPC_Frame) & his d_frames[i];
double t = Sampled_indexToX (thee, i);
Sound_into_Sound (sound.peek(), sframe.peek(), t - windowDuration / 2);
Vector_subtractMean (sframe.peek());
Sounds_multiply (sframe.peek(), window.peek());
try {
LPC_Frames_and_Sound_huber (lpc, sframe.peek(), lpcto, & struct_huber);
} catch (MelderError) {
frameErrorCount++;
}
iter += struct_huber.iter;
if ( (i % 10) == 1) {
Melder_progress ( (double) i / nFrames, L"LPC analysis of frame ",
Melder_integer (i), L" out of ", Melder_integer (nFrames), L".");
}
}
if (frameErrorCount) Melder_warning (L"Results of ", Melder_integer (frameErrorCount),
L" frame(s) out of ", Melder_integer (nFrames), L" could not be optimised.");
MelderInfo_writeLine4 (L"Number of iterations: ", Melder_integer (iter),
L"\n Average per frame: ", Melder_double (((double) iter) / nFrames));
huber_struct_destroy (&struct_huber);
return him.transfer();
} catch (MelderError) {
huber_struct_destroy (&struct_huber);
Melder_throw (me, ": no robust LPC created.");
}
}
