hmm

sim/genFn.R

@@ -24,10 +24,14 @@ genFn <- function(d){
                 d$p$pure[i,]  <- d$p$d[i] * idx + d$p$start
             }else{
                 d$t$pure[i,]  <- d$t$pure[i-1, d$M] + d$t$between + d$t$d[i] * idx
-                d$p$pure[i,]  <- d$p$pure[i-1, d$M] + d$t$between *d$p$mmean  + d$p$D + d$p$d * idx
+                d$p$pure[i,]  <- d$p$pure[i-1, d$M] + ## Anfangsdruck = Enddruck letzter SZ
+                  d$t$between *d$m$mean  + ## Druck fällt während Antrieb weiter
+                    d$p$D + ## ein Hub nach oben
+                      d$p$mmean * (d$t$T + d$t$between) +
+                      d$p$d * idx ## drockabfallvektor
             }
                         
-            d$p$rand[i,]      <- round(d$p$pure[i,] + rnorm(d$M, 0, d$p$randSd),d$p$round)
+            d$p$rand[i,]      <- round(d$p$pure[i,] + rnorm(d$M, 0, d$p$randSd), d$p$round)
             d$t$rand[i,]      <- d$t$pure[i,] + rnorm(d$M, 0, d$t$randSd)
 
             temp.lm           <- lm( d$p$rand[i,] ~ d$t$rand[i,])

sim/plotSz.R

@@ -9,8 +9,10 @@ plotSz <- function(d){
       geom_point(aes(x = d.t.rand, y = d.p.rand), color='red') +
         geom_abline(aes(intercept = d$p$randIntercept, slope = d$m$rand), color='red') +
           geom_vline(xintercept = d$t$rand0, linetype = "longdash") +
-            geom_point(aes(y = d$p$randMean, x = d$t$randMean), color="black", size=4) +
-              ggtitle("dp vs. dt")
+            geom_vline(xintercept = d$t$pure[,d$M], linetype = "dotted", color="blue") +
+              geom_hline(yintercept = d$p$pure[,d$M], linetype = "dotted", color="blue") +
+              geom_point(aes(y = d$p$randMean, x = d$t$randMean), color="black", size=4) +
+                ggtitle("dp vs. dt")
 
 
 }