Abstract
In this paper, we propose a numerical method for the voice generation process based on some mathematical model. As the model, we use Webster's horn equation for the 1D case. We then discretize this equation by FEM and calculate the frequency response function. We consider the complex eigenvalue problem corresponding to Webster's horn equation and give a variational formula for the complex eigenvalue with respect to the variation of vocal tract shape. We numerically confirm that the complex eigenvalue is closely related to the frequency response function, and propose a vocal tract shape design algorithm using variational formula and confirm its efficiency numerically.
