Model-based classification of nonstationary vocal fold vibrations

Classification of vocal fold vibrations is an essential task of the objective assessment of voice disorders. For historical reasons, the conventional clinical examination of vocal fold vibrations is done during stationary, sustained phonation. However, the conclusions drawn from a stationary phonation are restricted to the observed steady-state vocal fold vibrations and cannot be generalized to voice mechanisms during running speech. This study addresses the approach of classifying real-time recordings of vocal fold oscillations during a nonstationary phonation paradigm in the form of a pitch raise. The classification is based on asymmetry measures derived from a time-dependent biomechanical two-mass model of the vocal folds which is adapted to observed vocal fold motion curves with an optimization procedure. After verification of the algorithm performance the method was applied to clinical problems. Recordings of ten subjects with normal voice and ten dysphonic subjects have been evaluated during stationary as well as nonstationary phonation. In the case of nonstationary phonation the model-based classification into "normal" and "dysphonic" succeeds in all cases, while it fails in the case of sustained phonation. The nonstationary vocal fold vibrations contain additional information about vocal fold irregularities, which are needed for an objective interpretation and classification of voice disorders.     

Computation of physiological human vocal fold parameters by mathematical optimization of a biomechanical model

With the use of an endoscopic, high-speed camera, vocal fold dynamics may be observed clinically during phonation. However, observation and subjective judgment alone may be insufficient for clinical diagnosis and documentation of improved vocal function, especially when the laryngeal disease lacks any clear morphological presentation. In this study, biomechanical parameters of the vocal folds are computed by adjusting the corresponding parameters of a three-dimensional model until the dynamics of both systems are similar. First, a mathematical optimization method is presented. Next, model parameters (such as pressure, tension and masses) are adjusted to reproduce vocal fold dynamics, and the deduced parameters are physiologically interpreted. Various combinations of global and local optimization techniques are attempted. Evaluation of the optimization procedure is performed using 50 synthetically generated data sets. The results show sufficient reliability, including 0.07 normalized error, 96% correlation, and 91% accuracy. The technique is also demonstrated on data from human hemilarynx experiments, in which a low normalized error (0.16) and high correlation (84%) values were achieved. In the future, this technique may be applied to clinical high-speed images, yielding objective measures with which to document improved vocal function of patients with voice disorders.     

Three-dimensional biomechanical properties of human vocal folds: parameter optimization of a numerical model to match in vitro dynamics

The human voice signal originates from the vibrations of the two vocal folds within the larynx. The interactions of several intrinsic laryngeal muscles adduct and shape the vocal folds to facilitate vibration in response to airflow. Three-dimensional vocal fold dynamics are extracted from in vitro hemilarynx experiments and fitted by a numerical three-dimensional-multi-mass-model (3DM) using an optimization procedure. In this work, the 3DM dynamics are optimized over 24 experimental data sets to estimate biomechanical vocal fold properties during phonation. Accuracy of the optimization is verified by low normalized error (0.13 ± 0.02), high correlation (83% ± 2%), and reproducible subglottal pressure values. The optimized, 3DM parameters yielded biomechanical variations in tissue properties along the vocal fold surface, including variations in both the local mass and stiffness of vocal folds. That is, both mass and stiffness increased along the superior-to-inferior direction. These variations were statistically analyzed under different experimental conditions (e.g., an increase in tension as a function of vocal fold elongation and an increase in stiffness and a decrease in mass as a function of glottal airflow). The study showed that physiologically relevant vocal fold tissue properties, which cannot be directly measured during in vivo human phonation, can be captured using this 3D-modeling technique.     

Visualization and Quantification of the Medial Surface Dynamics of an Excised Human Vocal Fold During Phonation

