T. Popoviciu,Unele aspecte ale problemei preciziei în calculele numerice,Bul. mat. al Soc. Şt. Mat. Fiz. din R.P.R.,1(49) (1957) no. 4, pp. 473-478 (in Romanian).
Despre acest articol
Journal
Bull. Math. de la Soc. Sci, Math. Phys. de la R.P.R
Publisher Name
DOI
Not available yet.
Print ISSN
Not available yet.
Online ISSN
Not available yet.
Google Scholar Profile
Referințe
[1] A. Cauchy, Sur les fonctions interpolaires, C. R. Acad. Sci. Paris 11 (1840), 775.
[2] D. V. Ionescu, Formule de cubatură în care domeniul de integrare este un triunghi oarecare, Stud. Cerc. Științ., Filiala Cluj VI (1955), 7.
[3] E. Moldovan, O generalizare a noțiunii de convexitate, Stud. Cerc. Științ., Filiala Cluj VI (1955), 65.
[4] М. И. Морозов, О некоторых вопросах равномерного приближения непрерывных функций посредством функций интерполяционных классов, Изв. АН СССР 16 (1952), 75.
[5] T. Popoviciu, Asupra formei restului în unele formule de aproximare ale analizei, Lucrările Sesiunii generale științifice, 2–12 iunie 1950, Ed. Acad. R.P.R., București, 1951, p. 183.
[6] T. Popoviciu, Considerații teoretice asupra utilizării practice a unor formule de interpolare, Bul. Științ. Acad. R.P.R., Secț. Șt. Mat. Fiz. XII, 2 (1951), 441.
[7] T. Popoviciu, Folytonos függvények középértéklételeiről, A Magyar Tud. Akad. II Oszt. Közl. IV (1954), 353.
[8] T. Popoviciu, Despre precizia calculului numeric în interpolarea prin polinomul lui Newton pe noduri echidistante, Stud. Cerc. Științ., Filiala Cluj VI (1956), 20.
[9] T. Popoviciu, Despre precizia calculului numeric în interpolarea prin polinoame, Bul. Științ. Acad. R.P.R., Secț. Șt. Mat. Fiz. VII, 4 (1955), 953.
[10] J. Radon, Restausdrücke bei Interpolation und Quadraturformeln durch bestimmte Integralen, Monatsh. Math. Phys. 42 (1935), 389.
[11] L. Tornheim, On n-parameter Families of Functions and Associated Convex Functions, Trans. Amer. Math. Soc. 69 (1950), 457.
Paper (preprint) in HTML form
1957 b -Popoviciu- Bull. Math. - Unele aspecte ale problemei preciziei in calculele numerice
UNELE ASPECTE ALE PROBLEMEI PRECIZIEI ÎN CALGULELE NUMERICE
DETIBERIU POPOVICIU (Claj)
(Comunicare făcută la Congresul matematicienilor romîni, București, mai-iunie 1956)
Din imboldul necesităților practice, analiza numerică a luat o mare dezvotare în ultimele decenii. Problemele analizei numerice sînt foarte variate, dar toate se pot caracteriza prin aceea că prin ele se urmăreşte obținerea rezultatelor pînă la ultimele calcule numerice.
În cele ce urmează vom treee în revistă, foarte rezumativ, cîteva aspecte cu totul particulare ale unor probleme de analiză numerică, de altfel de mare importanță teoretică şi practică, aspecte care întră în cadrul preocupărilor Secției de matematică a Filialei din Cluj a Academiei R. P. Romine.
2. În multe probleme de analiză numerică se cere să se obţină un număr care depinde de o funcție ff, deci se cere valoarea unei funcționale A[f]A[f] de argument f. Pentru fixarea ideilor și pentru simplificare vom presupune deocamdată că ff este o funcție reală de o variabilă reală. Funcționala A[f]A[f] poate avea diferite structuri. Cazurile cele mai simple le întîlnim în calculul valorii funcției ff sau a derivatei sale de un ordin oarecare pe un punct dat, precum și în calculul integralei lui / într-un interval dat etc. Aceste din urmă cazuri particulare se întînesc frecvent, de exemplu la diferite metode de integrare numerică a ecuaţiilor diferențiale sau a altor ecuații funcționale.
Din cauza prea marei complicații sau a insuficientei cumoaşteri a funcției ff, în general, nu putem obține exact numărul A[f]A[f]. Trebuie atunci să ne mulțumim eu o aproximatie a acestui număr. Ceea ce face ca asemenea problemo de calcul aproximativ să fie importante este că în practică o aproximație convenabilă a numărului A[f]A[f] este totdeauna suficientă.
Există diferite procedee care permit să obținem o valoare aproximativă pentru AA [f]. Un procedeu destul de general se poate schematiza prin egalitatea aproximativă
deci se poate lua ca o aproximație pentru A[f]A[f] valoarea B[f]B[f] a unei alte functionale date, considerată «mai simplă» decît A[f]A[f], în sensul că efectiv calculul numărului B[f]B[f] se poate executa pînă la capăt și în anumite condiții impuse dinainte.
Un procedeu general pentru obținerea unei aproximații pentru funcționala A[f]A[f] constă în a înlocui funcția ff printr-o funcție de aproximare varphi\varphi a sa,
{:(2)t~~varphi.:}\begin{equation*}
t \approx \varphi . \tag{2}
\end{equation*}
Se poate lua atunci numărul A[varphi]A[\varphi] ca 0 aproximație pentru A[f]A[f], deci
De altfel egalitatea (3) este de forma (1) dacă funcția de aproximare varphi\varphi depinde univoc de funcţia ff, deci dacă varphi\varphi este valoarea unui operator
În cele ce urmează ne vom ocupa numai de astfel de formule de aproximare.
3. În utilizarea practică a unei aproximații de forma (1) este deosebit de important să se aprecieze croarea care se comite prin această aproximare. In căutarea unei astfel de aprecieri o delimitare convenabilă a restului
a formulei (1) este de mare importanță. In sfîrşit, pentru a obține o astfel de delimitare un studiu aprofundat al structurii restului R[f]R[f] este foarte util.
Dacă, de exemplu, R[f]R[f] este o funcțională liniară (aditivă și omogenă), în anumite condiții destul de generale se poate exprima R[f]R[f] eu ajutorul unei integrale (integrala lui Stieltjes dacă sîntem în condițiile teoremei lui F. Riesz). O astfel de exprimare a restului nu este suficientă, deoarece în delimitările practice ale restului intervin proprietăți ca acele care se referă la gradul de exactitate al formulei de aproximare (sau al restului ei). Gradul de exactitate (față de un polinom) a unei functionale liniare R[f]R[f] se defineşte ca fiind cel mai mare număr întreg nn (dacă există), astfel ca R[f]=0R[f]=0, pentru orice polinom ff de gradul nn. In acest caz restul se poate exprima printr-o funcţională liniară de derivata a (n+1)(n+1)-a a lui ff, dacà bineînțeles anumite condiți de existentă a acestor expresii subzistă. Diverse expresii de acest tip au fost studiate de mulți autori, mai cu seamă pentru restul formulelor de derivare şi de cuadratură numerică.
4. In cazul unei functionale liniare R[f]R[f], teoria funcților convexe de ordin superior ne permite să găsim o formă remarcabilă a lui R[f]R[f]. Dacă R[f]R[f] este de gradul de exactitate nn şi dacă ff este o funcție continuă într-un interval, avem
In această relatie, valabilă in condiții foarte generale, unde alpha,beta\alpha, \beta sint constante independente de funcția ff, am întrebuintat simbolul cunoscut al diferentelor divizate, iar xi_(1),xi_(2),dots,xi_(n+2)\xi_{1}, \xi_{2}, \ldots, \xi_{n+2} si xi_(1),xi_(2),dotsxi_(n+2)\xi_{1}, \xi_{2}, \ldots \xi_{n+2} sînt grupuri de cîte n+2n+2 puncte distincte care depind în general de funcția t,[5]t,[5].
Este deosebit de interesant cazul cînd unul din numerele alpha,beta\alpha, \beta se poate lua egal cu zero. Acest caz are loc dacă şi numai dacă avem R[f]!=0R[f] \neq 0 pentru orice functie (continuă şi) convexă de ordinul nn.
Forma (7) a restului precizează mult structura lui, legîndu-l de proprietățile diferențiale ale funcţiei ff. Legătura dintre diferențele divizate şi dintre proprieta-
țile diferențiale ale unei funcţii este asigurată de o scrio de proprietăți, aşa-zise de medie, dintre care cele mai simple sînt cunoscute deja de la Cauchy [1] iar altele au fost stabilite mai recent [7].
5. Nu insistăm mai pe larg asupra structurii restului R[f]R[f]; vom face numai cîteva observații.
a) Structura (7) a restului R[f]R[f], despre care presupunem că este o funcțională liniară, se poate studia într-un cadru mai general cu ajutorul unci convenabile extinderi a noțiunii de grad de exactitate. Această generalizare constă în a presupune că R[f]R[f] se anulează pe învelitoarea liniară a unui sistem format dintr-un număr finit de funcții liniar independente. Bineînțeles se întrebuințează generalizarea corespunzătoare a noțiunii de diferență divizată. Pe lîngă această condiție restul a fost deja obținut sub formă de integrală de către J. R a d o n [10].
b) Pînă în prézent nu se cunoaşte nimic analog pentru cazul cînd R[f]R[f] nu mai este o funcțională liniară. Asemenea funcționale ca resturi par a fi totuşi importante. De exemplu metoda de integrare numerică a lui Runge-Kutta a ecuaților diferentiale revine tocmai la un procedeu de aproximare de tipul general semnalat, în care însă R[f]R[f] nu este o funcțională liniară. In acest caz /// este o funcțe de două sau mai multe variabile independente. Faptul că restul în metoda lui Rung eKutta a fost puțin studiat pînă acum se datoreşte structurii sale complicate.
Se pare că problema structurii restului R[f]R[f], în cazul cînd R[f]R[f] nu mai este o funcțională liniară de jj, este în legătură cu generalizările noțiunii de convexitate față de un sistem interpolator de funcții. In accastă din urmă direcție trebuie să amintim lucrările lui M. I. Morozov [4], L. Tornheim [11], E. Moldovan [3].
c) În exemplul citat al metodei lui Runge-Kutta restul depinde de o funcție de mai multe variabile. Restul formulelor de aproximare de forma (1) sau (3), în care ff este o funcție de mai multe variabile reale, a fost mult mai puțin studiat pînă acum. Bineînțeles, se cunose numeroase rezultate cu privire la restul formulelor de derivare şi de cuadratură (cubatură) numerică. La noi în țară, în această direcție, trebuie să amintim rezultatele lui D. V. Ionescu [2]. Dar chiar în cazul cînd R[f]R[f] este o funcțională liniară, o formulă analoagă cu (7) nu a fost încă studiată în gencral. Facem, bincințeles, abstracție de cazurile simple cînd rezultatele se obțin plecind de la cazul unei singure variabile și procedînd apoi prin suprajunere succesivă a variabilelor.
6. Pentru a obține o aproximație a numărului A[f]A[f] trebuie să calculăm efectiv numărul B[f]B[f]. Să presupunem că acest număr este dat de formula (5) prin intermediul funcției de aproximare (4) a lui ff. Pentru calculul efectiv al lui B[f]B[f] trebuie efectuate un număr oarecare de operații care constau în general dintr-un număr finit de operații aritmetice elementare (adunare, scădere, înmulțire şi împărţire).
Pentru executarea calculelor, indiferent dacă ele se fac direct, eu ajutorul unor instrumente de calcul sau cu maşina, se urmează totdeauna un program care, în particular, precizează ordinea în care sînt executate diferitele operații.
La efectuarea fiecărei operații (elementare) și în special la înmulțiri şi împărțiri, se comit în general erori din cauză că în practică folosim numai anumite numere, aşa-zisele numere practice care se exprimă de obicei prin fracții zecimale limitate sau chiar limitate la un anumit număr de cifre zecimale. In felul acesta în cursul efectuării calculului lui B[f]B[f] se acumuleagă mai multe erori de calcul care vor
afceta rezultatul tinal urmărit. Aproximatia astfel calculată a lui B[f]B[f] va fi în definitiv aproximația efectiv obținută a numărului A[f]A[f].
7. În practica calculului numeric, stabilirea, în fiecare caz, a unui program de calcul anumit joacă un mare rol. Din punct de vedere teoretic asemenea programe se vor putca studia, cel puțin în cazurile concrete cele mai simple, cu ajutorul diferitcor expresii care reprezinta pe B[f]B[f].
00, Pentru a vedea întrucîtva varietatea maic de probleme care se pun în leğtură cu aceste considerații vom expune un caz concret simplu.
tom Să prosupunem că (4) este polinomul de interpolare al lui Lagrange varphi(x)=L(x_(1),x_(2),dots,x_(n+1);f∣x)\varphi(x)=L\left(x_{1}, x_{2}, \ldots, x_{n+1} ; f \mid x\right),
relativ la functia ff si pe nodurile de interpolare (distincte) x_(1),x_(2),dots,x_(n+1)x_{1}, x_{2}, \ldots, x_{n+1}, deci este polinomul de gradul nn care coincide cu funcția ff pe nodurile x_(i)x_{i}.
Să presupunem apoi că este vorba de problema interpolării cu ajutorul polinomului (8) al lui Lagrange, deci de calculul valorii f(x_(0))f\left(x_{0}\right) a funcţiei ff pe punctul x_(0)x_{0} (presupus diferit de noduri), cunoscind valorile funcției pe nodurile de interpolare. Sintem deci in cazul cînd A[t]=f(x_(0))A[t]=f\left(x_{0}\right) şi
Stud Calculul valorii (9) pe punctul x_(0)x_{0} a polinomului (8) se poate face în diferite feluri, după diferitele programe de calcul care pot fi indicate şi studiate prin expresiile explicite ale polinomului lui Lagrange.
Bre 8. Fie formula
unde indicele ii la simbolul Pi(^((i)):}\Pi\left({ }^{(i)}\right. arată că valoarea ii a lui jj este exclusă. Prin formula (10) se poate indica, de exemplu, că trebuie efectuată o sumă de n+1n+1 termem în prealabil calculați. Fiecare termen al sumei se calculcază în general prin 2n-12 n-1 îumulțiri şi o împărțire. Este de observat că profitînd de rezultatele parțiale, numărul înmulțirilor se poate reduce. Astfel pentru obținerea produselor
nu este nevoie de (n+1)(n-1)=n^(2)-1(n+1)(n-1)=n^{2}-1 înmulțiri, ci, în general, numai de 3n=33 n=3 înmulțiri. O asemenea reducere a numărului unor operațiuni de efectuat prezintă un mare interes în practică, bineînțeles în legătură și cu alte condiții pe care, în general, un program de calcul trebuie să le îndeplinească.
O altă observație, care se poate face relativ la calculul precedent, este că în obținerea ficeărui termen înmulţirea prin f(x_(i))f\left(x_{i}\right) poate să fie ulterioară sau anterioatč operatiei de împărțire, care trebuie in general efectuată pentru gasirea accstui termen.
Primul mod de a proceda, care constă în a calcula, în prealabil, coeficienții lui f(x_(i)),(i=1,2,dots,n+1)f\left(x_{i}\right),(i=1,2, \ldots, n+1) din formula (10) este avantajos dacă trebuie să aplicăm formula (10) la maj multe funcţii cu aceleaşi nodựi x_(i)x_{i} şi acelaşi punct x_(0)x_{0}.
Al doilea procedeu, care constă în lăsarea la urmă a cfectuării împărțirii în calculul termenilor (eventual numai a unora dintre termenii) sumei (10), poate prezenta avantaje în cazul unor date particulare simple, permițînd să se profite de anumite simplificări în cursul calculului, ca de exemplu de simplificări de fracții ordinare etc. si care nu sînt avantajoase în cazul primului mod de a organiza calculele. Trebuie să subliniem însă că asemenea simplificări nu pot fi prevăzute teoretic si nu intră în general în cadrul preocupărilor teoretice ale calculului iumeric.
9. Se cunose avantajele pe care le prezintă, în problema interpolării prin polinoame, formula lui Nowton
{:[L(x_(1),x_(2),dots,x_(n+1);f∣x_(0))=],[(12)=sum_(i=0)^(n)(x_(0)-x_(1))(x_(0)-x_(2))dots(x_(0)-x_(i))[x_(1),x_(2),dots,x_(i+1);f∣:}]:}\begin{gather*}
L\left(x_{1}, x_{2}, \ldots, x_{n+1} ; f \mid x_{0}\right)= \\
=\sum_{i=0}^{n}\left(x_{0}-x_{1}\right)\left(x_{0}-x_{2}\right) \ldots\left(x_{0}-x_{i}\right)\left[x_{1}, x_{2}, \ldots, x_{i+1} ; f \mid\right. \tag{12}
\end{gather*}
(primul termen find f(x_(1))f\left(x_{1}\right) ), unde [x_(1),x_(2),dots,x_(i+1);f]\left[x_{1}, x_{2}, \ldots, x_{i+1} ; f\right] este diferența divizată a functiei ff pe nodurile x_(1),x_(2),dots,x_(i+1)x_{1}, x_{2}, \ldots, x_{i+1}. Programul de calcul pe care-1 maca de obiceiu formula (12) constă în obținerea succesivă a numerelor y_(1),y_(2),dots,y_(n)y_{1}, y_{2}, \ldots, y_{n}, unde
care comportă în general nn înmulțiri, deci mai puține (pentru nn suficient de mare) decit programul indicat pentru formula (10). Nu trebuie însă să uităm că aici au fost calculate în prealabil diferențele divizate [x_(1),x_(2),dots,x_(i+1);f],(i=1,2,dots,n)\left[x_{1}, x_{2}, \ldots, x_{i+1} ; f\right],(i=1,2, \ldots, n). Calculul acestor diferențe divizate, de exemplu prin formarea tabloului diferențelor divizate, comportă în general erori de natura celor semnalate din cauza împărțirilor cu diferențe de noduri care intervin. Este cazul atunci, aşa cum am făcut cu altă ocazie [9], să se examineze influența erorilor de calcul, comise la formarea tabloului diferentelor divizate, asupra programului de calcul indicat pentru formula (12).
10. În încheiere și relativ la programul do calcul semnalat pentru formula (12) mai putem face următoarele observații.
a) Căutînd să reducem la minimum clorile de calcul care se acumulează prin acest program, găsim o justificare a întrebuințării diforitelor formule particulare de interpolare, după poritia punctului x_(0)x_{0} printre noduri. In felul acesta şi pe această cale se justifică, după cum am arătat altă dată [6], apicarea in practica a formulei lui Newton cu diferente ascendente sau descendente, a formulelor lui Euler, Stirling, Bessel, etc., după cum punctul x_(0)x_{0} este în apropierea extremităților, în apropierea centrului ete., a tabloului diferențelor (nodurile sînt presupuse echidistante și aranjate în ordine crescătoare sau descrescătoare).
b) In cazul nodurilor echidistante, în loc de diferențe divizate se întrebuințează diferențe obişnuite a căror formare nu comportă împărțiri. Se poate dar ádmite în general că tabloul diferențelor nu comportă erori de calcul. Am arătat într-o altă lucrare [8] că pentru programul (13) modificat în acest sens se poate preciza o delimitare a erorilor de calcul acumulate, întocmind şi tablouri suficiente în practică pentru găsirea rapidă a unor astfel de delimitări.
BIBLIOGRAFIE
A. Cauch y, Sur les fonctions interpolaires, C. R. de l'Acad. Sci. Paris 11 (1840), 775.
D. V. I one scu, Formule de cubalură in care domentul de integrare esle un triunghi ourecare, Stud. și Cerc. Ştiinţ., Filiala Cluj, VI (1955), 7.
E. Moldovan, O generalizare a notiunii de convexitate, Stud. și Cerc. Ştiin}. Filiala Cituj VI (1955), 65.
М. И. Морозов, О некоторых вопросах равномерного приближения непрерывных функ. ций посредством функций интерполяционных классов, Изв. АН СССР 16 (1952), 75
T. Popoviciu, Asupra formei restului in tunele formule de aproximare ale analizei. Lucrările Sesiunii generale știintifice din 2-12 iunie 1950, Ed. Acad. R.P.R., Bucureşti 1951, p. 183.
Consideratii teoretice asupra udilizării practice a unor formule de interpolare, Bul. știinţ. Acad. R.P.R., Secţ. Şt. Mat. l'iz, XIJ, 2 (1951), 441.
Folytonos fuggvények középértéklételeiröl, A Magyar Tud. Akad. II Oszt. Közl. IV (1954), 353.
Despre precizia calculului numeric în interpolarea prin polinomul lui Newlon pe noduri echidistante, Stud. și Cerc. Ştiint. Filiala Cluj VI (1956), 20.
Despre precizia calculului numeric în interpolarea prin polinoame, Bul. Stiinl.. Acad. R.P.R., Secț. Şt. Mat. Fiz., VII, 4 (1955), 953.
J. Radon, Restausdrücke bei Interpolation und Quadraturformeln durch bestimmle Integralen, Monatshefte lür Math. u. Physik 42 (1935), 389.
L. Tornheim, On n-parameter Families of Functions and Associated Convex Frunctions, Trans. Amer. Math. Soc. 69 (1950), 457.