Flute acoustics: measurement, modelling and design - School of ...

262 APPENDIX B. PROGRAM LISTINGS
maxf = p->x;
maxZ = p->y;
}
}
return maxf;
}
void harmonicity(Vector minv) {
int i;
Vector harmv;
/* for each Minimum in vector... */
for(i = 0; i < sizeVector(minv); i++) {
/* find all the harmonic minima of the particular Minimum */
harmv = harmonics(minv, i);
/* record the number of harmonics the Minimum has */
((Minimum)elementAt(minv, i))->numharm =
(double)sizeVector(harmv);
/* calculate and record the weighted average impedance of
these harmonic minima. */
((Minimum)elementAt(minv, i))->meanharmZ = harmAvgZ(harmv);
}
return;
}
Vector harmonics(Vector minv, int pos) {
int i;
int harmonic;
double window, left_bound, right_bound;
double df;
int curharmonic;
int found;
/* get the ratios of minima frequencies following
the Minimum in question */
double* ha = harmonicsArray(minv, pos);
int arraysize = sizeVector(minv) - pos - 1;
Vector harmv = createVector();
Harmonic h;
i = 1;
/* search for possible integer harmonics until
last Minimum is reached */
for(harmonic = 2; ; harmonic++) {
i--;
/* calculate a window in which harmonic will be defined.
The widest window allowable is 0.5. */
window = (HARMONIC_WINDOW/100.0)*harmonic;
if(window > 0.5)
window = 0.5;
left_bound = harmonic - window;
right_bound = harmonic + window;
found = 0;
/* search for any Minimum which exists in the defined window.
Compare this to previously found minima in the same