Gibbs-Boltzmann Zufallszahlen < Mathe-Software < Mathe < Vorhilfe
|
Status: |
(Frage) reagiert/warte auf Reaktion | Datum: | 23:22 Fr 18.12.2009 | Autor: | Nicole1007 |
Aufgabe | Hallo,
ich brauche dringend eure hilfe!! ich soll Gibbs-BOlzmann verteilte Zufalllszahlen in Matlab generieren!
|
Gibt es in Matlab einen Befehl? wie zum bsp. randn für die normalverteilung?
Oder hat jemand einen Tipp oder eine Idee wie ich das machen könnte?
Danke schon mal.
Gruß nicole
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
|
|
> Hallo,
> ich brauche dringend eure hilfe!! ich soll Gibbs-Boltzmann
> verteilte Zufalllszahlen in Matlab generieren!
>
> Gibt es in Matlab einen Befehl? wie zum bsp. randn für die
> normalverteilung?
> Oder hat jemand einen Tipp oder eine Idee wie ich das
> machen könnte?
> Danke schon mal.
> Gruß nicole
Hallo Nicole,
es wäre nützlich zu wissen, auf welche Formel(n)
du dich genau stützt und welche Daten jeweils
vorliegen.
LG Al-Chwarizmi
|
|
|
|
|
Aufgabe | ich habe die Verteilung f(z) = k*exp(-ß*H(z)) gegeben wobei H = 0.5*p*M'*p+v(x), wobei v(x) konstant und k die normierung des Integrals ist. und dass p ist 6-dimenstional. |
Nun soll ich verschiedene zufällige p aus der Verteilung mit hilfe von Matlab ziehen.
Ich dacht mir es müsste doch so wie bei der Normalverteilung eine solche funktion in Matlab geben aber ich habe im netz nichts gefunden und kenne selber keine.
Weiß jemand einen solchen befehl oer hat eine Idee wie man das Problem auf die Normalverteilung zurückführen kann?
Danke
Gruß nicole
|
|
|
|
|
> ich habe die Verteilung f(z) = k*exp(-ß*H(z)) gegeben
> wobei H = 0.5*p*M'*p+v(x), wobei v(x) konstant und k die
> normierung des Integrals ist. und dass p ist
> 6-dimenstional.
> Nun soll ich verschiedene zufällige p aus der Verteilung
> mit hilfe von Matlab ziehen.
> Ich dacht mir es müsste doch so wie bei der
> Normalverteilung eine solche funktion in Matlab geben aber
> ich habe im netz nichts gefunden und kenne selber keine.
> Weiß jemand einen solchen befehl oer hat eine Idee wie
> man das Problem auf die Normalverteilung zurückführen
> kann?
> Danke
> Gruß nicole
... damit kann ich leider praktisch nichts anfangen ...
wofür steht z [mm] (z\in\IR [/mm] oder [mm] z\in\IZ [/mm] oder [mm] z\in\{1,2,....,n\}) [/mm] ?
was ist unter H(z) bzw. H(x) zu verstehen, warum
verwendest du ein v(x), obwohl es konstant sein
soll ...
|
|
|
|
|
Hallo Nicole,
ich habe fast keine Ahnung von all dem, was Du da schreibst. Aber ohne Mühe finde ich verschiedene Algorithmen zur Erzeugung normalverteilter Zufallszahlen, auch im Zusammenhang mit der Gibbs-Boltzmann-Verteilung, zum Befehl rand(x) in Matlab, und verstehe so Dein Problem immer weniger. Es ist doch egal, ob v(x) konstant und p 6-dimensional ist, alles was Du brauchst, ist a) eine Erzeugung von Zufallszahlen mit bekannter Verteilung und b) einen Algorithmus, der diese Verteilung in eine andere überführt. Das kann doch nicht so schwierig sein.
Wo also ist das Problem?
Nebenbei: alle Links sind ein beliebiges Beispiel aus einer entsprechenden google-Suche, keine davon mit weniger als 13.000 Ergebnissen. Ist Dein google kaputt?
lg
reverend
|
|
|
|