source: src/spectral/spectral_centroid.c @ 3b3b03e

feature/autosinkfeature/cnnfeature/cnn_orgfeature/constantqfeature/crepefeature/crepe_orgfeature/pitchshiftfeature/pydocstringsfeature/timestretchfix/ffmpeg5pitchshiftsamplertimestretchyinfft+
Last change on this file since 3b3b03e was fc6c831, checked in by Paul Brossier <piem@piem.org>, 16 years ago

moved spectral_centroid to new file

  • Property mode set to 100644
File size: 1.2 KB
Line 
1/*
2   Copyright (C) 2007 Paul Brossier
3
4   This program is free software; you can redistribute it and/or modify
5   it under the terms of the GNU General Public License as published by
6   the Free Software Foundation; either version 2 of the License, or
7   (at your option) any later version.
8
9   This program is distributed in the hope that it will be useful,
10   but WITHOUT ANY WARRANTY; without even the implied warranty of
11   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12   GNU General Public License for more details.
13
14   You should have received a copy of the GNU General Public License
15   along with this program; if not, write to the Free Software
16   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
17
18*/
19
20#include "aubio_priv.h"
21#include "cvec.h"
22#include "spectral/spectral_centroid.h"
23
24smpl_t aubio_spectral_centroid(cvec_t * spectrum, smpl_t samplerate) {
25  uint_t i=0, j;
26  smpl_t sum = 0., sc = 0.;
27  for ( j = 0; j < spectrum->length; j++ ) {
28    sum += spectrum->norm[i][j];
29  }
30  if (sum == 0.) return 0.;
31  for ( j = 0; j < spectrum->length; j++ ) {
32    sc += (smpl_t)j * spectrum->norm[i][j];
33  }
34  return sc / sum * samplerate / (smpl_t)(spectrum->length);
35}
36
37
Note: See TracBrowser for help on using the repository browser.