Zakon velikih števil

izrek, ki opisuje rezultat izvajanja istega poskusa zelo velikokrat

Zákon velíkih števíl je v verjetnostnem računu in statistiki osnovni limitni izrek, ki opisuje rezultat izvajanja istega poskusa zelo velikokrat. Po zakonu mora biti srednja vrednost rezultatov, pridobljenih iz velikega števila poskusov, blizu pričakovane vrednosti, njena vrednost pa se vedno bolj približuje pričakovani vrednosti, če se izvaja vedno več poskusov.[1]

Prikaz zakona velikih števil z določeno izvedbo metov ene igralne kocke. Ko se število metov v tej izvedbi veča, se srednje vrednosti vseh rezultatov približujejo vrednosti 3,5. Čeprav bo vsaka izvedba z malih številom metov (na levi) kazala razločno obliko, bo oblika večjega števila metov (na desni) skrajno podobna.

Zakon velikih števil je pomemben, saj zagotavlja stabilne dolgoročne rezultate za srednje vrednosti poljubnih naključnih dogodkov.[1][2] Čeprav bo mogoče na primer kazina v enem zasuku kolesa rulete izgubila denar, se bodo z velikim številom zasukov vrednosti njenih zaslužkov približevale napovedljivemu odstotku. Vsak zmagovalni niz igralca bodo premagali parametri igre. Pomembno si je zapomniti, da zakon velja le, (kakor nakazuje tudi njegovo ime), kadar se obravnava veliko število opazovanj. Ne obstaja načelo, da bo malo število opazovanj sovpadalo s pričakovano vrednostjo, ali, da bo zmagovalni niz ene vrednosti takoj »uravnotežen« z drugimi (glej hazarderska zmota).

Zgledi uredi

Met igralne kocke uredi

Vsak met poštene (šeststrane) igralne kocke bo dal zalogo vrednosti števil  , vsako z enako verjetnostjo   in funkcijo verjetnosti  . Tako je pričakovana vrednost povprečne vrednosti metov enaka:

 

Pričakovani standardni odklon je enak:

 

Po zakonu velikih števil, če se vrže veliko metov poštene igralne kocke, bo srednja vrednost njihovih vrednosti (včasih imenovana vzorčna sredina) zelo verjetno blizu vrednosti 3,5 z večjo točnostjo pri večjem številu metov.

Iz zakona velikih števil sledi, da bo empirična verjetnost ugodnega izida v nizu Bernoullijevih poskusov konvergirala k teoretični verjetnosti. Za Bernoullijevo slučajno spremenljivko je pričakovana vrednost teoretična verjetnost ugodnega izida, srednja vrednost   takšnih spremenljivk, kjer se privzame, da so neodvisne in enakomerno porazdeljene (n.e.p.), pa je ravno relativna frekvenca (empirična verjetnost).

Met kovanca uredi

Met poštenega kovanca je na primer Bernoullijev poskus. Tu je  . Ko se pošteni kovanec vrže enkrat, bo teoretična verjetnost, da bo padla številka, enaka  . Zato bo po zakonu velikih števil razmerje padlih številk v »velikem« številu metov kovanca v grobem enako  . Še posebej bo razmerje padlih številk po   metih skoraj gotovo konvergiralo k  , ko se bo   približeval neskončnosti.

Čeprav se razmerje padlih številk (cifer) (in grbov (glav, mož)) približuje  , bo absolutna razlika števila padlih številk in grbov skoraj gotovo postala velika, ko bo število metov postalo veliko. Zato se verjetnost, da je absolutna razlika majhno število, približuje nič, ko število metov postane veliko. Tudi razmerje med absolutno razliko in številom metov se bo skoraj gotovo približevalo ničli:

 

Pričakovana razlika intuitivno narašča, vendar v manjši meri kot število metov.

Metode Monte Carlo uredi

Drugi dober zgled je zakon velikih števil metode Monte Carlo. Te metode so širok razred izračunavalnih algoritmov, ki se za pridobitev numeričnih rezultatov zanašajo na naključno vzorčenje. Večje je število ponovitev, boljša bo aproksimacija. Razlog, da je ta metoda pomembna, je v glavnem v tem, da je včasih težko ali nemogoče uporabiti druge pristope.[3]

Omejitev uredi

