Einzelnen Beitrag anzeigen
  #94  
Alt 27.04.12, 21:39
Benutzerbild von richy
richy richy ist offline
Singularität
 
Registriert seit: 01.05.2007
Ort: karlsruhe
Beitr?ge: 4.170
Standard AW: Math Verhulst 1989

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.

Ge?ndert von richy (27.04.12 um 21:42 Uhr)
Mit Zitat antworten