SUMMARY: The purpose of this investigation was to investigate physical mechanisms of vocal fold vibration during normal phonation through quantification of the medial surface dynamics of the fold. An excised hemilarynx setup was used. The dynamics of 30 microsutures mounted on the medial surface of a human vocal fold were analyzed across 18 phonatory conditions. The vibrations were recorded with a digital high-speed camera at a frequency of 4,000 Hz. The positions of the sutures were extracted and converted to three-dimensional coordinates using a linear approximation technique. The data were reduced to principal eigenfuctions, which captured over 90% of the variance of the data, and suggested mechanisms of sustained vocal fold oscillation. The vibrations were imaged as the following phonatory conditions were manipulated: glottal airflow, an adductory force applied to the muscular process, and an elongation force applied to the thyroid cartilage. Over the range of variables studied, only the variation in glottal airflow yielded significant changes in subglottal pressure and fundamental frequency. All recordings showed high correlation for the distribution of the dynamics across the medial surface of the vocal fold. The distribution of the different displacement directions and velocities showed the highest variations around the superior region of the medial surface. Although the computed vibration patterns of the two largest empirical eigenfunctions were consistent with previous experimental observations, the relative prominence of the two eigenfunctions changed as a function of glottal airflow, impacting theories of vocal efficiency and vocal economy.     

Multiparametric Analysis of Vocal Fold Vibrations in Healthy and Disordered Voices in High-Speed Imaging

Objectives

The aim of this study was to look for visual subjective and objective parameters of vocal fold dynamics being capable of differentiating healthy from pathologic voices in daily clinical practice applying endoscopic high-speed digital imaging (HSI).

Study Design and Methods

Four hundred ninety-six datasets containing 80 healthy and 416 pathologic subjects (232 functional dysphonia (FD), 13 bilateral, and 171 unilateral vocal fold nerve paralysis) were analyzed retrospectively. Videos at 4000 Hz (256 × 256 pixel) were recorded during sustained phonation. Subjective parameters were visually evaluated and complemented by an analysis of objective parameters. Visual subjective parameters were mucosal wave, glottal closure type, glottal closure insufficiency (GI), asymmetries of the vocal folds, and phonovibrogram (PVG) symmetry. After image segmentation, objective parameters were computed: closed quotient, perturbation measures (PMs) of glottal area, and left-right asymmetry values.

Results

HSI evaluation enabled to distinguish healthy from pathologic voices. For visual subjective parameters, GI, symmetrical behavior, and PVG symmetry exhibited statistical significant differences. For 95% of the data, objective parameters could be computed. Among objective parameters, closed quotient, jitter, shimmer, harmonic-to-noise ratio, and signal-to-noise ratio for the glottal area function differentiated statistically significant normal from pathologic voices. Applying linear discriminant analysis by combining visual subjective and objective parameters, accurate classifications were made for 63.2% of the female and 87.5% of the male group for the three-class problem (healthy, FD, and unilateral vocal fold nerve paralysis).

Conclusion

Actual acoustically applied PMs can be transferred to clinical beneficial HSI analysis. Combining visual subjective and objective basic parameters succeeds in differentiating pathologic from healthy voices. The presented evaluation can easily be included into everyday clinical practice. However, further research is needed to broaden our understanding of the variability within and across healthy and pathologic vocal fold vibrations for diagnosing voice disorders and therapy control.     

Spatio-temporal analysis of irregular vocal fold oscillations: biphonation due to desynchronization of spatial modes.

This report is on direct observation and modal analysis of irregular spatio-temporal vibration patterns of vocal fold pathologies in vivo. The observed oscillation patterns are described quantitatively with multiline kymograms, spectral analysis, and spatio-temporal plots. The complex spatio-temporal vibration patterns are decomposed by empirical orthogonal functions into independent vibratory modes. It is shown quantitatively that biphonation can be induced either by left-right asymmetry or by desynchronized anterior-posterior vibratory modes, and the term "AP (anterior-posterior) biphonation" is introduced. The presented phonation examples show that for normal phonation the first two modes sufficiently explain the glottal dynamics. The spatio-temporal oscillation pattern associated with biphonation due to left-right asymmetry can be explained by the first three modes. Higher-order modes are required to describe the pattern for biphonation induced by anterior-posterior vibrations. Spatial irregularity is quantified by an entropy measure, which is significantly higher for irregular phonation than for normal phonation. Two asymmetry measures are introduced: the left-right asymmetry and the anterior-posterior asymmetry, as the ratios of the fundamental frequencies of left and right vocal fold and of anterior-posterior modes, respectively. These quantities clearly differentiate between left-right biphonation and anterior-posterior biphonation. This paper proposes methods to analyze quantitatively irregular vocal fold contour patterns in vivo and complements previous findings of desynchronization of vibration modes in computer modes and in in vitro experiments.     

