# 5.
w = read.emusegs("H:/segfeb/sb.txt")

# 6a
dur = mudur(vowlax)

# 6b
dur[380:395]

# 6c
median(dur)

# 6d
n = nrow(vowlax)
dur[1:(n-20)]

# 7
v50 = dcut(vowlax.fdat, .5, prop=T)

# 7a
v50[100:130, c(1, 3)]


# 7b
v50[c(10, 15, 120),c(1,3)]

# 7c
b = bark(v50)
sum(b[,3] - b[,2] > 3.5)

# 8a
v50[vowlax.l=="E",2] > 1950

# 8b
mean(v50[vowlax.l=="a",2])


# 8c
sd(v50[vowlax.l %in% c("I", "E"),2])

# 9a
o = dcut(vowlax.grund, 1, prop=T)
o[vowlax.l=="a" & vowlax.spkr=="67"]

# 9b
vowlax.fund[vowlax.l %in% c("I", "E"),]

# 9c
o = dcut(vowlax.fdat, 0, prop=T)
o[vowlax.left=="d" & vowlax.spkr=="67" & vowlax.l %in% c("E", "a"),c(1,3)]

# 10
m = dcut(vowlax.fdat[,1:2], .5, prop=T)
temp = vowlax.l=="I"
eplot(m[temp,], vowlax.spkr[temp], form=T, dopoints=T)

# 11 a
fs = 2 * max(trackfreq(sib.dft))
fs


# 11b
N = 2 * (ncol(sib.dft) - 1)
N/(fs/1000)


# 12
s5 = dcut(sib.dft, .5, prop=T)
which = sib.w=="besucht"
s5[which,2000:4000]

# 13
plot(s5[,1000:7500], sib.l, fun="mean", power=T)

# 14
Lippenrundung verursacht eine groessere Vokaltraktlaenge
und daher eine allgemeine Frequenz-Senkung der
Energie im Spektrum: das wird wahrscheinlich der
Grund sein, weshalb die Frequenz des spektralen Gipfels 
vor gerundeten Vokalen bei ca. 4800 Hz niedriger
ist wie vor ungerundeten Vokalen (ca. 6000 Hz).