Kuliah Teknik Simulasi untuk minggu ke-3

Pembangkitan Bilangan Acak

Pseudo-Random Generation

Random yang semu

Karena barisan bilangan random dihasilkan bilangan suatu rumus :

x_n = ax_n-1 Modulo M

Misal :

a = 2, M = 5, x_o = 3

x₁ = 2(3) modulo 5 = 1

x₂ = 2(1) modulo 5 = 2

x₃ = 2(2) modulo 5 = 4

x₄ = 2(4) modulo 5 = 3

Bilangan Biner :

I = 1

II = 3

III = 7

Rumus untuk pseudo-random dapat diberikan seperti:

x_n = ax_n-1 Modulo M

x ~ u(a,b)

function uniform (a,b:real):real;

begin

uniform:=a+(b-a)x random;

end;

invers transform

F(x)= u, u ~ u(0,1)

x =...?

u ~ u(0,1)

f(x) = 1 - e^-x

F(x) = u

1 - e^-x = u

e^-x = 1 – u

Perintah ln di dalam Pascal

Function expo (mean:real):real;

Begin

expo = -loge(random)*mean;

end;

jika merandom data pada distribusi uniform dengan bantuan software minitab :

MTB > Random 1000 c1;

SUBC > Uniform 0,0 1,0.

MTB > Random 1000 c2;

SUBC > Uniform 0,0 1,0.

MTB > Plot c1*c2;

SUBC > symbol.

untuk lebih lengkapnya saya tampilkan dalam bentuk download mengenai distribusi uniform, eksponensial, dan normal. silahkan klik di sini :

http://rapidshare.com/files/279287706/Beberapa_Distribusi_Khusus.doc.html