Ok, ich lasse es mal ganz langsam angehen (Bin mit Laptop unterwegs). Das Thema Informationsverlust habe ich bereits auf meiner HP dokumentiert.
http://home.arcor.de/richardon/richy...nfoverlust.htm
Simuliere ich erstmal r=4 also chaotische Vorwaertsiteration und Rueckwaertsiteration mit richtigem Vorzeichen :
Das war einfach dieser Quellcode :
Zitat:
restart;
N:=100; # Laenge Infostring
N_s:=100; # Anzahl Anfangswerte 0..1
r:=1+sqrt(9);
ds:=1/N_s; s:=0;
for j from 1 to N_s do
#fuer alle Anfangswerte
++++++++++++++++++
k:=0; s:=ds*j;
sp[j]:=s; # Anfangswert von j
for i from 1 to N do
#Fuer alle N Iterationen
+++++++++++++++++++
s_alt:=s;
# VORWAERTSITERATION
*******************
s:=evalf(r*s*(1-s));
# SOFORTIGE BERCHNUNG RUECKITERATIONEN S0 und S1
*******************************************
s0:=evalf(1/2/r*(r-(r^2-4*r*s)^(1/2)));
s1:=evalf(1/2/r*(r+(r^2-4*r*s)^(1/2)));
Speichern fuer jeden Anfangswert und jede Iterationsnummer
*********************************************
(s0-s_alt)^2 < (s1-s_alt)^2 then inf[j][k]:=0;
else inf[j][k]:=1; fi; k:=k+1;
od:
od:
with(plots):
p := seq(seq(plots[polygonplot]([[i, j], [i+1, j], [i+1, j+1], [i, j+1]], color = COLOR(RGB, inf[i][j],0,0)), i = 1 .. N_s-1), j = 1 .. N-1):
plots[display]([p]);
|
Ok und jetzt will ich wissen wie sich die Rueckwaerstiterierte verhaelt. Konkret welcher Loesungszweig sich ergibt.