DIGITALNA ARHIVA ŠUMARSKOG LISTA
prilagođeno pretraživanje po punom tekstu




ŠUMARSKI LIST 1-2/1991 str. 58     <-- 58 -->        PDF

Do zadovoljavajućih tablica se može doći i primjenom kompjutera kako
je opisano u ovom radu.


NAČINI IZRADE TABLICA


Za simuliranje će se koristiti troznamenkasti slučajni broj od 001 do
999, a slučajnom broju 001 odgovara kumulativna frekvencija 0.001 odnosno
slučajnom broju 999 odgovara 0.999. Za simuliranje nam je potrebna varijabla
»u« koja je funkcija te kumulativne frekvencije, dakle potrebne su
nam tablice kod kojih je ulaz kumulativna frekvencija, a izlaz varijabla »u«.


Budući da takve tablice već postoje, najjednostavnije bi bilo da se cijele
unesu u kompjuter, međutim to se može izbjeći primjenom slijedećeg kratkog
programa u GWBASIC-u.


10 PI = 3.1415927: E = EXP(1):DIM U(1000)


20 DEF FN 0(U)= 1/SQR (2*PI)/E * (U A 2/2)


30MU = 500:SP = 0


40 FOR ZZ = .00005 TO 10.00005 STEP .0001:PU = FN O (ZZ)*.0001 :


SP = SP + PU:IF SP > = .001 THEN


MU=MU+1:U(MU) = ZZ:SP=SP—.001


50 NEXT ZZ


60 FOR ZZ = 1 TO 499:U(ZZ) = — U(1000—ZZ):NEXT ZZ


U liniji 20 definirana funkcija FN O (U) računa vrijednost ordinate normalne
Gaussove distribucije za neku varijablu »u« po izrazu:


Zadatak se rješava u liniji 40 gdje se uz pomoć petlje vrši strojno integriranje
s korakom 0.0001 i s početnom vrijednošću u = 0.00005. Sumiranje
ovih po volji malih površinica vrši se algoritmom SP = SP + PU sve dok se
ne postigne vrijednost veća ili jednaka 0.001. Kad se taj uvjet zadovolji,
u tom trenutku zatečena vrijednost ZZ je varijabla »u« koja se smješta u
memorijsko mjesto U(MU). Vrijednost MU/1000 je u stvari vrijednost integrala
Gaussove distribucije dakle na MU-tom memorijskom mjestu vektora
U je smještena varijabla »u«, koja odgovara integralu MU/1000.


U liniji 50 vrši se popunjavanje memorijskih mjesta od 1 do 499 tj.
dobivaju se u-varijable za vrijednosti integrala od 0.001 do 0.499 čija je vrijednost
simetrično jednaka onima od 0.999 do 0.501 samo s predznakom
»—«.


Podaci iz jednodimenzionalne matrice U se mogu odgovarajućim programom
ili naredbama ispisati u tablicu, snimiti u neku datoteku ili pak
upotrijebiti u nekom drugom programu za simuliranje u-varijable.


Postupak u nekom programu za simuliranje je jednostavan: Primjenom
RND naredbe može se »izmisliti« neki troznamenkasti broj u rasponu 001
do 999 (kombinacija 000 se izbacuje jer je za P(ll) = 0 i P(U; = 1 u = —oo