|
elchico |
Forum-Fortgeschrittener
|
|
Beiträge: 87
|
|
|
|
Anmeldedatum: 04.05.17
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 03.10.2020, 07:59
Titel: Poissrnd um 1 wie randn
|
|
Hallo zusammen,
ich habe eine Normalverteilung um müh+/- dev mit der randn Funktion berechnet:
Nun möchte ich eine analoge Poisson Verteilung berechnen und dachte eigentlich, es sollte relativ einfach sein mit:
Problem ist, dass poissrnd nur integer ausgibt (zumindest bei mir ...). Meine Frage wäre jetzt, was ich anders machen muss, damit doubles rauskommen?
Danke und LG
Michi
|
|
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.495
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 03.10.2020, 10:42
Titel:
|
|
Hallo,
Zitat: |
Problem ist, dass poissrnd nur integer ausgibt |
So ist die Poisson-Verteilung nun mal definiert. Nicht von MATLAB, sondern von Herrn Poisson.
Zitat: |
Meine Frage wäre jetzt, was ich anders machen muss, damit doubles rauskommen? |
Genau genommen kommen Doubles heraus, die ganzzahlig sind.
Wir sollten also einen oder zwei Schritte zurückgehen. Was möchtest du denn letztlich erreichen?
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
elchico |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 87
|
|
|
|
Anmeldedatum: 04.05.17
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 03.10.2020, 10:58
Titel:
|
|
|
|
|
Hallo Harald,
Zitat: |
So ist die Poisson-Verteilung nun mal definiert. Nicht von MATLAB, sondern von Herrn Poisson. |
Ah, das hätte der Herr Poisson ruhig anders definieren können ...
Ich möchte für meine Simulationen Zufallszahlen generieren (zu Anfangs eben Normalverteilt), die ich dann auf meinen Mittelwert +/- Standardabweichung verschiebe, so wie das mit der randn Funktion funktioniert (siehe oben). Das funktioniert eigentlich ganz gut.
Nun wollte ich mein Programm um die Funktion erweitern, auch Poisson-verteilte Zufallswerte zu generieren (zum Beispiel eben um 1 herum), um meine Simulationen um diese Funktion zu erweitern, und diese Kurve dann um meinen Mittelwert +/- SD zu verschieben / strecken (analog meiner randn-Zeile). Das geht aber nur, wenn ich nicht nur ganzzahlige double-Werte bekomme, sondern eben auch Werte zwischen den Integern generiert werden.
Zwei Gedanken noch dazu:
1) Ich weiß, dass die Poisson Verteilung bei größerem lambda immer mehr einer Gauss-Verteilung gleicht, weswegen ich eben die Poisson-Verteilung für kleine Werte (zum Beispiel eben 1) berechnen möchte, und diese dann nur um meinen Mittelwert +/- SD zu verschieben / strecken,
2) Der physikalische Background zu meinen Simulationen zeigt, dass nicht nur ganzzahlige Werte möglich sind, weswegen es mir nicht ausreichen würde, eine Poisson Verteilung für lambda = 100 zu berechnen, weil eben dann nur ganzzahlige Werte berechnet werden würden (mal davon abgesehen, dass ich bei großem lamda, wie unter 1) geschrieben, quasi eine Gauss-Verteilung bekomme).
Ich hoffe, das war verständlich beschrieben, sonst versuche gern ich nochmal mein Glück
LG
Michi
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.495
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 03.10.2020, 11:26
Titel:
|
|
Hallo,
du kannst eine Poisson-Verteilung schon so verschieben, dass du einen bestimmten Mittelwert und eine bestimmte Standardabweichung bekommst, z.B.
Was sich aber nicht ändern wird ist, dass du immer diskrete Werte bekommen wirst. So wollte das Herr Poisson nun mal.
Wenn du kontinuierliche Werte haben willst, dann wäre meine Empfehlung, andere kontinuierliche Verteilungen zu versuchen. Lognormal hat beispielsweise ein annähernd ähnliches Aussehen, was die Dichte angeht - nur eben kontinuierlich.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
elchico |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 87
|
|
|
|
Anmeldedatum: 04.05.17
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 06.10.2020, 08:51
Titel:
|
|
Hallo Harald,
die Lognormal Verteilung probiere ich mal aus. Danke Dir.
LG
Michi
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2025
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|