Formula de interpolare a lui Brahmagupta este o formulă de interpolare de ordinul al doilea polinomial, găsită de matematicianul și astronomul indian Brahmagupta (598-668) la începutul secolului al VII-lea d.Hr. O descriere poetică a acestei formule în sanscrită se găsește în partea suplimentară a Khandakhodyaka, o lucrare finalizată de Brahmagupta în 665 [1] . Același cuplet se găsește în lucrarea sa anterioară Dhyana-graha-adhikara, a cărei dată exactă nu a fost stabilită. Cu toate acestea, interconectarea internă a lucrărilor sugerează că a fost creată mai devreme decât lucrarea principală a omului de știință, finalizată în 628, „ Brahma-sphuta-siddhanta ”, astfel încât poate fi atribuită crearea unei formule de interpolare de ordinul doi. până în primul sfert al secolului al VII-lea [1] . Brahmagupta a fost primul care a găsit și a folosit formula diferențelor finite de ordinul doi în istoria matematicii [2] [3] .
Formula lui Brahmagupta coincide cu formula de interpolare de ordinul doi a lui Newton , care a fost găsită (redescoperită) după mai bine de o mie de ani.
Ca astronom, Brahmagupta a fost interesat să obțină valori precise pentru sinus din numărul mic de valori tabulate cunoscute pentru această funcție. Astfel, el s-a confruntat cu sarcina de a găsi valoarea , în funcție de valorile funcției disponibile în tabel:
… | … | |||||||
… | … |
Cu condiția ca valorile funcției să fie calculate în puncte cu un pas constant , ( pentru toți ), Aryabhata a sugerat utilizarea primelor diferențe finite (tabulare) pentru calcule:
Matematicienii înainte de Brahmagupta au folosit formula de interpolare liniară evidentă
,unde .
Brahmagupta a înlocuit această formulă cu o funcție arc de diferențe finite, ceea ce face posibilă obținerea unor valori mai precise ale funcției interpolate în ordine.
În terminologia lui Brahmagupta, diferența se numește segmentul trecut (गत काण्ड), segmentul util se numește (भोग्य काण्ड). Lungimea segmentului până la punctul de interpolare în minute se numește ciot (विकल). Noua expresie care trebuie înlocuită se numește segmentul util corect (स्फुट भोग्य काण्ड). Calculul segmentului util corect este descris în cupletul [4] [1] :
Conform comentariului lui Bhuttopala ( secolul X), versurile sunt traduse după cum urmează [ 1 ] [ 5 ] : Dacă mai mult, atunci scădeți. Veți obține diferența utilă corectă [6] .
900 de minute (15 grade) este intervalul dintre argumentele valorilor din tabelul sinusului folosit de Brahmagupta.
În notația modernă, algoritmul de calcul Brahmagupta este exprimat prin formulele:
Aceasta este formula de interpolare de ordinul doi a lui Newton [7] [8] .
Nu se știe cum a obținut Brahmagupta această formulă [1] . În timpul nostru, astfel de formule sunt dovedite folosind extinderea funcțiilor în dreptul de a crește egalități într- o serie Taylor la un punct . Cu toate acestea, formula poate fi demonstrată și prin metode elementare: după înlocuire, formula Brahmagupta stabilește o parabolă care trece prin trei puncte . Pentru a obține această formulă, este suficient să găsim coeficienții acestei parabole prin rezolvarea unui sistem de trei ecuații liniare definite de aceste puncte.
Calculul computerizat arată că având un tabel de 7 valori ale sinusului la nodurile cu un pas de 15 grade, Brahmagupta ar putea calcula această funcție cu o eroare maximă de cel mult 0,0012 și o eroare medie de cel mult 0,00042.