# 16-Punkt Cosinuswelle, eine Schwingung
crplot()

# Dasselbe mit Amplitude 1.5
crplot(A=1.5)

# 16-Punkt Cosinuswelle, 2 Schwingungen
crplot(k=2)

# 16-Punkt Sinuswelle (Cosinuswelle mit Phase um 180 Grad verschoben)
crplot(p=-pi/2)

# 8 Punkt Cosinuswelle, eine Schwingung
crplot(N=8)

# Die Faltung
# 16 Punkt Cosinuswelle, 3 Schwingungen.
crplot(k=3, N=16)
# 16 Punkt Cosinuswelle, 13 Schwingungen.
crplot(k=13, N=16)


# Segmentliste aller [c, x] Frikative dieser Datenbank
dor = 

# Label-Vektor
dor.l =

# Spektra fuer die Segmentliste
dor.dft =


# Spektra zum zeitlichen Mittelpunkt
dor.dft5 = 


# Anzahl der Spektralkomponente
ncol(dor.dft)
ncol(dor.dft5)

# Dauer (Punkte) des Fensters, mit der die Spektra berechnet wurden
N = 2 * (ncol(dor.dft)-1) 

# Die Frequenzen, zu denen diese Spektralkomponente vorkommen 
trackfreq(dor.dft)

# Die Höchstfrequenz
max(trackfreq(dor.dft))

# Daher die  Abtastrate
fs =

# Hz-Abstand zwischen den Spektralkomponenten
fs/N

# ms-Dauer des Fensters, mit der die Spektra berechnet wurden
1000*N/fs

# Abbildungen
# Alle Spektra (zum zeitlichen Mittelpunkt)
plot(dor.dft5)

# Alle Spektra nach Etikettierung kodiert
plot(dor.dft5, dor.l)

# Mittelwert pro Kategorie (ensemble-averaged spectra)
plot(dor.dft5, dor.l, fun=mean, dbnorm=T)


# Spektra zum zeitlichen Mittelpunkt zwischen 1000-2000 Hz
dor.dft5[,1000:2000]

# das gleiche, Segmente 4, 5, 7 Frequenzen 2000-2500 Hz
dor.dft5[c(4,5,7), 2000:2500]

# Abbildung, Spektra, pro Kategorie gemittelt, 500 – 5000 Hz
plot(dor.dft5[,500:5000], dor.l, fun=mean, dbnorm=T)