Galton Brett

[Blacky]

Hi hat jemand ne Idee wie man ein Galton Brett mit Python macht. Wäre echt nett wenn jemand hier nen Quelltext oder Tipps geben könnte
mfg
Herbert

[numerix]

Hi hat jemand ne Idee wie man ein Galton Brett mit Python macht. Wäre echt nett wenn jemand hier nen Quelltext oder Tipps geben könnte
mfg
Herbert

Der Weg ist ein anderer: DU sagst, welche Ideen du bisher hattest und zeigst möglichst auch schon einen Code-Ansatz dazu. Dann bekommst du sicher Hinweise, ob und was man besser machen könnte ...

[Blacky]

Ok hier ist mein Ansatz: Wusste ich nicht das das so geht sry.
import random
def galton(r, l):
print "galton(", r, ",", l, ")"
if r < l:
return None
else:
if (l == 0) or (r == l):
return 1
else:
return (galton(r-1, l-1) + galton(r-1, l))
galton(0, 2)

[/me]

Ok hier ist mein Ansatz: Wusste ich nicht das das so geht sry.

Code: Alles auswählen

import random
def galton(r, l):
    print "galton(", r, ",", l, ")"
    if r < l:
        return None
    else:
        if (l == 0) or (r == l):
            return 1
        else:
            return (galton(r-1, l-1) + galton(r-1, l))
galton(0, 2)

Du importierst random, nutzt es aber nicht. Möchtest du jetzt die Kugeln zufällig durchlaufen lassen oder möchtest du die mathematische Wahrscheinlichkeit ermitteln? Für zweiteres habe ich mal eben folgenden Code erstellt:

Code: Alles auswählen

def get_probability(rows_total, slot_number):
    def fak(value):
        [... 4 Zeilen Code aus Lerngründen entfernt ...]
    
    def n_over_k(n, k):
        [... 1 Zeile Code aus Lerngründen entfernt ...]
    
    nok = n_over_k(rows_total, slot_number)
    return nok * (1.0 / 2 ** rows_total)

rows = 4
slots = [get_probability(rows, i) for i in range(rows + 1)]
print(slots)

Dein Ansatz klingt eher nach der ersten beschriebenen Vorgehensweise, aber du solltest noch mal deutlich machen was du genau möchtest und was das Ergebnis sein soll.

[derdon]

Code: Alles auswählen

def get_probability(rows_total, slot_number):
    def fak(value):
        [... 4 Zeilen Code aus Lerngründen entfernt ...]
    
    def n_over_k(n, k):
        [... 1 Zeile Code aus Lerngründen entfernt ...]
    
    nok = n_over_k(rows_total, slot_number)
    return nok * (1.0 / 2 ** rows_total)

rows = 4
slots = [get_probability(rows, i) for i in range(rows + 1)]
print(slots)

Warum 4 Zeilen für fak?

Code: Alles auswählen

from operator import mul
def fac(num):
    return reduce(mul, xrange(1, num+1))

reicht doch auch.

[numerix]

Warum 4 Zeilen für fak?

Code: Alles auswählen

from operator import mul
def fac(num):
    return reduce(mul, xrange(1, num+1))

reicht doch auch.

Viel zu lang ...

Code: Alles auswählen

f=lambda x:x*f(x-1)if x else 1

reicht doch auch.

[Masaru]

Warum 4 Zeilen für fak?

Code: Alles auswählen

from operator import mul
def fac(num):
    return reduce(mul, xrange(1, num+1))

reicht doch auch.

Viel zu lang ...

Code: Alles auswählen

f=lambda x:x*f(x-1)if x else 1

reicht doch auch.

Viel zu unverständlich und sieht noch dazu echt ugly aus.

Code: Alles auswählen

import this

>>Masaru<<

[numerix]

Viel zu unverständlich und sieht noch dazu echt ugly aus.

Sehe ich nicht so. Es ist eine saubere, simple Umsetzung der rekursiven Definition der Fakultät. Und der "reduce - mul - xrange"-Variante an Verständlichkeit ganz gewiss überlegen.

@Blacky: Was genau soll deine Simulation denn leisten? Soll der ganze Prozess für eine bestimmte Anzahl an Kugeln quasi parallel simuliert werden? Soll dann für jede Zeile angezeigt werden, wie viele Kugeln sich in den einzelnen "Bahnen" gerade befinden oder genügt der Endzustand. Oder soll jeweils immer nur der Lauf einer Kugel simuliert werden und gezeigt werden, in welches Kästchen sie fällt?

