7 Mart 2020 Cumartesi

Matlab da poisson dagılımı hesaplatan program

clc
clear
tekrar=10000;
n=100;
L=10;
s=0;
for i=1:tekrar
    x=poissrnd(L,1,n);
    ort(i)=mean(x);
    z(i)=(ort(i)-L)/sqrt(L/n); %standardize edilmiş hali
    if z(i)<=1.96 %değeri 1,96 dan küçük olanları saydırıyorum
        s=s+1;
    end
    ol(i)=s/i; %olasılığını hesapladım
end
ol
figure(1)
plot(ol)
figure(2)
t=-4:0.05:4;
hist(z,t)
%negatif binomdan elde ettik örnek artimatik ortalamsı istatistiğinin
%dağılımını elde ettik

matlab da negatif binom hesaplatan program

clc
clear
tekrar=1000;
n=100;
p=0.5;
k=5;
for i=1:tekrar
    x=nbinrnd(k,p,1,n);
    ort(i)=mean(x);
end
t=3:0.02:7
hist(ort,t)

merkezi limit teoremi ile gama fonksiyonu bulduran program

clc
clear
L=10;
k=5;
n=100;
tek=1000;
for r=1:tek
for j=1:n
for i=1:k
u=rand;
x(i)=-L*log(1-u);
end
y(j)=sum(x);
end
ort(r)=mean(y)
end
mcort=mean(ort)
mcvar=var(ort)
t=30:0.5:70;
hist(ort,t)

MATLAB DA BEKLENEN DEĞERİ VEREN PROGRAM

clc
clear
s=0;
p=0.5
n=100
tekrar=10000
for j=1:tekrar
for i=1:n
    u=rand;
    if u<=1-p
        x(i)=0;
    else
        x(i)=1;
    end
end
x
y(j)=sum(x);
ort=sum(y)/j
if y(j)<=60
    s=s+1
end
ol=s/tekrar

end
s=var(y)
t=0:0.0001:n;
y
figure(1)
plot(ort)
figure(2)
hist(y)
ol

MATLAB DA BİNOM DAĞILIMINI VEREN KOD

clc
clear
tekrar=10000
s=0;
p=0.5
n=100
for j=1:tekrar
    k=1:n
for i=1:n
    u=rand;
    if u<=1-p
        x(i)=0;
    else
        x(i)=1;
    end
end
x;
y(j)=sum(x);
if y(j)<=55
    s=s+1
end
ort(j)=sum(y)/j;
end
olasilik=s/tekrar
s=var(y)
y;
%plot(ort)
hist(y,k)

MATLAB DA NORMAL DAĞILIMI VEREN KOD

clc
clear
tekrar=10000;
a=5^5;
m=(2^35)-31;
x(1)=1;
toplamf=0;
c=-1.96
d=1.96
for i=2:tekrar
    x(i)=a*x(i-1)-floor((a*x(i-1)/m))*m;
    y(i)=x(i)/m;
    y(i)=rand;
    f=(d-c)*(1/sqrt(2*pi))*exp((-1/2)*(((d-c)*y(i)+c)^2));
    toplamf=toplamf+f;
    ola(i)=toplamf/i;
end
plot(ola)
ola=toplamf/tekrar
toplamf=0;