Material parameter computation for multi-layered vocal fold models

Today, the prevention and treatment of voice disorders is an ever-increasing health concern. Since many occupations rely on verbal communication, vocal health is necessary just to maintain one's livelihood. Commonly applied models to study vocal fold vibrations and air flow distributions are self sustained physical models of the larynx composed of artificial silicone vocal folds. Choosing appropriate mechanical parameters for these vocal fold models while considering simplifications due to manufacturing restrictions is difficult but crucial for achieving realistic behavior. In the present work, a combination of experimental and numerical approaches to compute material parameters for synthetic vocal fold models is presented. The material parameters are derived from deformation behaviors of excised human larynges. The resulting deformations are used as reference displacements for a tracking functional to be optimized. Material optimization was applied to three-dimensional vocal fold models based on isotropic and transverse-isotropic material laws, considering both a layered model with homogeneous material properties on each layer and an inhomogeneous model. The best results exhibited a transversal-isotropic inhomogeneous (i.e., not producible) model. For the homogeneous model (three layers), the transversal-isotropic material parameters were also computed for each layer yielding deformations similar to the measured human vocal fold deformations.     

Nonlinear dynamics of the voice: Signal analysis and biomechanical modeling

Irregularities in voiced speech are often observed as a consequence of vocal fold lesions, paralyses, and other pathological conditions. Many of these instabilities are related to the intrinsic nonlinearities in the vibrations of the vocal folds. In this paper, bifurcations in voice signals are analyzed using narrow-band spectrograms. We study sustained phonation of patients with laryngeal paralysis and data from an excised larynx experiment. These spectrograms are compared with computer simulations of an asymmetric 2-mass model of the vocal folds. (c) 1995 American Institute of Physics.     

Physical mechanisms of phonation onset: a linear stability analysis of an aeroelastic continuum model of phonation

In an investigation of phonation onset, a linear stability analysis was performed on a two-dimensional, aeroelastic, continuum model of phonation. The model consisted of a vocal fold-shaped constriction situated in a rigid pipe coupled to a potential flow which separated at the superior edge of the vocal fold. The vocal fold constriction was modeled as a plane-strain linear elastic layer. The dominant eigenvalues and eigenmodes of the fluid-structure-interaction system were investigated as a function of glottal airflow. To investigate specific aerodynamic mechanisms of phonation onset, individual components of the glottal airflow (e.g., flow-induced stiffness, inertia, and damping) were systematically added to the driving force. The investigations suggested that flow-induced stiffness was the primary mechanism of phonation onset, involving the synchronization of two structural eigenmodes. Only under conditions of negligible structural damping and a restricted set of vocal fold geometries did flow-induced damping become the primary mechanism of phonation onset. However, for moderate to high structural damping and a more generalized set of vocal fold geometries, flow-induced stiffness remained the primary mechanism of phonation onset.     

Mechanisms of irregular vibration in a physical model of the vocal folds

Previous investigations have shown that one mechanism of irregular vocal fold vibration may be a desynchronization of two or more vibratory modes of the vocal fold tissues. In the current investigation, mechanisms of irregular vibration were further examined using a self-oscillating, physical model of vocal fold vibration, a hemi-model methodology, and high-speed, stereoscopic, digital imaging. Using the method of empirical eigen-functions, a spatiotemporal analysis revealed mechanisms of irregular vibration in subharmonic phonation and biphonation, which were not disclosed in a standard acoustic spectrum.     

Vibration in a self-oscillating vocal fold model with left-right asymmetry in body-layer stiffness

This study compares the phonatory behavior of an asymmetric vocal fold model to that of each individual vocal fold model in a hemi-configuration. Although phonation frequencies of the two folds in hemi-configurations had a ratio close to 1:3, a subharmonic synchronization between the two folds was not observed in the asymmetric model. Instead, the vibratory behavior was dominated by the dynamics of one fold only, and the other fold was enslaved to vibrate at the same frequency. Increasing subglottal pressure induced a shift in relative dominance between the two folds, leading to abrupt changes in both vibratory pattern and frequency.     

Vocal fold and ventricular fold vibration in period-doubling phonation: physiological description and aerodynamic modeling

