wie Zufallszahlen zu generieren, die Skew Normalverteilung in Matlab folgen

Question

Ich habe Wahrscheinlichkeitsdichtefunktion von Skew Normalverteilung.Ich möchte Zufallszahl, die der Skew Normalverteilung in Matlab folgen.

Answer 1

kann nicht für ihre Leistung/Angemessenheit bürgen, aber http://azzalini.stat.unipd.it/SN/ sagt die folgenden, und hat einen Link zu einer ZIP-Datei von MATLAB-Funktionen:

Die Bibliothek in Matlab von Nicola Sartori portiert wurde. Bisher bezieht sich dies auf update 0.21; daher sind Einrichtungen für die Skew-T-Verteilung nicht enthalten. Ein Teil der Einrichtungen für die Skew-t-Verteilung ist jedoch über a set of Matlab functions zur Verfügung, die geschrieben wurde, und zur Verfügung gestellt von Enrique Batiz (Enrique.Batiz [at] postgrad.mbs.ac.uk)

Auch this code sehen Das ist in Visual Basic, sollte aber leicht tragbar sein. Relevanter Auszug, der unten gezeigt wird. Dies verwendet RandNorm (auch in der verknüpften Webseite), die ein Paar Zahlen aus einer Einheitsnormalverteilung ist, und in MATLAB sollten Sie randn(2,1) verwenden können.

Function RandSkew(fAlpha As Single, _ 
        Optional fLocation As Single = 0, _ 
        Optional fScale As Single = 1, _ 
        Optional bVolatile As Boolean = False) As Single 

    ' shg 2008-0919 
    ' http://azzalini.stat.unipd.it/SN/faq.html 

    ' Returns a random variable with skewed distribution 
    '  fAlpha  = skew 
    '  fLocation = location 
    '  fScale > 0 = scale 

    Dim sigma As Single 
    Dim afRN() As Single 
    Dim u0  As Single 
    Dim v  As Single 
    Dim u1  As Single 

    If bVolatile Then Application.Volatile 
    Randomize (Timer) 

    sigma = fAlpha/Sqr(1 + fAlpha^2) 

    afRN = RandNorm() 
    u0 = afRN(1) 
    v = afRN(2) 
    u1 = sigma * u0 + Sqr(1 - sigma^2) * v 

    RandSkew = IIf(u0 >= 0, u1, -u1) * fScale + fLocation 
End Function