WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Poissrnd um 1 wie randn

 

elchico
Forum-Fortgeschrittener

Forum-Fortgeschrittener


Beiträge: 87
Anmeldedatum: 04.05.17
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.10.2020, 08:59     Titel: Poissrnd um 1 wie randn
  Antworten mit Zitat      
Hallo zusammen,

ich habe eine Normalverteilung um müh+/- dev mit der randn Funktion berechnet:

Code:
arrayValues = devValue*randn(numValues,1) + meanValue;


Nun möchte ich eine analoge Poisson Verteilung berechnen und dachte eigentlich, es sollte relativ einfach sein mit:
Code:
arrayValues = devValue*poissrnd(1,numValues,1) + meanValue;


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
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 03.10.2020, 11:42     Titel:
  Antworten mit Zitat      
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 ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
elchico
Themenstarter

Forum-Fortgeschrittener

Forum-Fortgeschrittener


Beiträge: 87
Anmeldedatum: 04.05.17
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.10.2020, 11:58     Titel:
  Antworten mit Zitat      
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 Wink
LG
Michi
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 03.10.2020, 12:26     Titel:
  Antworten mit Zitat      
Hallo,

du kannst eine Poisson-Verteilung schon so verschieben, dass du einen bestimmten Mittelwert und eine bestimmte Standardabweichung bekommst, z.B.
Code:
numValues = 10000;
devValue = 5;
meanValue = 3;

lambda = 8;
[m, v] = poisstat(lambda)

arrayValues = devValue/sqrt(v)*(poissrnd(lambda,numValues,1) - m) + meanValue;
mean(arrayValues)
std(arrayValues)


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 ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
elchico
Themenstarter

Forum-Fortgeschrittener

Forum-Fortgeschrittener


Beiträge: 87
Anmeldedatum: 04.05.17
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.10.2020, 09:51     Titel:
  Antworten mit Zitat      
Hallo Harald,

die Lognormal Verteilung probiere ich mal aus. Danke Dir.

LG
Michi
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 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.