FormantFilter Matrix_to_FormantFilter (I) {
iam (Matrix);
try {
autoFormantFilter thee = FormantFilter_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 FormantFilter.");
}
}
FormantFilter Sound_and_Pitch_to_FormantFilter (Sound me, Pitch thee, double analysisWidth, double dt,
double f1_hz, double fmax_hz, double df_hz, double relative_bw) {
try {
double t1, windowDuration = 2 * analysisWidth; /* gaussian window */
double nyquist = 0.5 / my dx, samplingFrequency = 2 * nyquist, fmin_hz = 0;
long nt, f0_undefined = 0;
if (my xmin > thy xmin || my xmax > thy xmax) Melder_throw
("The domain of the Sound is not included in the domain of the Pitch.");
double f0_median = Pitch_getQuantile (thee, thy xmin, thy xmax, 0.5, kPitch_unit_HERTZ);
if (f0_median == NUMundefined || f0_median == 0) {
f0_median = 100;
Melder_warning (L"Pitch values undefined. Bandwith fixed to 100 Hz. ");
}
if (f1_hz <= 0) {
f1_hz = 100;
}
if (fmax_hz <= 0) {
fmax_hz = nyquist;
}
if (df_hz <= 0) {
df_hz = f0_median / 2;
}
if (relative_bw <= 0) {
relative_bw = 1.1;
}
fmax_hz = MIN (fmax_hz, nyquist);
long nf = floor ( (fmax_hz - f1_hz) / df_hz + 0.5);
Sampled_shortTermAnalysis (me, windowDuration, dt, &nt, &t1);
autoFormantFilter him = FormantFilter_create (my xmin, my xmax, nt, dt, t1,
fmin_hz, fmax_hz, nf, df_hz, f1_hz);
// Temporary objects
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoMelderProgress progress (L"Sound & Pitch: To FormantFilter");
for (long i = 1; i <= nt; i++) {
double t = Sampled_indexToX (him.peek(), i);
double b, f0 = Pitch_getValueAtTime (thee, t, kPitch_unit_HERTZ, 0);
if (f0 == NUMundefined || f0 == 0) {
f0_undefined++; f0 = f0_median;
}
b = relative_bw * f0;
Sound_into_Sound (me, sframe.peek(), t - windowDuration / 2);
Sounds_multiply (sframe.peek(), window.peek());
Sound_into_FormantFilter_frame (sframe.peek(), him.peek(), i, b);
if ( (i % 10) == 1) {
Melder_progress ( (double) i / nt, L"Frame ", Melder_integer (i), L" out of ",
Melder_integer (nt), L".");
}
}
double ref = FilterBank_DBREF * gaussian_window_squared_correction (window -> nx);
NUMdmatrix_to_dBs (his z, 1, his ny, 1, his nx, ref, FilterBank_DBFAC, FilterBank_DBFLOOR);
return him.transfer();
} catch (MelderError) {
Melder_throw ("FormantFilter not created from Pitch & FormantFilter.");
}
}