Srednja vrednost rezultatov, pridobljenih z velikim številom poskusov, v nekaterim primerih morda ne bo konvergirala. Srednja vrednost   rezultatov iz Cauchyjeve porazdelitve ali nekaterih Paretovih porazdelitev (s pomembnostjo  ) na primer ne bodo konvergirale, ko bo   narastel čez vse meje – razlog je porazdelitev s težkimi repi. Cauchyjeva in Paretova porazdelitev predstavljata dva primera – Cauchyjeva porazdelitev nima pričakovanja,[4] pričakovanje Paretove porazdelitve ( ) pa je neskončno.[5] V drugem primeru so naključna števila enaka tangenti kota, enakomerno porazdeljenega med −90° in +90°. Mediana je enaka nič, pričakovana vrednost pa ne obstaja – povprečje   takšnih spremenljivk ima res enako porazdelitev kot ena takšna spremenljivka. Verjetnostno ne konvergira proti nič (ali h katerikoli drugi vrednosti), ko gre   proti neskončnosti.

Zgodovina uredi

 
Difuzija je zgled zakona velikih števil. Na začetku obstajajo raztopljene molekule na levi strani pregrade (škrlatna črta) in nobena na desni strani. Pregrada se odstrani in raztopina difundira po celotni posodi.
Na vrhu: z eno molekulo je gibanje videti precej naključno.
V sredini: z več molekulami je razvidna jasna tendenca, kjer raztopina polni posodo vedno bolj enakomerno z občasnimi naključnimi fluktuacijami.
Spodaj: z velikanskim številom raztopljenih molekul (preveč, da bi se videle), nakjučnost bistveno izgine – raztopina, kot je videti, se giblje gladko in sistematično iz območij z visoko koncentracijo v območja z nizko koncentracijo. V stvarnih razmerah lahko kemiki opišejo difuzijo kot deterministični makroskopski fenomen (glej difuzijski zakon], navkljub njeni podvrženi naključni naravi.

Italijanski matematik Gerolamo Cardano (1501–1576) je brez dokaza navedel, da se točnosti empiričnih statistik izboljšujejo z večjim številom poskusov.[6][7] To je bilo potem formalizirano kot zakon velikih števil. Cardanova matematika je pripadala obdobju v katerem se je izraz slutil s pomočjo formul. Zakona eksplicitno na ta način ni zmogel zapisati, vendar se bo po njem dogodek zgodil z vrednostjo, ki je blizu  , če je verjetnost zanj   in   veliko število ponavljanj.[7]

Posebno obliko zakona velikih števil (za dvojiško slučajno spremenljivko) je prvi dokazal Jakob Bernoulli.[8] Za razvoj dovolj strogega matematičnega dokaza je potreboval 20 let. Objavil ga je v četrtem delu svojega dela Umetnost domnevanja (Ars Conjectandi) leta 1713. To je imenoval »zlati izrek«, vendar je postal splošno znan kot »Bernoullijev izrek«.[7] Ne sme se zamenjevati z Bernoullijevim načelom, imenovanim po njegovem nečaku Danielu Bernoulliju. Verjetnost po tem ni bila le matematični abstraktni koncept, ampak je bila količina, ki se je lahko z naraščajočim številom vzorcev ocenila s povečanim zaupanjem.[7] Leta 1837 je Siméon-Denis Poisson izrek naprej opisal pod imenom »zakon velikih števil« (»la loi des grands nombres«).[a][9] Tako je bil zakon znan pod obema imenoma, vendar se »zakon velikih števil« rabi pogosteje.

Po objavi Bernoullijevih in Poissonovih prizadevanj so tudi drugi matematiki prispevali k prečiščenju zakona, med drugim Čebišov,[10] Markov, Borel, Cantelli in Kolmogorov ter Hinčin. Čebišov je podal splošno formulacijo zakona velikih števil – če so pričakovane vrednosti niza slučajnih spremenljivk in kvadrati teh pričakovanj v celoti končni, bo aritmetična sredina z njihovo rastjo zelo verjetno konvergirala k aritmetični sredini njihovih pričakovanj. Markov je pokazal, da se lahko zakon uporabi za slučajno spremenljivko, ki pod določenim drugim šibkejšim privzetkom nima končne variance, Hinčin pa je leta 1929 pokazal, da, če je vrsta sestavljena iz neodvisno enakomerno porazdeljenih slučajnih spremenljivk, je dovolj, da za veljavnost šibkega zakona velikih števil pričakovana vrednost obstaja.[11][12] Ta nadaljnja raziskovanja so dala dve pomembni obliki zakona velikih števil. Ena se imenuje »šibki« zakon, druga pa »krepki« zakon, glede na dva različna načina konvergence kumulativne vzorčne sredine k pričakovani vrednosti – še posebej, kakor je pojasnjeno spodaj – iz krepke oblike sledi šibka oblika zakona.[11][13]

