void Spectrum_drawPhases (Spectrum me, Graphics g, double fmin, double fmax, double phase_min, double phase_max, int unwrap, int garnish) {
autoMatrix thee;
int reverse_sign = my z[1][1] < 0;
if (unwrap) {
thee = Spectrum_unwrap (me);
} else {
thee = Matrix_create (my xmin, my xmax, my nx, my dx, my x1, 1.0, 2.0, 2, 1.0, 1.0);
for (long i = 1; i <= my nx; i ++) {
thy z[2][i] = PPVPHA (my z[1][i], my z[2][i], reverse_sign);
}
}
Matrix_drawRows (thee.get(), g, fmin, fmax, 1.9, 2.1, phase_min, phase_max);
if (garnish) {
}
}
Cepstrum Spectrum_to_Cepstrum (Spectrum me) {
try {
autoMatrix unwrap = Spectrum_unwrap (me);
autoSpectrum sx = Data_copy (me);
// Copy magnitude-squared and unwrapped phase.
for (long i = 1; i <= my nx; i ++) {
double xa = unwrap -> z[1][i];
sx -> z[1][i] = xa > 0 ? 0.5 * log (xa) : -300;
sx -> z[2][i] = unwrap -> z[2][i];
}
// Compute complex cepstrum x.
autoSound x = Spectrum_to_Sound (sx.peek());
autoCepstrum thee = Cepstrum_create (0, x -> xmax - x -> xmin, x -> nx);
NUMvector_copyElements (x -> z[1], thy z[1], 1, x -> nx);
return thee.transfer();
} catch (MelderError) {
Melder_throw (me, ": no Cepstrum created.");
}
}
