T.Popoviciu,Consideraţii teoretice asupra utilizării practice a unor formule de interpolare, Bul. şt. Acad. R.P.R.,3(1951) no. 4, pp. 441-449 (in Romanian) .
Despre acest articol
Journal
Buletinul Stiintific al Academiei R.P.R.
Publisher Name
Academia Republicii S.R.
DOI
Not available yet.
Print ISSN
Not available yet.
Online ISSN
Not available yet.
Google Scholar Profile
Referințe
Lucrare in format HTML
1951 a -Popoviciu- Bul. St. Sect. St. Mat. Fiz. - Consideratii teoretice asupra utilizarii practice
CONSIDERAŢII TEORETICE ASUPRA UTILIZARIII PRACTICE A UNOR FORMULE DE INTERPOLARE
DETIBERIU POPOVICIUMEMBRU CORESPONDENT AL ACADEMIEI R.P.R.
Comunicare prezentată in şedinţa din 10 Iulie 1951
Polinomul lui Lagrange relativ la funcția f(x)f(x) și la nodurile de interpolare
acestei formule.
Aceşti coeficienţi se vor găsi pe latura superioară (descendentă) a tabloului de diferențe divizate corespunzător nodurilor luate in ordinea (4). Insă, pentru formarea acestui tablou, va fi avantajos ca nodurile (4) să fie în ordinea crescătoare a valorilor lor numerice, căci atunci toate împărţirile ce se fac cu diferenţa de noduri în vederea formării tabloului, se fac cu numere pozitive pe baza relaţiei de recurenţă
In felul acesta sunt eliminate erorile care ar putea proveni din neatentie la semnele diferenţelor de noduri.
2. Pentru simplificarea notaţiilor, să presupunem că nodurile, pe care le socotim distincte, sunt luate în ordinea crescătoare, deci că
Atunci, pe baza observatiei precedente, tabloul (3) este cel mai avantajos dintre toate cele n+1n+1 tablouri analoage obţinute permutând în toate modurile cele n+1n+1 noduri.
Observăm că nu numai formula de interpolare (2) se bucură de proprietatea că toţi coeficienţii săi sunt cuprinşi in tabloul (3). Să găsim atunci toate formulele (5) care au toţi coeficientii lor in tabloul (3). Vom zice că aceste formule apartin tabloului (3). Pentru aceasta este necesar şi suficient ca pentru orice i=0,1,dots,ni=0,1, \ldots, n punctele x_(nu_(1)),x_(nu_(2)),dots,x_(nu_(i+1))x_{\nu_{1}}, x_{\nu_{2}}, \ldots, x_{\nu_{i+1}} să fie i+1i+1 puncte consecutive (luate într'o anumită ordine) ale şirului (1), deci ca permutarea
a numerelor 1,2dots n+11,2 \ldots n+1 să se bucure de proprietatea că orice sectiune
v_(1),v_(2),dots,v_(i+1),i=0,1,dots,nv_{1}, v_{2}, \ldots, v_{i+1}, i=0,1, \ldots, n
a şirului (9) să fie formată din i+1i+1 numere naturale consecutive luate într'o anumită ordine.
Dacă convenim a spune despre o astfel de permutare (9) că este o permutare consecutivă permutării inițiale 1,2,dots,n+11,2, \ldots, n+1, vedem că:
conditia necesară şi suficientă pentru ca formula de interpolare (5) să apartină tabloului (3) este ca permutarea (9) să fie o permutare consecutivă permutării initiale 1,2,dots,n+11,2, \ldots, n+1 [permutării cu care s'a format tabloul (3)].
Să numim elemente vecine unui element al tabloului (3), acele elemente care se găsesc in coloana imediat precedentă şi în coloana imediat următoare, situate imediat deasupra si imediat dedesubtul elementului considerat. Elementele vecine elementului [x_(i),x,dots,x_(i+j);f]\left[x_{i}, x, \ldots, x_{i+j} ; f\right] sunt atunci figurate in schema
Un element are deci, în general, 4 elemente vecine. Fac exceptie numai elementele de pe laturile tabloului triunghiular şi anume : elementele [x_(1);f][x_(n+1);f]\left[x_{1} ; f\right]\left[x_{n+1} ; f\right] care au un singur element vecin, [x_(2);f],[x_(3);f],dots,[x_(n);f],[x_(1),x_(2),dots,x_(n+1);f]\left[x_{2} ; f\right],\left[x_{3} ; f\right], \ldots,\left[x_{n} ; f\right],\left[x_{1}, x_{2}, \ldots, x_{n+1} ; f\right] numai două elemente vecine, iar[x_(1),x_(2);f],[x_(1),x_(2),x_(3);f],dots,[x_(1),x_(2),dots,x_(n);f][x_(n),x_(n+1);f],[x_(n-1),x_(n),x_(n+1);f],[x_(2),x_(3),dots,x_(n+1);f]\operatorname{iar}\left[x_{1}, x_{2} ; f\right],\left[x_{1}, x_{2}, x_{3} ; f\right], \ldots,\left[x_{1}, x_{2}, \ldots, x_{n} ; f\right] \left[x_{n}, x_{n+1} ; f\right],\left[x_{n-1}, x_{n}, x_{n+1} ; f\right],\left[x_{2}, x_{3}, \ldots, x_{n+1} ; f\right] numai 3 elemente vecine. Este uşor de văzut că dacă n=1n=1 exista numai elemente cu unul sau două elemente vecine, iar dacă n=2n=2, numai elemente cu unu, două sau trei elemente vecine.
Coeficientii (8) ai formulelor (5) care aparțin tabloului (3) sunt aiunci caracterizati de proprietatea că doi coeficienţi consecutivi oarecare sunt elemente vecine in tabloul (3).
3. Numărul formulelor de interpolare (5), aparţinând tabloului (3) se poate usor determina. O formulă (5) se obţine alegând coeficienţii câte unul din fiecare coloană, doi coeficienti consecutivi find vecini. Numărul N_(n)N_{n} al acestor formule este deci egal cu numărul drumurilor ce se pot descrie plecând din prima coloană până în ultima, trecând numai prin elemente vecine. Toate aceste drumuri se termină cu elementul [x_(1),x_(2)dots,x_(n+1);f]\left[x_{1}, x_{2} \ldots, x_{n+1} ; f\right]. Dar elementele [x_(1),x_(2),dots,x_(n);f],[x_(2),x_(3),dots,x_(n+1);f]\left[x_{1}, x_{2}, \ldots, x_{n} ; f\right],\left[x_{2}, x_{3}, \ldots, x_{n+1} ; f\right] sunt vecine cu elementul din ultima coloană, iar numărul drumurilor care trec prin fiecare din aceste elemente este evident egal cu N_(n-1)N_{n-1}. Rezultă deci că N_(n)=2N_(n-1)N_{n}=2 N_{n-1}. Insă, N_(1)=2N_{1}=2 deci N_(n)=2^(n)N_{n}=2^{n}.
Putem de asemenea determina numărul formulelor (5) aparţinând tabloului (3) în care primul nod este x_(i)x_{i}. Dacă notăm acest număr cu C_(n)^(i-1)C_{n}^{i-1}, avem evident C_(n)^(o)=1,C_(n)^(n)=1C_{n}^{o}=1, C_{n}^{n}=1 si pe baza recurenței precedente
Se deduce cu uşurinţă de aici că C_(n)^(i-1)=((n)/(i-l))C_{n}^{i-1}=\binom{n}{i-l}.
Mai general, numărul formulelor (5) apartinând tabloului (3) în care primele jj noduri sunt x_(i),x_(i+1),dots,x_(i+j-1)x_{i}, x_{i+1}, \ldots, x_{i+j-1} intr'o anumită ordine oarecare, este egal ou 2^(j-1)((n-j+1)/(i-1))2^{j-1}\binom{n-j+1}{i-1}.
4. Să presupunem că (n+1)(n+1) ! se divide cu 2^(n)2^{n} şi fie atunci (n+1)!=2^(n)p(n+1)!=2^{n} p. Se poate pune însă întrebarea dacă nu pot fi găsite pp tablouri de diferenţe divizate, obținute prin permutări convenabile ale nodurilor (1), astfel ca ele să contină toate cele (n-1)(n-1) ! formule (5). In acest caz ar trebui ca fiecare din cele (n+1)(n+1) ! formule (5), să aparţină unuia și numai unui singur tablou dintre cele pp tablouri considerate. Vom arăta că acest lucru este imposibil afară de cazul banal n=1\mathrm{n}=1.
In primul rând, trebue să observăm că (n+1)(n+1) ! se divide cu 2^(n)2^{n} dacă şi numai dacă n+1n+1 este o putere a lui 2 . Intr'adevăr, se știe din teoria numerelor că factorul prim 2 figurează în ( n+1n+1 )! exact la puterea ( n > 1n>1 )
n+1-sum_(i=0)^(k)a_(0)n+1-\sum_{i=0}^{k} a_{0}
unde a_(0)a_(1)dotsa_(k)a_{0} a_{1} \ldots a_{k} este numărul n+1n+1 scris în sistemul diadic. Trebue deci să avem sum_(i=0)^(k)a_(i)=1\sum_{i=0}^{k} a_{i}=1, de unde a_(0)=1,a_(1)=a_(2)=dots=a_(k)=0a_{0}=1, a_{1}=a_{2}=\ldots=a_{k}=0, adică
Mai exact, se vede că ( n+1n+1 )! se divide cu 2^(n)2^{n} numai dacă avem (10) și atunci câtul pp este un număr impar.
Să presupunem acum că avem (10) si că n > 1n>1, deci k > 1k>1. Fie atunci pp tablouri de diferențe divizate. Dacă toate cele ( n+1n+1 )! formule (5) ar apartine la aceste pp tablouri, ar trebui în particular să aparţină acestor tablouri exact nn ! formule cu primul nod x_(1)x_{1}. Insă, oricărui tablou de diferențe divizate aparține cel puţin o formulă cu primul nod x_(1)x_{1} și am văzut că numărul formulelor cu primul nod x_(1)x_{1} care aparţine unui tablou nu poate fi decât
1,((n)/(1)),((n)/(2)),dots,((n)/(n-1))" sau "((n)/(n))1,\binom{n}{1},\binom{n}{2}, \ldots,\binom{n}{n-1} \text { sau }\binom{n}{n}
Fie atunci a numărul tablourilor, dintre cele pp considerate, cărora apartin ((n)/(i))\binom{n}{i} formule cu primul nod x_(1)x_{1}. Trebue să avem
care este imposibil, deoarece în ipoteza (10) se demonstrează că toate numerele ((n)/(1))-1,i=1,2dots,n-1\binom{n}{1}-1, i=1,2 \ldots, n-1 sunt pare, iar numărul n!-pn!-p este impar.
5. Să presupunem că nodurile (1) sunt în ordinea crescătoare, deci că x_(1) < x_(2) < dots < x_(n+1)x_{1}<x_{2}<\ldots<x_{n+1} si fie xx un punct fix al axei reale diferit de aceste noduri. Ne propunem atunci să determinăm permutarea (4) astfel ca numerele
{:(12)|(x-x_(v_(1)))(x-x_(v_(2)))dots(x-x_(v_(i)))|","i=1","2","dots","n:}\begin{equation*}
\left|\left(x-x_{v_{1}}\right)\left(x-x_{v_{2}}\right) \ldots\left(x-x_{v_{i}}\right)\right|, i=1,2, \ldots, n \tag{12}
\end{equation*}
să fie cele mai mici posibile.
In utilizarea formulei (5) a lui Newton, această proprietate se pune atunci când voim ca multiplicatorul l(x)l(x) al restului (6) să fie cel mai mic posibil în valoare absolută, atunci când oprim această formulă succesiv la 1,2,dots,n+11,2, \ldots, n+1 termeni. Intr'adevăr, avem, pe baza formulelor (5), (6), (7),
să fie nedescrescător.
Prin această proprietate, permutarea (4) nu este neapărat determinată în mod unic, dar vom demonstra că:
x fiind un punct dat al axei reale, orice permutare (4) pentru care şirul (13) este nedescrescător este o permutare consecutivă permutării inițiale 1,2,dots,n+11,2, \ldots, \mathrm{n}+1.
Proprietatea aceasta subsistă indiferent de faptul că xx coincide sau mu eu un nod.
Pentru a demonstra proprietatea va fi destul să arătăm că dacă v_(1),v_(2),dots,v_(i)v_{1}, v_{2}, \ldots, v_{i} sunt numere naturale consecutive intr'o anumită ordine şi numerele naturale v_(1),v_(2),dots,v_(i),v_(i+1)v_{1}, v_{2}, \ldots, v_{i}, v_{i+1} se bucură de aceastǎ proprietate.
( x <= (x_(1)+x_(2))/(2)x \leqq \frac{x_{1}+x_{2}}{2} pentru r=1r=1 si (x_(n)+x_(n+1))/(2) <= x\frac{x_{n}+x_{n+1}}{2} \leqq x pentru r=n+1r=n+1 )
Să presupunem acum că v_(i_(1)),v_(i_(2)),dots,v_(i)v_{i_{1}}, v_{i_{2}}, \ldots, v_{i} sunt ii numere naturale consecutive si fie s,s+1,dots,s+i-1s, s+1, \ldots, s+i-1 aceste numere. Avem atunci
(x_(s-1)+x_(s+i-1))/(2) <= x <= (x_(s)+x_(s+i))/(2)\frac{x_{s-1}+x_{s+i-1}}{2} \leqq x \leqq \frac{x_{s}+x_{s+i}}{2}
prin urmare cu atât mai mult
x_(s-1) < x < x_(s+i)x_{s-1}<x<x_{s+i}
şi pe baza observaţiei de mai sus
min quad(|x-x_(j)|)=min(|x-x_(s-1)|,|x-x_(s+i)|)\min \quad\left(\left|x-x_{j}\right|\right)=\min \left(\left|x-x_{s-1}\right|,\left|x-x_{s+i}\right|\right)
j=1,2,dots,s-1,s+i,s+i+1,dots,n+1j=1,2, \ldots, s-1, s+i, s+i+1, \ldots, n+1
adică v_(i+1)=s-1v_{i+1}=s-1 sau s+is+i, tocmai ceea ce trebuia demonstrat.
In cazul când s=1s=1, sau s+i-1=n+1s+i-1=n+1, demonstratia se modifică puţin simplificându-se și este inutil să o mai reproducem.
Dacă permutarea (4) se bucură de proprietatea că este posibil de a găsi un xx, astfel ca şirul (13) să fie nedescrescător, vom zice că această permutare sau formula (5) corespunzătoare verifică proprietatea de minimum.
6. Formulele de interpolare (5), pentru care şirul (13) devine nedescrescător pentru un xx convenabil, aparțin deci toate tabloului (3) format cu nodurile luate în ordinea crescătoare.
Nu toate cele 2^(n)2^{n} formule (5) aparţinând tabloului (3) verifică proprietatea de minimum studiată.
Se poate afirma cu siguranţă că cel puțin 2n2 n dintre aceste formule verifică proprietatea de minimum. Şi anume, dacă luăm succesiv pentru xx mijloacele intervalelor (x_(i),x_(i+1)),i=1,2,dots,n\left(x_{i}, x_{i+1}\right), i=1,2, \ldots, n, vedem că permutările
1,2,dots,n+1quad;n+1,n dots,2,11,2, \ldots, n+1 \quad ; n+1, n \ldots, 2,1
precum și cel puţin câte una dintre permutările de formele
i+1,i+2,dotsquad;quad i+1,i,dotsi+1, i+2, \ldots \quad ; \quad i+1, i, \ldots
verifică proprietatea de minimum pentru i=1,2,dots,n-1i=1,2, \ldots, n-1.
Pentru n=1,2n=1,2 avem 2n=2^(n)2 n=2^{n}, deci, în aceste cazuri, cele 2 , respectiv 4, formule aparţinând tabloului (3), verifică proprietatea de minimum.
Dacă n=3n=3 avem 2^(n)-2n=22^{n}-2 n=2. Luând succesiv pentru xx mijloacele intervalelor (x_(1),x_(2)),(x_(1),x_(3)),(x_(2),x_(4)),(x_(3),x_(4))\left(x_{1}, x_{2}\right),\left(x_{1}, x_{3}\right),\left(x_{2}, x_{4}\right),\left(x_{3}, x_{4}\right), se verifică că 6 din cele 8 permutări consecutive permutării 1,2,3,41,2,3,4 verifică întotdeauna proprietatea de minimum. Permutările rămase, care sunt 2,3,4,12,3,4,1 şi 3,2,1,43,2,1,4, verifică proprietatea de minimum dacă avem respectiv
(x_(1)+x_(4))/(2) <= x <= (x_(2)+x_(3))/(2)" si "(x_(2)+x_(3))/(2) <= x <= (x_(1)+x_(4))/(2)\frac{x_{1}+x_{4}}{2} \leqq x \leqq \frac{x_{2}+x_{3}}{2} \text { si } \frac{x_{2}+x_{3}}{2} \leqq x \leqq \frac{x_{1}+x_{4}}{2}
Se vede că, dacă x_(1)+x_(4)!=x_(2)+x_(3),7x_{1}+x_{4} \neq x_{2}+x_{3}, 7 dintre cele 8 formule apartinand tabloului (3) verifică proprietatea de minimum, iar dacă x_(1)+x_(4)=x_(2)+x_(3)x_{1}+x_{4}=x_{2}+x_{3}, toate cele 8 formule verifică proprictatea de minimum.
Dacă n >= 4n \geqq 4, nu toate cele 2^(4)2^{4} permutări consecutive lui 1,2,dots,n+11,2, \ldots, n+1 verifică proprietatea de minimum. Astfel, este imposibil ca cele patru permutări (în care este suficient să scriem primele 5 elemente)
să verifice proprietatea. Intr'adevăr, pentru cazul contrar ar trebui să avem x_(1)+x_(4)=x_(2)+x_(3)x_{1}+x_{4}=x_{2}+x_{3} si x_(1)+x_(5)=x_(2)+x_(3)x_{1}+x_{5}=x_{2}+x_{3} deci x_(4)=x_(5)x_{4}=x_{5}, ceea ce este in contradictie cu ipoteza că nodurile sunt distincte.
O analiză mai amănunțită ne arată că în cazul n=4n=4, dintre cele 16 formule (5) aparţinând tabloului (3), cel puţin 11 şi cel mult 14 verifică proprietatea de minimum.
Ar fi interesant de văzut, pentru nn oarecare, care este numărul minim ( >= 2n)(\geq 2 n) și care este numărul maxim ( < 2^(n)<2^{n} dacă n > 3n>3 ) de permutări consecutive permutării initiale care pot verifica proprietatea de minimum.
Pentru o configurație dată de noduri număıul formulelor care verifică proprietatea de minimum se poate determina uşor. Pentru determinarea acestui număr nu este nevoie să examinăm monotonia șirului (13) pentru toate valorile lui xx, ci din motive de continuitate, numai acele valori ale lui xx pentru care șirul (13) este nedescrescător fără a fi crescător. Aceasta revine la a examina valorile (în număr finit) ale lui xx, pentru care avem
Este destul deci să examinăm numai punctele xx care coincid cu mijloacele segmentelor determinate de două noduri.
In particular, dacă nodurile (1) sunt echidistante, punctele xx care intervin sunt nodurile înseşi precum şi mijloacele segmentelor determinate de câte două noduri consecutive.
O enumerare detaliată, pe care nu o mai reproducem, ne conduce la rezultatul că dacă nodurile sunt echidistante, numărul formulelor (5) aparţinând tabloului (3) şi care verifică proprietatea de minim este egal cu
unde [alpha][\alpha] este intregul cuprins în alpha\alpha.
7. Pentru aplicarea formulei de interpolare, în punctul xx se va alege acea formulă (5) care verifică proprietatea că şirul (13) este nedescrescător.
Să presupunem, în particular, că nodurile sunt simetrice fată de un punct al axei reale. Acest punct poate fi ales originea 0. Trebue atunci să distingem două cazuri, după cum nodurile sunt in număr impar sau par. Dacă presupunem că se interpolează in vecinătatea orjginii, trebue de asemenea să distingem două cazuri după cum punctul xx este la stânga sau la dreapta lui 0 .
Cazul I. Număr impar de noduri. Nodurile în ordine crescătoare se pot scrie atunci ( n=2kn=2 k )
după cum xx este la stânga sau la dreapta lui 0 .
Dacă xx este suficient de aproape de origine, acestea sunt de altfel singurele permutari care verifica proprietatea de minimum.
Formulele de interpolare corespunzătoare sunt generalizări cunoscute cehidis a căcând media aritmetică a acestor educ când nodurile sunt ralizare cunoscută a formulei de interpolare ormule, se deduce o genescriem aici explicit aceste formule.
Cazul 11. Număr par de noduri scrie atunci ( n=2k-1n=2 k-1 )
după cum xx este la stânga sau la dreapta lui 0 . Aceste permutări sunt de altfel singurele care verifică proprietatea de minimum dacă xx este suficient de aproape de origine.
Formulele de interpolare corespunzătoare sunt generalizări cunoscute ale altor formule ale lui Euler, care de altfel se reduc dacă nodurile sunt echidistante. Făcând media aritmetică a acestor două formule, obtinem o generalizare cunoscută a formulei de interpolare a lui Bessel. Inutil să mai scriem explicit aceste formule.
Proprietatea de minimum pe care am stabilit-o, constitue o justificare a utilizarii acestor formule, în cazul când se interpolează în mijlocul tabloului diferentelor divizate.
ТЕОРЕТИЧЕСКИЕ СООБРАЖЕНИЯ ОТНОСИТЕЛЬНО ПРАКТИЧЕСКОГО ПРИМЕНЕНИЯ НЕКОТОРЫХ ФОРМУЛ ИНТЕРПОЛЯЦИИ
(КРАТКОЕ СОДЕРЖАНИЕ)
Всякой нерестановке данных узлов (1) интерполяции соответствует одна формула интерполяции (5) Ньютона. Коэффициенты (8) этой формулы имеются в таблице разделенных разностей. Наиболее выгодной является таблица (3), составленная из узлов, расположенных в порядке возрастающих величин. Докавывается, что формулы интерполяции (5), для которых остаток лучше всего ограничен для всех частей этих формул, имеют свои коәффициенты, заключенные в составленных таким способом таблицах (3).
Эти замечания являются оправданием применения классических формул Эйелра, Стирлинга и Бесселя.
CONSIDÉRATIONS THEORIQUES SUR L'UTILISATION PRATIQUE DE GERTAINES FORMULES D'INTERPOLATION
(RÉSUMÉ)
A toute permutation de noeuds (1) d'interpolation donnés, correspond une formule d'interpolation (5) de Newton. Les coefficients (8) de cette formule figurent dans le tableau des différences divisées. Le tableau le plus avantageux est le tableau (3) construit avec les noeuds pris dans l'ordre de grandeur croissante. On démontre que les formules d'interpolation (5) pour lesquelles le reste est le meilleur, pour toute section de ces formules, ont leurs coefficients compris dans le tableau ainsi construit.
Ces remarques constituent une justification de l'utilisation des formules classique d'Euler, de Stirling et de Bessel.
^(1)){ }^{1)} Termenul pentru i=0i=0 se reduce la [x_(1);f]=f(x_(3))\left[x_{1} ; f\right]=f\left(x_{3}\right). Ne folosim de notatiile {:x_(1),x_(2),dots,x_(i);f]\left.x_{1}, x_{2}, \ldots, x_{i} ; f\right] pentru diferența divizată a funcției f(x)f(x) pe nodurile x_(1),x_(2),dots,x^(i)x_{1}, x_{2}, \ldots, x^{i}.