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 |