# 3.14.harmonic4.draw # # Praat script that creates figure 3.14 of # Boersma (1998): Functional Phonology. # http://fonsg3.hum.uva.nl/paul/ # Erase all Times 10 right = 4 Viewport... 0 2.9 0 3 Create Sound... implicit 0 4 500 2*arcsin(1/2*sqrt(x^2/(1+x^2))) Draw... 0 4 0 pi no Remove Create Sound... 2ndOrder 0 4 500 2*arcsin(1/2*sqrt(x^2/(1+1/2*x^2))) Draw... 0 4 0 pi no Remove Create Sound... explicit 0 2 500 2*arcsin(1/2*x) Draw... 0 4 0 pi no # "exact": Draw line... 0 0 pi pi Draw line... pi pi 4 2.28319 Draw inner box One mark left... 0 yes yes no One mark left... 1 yes yes no One mark left... 2 yes yes no One mark left... 3 no yes no One mark left... pi no yes no %\pi One mark bottom... 0 yes yes no One mark bottom... 1 yes yes no One mark bottom... 2 yes yes no One mark bottom... 3 no yes no One mark bottom... pi no yes yes %\pi One mark bottom... 4 yes yes no Text bottom... yes True frequency %\om\De%t y = 2*pi - right One mark right... y no yes no "exact" y = 2*arcsin(1/2*sqrt(right^2/(1+right^2))) One mark right... y no yes no "implicit" y = 2*arcsin(1/2*sqrt(right^2/(1+1/2*right^2))) One mark right... y no yes no "2nd-order" One mark top... 2 no yes no "explicit" Viewport... 0.2 2.9 0 3 Text left... yes Simulated frequency %\om\'p\De%t Viewport... 3.2 5.8 0 3 Text left... yes Amplitude Viewport... 2.9 5.8 0 3 Create Sound... implicit 0 4 500 sqrt((1+x^2)/(1+3/4*x^2)) Draw... 0 4 0 4 no Remove Create Sound... 2ndOrder 0 4 500 sqrt((1+1/2*x^2)/(1+1/4*x^2)) Draw... 0 4 0 4 no Remove Create Sound... exact 0 4 500 abs(1/(cos(1/2*x))) Draw clipped... 0 4 0 4 no Remove Create Sound... explicit 0 2 500 1/sqrt(1-1/4*x^2) Draw clipped... 0 4 0 4 no Remove Draw inner box Marks bottom every... 1 1 yes yes no Marks left... 2 yes yes no One mark left... 1 yes yes yes Text bottom... yes True frequency %\om\De%t y = sqrt((1+right^2)/(1+3/4*right^2)) One mark right... y no yes no "implicit" y = sqrt((1+1/2*right^2)/(1+1/4*right^2)) One mark right... y no yes no "2nd-order" One mark top... 1.95 no yes no "explicit" Text special... pi Centre 2.9 Half Times 10 0 "exact" Viewport... 0.2 6.15 0.2 3