Groups | Search | Server Info | Login | Register


Groups > de.sci.electronics > #353477

Re: random Statistik

From Helmut Schellong <var@schellong.biz>
Newsgroups de.sci.electronics
Subject Re: random Statistik
Date 2024-06-17 15:05 +0200
Message-ID <v4pcan$1pb73$1@solani.org> (permalink)
References (12 earlier) <666D7563.3B0069E4@proton.me> <v4k3ru$1me5m$1@solani.org> <666EADAE.44F26696@proton.me> <v4njn7$1oabj$1@solani.org> <66700BA8.4C2E8109@proton.me>

Show all headers | View raw


Carla Schneider wrote:
> Helmut Schellong wrote:

>>>> Ich schrieb bereits in [08.06.2024, 15:32] vom Pendeln um den Mittelwert.
>>>
>>> Das entspricht einem eindimensionalen
>>> https://en.wikipedia.org/wiki/Random_walk
>>
>> Wird wohl so sein.
>> Jedenfalls nehmen die Frequenz und die Amplitude des Pendelns fortlaufend ab.
> 
> Eigentlich nimmt die absolute Amplitude zu mit Wurzel(Anzahl).
> Was abnimmt ist die relative Amplitude als Wurzel(Anzahl)/Anzahl.
> 
>> Sie werden erst im Unendlichen 0.
> 
> Die absolute Amplitude wird auch im Unendlichen Unendlich...

1 4 4 0.000000 -4
2 1 5 2.000000 -5
3 5 10 1.666667 -7
4 5 15 2.500000 -8
5 5 20 3.000000 -8
6 4 24 3.333333 -9
7 2 26 3.428571 -9
8 4 30 3.250000 -9
9 3 33 3.333333 -9
10 1 34 3.300000 -9
11 5 39 3.090909 -10
...
99993 5 350380 3.503995 1020
99994 1 350381 3.504010 1020
99995 1 350382 3.503985 1020
99996 2 350384 3.503960 1020
99997 2 350386 3.503945 1020
99998 1 350387 3.503930 1020
99999 4 350391 3.503905 1020
100000 1 350392 3.503910 1020

Ich rede in der Regel stets von den Zahlenfolgen, wie vorstehend.
Es ist unübersehbar, daß der Mittelwert um den Erwartungswert 3.500000 pendelt.
Es resultieren positive und negative Differenzen.
Deshalb pendelt die aufaddierte Summe der Differenzen zwischen positiv und negativ.
So starke Abweichungen wie zu Beginn gibt es bei späteren Zeilen nicht mehr.
Und zwar nie, weil die Differenzen immer kleiner werden - sie gehen gegen Null.
Auch die Pendelfrequenz geht deshalb gegen Null.
Es ist erkennbar, daß sich die dritte und vierte Nachkommastelle MW unten fast nicht mehr ändern.

>> Genau deshalb wird die Anzahl der Zufallszahlen irgendwann extrapolierend irrelevant.

>>> Wenn ich die Zufallszahlen in ein File schreiben wollte damit andere
>>> sie testen koennen wuerde ich das in binaerer Form tun, also 3 bit pro Zahl (von 1-6),
>>> das koennte man in 10 Zahlen pro 4 Byte tun. Dann wuerden 10Mio Zahlen 4Mbyte brauchen.
>>> Man koennt auch die Bytes des SPRITZ Algorithmus so ausgeben, das waere dann 1 Byte pro Zahl.
>>>
>>> In C koennte man das einfach tun, mit fopen,fwrite,fclose, aber ginge das auch in deinem
>>> shell-script ?
>>> Ich meine jetzt nicht dass du das machen sollst, sondern es geht mir nur um die
>>> Praktikabilitaet soetwas als Shellscript zu programmieren.
>>
>> Ja, das wird nur mit meiner Shell bish funktionieren.
>>      http://www.schellong.de/htm/bishmnk.htm#bishK
>> Die bish ist komplett binär.
>> Variablen können alle Byte-Werte von 0..255 enthalten, >300 Mio. davon, falls erforderlich.
>> Dateien werden geöffnet per  > < >>  und geschlossen per ><<< .
>>
>>         catv  [quelle(n) ...]  [=ziel]
>>               [[offset,]max,]objekt ...     objekt: name|handle
>>               [:+-][offset][:+-]
>>               [:.]objekt[:.]
>>               /konstante
>>               /abc%x%255xyz
>>               .konstante
>>               , ,. ,r
>> Vorstehendes Kommando kann Dateien und Variablen lesen und schreiben
>> und Konstanten lesen und konvertieren.
>> Es können beliebig viele Variablen und Dateien, mit jeweils Offset und Anzahl Byte, und Konstanten
>> gemischt verkettet in ein beliebiges Ziel ebenfalls mit Offset und Anzahl geschrieben werden.
> 
> Er muesste also die Zahlen die bei dir ja Dezimalszahlen als Strings sind,
> erst mal in Binaerzahlen als Bytes umwandeln und die dann in eine Datei schreiben...
> 
> 
>>
>> Mittels  set sa:025  wird die Variable  sa  erzeugt und mit 25 Ziffern '0' initialisiert.
>> Mittels  set sa:.25  wird die Variable  sa  erzeugt und mit 25 Zeichen '.' initialisiert.
>> Mittels  set sa:46.25  wird die Variable  sa  erzeugt und mit 25 Zeichen '.' initialisiert.
>> Mittels  set sa:0.25000000  wird die Variable  sa  erzeugt und mit 25 Mio. 0-Byte initialisiert.
> 
> sa ist also eigentlich ein string, das ist fuer eine shell-Varible auch vernuenftig weil
> es der normale Einsatzzweck ist, z.B. der Library search path...
> 
> Zum Rechnen ist es eher unguenstig, weil bei jeder Rechnung die Strings erst mal in
> Zahlen konvertiert werden muessen, und danach wieder zurueck in Strings.

       typeset  [ -+i[4|8]ISfF{s[obase]|obase} ]  [Name[=wert]]...
                      Dieses Kommando ist in Verbindung mit arithmetischen
                      Ausdrücken besonders interessant.
                      Option -f markiert Gleitkomma-Variablen als binär-float.
                      Option -i markiert Variablen als binär-integer.
                      Sie werden solchermaßen besetzt binär gespeichert
                      und nicht als dezimale ascii-Zeichenkette. ... ...

