sigpr_frame.cc (speech_tools-2.4-release) | : | sigpr_frame.cc (speech_tools-2.5.0-release) | ||
---|---|---|---|---|
skipping to change at line 726 | skipping to change at line 726 | |||
PI_over_liftering_parameter = PI; // since sin(n.PI) == 0 | PI_over_liftering_parameter = PI; // since sin(n.PI) == 0 | |||
// if we are not including cepstral coeff 0 (c0) then we need | // if we are not including cepstral coeff 0 (c0) then we need | |||
// to do a cosine transform 1 longer than otherwise | // to do a cosine transform 1 longer than otherwise | |||
cos_xform_order = include_c0 ? mfcc_vec.length() : mfcc_vec.length() + 1; | cos_xform_order = include_c0 ? mfcc_vec.length() : mfcc_vec.length() + 1; | |||
const_factor = sqrt(2 / (float)(fbank_vec.length())); | const_factor = sqrt(2 / (float)(fbank_vec.length())); | |||
for(i=0;i<mfcc_vec.length();i++) | for(i=0;i<mfcc_vec.length();i++) | |||
{ | { | |||
actual_mfcc_index = include_c0 ? i : i+1; | actual_mfcc_index = include_c0 ? i : i+1; | |||
pi_i_over_N = | pi_i_over_N = | |||
PI * (float)(actual_mfcc_index) / (float)(fbank_vec.length()); | PI * (float)(actual_mfcc_index) / (float)(fbank_vec.length()); | |||
for(j=0;j<fbank_vec.length();j++) | for(j=0;j<fbank_vec.length();j++) | |||
// j + 0.5 is because we want (j+1) - 0.5 | // j + 0.5 is because we want (j+1) - 0.5 | |||
mfcc_vec[i] += fbank_vec(j) * cos(pi_i_over_N * ((float)j + 0.5)); | mfcc_vec[i] += fbank_vec(j) * cos(pi_i_over_N * ((float)j + 0.5)); | |||
mfcc_vec[i] *= const_factor; | mfcc_vec[i] *= const_factor; | |||
End of changes. 1 change blocks. | ||||
1 lines changed or deleted | 0 lines changed or added |