Normalverteilung Dichtefunktion

Wenn du dir nicht sicher bist, in welchem der anderen Foren du die Frage stellen sollst, dann bist du hier im Forum für allgemeine Fragen sicher richtig.
Antworten
nono191
User
Beiträge: 38
Registriert: Dienstag 6. Oktober 2015, 13:01

Mein Problem ist, dass ich ein Programm schreiben will dass die Dichte berechnet, jedoch gibt es immer falsche Ergebnisse und ich konnte keinen Fehler finden.

Code: Alles auswählen

import math
import scipy.stats as stats


class NormalDistribution:

    def __init__(self, mu, sigma_square):

        self.mu = mu
        self.sigma = math.sqrt(sigma_square)

    def pdf(self, x):
        return 1 / (self.sigma * math.sqrt(2 * math.pi)) * math.exp(-(math.pow(x - self.mu, 2) / (2 * math.pow(self.sigma, 2))))


print(NormalDistribution(0, 3).pdf(1))
print(stats.norm.pdf(1, 0, 3))

Code: Alles auswählen

Output:
0.19496965572274114 # Ausgabe des Programms
0.12579440923099774 # Korrekte Ausgabe
Ich wäre sehr dankbar wenn mir jemand auf die Sprünge helfen könnte.

LG nono
nono191
User
Beiträge: 38
Registriert: Dienstag 6. Oktober 2015, 13:01

P.S Die Formel ist Bild
Sirius3
User
Beiträge: 17703
Registriert: Sonntag 21. Oktober 2012, 17:20

@nono191: wie kommst Du auf die Idee, dass Du sigma^2 angeben mußt?

Code: Alles auswählen

import math

def pdf(x, mu, sigma):
    return math.exp(-(x - mu)**2 / (2 * sigma**2)) / (sigma * (2 * math.pi)**0.5)

print(pdf(1, 0, 3))
nono191
User
Beiträge: 38
Registriert: Dienstag 6. Oktober 2015, 13:01

Hallo Sirius,
Erstmal danke für deine Antwort.

Ich habe gedacht, dass man die Varianz angeben muss.

Danke
nono
Antworten