P R E S E K List za mlade matematike, fizike, astronome in računalnikarje ISSN 0351-6652 Letnik 26 (1998/1999) Številka 4 Strani 234-240 Marino Pavletic: VERIŽNI PRIBLIŽKI Ključne besede: matematika, teorija števil, verižni ulomki, približki. Elektronska verzija: http://www.presek.si/26/1376-Pavletic.pdf © 1999 Društvo matematikov, fizikov in astronomov Slovenije © 2010 DMFA - založništvo 234 Matematika VERIŽNI PRIBLIŽKI V matematičnih učbenikih kar mrgoli nalog, kot je tale; Poenostavi: 2 + —-— . 1+ 1 «+Š Brez posebnega truda ugotovimo, daje rezultat naloge Ko opazujemo nalogo in rezultat, nam pride na misel vprašanje, kako je pisec nalogo sestavil. Pogosta zvijača pri sestavljanju matematičnih nalog je, da začnemo pri koncu in rezultat preoblikujemo tako dolgo, dokler se nam izraz ne zdi dovolj zapleten. Takega zastavimo kot nalogo za učence. Začnimo torej z ulomkom x — |g- Kot vidimo, ga lahko razdelimo na celi in ulomljeni del 54 + le; 19 "2 + 19' Ulomljeni del zdaj spremenimo v dvojni ulomek 16 Pod glavno ulomkovo črto smo dobili ulomek ki ga spet lahko razdelimo na celi in ulomljeni del = 1 -I- fg. Prvotni ulomek je torej enak 1=2+- 1 1 + A 16 Postopek lahko nadaljujemo. Zadnji dobljeni ulomljeni del ^ spet spremenimo v dvojni ulomek in dobimo s = 2+-L 3 Če še delimo na celi in ulomljeni del 5+ lahko z: izrazimo kol a = 2 + --. (1) 1 + To pa je izraz, zapisan na začetku. Matematika 235 Sestavljalcn naloge se je zdel ta izraz očitno že dovolj zapleten, mi pa se vprašajmo, ali bi lahko s postopkom nadaljevali. Kot vidimo, smo uporabljali izmenično dva koraka: A, ulomek (večji od 1) smo razdelili na celi in ulomljeni del (ulomek med 0 in 1), B, dobljeni ulomljeni del srno spremenili v dvojni ulomek oziroma poiskali smo njegovo recipročno vrednost. Ker je le-ta večja od 1, smo lahko spet uporabili korak A. Zadnji ulomljeni del, ki smo ga dobili v našem primeru, je enak i in njegova recipročna vrednost je 3, torej celo število. Ce poskusimo spet uporabiti korak A, dobimo neceli del enak 0, kar pomeni, da postopka ne moremo več nadaljevati. Z uporabo enakosti 3 = 2 + | lahko izsilimo še en korak, ki nas pripelje do pravilnega, a nič kaj lepega zapisa z = 2 + ---K.--. (2) 5+ ■ 1 Pravi matematik zdajle že premišljuje, kako bi se dalo nova spoznanja posplošiti. Tudi nam po glavi rojijo vprašanja: AH lahko podoben postopek uporabimo tudi za druge ulomke? Ali morda tudi za ostala števila? Kako bi lahko dobljeni rezultat uporabili v praksi? Da potešimo radovednost, povejmo nekaj osnovnih ugotovitev, do katerih so se prikopali veliki matematiki pred nami. Ugotovili so, da lahko postopek izvajamo tudi na splošno. S ponavljanjem korakov A in B lahko vsak ulomek u pretvorimo v obliko ti = b0 + ------^---. (3) b i + i>2 + 1 bn Pravimo, da smo število u razvili v verižni ulomek. Zaradi načina dela (korak B) so vsi števci enaki 1 in vrednost verižnega ulomka določajo števila 60,61,___,bn, ki jim pravimo verižni koeficienti. Ker je zapis verižnega ulomka v običajni matematični obliki (3) precej nepregleden, ga na kratko zapišemo tako, da navedemo njegove verižne koeficiente: u = = [6q,&i,..., 6„]. Torej bomo število, s katerim smo začeli to zgodbo, pisali kot ^ = [2,1,5,3] (ali pa, če upoštevamo drugo, gršo možnost, f§ = [2,1,5,2,1]). Velja izrek, ki ga navajamo brez dokaza: Izrek. Poljubno pozitivno racionalno število lahko razvijemo v verižni ulomek na natanko dva načina. Oba načina smo že spoznali, to sta obliki (1) in (2). Običajno uporabljamo "lepši" način - tisti, ki se ne konča z j. Cc je verižni ulomek preprost, ga lahko spremenimo nazaj v običajni ulomek že z nekaj preprostimi računskimi operacijami - računati moramo od spodaj navzgor. Težave nastopijo, če je verižni ulomek zapleteuejši. Za ta primer so matematiki izdelali posebno shemo, ki nam omogoča hitrejše računanje. Oglejmo si jo na našem zgledu: verižni ulomek [2,1,5,3] bomo poskusili preoblikovati v . Najprej napišimo shemo 2 15 3 0 1 1 0 Števila 2, 1, 5, 3 v prvi vrstici so verižni koeficienti našega števila, ničle in enke na levi strani pa so začetne vrednosti, ki so potrebne za pravilno delovanje sheme (njihova razporeditev je enaka za vse ulomke). Zdaj moramo napolniti desni del spodnjih dveh vrstic. Ravnali se bomo po principu a+b-c. Ce je v zgornji (prvi) vrstici število c, v spodnji (drugi ali tretji) pa sta v prejšnjih stolpcih števili a in b, bomo v prazen prostorček pod c vpisali vrednost a + b-c. Enak postopek uporabljamo za drugo in za tretjo vrstico, paziti je treba le, da število c obakrat, vzamemo iz prve vrstice. Tako bomo pod število 2 vpisali v drugo vrstico 0 +1 ■ 2 = 2, v tretjo pa 1 + 0 ■ 2 = 1 2 5 3 0 1 2 1 0 1 Enako računamo do konca 2 1 5 3 0 1 2 3 17 54 1 0 1 1 6 19 V zadnjem stolpcu že vidimo števili 54 in 19, torej se rezultat skriva v zadnjem stolpcu. Zapisan jc že v obliki ulomka. manjka le ulomkova črta. Tudi ostale pare števil si lahko predstavljamo kot ulomke - dodajmo še ulomkove črte 2 1 5 3 0 1 2 ,3 17 54 1 0 1 1 6 19 Poleg rezultata smo dobili še tri ulomke: = f = 2, = j = 3 in X2 = Pomožnim ulomkom, ki jih dobimo z opisano shemo, bomo rekli verižni približki danega števila. Verižne približke dobimo pravzaprav tako, da zadnjih nekaj verižnih koeficientov v razvoju [bo, bi,,,., &n] izpustimo. V našem primeru je 2 — [2], 3= [2,1], g7-[2,1,5], ^ — [2,1,5,3]. Verižni približki so za matematiko zanimivi, ker so zelo dobri približki. Med ulornki. ki imajo imenovalec manjši ali enak 6, je nlomek najboljši približek za vrednost To velja tudi splošneje, z eno samo majhno izjemo: prvi verižni približek je celi del danega števila in zato ni nujno najboljši celoštevilski približek. Tudi v našein primeru prvi približek (j"o = 2) ni najboljši celoštevilski približek - drugi približek (i] — 3) je boljši. Izrek. Vsak verižni približek danega števila, razen morda prvega, je boljši približek tega števila kot katerikoli drug ulotnek z manjšim ali enakim imenovalcem. Vsak nlomek lahko zapišemo tudi v decimalni obliki - tako trdi teorija. V praksi pa se pojavijo težave. Decimalni zapis ima skoraj vedno neskončno mnogo decimalk in iz čisto praktičnih razlogov ga ponavadi zaokrožimo na nekaj mest. Poglejmo si naslednjo nalogo. Pri računanju z računalom se na za-slončku izpiše število y = 2.358208955. To seveda ni pravi rezultat računa, ki smo ga želeli izračunati, jc le decimalni približek. Ali bi se dalo ta decimalni približek nekako nadomestiti z ulomkom? Pomagali si bomo s postopkom razvijanja v verižni ulomek. Koraka A in B namreč lahko izvajamo tudi, če je število zapisano v decimalni obliki. Pa poglejmo: y = 2.358208955 = 2 + 0.358208955 = 2 + 1 2.791666668 =2+-—„ * _____— 2 + 1 2 + 0.791666668 " 2 , 1 : 1,263157892 Vidimo, da moramo na računalu samo izmenično uporabljati tipki FRAC (— ulomJjeni del)1 in 1 /x (= recipročna vrednost). Pri tem ne pozabimo zapisovati dobljenih celih delov, to so namreč verižni koeficienti. Po nekaj korakih dobimo naslednji razvoj y = {2,2,1, 3,1,4, 997446,1,1, 6, 3,...] . Iz koeficientov dobimo naslednje verižne približke [2] =2 [2,21 = 1 [2,2,1] = l [2,2,1,3] = jj [2,2,1,3,1]=^ 1 58 [2,2,1,3,1,4] = [2 2 1 3 1 4 9974461 = 157596501 ■' ' ' 'J 66828896 Žepno računalo nam pokaže, da se že ulomek ^ ujema s številom y na vsa mesta, kar jih je na zaslončku, To velja tudi za zadnji ulomek, ki pa ima neprimerno večji števec in imenovalec. Ulomek -J^ je torej precej lepši približek za rezultat računa, ki smo ga želeli izračunati. 1 Na nekaterih računalih te tipke ni. Pomagamo si tako, da celi del kar odštejemo. Ker ostalih decimalk rezultata ne poznamo, ne moremo vedeti, če se ga sploh da pisati kot ulomek; zgornji račun pa nam kaže, da lahko za praktične potrebe privzamemo, da je rezultat skoraj enak j0, Celo več: ujemanje je tako lepo, daje možno, daje pravi rezultat enak naslednje verižne koeficiente pa smo dobili le zaradi napak pri zaokrožanju. Tisti, ki se rnu zdi, daje število 997446 preveliko za napako pri zaokrožanju, naj se spomni, da smo ga dobili kot recipročno vrednost števila 0,00000100256. Pojav tako velikega verižnega koeficienta nas navadno opozarja, da smo prišli do meje natančnosti, ki jo lahko dosežemo z računalom: naslednji koeficienti so praviloma napačni (zaradi napak pri zaokrožanju). Poskusimo podobno še s številom z — 1.414213562. Hitro dobimo razvoj ^ = [1,2,2,2, 2,2,2,...] in ustrezne verižne približke 3 7 41 99 ' 2' 5' 29' 70; "' Število 2 je desetmestni približek Števila \/2. Števila \/2 ne moremo zapisati z ulomkom — pravimo, da je iracionalno. Zato se verižni koeficienti nadaljujejo v neskončnost. Če bi bilo koeficientov končno mnogo, bi bil zadnji verižni približek neki ulomek, enak številu \/2, kar ni mogoče. Tudi na splošno velja: Izrek. V verižni ulomek lahko razvijemo poljubno pozitivno realno število. Razvoj racionalnih števil je končen, razvoj iracionalnih števil pa neskončen. Opomba: Zgled «/2 = [1,2,2,2,2,...] potrebuje še malo razlage. Kdor je zadosti vztrajen, bo prišel po nekaj korakih računanja z računalom tudi do koeficientov, različnih od 2. To je posledica napak pri zaokrožanju, pri pravem razvoju se nadaljujejo v neskončnost same dvojke. Bralca vabim, da. poskusi napisati računalniški program za (približno) pretvarjanje decimalnih števil v ulonike. Program naj izpiše vse verižne približke, ki jih je pri danem številu smiselno izpisati. To pomeni, da mora upoštevati, ria koliko mest natančni so podatki, in nehati z ¡¡spisom, ko pridejo na vrsto verižni koeficienti, pri katerih se napake pri zaokrožanju preveč poznajo. Za konec si oglejmo še verižne približke za število 7r, Ce izhajamo iz približka ~ = 3.141592654 in računamo na deset mest (z računalom), dobimo razvoj tt=[3, 7,15,1,292,1,1,1,2,1,1,4, ...j. Temu razvoju pripadajo verižni približki 22 333 355 104348 ' 7' Tog1 m' 33215' Nadaljnji verižni približki se razlikujejo od prave vrednosti šele na desetem decimalnem mestu. Sklepamo, da naslednji verižni koeficienti zaradi napak pri zaokrožanju najbrž niso natančni. Dober računalniški program bo torej izpisal samo navedene približke. Marino Pavletič