piątek, 15 sierpnia 2014

O wyciąganiu pierwiastków

W tym odcinku o szybkim przybliżonym obliczaniu pierwiastków dla liczb niebędących kwadratami liczb naturalnych.

<Wyjaśnienie, jak to działa, jak kto chce tylko używać sposobu w praktyce, to nie jest to koniecznie potrzebne>

Wiadomo, że pochodna danej funkcji f(x) ma wzór:
f'(x)=lim f(x+Δx)-f(x)Δx
           x→0

Usuwając ograniczenie x→0 uzyskuje się w przybliżeniu:


f'(x)=f(x+Δx)-f(x)Δx, czyli:

f(x+Δx)=f(x)+f'(x)Δx,

co daje tym lepsze przybliżenie, im mniejsze jest Δx w stosunku do x.

Dla funkcji x można zapisać:

x+Δxx + Δx 12x

<Koniec wyjaśnienia>


Gdy obliczamy pierwiastek liczby niebędącej kwadratem liczby naturalnych, wyszukujemy najbliższy jej kwadrat liczby naturalnej i podstawiamy do równania:

x+Δxx12x


 Przykładowo mamy więc:


 5 = 4+14 + 1 124=2+14 = 94,
(94)2= 8116 =5+ 116, czyli błąd względny wynosi 116:5=1,25%,


  10 = 9+19 + 1 129=3+16 = 196,
(196)2= 36136 =10+ 136, czyli błąd względny wynosi 136:10=0,28%,


  3 = 4+(-1)4 +(- 1) 124=2-14 = 74,
(74)2= 4916 =3+ 116, czyli błąd względny wynosi 116:3=2,08%.


Na przykładzie 3° zauważmy jak można w prosty sposób zwiększać dokładność przybliżenia, podstawiając do wzoru za  x nie najbliższy kwadrat liczby naturalnej, ale liczbę obliczoną w poprzednim kroku i tak ad infinitum (choć w pewnym momencie trzeba mieć kartkę i długopis, żeby to zrobić):

3 = 491611674 - 116 1274= 9756,

(9756)2= 94093136 =3+ 13136, czyli błąd względny wynosi 13136:3=0,01%.

Temat można jeszcze pogłębiać, ale chyba już starczy.

Na koniec mała uwaga – w przykładzie 3° od pierwszego przybliżenia 74 do drugiego 9756 można dojść również i tak, że nowy mianownik jest iloczynem licznika i mianownika starego ułamka oraz liczby 2 a nowy licznik jest dwukrotnością kwadratu starego licznika pomniejszoną o 1:
56=7x4x2                   97=(7x7x2)-1
Na tej samej zasadzie można znaleźć kolejne przybliżenie 3:
10864=97x56x2         18817=(97x97x2)-1
(1881710864)2=3+ 1118026496

Proponuję posprawdzać dla innych liczb niż 3, także takich, gdy najbliższy kwadrat jest oddalony o więcej niż 1.

Brak komentarzy:

Prześlij komentarz