Ekonomi ve İstatistik Portalı » Simülasyon Teknikleri

Matlab da Gamma Dağılımından Sayı Üretme Algoritması

OkanSarioglu — Fri, 26 Jun 2009 20:39:24 +0000

function SONUC=gamma(a,b,n)

for i=1:n
B=1;
for j=1:b
u=unifrnd(0,1);
B=B*u;
end
N(i)=(-1/a)*log(B);
end
SONUC=[N];
hist(N)
end

Matlab da İntegral Hesaplayan Algoritma

OkanSarioglu — Fri, 26 Jun 2009 20:36:44 +0000

function SONUC=integral(a,b,n)
fmax = 1/3;
k=0;
for i=1:n
Ux=a+(b-a)*unifrnd(0,1);
fux =(1/3)*(x^2);
uy=fmax*unifrnd(0,1);
if uy<=fux
k=k+1;
end
end
Sonuc=(k/n)*(b-a)*fmax;
end

Buradaki fmax değeri integralini alacağınız fonksiyonun maksimum değeridir. Buradaki fonksiyon fux=(1/3)*(x^2) ‘ dir.

Matlab da Üstel Dağılımdan Sayı Üretme Algoritması

OkanSarioglu — Fri, 26 Jun 2009 20:31:10 +0000

function SONUC=ustel (m,xo,n)
a=1/(m-xo);
for i=1:n
u=unifrnd (0,1);
x(i)=xo-(1/a)*log(u);
end
SONUC=x;
end

Matlab da Poisson Dağılımından Sayı Üretme Algoritması

OkanSarioglu — Fri, 26 Jun 2009 20:30:19 +0000

function SONUC=poisson(n,alfa)
a=exp(-alfa);
for i=1:n
b=1;
j=0;
while b>=a
u=unifrnd(0,1);
while u==0
u=unifrnd(0,1);
end
b=b*u;
j=j+1;
end
X(i) = j-1;
end
SONUC=X’;
end

Matlab da Normal Dağılımdan Sayı Üretme Algoritması

OkanSarioglu — Fri, 26 Jun 2009 20:27:56 +0000

function SONUC=normal(a,b,n)
%n normal dagilimdan kaç sayi istiyoruz.
for i=1:n
toplam=0;
for j=1:12
u1=unifrnd(0,1);
toplam=toplam+u1;
end
X(i)=a+b*(toplam-6);
end
SONUC=[X];
end