MATLAB kommandoer for noen diskrete fordelinger

 

 

% GEOMETRISK FORDELING - et eksempel

%

% Vi kaster en terning til vi f?r sekser f?rste gang,

% og lar X v?re antall kast vi m? gj?re.

% Hvis vi bruker notasjonen i boka til Rice (side 40),

% er X geometrisk fordelt med p=1/6.

%

% Vi finner punktsannsynligheten (for k=1,2, ...,30)

% ved kommandoene (merk at punktsannsynligheten ogs? er

% definert for k=31,32, ....)

 

p=1/6;

k=1:30;

pk=(1-p).^(k-1).*p;

 

% Alternativt kan vi finne dem ved kommandoen

 

pk=geopdf(k-1,p);

 

% Merk at det f?rste argumentet her er k-1.

% Det kommer av at Matlab har en litt annen definisjon

% av den geometriske fordelingen enn boka til Rice.

 

% Vi kan lage et stolpediagram av punktsannsynligheten ved kommandoen

 

bar(k,pk,0.3)

 

 

 

% NEGATIV BINOMISK FORDELING - et eksempel

%

% Vi kaster en terning til vi f?r sekser for tredje gang,

% og lar X v?re antall kast vi m? gj?re.

% Hvis vi bruker notasjonen i boka til Rice (side 41),

% er X negativt binomisk fordelt med p=1/6 og r=3

%

% Vi finner punktsannsynligheten (for k=1,2,...,60)

% ved kommandoene (merk at punktsannsynligheten ogs? er

% definert for k=61,62,....)

 

p=1/6;

k=1:60;

pk=nbinpdf(k-3,3,p);

 

% Merk at det f?rste argumentet her er k-3.

% Det kommer av at Matlab har en litt annen definisjon

% av negativ binomisk fordelingen enn boka til Rice.

 

% Vi kan lage et stolpediagram av punktsannsynligheten ved kommandoen

 

bar(k,pk,0.3)

 

 

% HYPERGEOMETRISK FORDELING - et eksempel

%

% Vi tipper en lottorekke og lar X v?re antall vinnertall

% vi tipper riktig.

% Hvis vi bruker notasjonen i boka til Rice (side 42),

% er X hypergeometrisk fordelt med n=34, r=7 og m=7

%

% Vi finner punktsannsynligheten (for k=0,1,...,7)

% og lager en "tabell over den ved kommandoene

 

n=34;

r=7;

m=7;

k=0:7;

pk=hygepdf(k,n,r,m);

format short g

[k',pk']