1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
// pitchtracker.h
/******************************************************************************
*
* MiXViews - an X window system based sound & data editor/processor
*
* Copyright (c) 1993, 1994 Regents of the University of California
*
* Author: Douglas Scott
* Date: December 13, 1994
*
* Permission to use, copy and modify this software and its documentation
* for research and/or educational purposes and without fee is hereby granted,
* provided that the above copyright notice appear in all copies and that
* both that copyright notice and this permission notice appear in
* supporting documentation. The author reserves the right to distribute this
* software and its documentation. The University of California and the author
* make no representations about the suitability of this software for any
* purpose, and in no event shall University of California be liable for any
* damage, loss of data, or profits resulting from its use.
* It is provided "as is" without express or implied warranty.
*
******************************************************************************/
// PitchTracker determines the frequency of a frame's worth of sound and puts
// the result into a PCHData object (via the 'frame' pointer).
#ifndef PITCHTRACKER_H
#ifdef __GNUG_
#pragma interface
#endif
#define PITCHTRACKER_H
#include "arrayfunction.h"
#include "processfun.h"
#include "range.h"
#define maxFrameSize 1010 /* maximum sample frame size */
#define MAXMM ((maxFrameSize / 10 + 1) / 2)
#define TWO_PI 6.28318530
class PitchTracker : public ArrayFunction {
typedef ArrayFunction Super;
friend class PitchTrackRequester;
typedef float P3DA5[50][5][MAXMM];
typedef float P2DA5[50][5];
typedef float P3DA6[50][6][MAXMM];
typedef float P2DA6[50][6];
typedef float PA5MAX[5][MAXMM];
typedef float PA6MAX[6][MAXMM];
public:
PitchTracker(Data* data, const Range& freqs, int framesize, int frmoffset);
PitchTracker(Data* data);
redefined const char* message() { return "Running pitchtrack: "; }
redefined int operator () (double *input, Data *frame);
void filter(double *array, int size=0);
protected:
redefined Requester* createRequester();
redefined void initialize();
redefined int doApply(Data *);
redefined int maxFrameLength() { return maxFrameSize; }
Modifier* createFilter(Data *);
double getPitch(double *array);
double getRMS(double *array);
private:
double search(double qsum, float *g, float *h);
void filfun(double qsum, float *fun, float *g, float *h);
double nyquist() { return sampRate()/2.0; }
void trigpo(double omega, PA5MAX phi, PA6MAX psi,
float *gampsi, float *gamphi, int n);
private:
float freq[50];
P3DA5 tphi;
P3DA6 tpsi;
P2DA5 tgamph;
P2DA6 tgamps;
Range freqRange;
int JMAX;
int MM;
enum FilterType { BandPass = 0x1, LowPass = 0x2 };
ChoiceValue filterType;
};
// this low-pass filters sound in preparation for pitch tracking
class PitchTrackFilter : public SimpleFunction {
public:
PitchTrackFilter(Data *d);
virtual ~PitchTrackFilter() {}
redefined const char* message() {
return "Running low pass filter...";
}
redefined double operator () (double input);
protected:
redefined void initialize();
private:
double w1,w11,w12,w2,w21,w22,w3,w31,w32,w4,w41,w42;
};
#endif
|
