pitch_internal.cc (2162B)
1 /* 2 * Copyright (c) 2012 The WebRTC project authors. All Rights Reserved. 3 * 4 * Use of this source code is governed by a BSD-style license 5 * that can be found in the LICENSE file in the root of the source 6 * tree. An additional intellectual property rights grant can be found 7 * in the file PATENTS. All contributing project authors may 8 * be found in the AUTHORS file in the root of the source tree. 9 */ 10 11 #include "modules/audio_processing/vad/pitch_internal.h" 12 13 #include <cmath> 14 15 namespace webrtc { 16 17 // A 4-to-3 linear interpolation. 18 // The interpolation constants are derived as following: 19 // Input pitch parameters are updated every 7.5 ms. Within a 30-ms interval 20 // we are interested in pitch parameters of 0-5 ms, 10-15ms and 20-25ms. This is 21 // like interpolating 4-to-6 and keep the odd samples. 22 // The reason behind this is that LPC coefficients are computed for the first 23 // half of each 10ms interval. 24 static void PitchInterpolation(double old_val, const double* in, double* out) { 25 out[0] = 1. / 6. * old_val + 5. / 6. * in[0]; 26 out[1] = 5. / 6. * in[1] + 1. / 6. * in[2]; 27 out[2] = 0.5 * in[2] + 0.5 * in[3]; 28 } 29 30 void GetSubframesPitchParameters(int sampling_rate_hz, 31 double* gains, 32 double* lags, 33 int num_in_frames, 34 int num_out_frames, 35 double* log_old_gain, 36 double* old_lag, 37 double* log_pitch_gain, 38 double* pitch_lag_hz) { 39 // Gain interpolation is in log-domain, also returned in log-domain. 40 for (int n = 0; n < num_in_frames; n++) 41 gains[n] = log(gains[n] + 1e-12); 42 43 // Interpolate lags and gains. 44 PitchInterpolation(*log_old_gain, gains, log_pitch_gain); 45 *log_old_gain = gains[num_in_frames - 1]; 46 PitchInterpolation(*old_lag, lags, pitch_lag_hz); 47 *old_lag = lags[num_in_frames - 1]; 48 49 // Convert pitch-lags to Hertz. 50 for (int n = 0; n < num_out_frames; n++) { 51 pitch_lag_hz[n] = (sampling_rate_hz) / (pitch_lag_hz[n]); 52 } 53 } 54 55 } // namespace webrtc