From 746bf195c3f03a51503677262c0df258ae4f0020 Mon Sep 17 00:00:00 2001
From: wactbprot <thsteinbock@web.de>
Date: Sun, 18 Aug 2013 19:41:07 +0200
Subject: [PATCH] und jetzt alles als list() und ab in eine db
---
sim/theo.R | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/sim/theo.R b/sim/theo.R
index 2bae4cc..ccc4719 100644
--- a/sim/theo.R
+++ b/sim/theo.R
@@ -1,12 +1,13 @@
library(ggplot2)
### --------------------------------------------------------
+set <- list()
p.fill <- 44.6 # mbar
M <- 100 # Anz. Messpunkte pro SZ
N <- 5 # Anz. der SZ
p.mean <- 0.0 # mbar
m.mean <- -1.8e-3 # mbar/s
-t.between <- 5.0
-t.rand.sd.rel <- 1e-3
+t.between <- 5.0 # s
+t.rand.sd.rel <- 1e-3 #
p.rand.sd.rel <- 1e-8
p.mean.sz <- rep(p.mean, N) # Druckmittelwert der SZ
m.sz <- rep(m.mean, N) # Ausgangssteigung der SZ
@@ -62,7 +63,7 @@ for(i in 1:N){
c.i <- p.rand.mean - m.rand * t.rand.mean
t.rand.0 <- (mean(p.rand.mean) - c.i)/ m.rand
-Delta.t <- diff(t.rand.0)
+Delta.t.sd <- sd(diff(t.rand.0))
pt.df <- data.frame(p.new = as.vector(p.new),
t.new = as.vector(t.new),
@@ -72,5 +73,6 @@ ggplot(pt.df) +
geom_point(aes(x = t.new, y = p.new), color='black') +
geom_point(aes(x = t.rand, y = p.rand), color='red') +
geom_abline(aes(intercept = c.rand, slope = m.rand), color='red') +
- geom_vline(xintercept = t.rand.0) +
+ geom_vline(xintercept = t.rand.0, linetype = "longdash") +
+ geom_point(aes(y = p.rand.mean, x = t.rand.mean), color="black", size=4) +
ggtitle("dp vs. dt")
--
GitLab