# 5.07_lungsPitch.draw
#
# Praat script that creates figure 5.10 of
#    Boersma (1998): Functional Phonology, section 5.7.
# http://fonsg3.hum.uva.nl/paul/
#
# Expects the file 5.07_lungs.data in the same directory;
# that file can be created by 5.07_lungs.make

;Read from file... 5.07_lungs.data
# Run this script after 5.07_lungsIntensity.draw, so you'll have a Sound called "smooth".
Erase all
Times
10

# Pitch contour.

select Sound a_speaker
; Compute signal at 1 metre instead of 20 cm.
Formula... self/5
To Pitch (ac)... 0.001 100 15 no 0.06667 0.75 0.01 0.35 0.14 600
Viewport... 0 6 0 2
Viewport text... Right Half 0 (a)
Viewport... 0 6.3 0 2
Text top... no Pitch contour
Draw... 0 1 100 300 no
Marks bottom every... 1e-3 100 yes yes no
Marks left every... 1 50 yes yes no
Draw inner box
Text bottom... yes Time (ms)
Text left... yes Pitch (Hz)

# Zoom in on the three periodic parts.

select Sound a_speaker

Viewport... 0 2.7 2 3.5
Viewport text... Right Half 0 (b)
Viewport... 0 3 2 3.5
Text top... 0 Acoustic signal at 1 metre
Draw... 0.19 0.21 -1 1 no
Draw inner box
Text bottom... yes Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... yes Sound pressure (Pa)
Marks left every... 1 1 yes yes yes

Viewport... 0 2.7 3.5 5
Viewport text... Right Half 0 (c)
Viewport... 0 3 3.5 5
Text top... 0 Acoustic signal at 1 metre
Draw... 0.39 0.41 -1 1 no
Draw inner box
Text bottom... 1 Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... 1 Sound pressure (Pa)
Marks left every... 1 1 yes yes yes

Viewport... 0 2.7 5 6.5
Viewport text... Right Half 0 (d)
Viewport... 0 3 5 6.5
Text top... no Acoustic signal at 1 metre
Draw... 0.49 0.51 -1 1 no
Draw inner box
Text bottom... yes Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... yes Sound pressure (Pa)
Marks left every... 1 1 yes yes yes

# Vocal cord motions.

Viewport... 3 6 2 3.5
Viewport text... Right Half 0 (e)
Viewport... 3.3 6.3 2 3.5
Text top... no Regular periodic vocal-fold motion
select Sound width37
Draw... 0.19 0.21 -1e-3 2e-3 no
Dotted line
select Sound width36
Draw... 0.19 0.21 -1e-3 2e-3 no
Plain line
Draw inner box
Text bottom... yes Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... 1 Width (mm)
One mark left... -1e-3 no yes no -1
One mark left... 0 no yes yes 0
One mark left... 2e-3 no yes no 2

Viewport... 3 6 3.5 5
Viewport text... Right Half 0 (f)
Viewport... 3.3 6.3 3.5 5
Text top... no Aperiodic vocal-fold motion
select Sound width37
Draw... 0.39 0.41 -1e-3 2e-3 no
Dotted line
select Sound width36
Draw... 0.39 0.41 -1e-3 2e-3 no
Plain line
Draw inner box
Text bottom... yes Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... yes Width (mm)
One mark left... -1e-3 no yes no -1
One mark left... 0 no yes yes 0
One mark left... 2e-3 no yes no 2

Viewport... 3 6 5 6.5
Viewport text... Right Half 0 (g)
Viewport... 3.3 6.3 5 6.5
Text top... no Irregular periodic vocal-fold motion
select Sound width37
Draw... 0.49 0.51 -1e-3 2e-3 no
Dotted line
select Sound width36
Draw... 0.49 0.51 -1e-3 2e-3 no
Plain line
Draw inner box
Text bottom... yes Time (ms)
Marks bottom every... 0.001 5 yes no no
Marks bottom every... 0.001 1 no yes no
Text left... 1 Width (mm)
One mark left... -1e-3 no yes no -1
One mark left... 0 no yes yes 0
One mark left... 2e-3 no yes no 2

# Pitch as a function of lung pressure.

Viewport... 0 6 6.5 9
Viewport text... Right Half 0 (h)
Viewport... 0 6.3 6.5 9
Text top... no Fundamental frequency as a `function' of subglottal pressure
select Pitch a_speaker
To Matrix
Create Matrix... poly 0 1 1000 1e-3 0.5e-3 1 2 2 1 1 0
Formula... if row=1 then Sound_smooth(x) else Matrix_a_speaker(x,nearest 1) fi
Formula... if row=2 and self<100 then 100 else self fi
To Polygon
Draw... 0 6000 100 300
Draw inner box
Text bottom... 1 Subglottal pressure (kPa)
Marks bottom every... 1000 1 yes yes yes
Text left... 1 Pitch (Hz)
Marks left every... 1 50 yes yes yes

Viewport... 0 6 0.15 9