ich habe einige Daten, von denen ich den Frequenzbereich ermitteln möchte. Die folgende Funktion habe ich dafür gefunden, der Code funktioniert soweit auch.
Code: Alles auswählen
d <- c(-0.538, -0.785, -1.123, -1.062, -1.031, -0.877, -0.846, -0.846, -0.877,-0.692,-0.508,-0.508,-0.538,-0.569,-0.538,-0.508,-0.600,-0.754,
-0.815,-0.877,-1.000,-1.062,-1.092,-1.123,-1.123,-1.369,-1.585,-1.154,-1.000,-1.062,-1.031,-1.092,-0.969,-0.938,-1.000,-0.938,-0.877,-0.723,
-0.538,-0.385,-0.415,-0.415,-0.600,-0.446,-0.477,-0.600,-0.785,-0.877,-1.000,-1.215,-1.369,-1.369,-1.369,-1.277,-1.862,-1.400,-0.938,-0.662,
-0.815,-1.215,-1.062,-1.062,-0.846,-0.846,-0.908,-0.723,-0.508,-0.385,-0.477,-0.508,-0.508,-0.477,-0.508,-0.569,-0.754,-0.969,-1.092,-1.092)
e <- c(0.344,0.313,0.313,0.219,0.063,0.000,-0.031,-0.094,-0.063,0.000,0.031,0.156,0.094,0.156,0.125,0.156,0.313,0.719,0.500,0.438,0.156, 0.125,0.031,-0.031,0.156,0.469,0.469,0.344,0.219,0.094,0.125,0.281,0.219,0.219,0.344,0.281,0.281,0.219,0.031,0.063,0.188,0.313,0.375,0.313,
0.219,0.063,0.094,0.219,0.406,0.594,0.500,0.313,0.344,0.281,0.219,0.281,0.281,0.281,0.250,0.063,0.063,0.000,-0.031,0.000,0.000,0.094,0.094,
0.125,0.219,0.313,0.344,0.563,0.500,0.094,-0.125,0.125,0.594,0.219)
res <- sqrt(d^2 + e^2)
xraw.spec <- spec.pgram(res, taper = 0)
plot(xraw.spec, log = "no")
Vielen Dank vorab für die Hilfe
Marge