Ohne eine klare Beschreibung, was du unter "Galton Brett mit Python machen" verstehst, wird man dir nicht helfen können.

[/me]

Warum 4 Zeilen für fak?

Code: Alles auswählen

from operator import mul
def fac(num):
    return reduce(mul, xrange(1, num+1))

reicht doch auch.

Nun, die 4 Zeilen liegen wohl in erster Linie daran, dass 25 Jahre mit anderen Programmiersprachen so schwer aus dem Kopf gehen.

reduce() war mir beim Schreiben geistig einfach nicht präsent.

[Masaru]

.... Und der "reduce - mul - xrange"-Variante an Verständlichkeit ganz gewiss überlegen. ....

Das sehe ich wiederrum nicht so ... Deine lambda-Variante mag zwar eine "Alternative" sein, aber ist sie mit Sicherheit nicht "überlegen": völliger Quatsch .

[/me]

Code: Alles auswählen

f=lambda x:x*f(x-1)if x else 1

Nettes Beispiel, aber ein Entwickler der die Fakultät ernsthaft via Rekursion berechnet gehört geteert und gefedert.

[Dav1d]

Eine tolle Mischung

Code: Alles auswählen

f = lambda x: reduce(__import__('operator').mul, xrange(1, x+1))

[BlackJack]

@numerix: Aber das ist so rekursiv. Würde ich in Python nicht verwenden.

Aber warum überhaupt die Fakultät implementieren:

Code: Alles auswählen

import gmpy
n_over_k = gmpy.bincoef

[/me]

Aber warum überhaupt die Fakultät implementieren:

Code: Alles auswählen

import gmpy
n_over_k = gmpy.bincoef

Ha, wieder was gelernt. Merci!

[numerix]

@numerix: Aber das ist so rekursiv. Würde ich in Python nicht verwenden.

Das sollte auch kein Beispiel für eine performante Implementation sein. Es ging in dem Zusammenhang ja um die Kürze des Codes und ich bin nach wie vor der Ansicht, dass das mal ein Beispiel ist, dass ein Code kurz UND verständlich sein kann (auch wenn Masaru das für "Quatsch" hält). Aus meiner Sicht ist das ein Stück eleganter Code (klingt natürlich doof, wenn man es selbst präsentiert, aber ich bin nicht mal sicher, ob ich das nicht mal woanders gesehen habe und davon inspiriert wurde).

Natürlich sieht eine performante Implementierung anders aus. Vor einiger Zeit habe ich mal eine äußerst interessante Seite dazu gefunden, wo jemand die Fakultät zu seinem Hobby gemacht hat. Das ist eine echte Fundgrube und es ist wirklich erstaunlich, wie effizient man die Fakultät (für große n) berechnen kann - die Algorithmen passen dann aber auch nicht mehr in eine Zeile ...

Aber warum überhaupt die Fakultät implementieren:

Weil es Spaß macht. Wem es keinen macht:

Code: Alles auswählen

from math import factorial

Edit: Natürlich braucht man für das Galtonbrett überhaupt keine Fakultät und auch keine Binomialkoeffizienten - da kommt man auch mit reiner Addition aus.

[derdon]

Ich finde die Variante von numerix elegant und leserlich. Dass sie kurz ist, wird wohl keiner bestreiten können. Meine dagegen ist unleserlicher; der einzige Vorteil liegt in Performance und Speicherverbrauch (sehe ich als einen Punkt).

[Masaru]

... auch wenn Masaru das für "Quatsch" hält...

Whop, ich habe lediglich gesagt, dass Deine Aussage die "lambda-Variante wäre überlegen" völliger Quatsch sei. Wer lesen kann ist klar im Vorteil .

[numerix]

Du bist aber giftig heute.

[derdon]

In Lesbarkeit ist sie das aber, in Performance nicht. numerix hat nicht gesagt, was er mit "überlegen" meinte.

[Masaru]

*popcornkau* doch, und zwar meine in Punkto "Verständlichkeit" ... was lambda jedoch nie sein wird. Falls Du mir zeigen kannst, wo etwas Gegenteiliges steht, dann ziehe ich meinen Hut und nehme alles mit tiefster Demütigung und Entschuldigung zurück.

@numerix: wann war ich denn mal nicht giftig? *grübel*