Skip to content
Snippets Groups Projects
Commit cf2dcb66 authored by wactbprot's avatar wactbprot
Browse files

array mit neuen p; hunger

parent 67b460ab
No related branches found
No related tags found
No related merge requests found
p.fill <- 44.6 # mbar
N <- 4
dp <- p.fill/10
tms <- c(
1366810108219.4746,
1366810145154.9,
1366810182083.075,
1366810219633.6438,
1366810256755.1125
) # ms
ts <- (tms - min(tms))/1000 # s
p <- c(
-0.016786,
0.008593,
0.034127,
0.058833,
0.08419
) # mbar
m <- c(-0.0000018029,
-0.0000018078,
-0.0000018123,
-0.0000018144,
-0.000001821
)*1000 # mbar/s
## p = m*t + b
## b = p - m*t
b <- p - m*ts # mbar
mp <- mean(p)
plot( ts, p,
xlim=c(-50,200),
ylim=c(-0.05,0.1),
xlab="t in s",
ylab="p in mbar")
grid()
abline(h=mp, col=2)
sapply(1:length(b),
function(i){
abline(a=b[i], b=m[i], col=3)
})
## finde t0s:
## mp = m * t0 - b
## t0 = (mp - b)/m
t0 <- (mp -b) /m
abline(v=t0, col=4)
dt <- diff(t0)
## ins Blaue:
## Korrektur:
## sin(winkel_p_anstieg)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment