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

APPENDIX B. PROGRAM LISTINGS 193
*/
The complex radiation impedance.
complex flangedZ(double f, double c, double rho, double a);
/*
Calculates the radiation impedance of a flanged pipe.
Refer to Dalmont (2001) JSV 244, 505-34.
Parameters:
f: the frequency.
c: the speed of sound.
rho: the density of air.
a: the radius of the tube in metres.
Returns:
The complex radiation impedance.
*/
TransferMatrix tubeMatrix(double f, double c, double rho, double L,
double a, double alphacorrection);
/*
Calculates the transfer matrix of a cylindrical tube.
The transfer matrix T relates pressure and flow at the input
to the pressure and flow at the output. Flow is positive
into the input and out of the output.
Parameters:
f: the frequency.
c: the speed of sound.
rho: the density of air.
L: the length of the tube in metres.
a: the radius of the tube in metres.
alphacorrection: the multiplicative attenuation coefficient
factor.
Returns:
The transfer matrix of the tube.
*/
TransferMatrix coneMatrix(double f, double c, double rho, double L,
double a1, double a2, double alphacorrection);
/*
Calculates the transfer matrix of a truncated cone.
The transfer matrix T relates pressure and flow at the input
to the pressure and flow at the output. Flow is positive
into the input and out of the output.
Parameters:
f: the frequency.
c: the speed of sound.
rho: the density of air.
L: the length of the pipe in metres.
a1: the radius of the pipe at input in metres.
a2: the radius of the pipe at output in metres.
alphacorrection: the multiplicative attenuation coefficient
factor.
Returns:
The transfer matrix of the pipe.
*/
TransferMatrix discontinuityMatrix(double f, double c, double rho,
double a1, double a2);
/*