Eigentlich war das als Antwort für den Thread Spritz-Ess-Abstand bei NovoRapid gedacht.
Aber es ist ein bischen off-topic und der Post so lang mit den Bildern, dachte ich lieber ein eigener Thread.
---
Ich stelle mir das mit dem SEA so vor:
Essen lässt den BZ ansteigen.
Insulin senkt den BZ.
Wenn man die beiden Wirkkurven von IE und BE addiert, erhält man eine neue Funktion. Nämlich die Änderung des Blutzuckers.
Ich hab da mal was vorbereitet.
(Irgendwelche Extra-Effekte wie Resistenz bei hohem BZ usw werden ignoriert und BE-Faktor ist einfach 1)
Obwohl das alles natürlich total vereinfacht ist, ist es doch halbwegs anschaulich und interessant was dabei rauskommt.
Im Idealfall sieht das zum Beispiel so aus, gegessen wird eine BE und gleichzeitig wird eine IE gespritzt:
[Blockierte Grafik: https://i.imgur.com/07612Jy.png]
'essen' geht nach oben und würde den BZ um 20 angeben.
'insulin' geht noch unten und würde den BZ um 20 senken
Blutzuckeränderung = essen + insulin
Hier ergibt das in der Summe sehr schön Null, also eine gerade Linie, weil die Kurven von essen und insulin sich genau aufheben.
Kurzer Einschub wie es aussieht wenn man zu wenig spritzt, zum Beispiel 2 BE essen aber nur 1 IE spritzen:
[Blockierte Grafik: https://i.imgur.com/NYKtFc1.png]
essen würden den BZ um 40 angeben, insulin aber nur um 20 senken.
In der Summe "fehlt am Ende etwas", die rote BZ-Linie geht nicht wieder auf 0 zurück.
(Man müsste also nachspritzen oder vielleicht wollte man absichtlich einen hohen BZ)
Problem ist dass es arger Zufall wäre, wenn die Wirkkurven von Essen und Insulin immer zusammen passen.
Sagen wir die BE wirken doppelt so schnell wie das Insulin. 1 BE essen, 1 IE spritzen, aber die BE wirkt doppelt so schnell wie das Insulin:
[Blockierte Grafik: https://i.imgur.com/2KJQ6qu.png]
Die orangene BE-Linie steigt nun steiler als die grüne Insulin-Kurve runtergeht.
Man sieht das nach ~4 Stunden der BZ wieder beim Anfangswert ankommt, aber zwischen hoch ging.
Maximalwert etwa bei 12.
Nun mit 1 IE spritzen, eine Zeiteinheit warten, dann 1 BE essen;
[Blockierte Grafik: https://i.imgur.com/67nHIjP.png]
Jetzt ist der BZ nur noch um umgefähr 3 gestiegen. (vorher 12!)
Allerdings ging der BZ vorher auch erstmal runter, wenn man also schon einen tiefen BZ-Wert vor dem Essen hatte so würde man erstmal unterzuckern.
Die Unterzucker wäre allerdings nur durch den zu langen SEA verursacht und würde dann wohl von selbst wieder weggehen sobald das essen wirkt.
Aber gut ist das ja trotzdem nicht.
Versuchen wir mal mit kürzerem SEA:
Nun mit 1 IE spritzen, eine halbe Zeiteinheit warten, dann 1 BE essen:
[Blockierte Grafik: https://i.imgur.com/DJFWGQi.png]
der BZ geht nun um etwa 8 nach oben.
Fazit:
SEA=0 -> BZ ging um 12 hoch
SEA=1 -> BZ ging nur um 3 hoch, dafür leider Unterzucker während der Mahlzeit
SEA=0,5 -> BZ geht um 8 hoch. Nicht so gut wie bei SEA=1, aber dafür ohne Unterzuckerung und besser als gar kein SEA.
---
Das Skript:
###################
#############
#############
SEA = 0
BE = 5
IE = 3
BEspeed = 0.5
IEspeed = 1
##########
#############
reset
# Scale font and line width (dpi) by changing the size! It will always display stretched.
set terminal svg size 400,300 enhanced fname 'arial' fsize 10 butt solid
set output 'out.svg'
set grid
set xtics 1
set ytics 10
set yrange [:]
set xrange [-3:]
set offset graph 0.1, graph 0.1, graph 0.1, graph 0.1
set ylabel "Änderung BZ"
set xlabel "Zeit"
xoffset = -2
#set xtics add ("essen & spritzen 0" "0") rotate by 90 right
set xtics add ("spritzen" SEA) rotate by 90 right
set xtics add ("essen" "0") rotate by 90 right
essen(x) = ((( tanh(x*BEspeed+xoffset)) *BE)+BE)*10
insulin(x) = ((( -tanh(x*IEspeed+xoffset-SEA)) *IE)-IE)*10
BZ(x) = essen(x) + insulin(x)
set title "BE:".BE ." IE:".IE . " SEA:".sprintf("%1.2f",SEA)."\n"." BE-Speed:".sprintf("%1.2f",BEspeed)."\nInsulin-Speed:".sprintf("%1.2f",IEspeed)
set style line 1 linewidth 3 lc 7 lt 3
set style line 2 linewidth 3 lc 3 lt 3
set style line 3 linewidth 4 lc 1 lt 3
set key outside bottom horizontal box
plot essen(x) ls 1, insulin(x) ls 2, BZ(x) ls 3 title "BZ(x)=essen(x)+insulin(x)"
###############################
einfach kopieren und auf http://gnuplot.respawned.com/ einfügen. (oder normales gnuplot installieren)
In den ersten Zeilen in "####-Kasten" kann man die Werte für BE, IE usw ändern und gucken was für eine Grafik dann dabei rauskommt.