Die Umwandlungsfunktionen wie z.B. atoi() und itoa() sind ultraschnell implementiert.
Sind in meinen C-Büchern.

>>> Er interpretiert set sa:025 dass sa eine oktalzahl zugewiesen wird
>>> und damit sa als integer. Ein paar Zeilen spaeter wird sa dann 0.0 zugewiesen
>>> womit sa eine floatingpoint zahl sein  soll...
>>
>> Variablen in einer Shell haben keinen Typ.
> 
> Haette ja sein koennen dass das bei deiner Shell anders ist.

Prinzipiell ja, ich hatte nach Einführung von Gleitkomma darüber gegrübelt.
Du hast aber sicher keine entsprechenden Deklarationen bei den Variablen gesehen.
Eben weil ich das nicht implementierte.
Der Gleitkommatyp ist nur am Inhalt erkennbar: Ist ein '.' oder 'e' vorhanden?
Der Arithmetik-Modul muß solches ja erkennen!
Es ist auch Ausgabe per  10^3  10^6  10^9  10^12  konfigurierbar - eine Engineer-Ausgabe.

>> Standardmäßig enthalten sie Zeichenketten.
>> sa:025 erzeugt Variable sa und initialisiert mit 25 '0'-Ziffern.
>> Danach wird 0.0 zugewiesen, um ihr den Gleitkommatyp zuzuordnen.
>> Nach 0.0 sind nur die 3 Zeichen 0.0 enthalten, jedoch 22 weitere Byte reserviert.
> 
> D.h. es ist eigentlich ein String der Laenge 22 in dem "0.0" steht.

Länge 25. 3 + 22 undefinierte Byte.

> In C dagegen macht chatgtp daraus eine "long double" .
> Auf einem Intel oder Amd PC ist das eine 10byte floating-point Zahl,
> mit 63 bit Fraction und 15 Bit exponent,das groesste was es auf diesen Prozessoren gibt.
> 
Ich hatte mal mit long double == 128 Bit auf HP-HW gearbeitet.
33 oder 34 signifikante Dezimal-Stellen!

Das aktuelle Skript wird langsam professioneller.
Bei den Neuigkeiten sind weitere Besonderheiten/Alleinstellungsmerkmale erkennbar:
.   expr "$cnt" :: '[^0-9]'
.   do  c$r$e=$i; c$r$e=0; done
.   do  s=c$r$e; echo $r$e:${{s}} %c; done