Occurrences of period-doubling are found in human phonation, in particular for pathological and some singing phonations such as Sardinian A Tenore Bassu vocal performance. The combined vibration of the vocal folds and the ventricular folds has been observed during the production of such low pitch bass-type sound. The present study aims to characterize the physiological correlates of this acoustical production and to provide a better understanding of the physical interaction between ventricular fold vibration and vocal fold self-sustained oscillation. The vibratory properties of the vocal folds and the ventricular folds during phonation produced by a professional singer are analyzed by means of acoustical and electroglottographic signals and by synchronized glottal images obtained by high-speed cinematography. The periodic variation in glottal cycle duration and the effect of ventricular fold closing on glottal closing time are demonstrated. Using the detected glottal and ventricular areas, the aerodynamic behavior of the laryngeal system is simulated using a simplified physical modeling previously validated in vitro using a larynx replica. An estimate of the ventricular aperture extracted from the in vivo data allows a theoretical prediction of the glottal aperture. The in vivo measurements of the glottal aperture are then compared to the simulated estimations.     

Bifurcations and chaos in register transitions of excised larynx experiments

Experimental data from an excised larynx are analyzed in the light of nonlinear dynamics. The excised larynx provides an experimental framework that enables artificial control and direct observation of the vocal fold vibrations. Of particular interest in this experiment is the coexistence of two distinct vibration patterns, which closely resemble chest and falsetto registers of the human voice. Abrupt transitions between the two registers are typically accompanied by irregular vibrations. Two approaches are presented for the modeling of the excised larynx experiment; one is the nonlinear predictive modeling of the experimental time series and the other is the biomechanical modeling (three-mass model) that takes into account basic mechanisms of the vocal fold vibrations. The two approaches show that the chest and falsetto vibrations correspond to two coexisting limit cycles, which jump to each other with a change in the bifurcation parameter. Irregular vibrations observed at the register jumps are due to chaos that exists near the two limit cycles. This provides an alternative mechanism to generate chaotic vibrations in excised larynx experiment, which is different from the conventionally known mechanisms such as strong asymmetry between the left and right vocal folds or excessively high subglottal pressure.     

Inverse filtering

Inverse filtering is a noninvasive method of producing a glottogram thought to reflect the vibratory motions of the vocal fold. The flow glottogram provides information related to the type of phonation, the sound pressure level, the regularity of vocal fold vibrations, and to the presence or absence of vocal fold closure. Limited data, speech samples, filter tuning, and lack of unanimity on waveform display and interpretation have contributed to slow application of the techniques to the clinical population. The author argues the technique has utility in both diagnosis and treatment     

Irregular vocal-fold vibration--high-speed observation and modeling

Direct observations of nonstationary asymmetric vocal-fold oscillations are reported. Complex time series of the left and the right vocal-fold vibrations are extracted from digital high-speed image sequences separately. The dynamics of the corresponding high-speed glottograms reveals transitions between low-dimensional attractors such as subharmonic and quasiperiodic oscillations. The spectral components of either oscillation are given by positive linear combinations of two fundamental frequencies. Their ratio is determined from the high-speed sequences and is used as a parameter of laryngeal asymmetry in model calculations. The parameters of a simplified asymmetric two-mass model of the larynx are preset by using experimental data. Its bifurcation structure is explored in order to fit simulations to the observed time series. Appropriate parameter settings allow the reproduction of time series and differentiated amplitude contours with quantitative agreement. In particular, several phase-locked episodes ranging from 4:5 to 2:3 rhythms are generated realistically with the model.     

Spatial arrangement of the structural elements ofvocal fold layers: An adjustment to the vibration process

It is well established that the multilayered structure of the vocal fold is highly adjusted to the requirements of the vibration process during phonation. There is also some partial data indicating that the spatial arrangement of each vocal fold layer corresponds to the functional requirements, and thus facilitate the phonation process. Nevertheless, all reports on the spatial arrangement of the vocal fold structures deal only with an individual element of the vocal fold histologic structure. The present study encompasses the spatial histologic analysis of all major elements of the vocal fold layers. It was demonstrated that the vocal fold epithelial cells, the connective and muscle fibers, and even the blood vessels run parallel to the vocal fold free edge, which indicates a high adjustment to the phonation requirements and the vibration process.