BarkFilter Matrix_to_BarkFilter (I) {
iam (Matrix);
try {
autoBarkFilter thee = BarkFilter_create (my xmin, my xmax, my nx, my dx, my x1,
my ymin, my ymax, my ny, my dy, my y1);
NUMmatrix_copyElements (my z, thy z, 1, my ny, 1, my nx);
return thee.transfer();
} catch (MelderError) {
Melder_throw (me, U": not converted to BarkFilter.");
}
}
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) {