Dve obliki zakona uredi

Obstajata dve različni različici zakona velikih števil, ki sta opisani spodaj. Imenujeta se krepki zakon velikih števil in šibki zakon velikih števil.[14][1] Zakona za primer kjer je   neskončno zaporedje neodvisnih in enakomerno porazdeljenih slučajnih spremenljivk integrabilnih po Lebesgu s pričakovano vrednostjo   v obeh oblikah pravita, da z navidezno gotovostjo vzorčna sredina:

 

konvergira k pričakovani vrednosti :

  (1. zakon)

(Integrabilnost po Lebesgu   pomeni, da pričakovana vrednost   obstaja glede na Lebesguov integral in je končna. Ne pomeni, da je povezana mera verjetnosti absolutno zvezna glede na Lebesguovo mero.)

Na podlagi (nepotrebnega – glej spodaj) privzetka o končni varianci   (za vsak  ) in brez korelacije med slučajnimi spremenljivkami, je varianca srednje vrednosti   slučajnih spremenljivk enaka:

 

Ta privzetek o končni varianci   ni potreben. Zaradi večje ali neskončne variance bo konvergenca počasnejša, zakon velik števil pa bo vseeno veljal. Ta privzetek se pogosto rabi, saj so dokazi lažji in krajši.

Obojestranska neodvisnost slučajnih spremenljivk se lahko zamenja z neodvistnostjo po parih v obeh različicah zakona.[15]

Razliko med krepko in šibko različico upošteva način konvergence, ki se obravnava. Za interpretacijo teh načinov glej konvergenca slučajnih spremenljivk.

Šibki zakon uredi

 
Simulacija, ki prikazuje zakon velikih števil. V vsakem koraku se vrže kovanec, ki je na eni strani rdeč, na drugi pa moder. V ustrezni stolpec se doda pika. Tortni diagram kaže razmerje med rdečo in modro barvo strani kovanca, vrženega do tedaj. V začetku se razmerje sicer zelo spreminja, potem pa se z večanjem števila metov približuje vrednosti 50 %.

Šibki zakon velikih števil (imenovan tudi Hinčinov zakon) pravi, da vzorčna sredina verjetnostno konvergira k pričakovani vrednosti:[16][17]

  (2. zakon)

Tako za poljubno pozitivno število   velja:

 

Šibki zakon ob opisu tega rezultata pravi, da bo za poljubno neničelno mejo, ne glede na to kako majhna je, z dovolj velikim vzorcem obstajala zelo velika verjetnost, da bo srednja vrednost opazovanj blizu pričakovane vrednosti – to je znotraj te meje.

Kakor je omenjeno prej, šibki zakon velja za primer neodvisnih in enakomerno porazdeljenih slučajnih spremenljivk, velja pa tudi v nekaterih drugih primerih. Varianca je lahko na primer za vsako slučajno spremenljivko v nizu različna, pri čemer je pričakovana vrednost konstantna. Če je varianca omejena, zakon velja, kakor je pokazal Čebišov že leta 1867. (Če se pričakovane vrednosti med nizom spreminjajo, se lahko zakon preprosto uporabi za povprečni odmik od ustrezne pričakovane vrednosti. Zakon potem pravi, da to verjetnostno konvergira k nič.) Dejansko dokaz Čebišova deluje tako dolgo, dokler gre varianca povprečja prvih   vrednosti proti nič, ko gre   proti neskončnosti.[12] Naj se na primer privzame, da za vsako slučajno spremenljivko v nizu velja Gaussova porazdelitev s srednjo vrednostjo enako 0, vrednostjo variance pa  , ki ni omejena. V vsakem koraku bo povprečje normalno porazdeljeno (kot povprečje množice normalno porazdeljenih sprememnljivk). Varianca vsote je enaka vsoti varianc, ki je asimptota k  . Varianca povprečja je tako asimptotična k   in gre k nič.

