During voiced speech, vocal folds interact with the 2 vocal tract acoustics. The resulting glottal source-3 resonator coupling has been observed using mathe-4 matical and physical models as well as in in vivo 5 phonation. We propose a computational time-domain 6 model of the full speech apparatus that contains a 7 feedback mechanism from the vocal tract acoustics 8 to the vocal fold oscillations. It is based on nu-9 merical solution of ordinary and partial differential 10 equations defined on vocal tract geometries that have 11 been obtained by magnetic resonance imaging. The 12 model is used to simulate rising and falling pitch glides 13 of [A, i] in the fundamental frequency (f o) interval 14 [145 Hz, 315 Hz]. The interval contains the first vo-15 cal tract resonance f R1 and the first formant F 1 of [i] 16 as well as the fractions of the first resonance f R1 /5, 17 f R1 /4, and f R1 /3 of [A]. The glide simulations reveal 18 a locking pattern in the f o trajectory approximately