Metropolisalgorithmus mit Matlab (mhsample)

Neue Frage »

GozyllaR Auf diesen Beitrag antworten »
Metropolisalgorithmus mit Matlab (mhsample)
Hallo,

ich habe ein Problem welches an die Stochastikexperten unter euch geht.

Ich habe eine Funktion. Diese beruht auf statistischen Daten und soll eine Wahrscheinlichkeitsdichtefunktion darstellen. Ob das mathematisch gesehen wirklich eine ist, sei mal dahingestellt.

Wie ich zu der Funktion komme:

code:
1:
2:
3:
4:
5:
6:
7:
8:
9:
x = [750 1250 1750 2250 2750 3250 3750 4250 4750 5250 5750 6250 6750 7250 7750 8250 8750 9250 9750]; 
y = [0.1 0.12 0.12 0.115 0.105 0.09 0.07 0.055 0.04 0.03 0.025 0.023 0.021 0.02 0.018 0.015 0.012 0.011 0.01]; 

p = polyfit(x,y,5); 
xFit = 500:1:10000; 
yFit = polyval(p,xFit); 

plot(xFit,yFit)


Nun möchte ich mit Matlab zufällig Zahlen in ein Array generieren, welche dieser Funktion folgen. Das bedeutet beispielsweise das Histogram des generierten Arrays soll so aussehen wie die Funktion, die ich oben erstellt habe.

Ich vermute ich komme mit der Funktion mhsample weiter. Leider bekomme ich es nicht auf die Reihe diese Funktion richtig zu benutzen. Da fehlt mir das mathematische Verständnis für den Metropolisalgorithmus.

Eventuell muss ich auch vorher zunächst ein Fitting einer Verteilung durchführen. Aber ob das notwendig ist oder wie das genau funktioniert weiß ich auch nicht.

Vielen Dank vorab für die Hilfe.
Grüße
Neue Frage »
Antworten »



Verwandte Themen

Die Beliebtesten »
Die Größten »
Die Neuesten »