Obstajajo tudi primeri šibkega zakona, ki velja četudi pričakovana vrednost ne obstaja.

Krepki zakon uredi

Krepki zakon velikih števil (imenovan tudi zakon Kolmogorova) pravi, da vzorčna sredina skoraj gotovo konvergira k pričakovani vrednosti:[18]

  (3. zakon)

Tako velja:

 

To pomeni, da je verjetnost, da bo srednja vrednost opazovanj konvergirala k pričakovani vrednosti, ko gre število poskusov   proti neskončnosti, enaka 1.

Dokaz je bolj zapleten od dokaza šibkega zakona.[19] Ta zakon opravičuje intuitivno interpretacijo pričakovane vrednosti (le za Lebesguov integral) slučajne spremenljivke, ko se večkrat vzorči kot »povprečje na dolgi rok«.

Skoraj gotova konvergenca se imenuje tudi krepka konvergenca slučajnih spremenljivk. Ta različica se imenuje krepki zakon, ker bodo slučajne spremenljivke, ki konvergirajo krepko (skoraj gotovo), zagotovo konvergirale šibko (verjetnostno, v verjetnosti). Vendar je za šibki zakon znano, da velja v določenih pogojih, v katerih krepki zakon ne velja, in je tedaj konvergenca le šibka (verjetnostno). Glej #Razlike med šibkim in krepkim zakonom.

Na sam krepki zakon velikih števil se lahko gleda kot na posebni primer točkovnega ergodičnega izreka

Krepki zakon velja za neodvisne enakomerno porazdeljene slučajne spremenljivke, ki imajo pričakovano vrednost (kakor šibki zakon). To je dokazal Kolmogorov leta 1930. Velja lahko tudi v drugih primerih. Kolmogorov je leta 1933 dokazal tudi, da, če so spremenljivke neodvisne in enakomerno porazdeljene, potem je, da bo povprečje skoraj gotovo konvergiralo k nečemu (ta velja za drugo obliko krepkega zakona), potrebno, da imajo pričakovano vrednost (in bo potem seveda povprečje skoraj gotovo konvergiralo k temu).[20]

Če so sumandi neodvisni ne pa tudi enakomerno porazdeljeni, velja:

 

tako, da ima vsak   končni drugi moment,[b] in:

 

Ta izjava je znana kot krepki zakon Kolmogorova, glej na primer Sen; Singer (1993), Izrek 2.3.10.

Zgled niza v katerem velja šibki zakon ne pa krepki zakon je kadar je   enak plus ali minus   (s pričetkom v dovolj velikem  , da je imenovalec pozitiven) z verjetnostjo enako   za vsakega.[20] Varianca   je potem enaka  . Krepkizakon Kolmogorova ne velja, ker je delna vsota v njegovem kriteriju do   asimptotična k  , kar je brez meje.

Če se slučajne spremenljivke zamenja z Gaussovimi spremenljivkami z enako varianco, namreč  , bo povprečje v vsaki točki tudi normalno porazdeljeno. Širina porazdelitve povprečja se bo približevala nič (standardni odklon k  ), za dano število   obstaja verjetnost, ki z naraščajočim   ne gre proti nič, povprečje pa včasih po  -tem poskusu gre nazaj k  . Ker širina porazdelitve povprečja ni enaka nič, mora imeti pozitivno spodnjo mejo  , kar pomeni, da obstaja verjetnost enaka vsaj  , da bo povprečje doseglo   po   poskusih. To se bo zgodilo z verjetnostjo enako   pred določenim  , kar je odvisno od  . Vendar tudi po   poskusih obstaja še vedno verjetnost enaka vsaj  , da se bo to zgodilo. (To zgleda nakazuje, da je  , povprečje pa bo doseglo   neskončno mnogokrat.)

Razlike med šibkim in krepkim zakonom uredi

Šibki zakon pravi, da bo za poljubni veliki   povprečje   verjetno blizu  . Tako pušča odprto vprašanje verjetnosti, da se bo   zgodilo neskončno mnogokrat, čeprav ne v rednih intervalih. (Ni nujno, da velja   za vse  ).

