Groups | Search | Server Info | Login | Register

Groups > de.sci.electronics > #353341

Re: random Statistik

Show all headers | View raw

Helmut Schellong wrote:
> 
> Helmut Schellong wrote:
> 
> Aktuelles Skript:
> ====================================================================
> ARITH_FMT='%.6Lf'
> set r:025 r0:025 z:025 n:020 e:020 s:025 m:025 m0:025 sa:025 dp:025 dm:025 d0:025
> s=0 m=0.0 m0=-1.0 sa=0.0 d=00 i=0000000000
> array 1 6  Z  $i $i $i $i $i $i
> 
> echo Anzahl Zufallszahl SummeZ Mittelwert Abweichungen
> 
> for n from 1 to 1000000 repeat
> do
>     let "r= random(6)+1" "d=r-r0, r0=r" "d>=0"
>     run? "let 'dp+=d'" "let 'dm-=d'"
>     let "s+=r" "m=s/$n.0" "sa+=m-3.5" "++Z[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
> ====================================================================
> 
> Das Skript enthält nun neben neuen Berechnungen Zahlenbeschriftungen.

Chatgpt:
-------
The script appears to be written in a Unix shell scripting language, 
specifically using syntax and constructs typical of bash or sh. 
--------

Laeuft allerdings bei mir nicht:

------
bash: line 4: array: command not found
Anzahl Zufallszahl SummeZ Mittelwert Abweichungen
bash: line 8: syntax error near unexpected token `from'
bash: line 8: `for n from 1 to 1000000 repeat'
------
Ich vermute mal das ist fuer deine selbstgeschriebene Spezial-shell.

Es ist auch keine gute Idee sowas als shell script zu programmieren...
Gluecklicherweise gibts chatgpt um das ganze in C zu uebersetzen:

---------
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

int main() {
    long double ARITH_FMT = 0.0;
    long long r = 0, r0 = 0, z = 0, n = 0, e = 0, s = 0, m = 0, m0 = 0;
    long double sa = 0.0, dp = 0.0, dm = 0.0;
    int d0 = 0, d = 0, i = 0;
    long long Z[7] = {0};  // Array to store counts of dice rolls from 1 to 6
    long long count = 1000000;

    printf("Anzahl Zufallszahl SummeZ Mittelwert Abweichungen\n");

    // Initialize random number generator
    srand(time(NULL));

    for (n = 1; n <= count; n++) {
        r = (rand() % 6) + 1;  // Generate random number between 1 and 6
        d = r - r0;
        r0 = r;

        if (d >= 0) {
            dp += d;
        } else {
            dm -= d;
        }

        s += r;
        m = s / (long double)n;
        sa += m - 3.5;
        Z[r]++;

        if (d == 0) {
            d0++;
        }

        printf("%lld %lld %lld %.6Lf %lld\n", n, r, s, m, llround(sa));
    }

    r = (count - 1) / 6;

    printf("%lld %lld %lld %lld %lld %lld  Anzahlen von 1 2 3 4 5 6\n",
           Z[1], Z[2], Z[3], Z[4], Z[5], Z[6]);

    printf("%lld %lld %lld %lld %lld %lld  %% Differenzen zum Erwartungswert\n",
           Z[1] - r, Z[2] - r, Z[3] - r, Z[4] - r, Z[5] - r, Z[6] - r);

    printf("%.6Lf %.6Lf %d  Summen der pos. und neg. Differenzen und Gleichheiten Z-Z0\n", dp, dm,
d0);

    return 0;
}
-----------

Wenn man aufs ausdrucken der count Zeilen verzichtet laeuft es sogar richtig schnell:


Anzahl Zufallszahl SummeZ Mittelwert Abweichungen
166214 166518 166529 166583 167060 167096  Anzahlen von 1 2 3 4 5 6
-452 -148 -137 -83 394 430  % Differenzen zum Erwartungswert
973651.000000 973646.000000 166457  Summen der pos. und neg. Differenzen und Gleichheiten Z-Z0


Mit count=1000000000 (100*100 Mio)  dauert es 4 Minuten:

bash-5.2# date; ./script2 ; date
Fri Jun 14 12:19:19 CEST 2024
Anzahl Zufallszahl SummeZ Mittelwert Abweichungen
1666727685 1666703808 1666627423 1666644804 1666640939 1666655341  Anzahlen von 1 2 3 4 5 6
61019 37142 -39243 -21862 -25727 -11325  % Differenzen zum Erwartungswert
9722203361.000000 9722203357.000000 1666717345  Summen der pos. und neg. Differenzen und
Gleichheiten Z-Z0
Fri Jun 14 12:23:34 CEST 2024


Man sieht schoen dass die Differenzen wie erwartet etwa 100 mal so gross geworden sind.
Wenn man eine Statistik ueber die Differenzen machen will sollte man das ganze aber oefter
rechnen, also statt 10Mrd lieber 10000 mal eine Mio.


Nochmal das 100fache waeren mit dann 400Minuten auch kein Problem und an einem Tag machbar.

Back to de.sci.electronics | Previous | Next — Previous 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