# 5.07_lungsIntensity.draw # # Praat script that creates figure 5.9 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 Erase all Times 10 Read from file... 5.07_lungs.data # Subglottal pressure and its smoothed version Viewport... 0 6 0 3 Text top... no Subglottal pressure + smoothed version (Gaussian smoothing window: 10 ms) Text left... yes Raw subglottal pressure (kPa) Text bottom... yes Time (s) Draw inner box select Sound pressure30 Draw... 0 1 0 8100 no Marks left every... 1000 2 yes yes yes Marks bottom every... 1 0.1 yes yes no ! Smooth with a Hanning-equivalent width of 10 milliseconds. select Sound pressure30 Copy... smooth To Spectrum Formula... self*exp(-(x*0.010)^2) To Sound Draw... 0 1 -1999.9 6100 no Marks right every... 1000 2 yes yes yes ! Draw the smoothing window in an inset circle. Paint circle... 1 0.9 4500 0.07 Draw circle... 0.9 4500 0.07 Create Matrix... gauss 0 1 601 1e-4 0.87 1 1 1 1 1 0 To Sound (slice)... 1 Formula... 4000+1200*exp(-(pi*(x-0.9)/0.01)^2) Draw... 0 1 -1999.9 6100 no plus Matrix gauss Remove Viewport... 0 6 0 2.5 Text right... yes Smoothed subglottal pressure (kPa) Viewport... 0 6 2 3 Viewport text... Right Half 0 (a) # Acoustic output at one metre and its intensity. Viewport... 0 6 3 6 Text top... no Acoustic output at 1 metre + intensity (Gaussian smoothing window: 20 ms) Draw inner box Viewport... 0 6 3.7 6 Text left... 1 Sound pressure (Pa) select Sound a_speaker ; Convert from 20 cm to 1 m Formula... self/5 Draw... 0 1 -1 1 no Text bottom... yes Time (s) Text left... yes Sound pressure (Pa) Marks left every... 1 1 yes yes no Marks bottom every... 1 0.1 yes yes no ! Smooth with a Hanning-equivalent width of 20 milliseconds. Viewport... 0 6 3 5 select Sound a_speaker ; 3.2 / 0.020 = 160 To Intensity... 160 1e-3 Draw... 0 1 60 90 no Text right... yes Intensity (dB) One mark right... 60 yes yes no One mark right... 70 yes yes yes One mark right... 80 yes yes yes One mark right... 90 yes yes no Viewport text... Right Bottom 0 (b) ! Draw window. Paint circle... 1 0.9 80 0.07 Draw circle... 0.9 80 0.07 Create Matrix... gauss 0 1 141 1e-4 0.893 1 1 1 1 1 0 To Sound (slice)... 1 Formula... 86+20*log10(exp(-(pi*(x-0.9)/0.02)^2)) Draw... 0 1 60 90 no plus Matrix gauss Remove select Intensity a_speaker Draw... 0 1 60 90 no # Parametrized curve: intensity as a function of lung pressure Viewport... 0 6 6 9 Text top... no Intensity as a function of lung pressure Down to Matrix To Sound (slice)... 1 plus Sound smooth To ParamCurve Draw... 0 0 1e-3 0 6000 57 82 no Draw inner box Text bottom... yes Subglottal pressure (kPa) Marks bottom every... 1000 1 yes yes yes Text left... yes Intensity at 1 metre (dB %%re% 2\.c10^^-5^ Pa) Marks left every... 1 5 yes yes yes Viewport text... Right Half 0 (c) Viewport... 0 6 0.15 9