98 lines
4.2 KiB
C++
98 lines
4.2 KiB
C++
/************************************************************************/
|
|
/* */
|
|
/* Centre for Speech Technology Research */
|
|
/* University of Edinburgh, UK */
|
|
/* Copyright (c) 1996,1997 */
|
|
/* All Rights Reserved. */
|
|
/* */
|
|
/* Permission is hereby granted, free of charge, to use and distribute */
|
|
/* this software and its documentation without restriction, including */
|
|
/* without limitation the rights to use, copy, modify, merge, publish, */
|
|
/* distribute, sublicense, and/or sell copies of this work, and to */
|
|
/* permit persons to whom this work is furnished to do so, subject to */
|
|
/* the following conditions: */
|
|
/* 1. The code must retain the above copyright notice, this list of */
|
|
/* conditions and the following disclaimer. */
|
|
/* 2. Any modifications must be clearly marked as such. */
|
|
/* 3. Original authors' names are not deleted. */
|
|
/* 4. The authors' names are not used to endorse or promote products */
|
|
/* derived from this software without specific prior written */
|
|
/* permission. */
|
|
/* */
|
|
/* THE UNIVERSITY OF EDINBURGH AND THE CONTRIBUTORS TO THIS WORK */
|
|
/* DISCLAIM ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING */
|
|
/* ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS, IN NO EVENT */
|
|
/* SHALL THE UNIVERSITY OF EDINBURGH NOR THE CONTRIBUTORS BE LIABLE */
|
|
/* FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES */
|
|
/* WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN */
|
|
/* AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, */
|
|
/* ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF */
|
|
/* THIS SOFTWARE. */
|
|
/* */
|
|
/*************************************************************************/
|
|
/* */
|
|
/* Author: Richard Caley (rjc@cstr.ed.ac.uk) */
|
|
/* -------------------------------------------------------------------- */
|
|
/* Some generally useful signal processing things. */
|
|
/* */
|
|
/*************************************************************************/
|
|
|
|
#include "EST_Wave.h"
|
|
#include "EST_TBuffer.h"
|
|
#include "sigpr/EST_misc_sigpr.h"
|
|
|
|
#include "EST_inline_utils.h"
|
|
|
|
static void short_set(EST_Wave &to, int ch,
|
|
const EST_TBuffer<double> &from,
|
|
double maxval)
|
|
{
|
|
for(int i=0; i<to.num_samples(); i++)
|
|
to.a_no_check(i,ch) = irint(from(i)/maxval*10000);
|
|
}
|
|
|
|
/* redundant - see new version in filter */
|
|
void EST_pre_emphasis(EST_Wave &signal, EST_Wave &psignal, float a)
|
|
{
|
|
int num_samples = signal.num_samples();
|
|
short last=0;
|
|
double maxval=0;
|
|
EST_TBuffer<double> fpdata(num_samples);
|
|
|
|
for(int i=0; i<num_samples; i++)
|
|
{
|
|
fpdata[i] = signal.a(i) - a*last;
|
|
last = signal.a(i);
|
|
if (absval(fpdata[i]) > maxval)
|
|
maxval = absval(fpdata[i]);
|
|
}
|
|
|
|
psignal.resize(num_samples,1,FALSE);
|
|
psignal.set_sample_rate(signal.sample_rate());
|
|
short_set(psignal, 0, fpdata, maxval);
|
|
}
|
|
|
|
/* redundant - see new version in filter */
|
|
void EST_post_deemphasis(EST_Wave &signal, EST_Wave &dsignal, float a)
|
|
{
|
|
int num_samples = signal.num_samples();
|
|
double last=0;
|
|
double maxval=0;
|
|
EST_TBuffer<double> fddata(num_samples);
|
|
|
|
for(int i=0; i<num_samples; i++)
|
|
{
|
|
fddata[i] = signal.a(i) + a*last;
|
|
last = fddata[i];
|
|
if (absval(fddata[i]) > maxval)
|
|
maxval = absval(fddata[i]);
|
|
}
|
|
|
|
dsignal.resize(num_samples,1,FALSE);
|
|
dsignal.set_sample_rate(signal.sample_rate());
|
|
short_set(dsignal, 0, fddata, maxval);
|
|
}
|
|
|
|
|
|
|