Krepki zakon kaže, da se to skoraj gotovo ne bo zgodilo. Iz njega še posebej izhaja, da z verjetnostjo enako 1 za poljubni   neenakost   velja za dovolj veliki  .[21]

Krepki zakon ne velja v naslednjih primerih, šibki pa velja:[22][23][24]

1. Naj je   eksponentno porazdeljena slučajna spremenljivka s parametrom 1. Slučajna spremenljivka   po Lebesguovem integralu nima pričakovane vrednosti, vendar se s pogojno konvergenco in interpretacijo integrala kot Dirichletov integral, ki je nepravi Riemannov integral, lahko zapiše:

 

2. Naj je   geometrična porazdelitev z verjetnostjo  . Slučajna spremenljivka   v običajnem smislu nima pričakovane vrednosti, ker neskončna vrsta ni absolutno konvergentna. S pogojno konvergenco pa se lahko zapiše:

 

3. Če je zbirna funkcija verjetnosti slučajne spremenljivke enaka:

 
 

potem nima pričakovane vrednosti, šibki izrek pa velja.[25][26]

Enolični zakon velikih števil uredi

Naj je   neka funkcija definirana za   in zvezna za  . Potem bo za poljubni   zaporedje   takšno zaporedje neodvisnih in enakomerno porazdeljenih slučajnih spremenljivk, da bo vzorčna sredina tega zaporedja konvergirala verjetnostno k  . To je točkovna konvergenca (v  ).

Enolični zakon velikih števil podaja pogoje pod katerimi se konvergence v   zgodijo enolično. Če je:[27][28]

  1.   kompaktna,
  2.   zvezna v vsaki   za skoraj vse spremenljivke   in merljiva funkcija spremenljivke   v vsaki  ,
  3. potem obstaja takšna prevladujoča funkcija  , da je   in
     

Potem je   zvezna v   in:

 

Ta rezultat je uporaben za izpeljavo konsistence velikega razreda cenilk (glej ekstremalna cenilka).

Borelov zakon velikih števil uredi

Borelov zakon velikih števil iz leta 1909, imenovan po Émileu Borelu, pravi, da, če se poskus ponavlja velikokrat, neodvisno pod enakimi pogoji, je razmerje kolikokrat se poljubni določeni dogodek pojavi, približno enako verjetnosti pojavitve dogodka v poljubnem posameznem poskusu – večje je število ponavljanj, boljši bo približek. Točneje, če   označuje iskani dogodek,   verjetnost njegove pojavitve in   število kolikokrat se dogodek   pojavi v prvih   poskusih, bo z verjetnostjo 1:[29]

 

Izrek strogo podaja intuitivno predstavo verjetnosti kot dolgoročni relativni frekvenci pojavitve dogodka. Je posebni primer enega od več splošnejših zakonov velikih števil v verjetnostnem računu.

Neenakost Markova (prva neenakost Čebišova): če je   nenegativna slučajna spremenljivka, za poljubno realno število   velja:

 

Neenakost Čebišova (druga neenakost Čebišova). Neenakost ocenjuje kakšna je verjetnost, da se slučajna spremenljivka veliko razlikuje od končne pričakovane vrednosti  . Naj je   poljubna realna slučajna spremenljivka s končno pričakovano vrednostjo   in končno neničelno varianco  . Neenakost Čebišova sledi iz neenakosti Markova za  . Za vsako realno število   velja:

 

Dokaz šibkega zakona uredi

Za dano neskončno zaporedje neodvisnih in enakomerno porazdeljenih slučajnih spremenljivk   s končno pričakovano vrednostjo   je treba določiti konvergenco vzorčne sredine:

 

Šibki izrek velikih števil pravi:

Izrek:   (2. zakon)

Dokaz z neenakostjo Čebišova s privzetkom končne variance uredi

Pri tem dokazu se privzame končna varianca   (za vse  ). Neodvisnost slučajnih spremenljivk narekuje, da med njimi ni korelacij. Velja tudi:

 

Skupna srednja vrednost   zaporedja je srednja vrednost vzorčne sredine:

 

Neenakost Čebišova za   da:

 

Od tod sledi naslednje:

 

Ko se   približuje neskončnosti, vrednost izraza teži k 1. Iz definicje verjetnostne konvergence sledi:

  (2. zakon)

Dokaz s konvergenco karakterističnih funkcij uredi