Aktuelles Skript:
-----------------------------------------------------------------------------------
set cnt:020
[ $# -ge 1 ] && cnt="$1"
expr "$cnt" :: '[^0-9]' && cnt=10000
[ "$cnt" -lt 1 ] && cnt=10000

ARITH_FMT='%.6Lf'
set r:025 r0:025 n:020 e:020 s:025 m:025 m0:025
set sa:025 dp:025 dm:025 d0:025
s=0 m=0.0 m0=-1.0 sa=0.0 d=00 skip=1 i=0000000000
array 1 6  Z  $i $i $i $i $i $i

for r in 1 2 3 4 5 6
do for e in 1 2 3 4 5 6
    do  c$r$e=$i; c$r$e=0; done
done

echo Anzahl Zufallszahl SummeZz Mittelwert AbweichungenMw

for n from 1 to $cnt repeat
do
    let "r= random(6)+1" "d=r-r0" "d>=0"
    run? "let 'dp+=d'" "let 'dm-=d'"
    if [ skip -le 0 ]
    then let "++c$r0$r, skip=1"
    else skip=0
    fi
    let "s+=r" "m=$s.0/$n.0" "sa+=m-3.5" "++Z[r], r0=r" "d==0" && let "++d0"
    echo "$n $r $s $m $((llround(sa)))"
done

let "r=($n-1)/6"
echo "%n$Z1 $Z2 $Z3 $Z4 $Z5 $Z6  " Anzahlen von 1 2 3 4 5 6
echo "$((Z1-r)) $((Z2-r)) $((Z3-r)) $((Z4-r)) $((Z5-r)) $((Z6-r))  " %
        Differenzen zum Erwartungswert
echo "$dp $dm $d0  " Summen der pos. und neg. Differenzen und Gleichheiten Z-Z0
echo "Anzahlen der 36 Zahlenfolgen 11 .. 25 .. 66:"
for r in 1 2 3 4 5 6
do for e in 1 2 3 4 5 6
    do  s=c$r$e; echo $r$e:${{s}} %c; done
    echo
done
-----------------------------------------------------------------------------------
99989 1 350368 3.504055 1020
99990 4 350372 3.504030 1020
99991 1 350373 3.504035 1020
99992 2 350375 3.504010 1020
99993 5 350380 3.503995 1020
99994 1 350381 3.504010 1020
99995 1 350382 3.503985 1020
99996 2 350384 3.503960 1020
99997 2 350386 3.503945 1020
99998 1 350387 3.503930 1020
99999 4 350391 3.503905 1020
100000 1 350392 3.503910 1020

16623 16560 16740 16655 16723 16699   Anzahlen von 1 2 3 4 5 6
-43 -106 74 -11 57 33   Differenzen zum Erwartungswert
97563 97562 16437   Summen der pos. und neg. Differenzen und Gleichheiten Z-Z0
Anzahlen der 36 Zahlenfolgen 11 .. 25 .. 66:
11:1343 12:1365 13:1334 14:1339 15:1397 16:1412
21:1467 22:1285 23:1396 24:1406 25:1382 26:1423
31:1407 32:1375 33:1423 34:1389 35:1428 36:1372
41:1404 42:1410 43:1422 44:1378 45:1402 46:1394
51:1391 52:1412 53:1390 54:1391 55:1381 56:1388
61:1421 62:1354 63:1381 64:1342 65:1380 66:1416
-----------------------------------------------------------------------------------


-- 
Mit freundlichen Grüßen
Helmut Schellong   var@schellong.biz
http://www.schellong.de/c.htm  http://www.schellong.de/c2x.htm  http://www.schellong.de/c_padding_bits.htm
http://www.schellong.de/htm/bishmnk.htm  http://www.schellong.de/htm/rpar.bish.html  http://www.schellong.de/htm/sieger.bish.html
http://www.schellong.de/htm/audio_proj.htm  http://www.schellong.de/htm/audio_unsinn.htm  http://www.schellong.de/htm/tuner.htm
http://www.schellong.de/htm/string.htm  http://www.schellong.de/htm/string.c.html  http://www.schellong.de/htm/deutsche_bahn.htm
http://www.schellong.de/htm/schaltungen.htm  http://www.schellong.de/htm/math87.htm  http://www.schellong.de/htm/dragon.c.html

Back to de.sci.electronics | Previous | NextPrevious in thread | Next in thread | Find similar


Thread

Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-10 11:24 +0200
  Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-10 17:05 +0200
    Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-11 11:48 +0200
      Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-11 13:00 +0200
        Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-11 15:50 +0200
          Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-11 16:26 +0200
            Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-11 17:51 +0200
      Re: random Statistik Rolf Bombach <rolfnospambombach@invalid.invalid> - 2024-06-11 14:31 +0200
      Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-11 15:45 +0200
      Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-12 05:04 +0200
        Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-12 08:43 +0200
          Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-12 16:36 +0200
            Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-12 21:41 +0200
              Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-13 21:01 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-14 10:13 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-14 12:36 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-15 00:44 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-15 13:05 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-15 15:10 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-16 11:17 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-16 22:59 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-17 12:10 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-17 15:05 +0200
                Re: random Statistik Rolf Bombach <rolfnospambombach@invalid.invalid> - 2024-06-15 20:34 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-15 23:12 +0200
                Re: random Statistik Rolf Bombach <rolfnospambombach@invalid.invalid> - 2024-06-17 15:16 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-17 16:40 +0200
                Re: random Statistik - Windows PowerShell Helmut Schellong <var@schellong.biz> - 2024-06-18 00:17 +0200
            Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-13 09:43 +0200
              Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-13 12:55 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-13 14:05 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-13 16:11 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-13 16:57 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-13 20:47 +0200
                Re: random Statistik Carla Schneider <carla_schn@proton.me> - 2024-06-18 10:59 +0200
                Re: random Statistik Helmut Schellong <var@schellong.biz> - 2024-06-18 17:37 +0200

csiph-web