Po Taylorjevem izreku za kompleksne funkcije se lahko karakteristična funkcija poljubne slučajne spremenljivke   s končno srednjo vrednostjo   zapiše kot:

 

Vse spremenljivke   imajo enako karakteristično funkcijo, zato se lahko to označi enostavno kot  .

Med osnovne značilnosti karakterističnih funkcij spadata:

 , če sta   in   neodvisni.

S tema praviloma se lahko izračuna karakteristično funkcijo  , izraženo s  :

 

Limita   je karakteristična funkcija konstantne slučajne spremenljivke   in zaradi tega po Lévyjevem izreku zveznosti   porazdelitveno konvergira k  :

 

  je konstanta, kar nakazuje, da sta porazdelitvena konvergenca k   in verjetnostna konvergenca k   enakovredni (glej konvergenca slučajnih spremenljivk.) Tako velja:

  (2. zakon)

To kaže, da vzorčna sredina konvergira k odvodu karakteristične funkcije v izhodišču vse dokler ta funkcija obstaja.

Posledice uredi

Zakon velikih števil zagotavlja pričakovanje neznane porazdelitve iz realizacije zaporedja in tudi vsako značilnost verjetnostne porazdelitve.[1] Z Borelovim zakonom velikih števil se lahko preprosto pridobi funkcija verjetnosti. Za vsak dogodek obravnavane funkcije verjetnosti se lahko aproksimira verjetnost pojavitve dogodka z razmerjem kolikokrat se poljubno določeni dogodek pojavi. Večje je število ponavljanj, boljši bo približek. Glede na zvezni primer   za mali pozitivni   tako za velik   velja:

 

S to metodo se lahko pokrije celotna os   z mrežo (velikosti  ) in tvori palčni graf, ki se imenuje histogram.

Glej tudi uredi

Opombe uredi

  1. Poisson je imenoval »zakon velikih števil« (la loi des grands nombres) v Poisson (1837), str. 7. Poskušal je podati dokaz zakona v dveh delih na straneh 139–143 in 277.
  2. Drugi centralni moment ali varianca  .

Sklici uredi

  1. 1,0 1,1 1,2 1,3 Dekking (2005), str. 181.
  2. Yao; Gao (2016).
  3. Kroese idr. (2016).
  4. Dekking (2005), str. 92.
  5. Dekking (2005), str. 63.
  6. Mlodinow (2008), str. 50.
  7. 7,0 7,1 7,2 7,3 Gorroochurn (2012), str. 16.
  8. Bernoulli (1713), § 4.
  9. Hacking (1983).
  10. Čebišov (1846).
  11. 11,0 11,1 Seneta (2013).
  12. 12,0 12,1 Prohorov, Jurij Vasiljevič, »Law of large numbers«, Encyclopedia of Mathematics (v angleščini)
  13. Stöcker (2006), str. 697
  14. Bhattacharya; Lin; Patrangenaru (2016).
  15. Etemadi (1981).
  16. Hinčin (1929).
  17. Loève (1977), § 1.4, str. 14.
  18. Loève (1977), § 17.3, str. 251.
  19. »The strong law of large numbers – What's new« (v angleščini). Terrytao.wordpress.com. Pridobljeno 9. junija 2012.
  20. 20,0 20,1 Prohorov, Jurij Vasiljevič. »Strong law of large numbers«. Encyclopedia of Mathematics (v angleščini).
  21. Ross (2009).
  22. Lehmann; Romano (2006).
  23. »A NOTE ON THE WEAK LAW OF LARGE NUMBERS FOR EXCHANGEABLE RANDOM VARIABLES« (PDF) (v angleščini). Dguvl Hun Hong and Sung Ho Lee. Arhivirano iz prvotnega spletišča (PDF) dne 1. julija 2016. Pridobljeno 28. junija 2014.
  24. »weak law of large numbers: proof using characteristic functions vs proof using truncation VARIABLES« (v angleščini).
  25. Mukherjee, Sayan. »Law of large numbers« (PDF) (v angleščini). Arhivirano iz prvotnega spletišča (PDF) dne 9. marca 2013. Pridobljeno 28. junija 2014.
  26. Geyer (2013).
  27. Newey; McFadden (1994), Lema 2.4.
  28. Jennrich (1969).
  29. Wen (1991).

Viri uredi

Zunanje povezave uredi