elemente de teoria erorilor

22
ELEMENTE DE TEORIA ERORILOR §1.1 Numere aproximative. Erori a) Sursele şi clasificarea erorilor . În rezolvarea numerică a unei probleme deosebim - în general - trei feluri de erori: 1) Erori inerente (iniţiale) . Aceste erori provin din: simplificarea modelului fizic pentru a putea fi descris printr-un model matematic; măsurători iniţiale; calcule anterioare; datele problemei; 2) Erori de metodă (de trunchiere) care se datorează preciziei insuficiente a metodei folosite. Majoritatea metodelor necesită un număr mare de operaţii aritmetice (adesea un număr infinit!) pentru a ajunge la soluţia exactă a problemei. Acest fapt necesită trunchieri (renunţarea la o mulţime - eventual infinită - de operaţii) şi aproximaţii ; 3) Erori de rotunjire (de calcul) care apar în datele de intrare, în calculele pe parcurs şi în datele de ieşire. Astfel, datele de start nu pot fi numere cu o infinitate de cifre deoarece nu putem opera cu acestea. Se înlătură multe cifre - eventual o infinitate - rotunjind numerele date, prin reţinerea unui număr finit de cifre. Pe parcursul calculelor numărul cifrelor creşte şi se impune de asemenea rotunjirea. Chiar rezultatul final se rotunjeşte reţinându - se un număr de cifre corespunzător preciziei dorite. Eroarea totală se compune din cele trei erori: Eroarea totală = erori inerente + erori de metodă + erori de rotunjire. Pentru o anumită metodă numerică de rezolvare a unei probleme, eroarea inerentă şi eroarea metodei constituie eroarea ireductibilă care nu poate fi influenţată de exactitatea efectuării calculelor. Exemplul 1.1.1 Să se calculeze volumul corpului format dintr-un cilindru circular drept şi un con circular drept cu înălţimi egale şi având baza comună, unghiul secţiunii axiale de la vârful conului fiind 2u, iar generatoarea a. Soluţie. Fie h înălţimea cilindrului, egală cu a conului şi fie R=OB raza bazei conului şi a cilindrului. Din triunghiul dreptunghic AOB ( AB =a, OA=h, OB=R) deducem R=a sin u, h=a cos u (cum se vede din figura1.1.1).Avem 5

Upload: laura-tarita

Post on 05-Dec-2014

89 views

Category:

Documents


11 download

DESCRIPTION

Elemente de Teoria Erorilor

TRANSCRIPT

Page 1: Elemente de Teoria Erorilor

ELEMENTE DE TEORIA ERORILOR

§1.1 Numere aproximative. Erori

a) Sursele şi clasificarea erorilor. În rezolvarea numerică a unei probleme deosebim - în general - trei feluri de erori:

1) Erori inerente (iniţiale) . Aceste erori provin din: simplificarea modelului fizic pentru a putea fi descris

printr-un model matematic; măsurători iniţiale; calcule anterioare; datele problemei;

2) Erori de metodă (de trunchiere) care se datorează preciziei insuficiente a metodei folosite. Majoritatea

metodelor necesită un număr mare de operaţii aritmetice (adesea un număr infinit!) pentru a ajunge la

soluţia exactă a problemei. Acest fapt necesită trunchieri (renunţarea la o mulţime - eventual infinită - de

operaţii) şi aproximaţii;

3) Erori de rotunjire (de calcul) care apar în datele de intrare, în calculele pe parcurs şi în datele de ieşire.

Astfel, datele de start nu pot fi numere cu o infinitate de cifre deoarece nu putem opera cu acestea. Se

înlătură multe cifre - eventual o infinitate - rotunjind numerele date, prin reţinerea unui număr finit de

cifre. Pe parcursul calculelor numărul cifrelor creşte şi se impune de asemenea rotunjirea. Chiar rezultatul

final se rotunjeşte reţinându - se un număr de cifre corespunzător preciziei dorite.

Eroarea totală se compune din cele trei erori:

Eroarea totală = erori inerente + erori de metodă + erori de rotunjire.

Pentru o anumită metodă numerică de rezolvare a unei probleme, eroarea inerentă şi eroarea metodei constituie

eroarea ireductibilă care nu poate fi influenţată de exactitatea efectuării calculelor.

Exemplul 1.1.1 Să se calculeze volumul corpului format dintr-un cilindru circular drept şi un con circular drept cu înălţimi egale şi având baza

comună, unghiul secţiunii axiale de la vârful conului fiind 2u, iar generatoarea a.

Soluţie.Fie h înălţimea cilindrului, egală cu a conului şi fie R=OB raza bazei conului şi a

cilindrului. Din triunghiul dreptunghic AOB ( AB =a, OA=h, OB=R) deducem R=a sin u,

h=a cos u (cum se vede din figura1.1.1).Avem

V=Vcorp=Vcil.+Vcon=

Aşadar, valoare volumului căutat este:

.

Considerăm datele iniţiale a şi u ca rezultat al unei măsurări. Fie şi rezultatele măsurării, a şi

u. Valoarea volumului cu aceste date se scrie

Eroare inerentă este .

Pentru calculul valorilor funcţiilor trigonometrice folosim sumele parţiale ale seriilor de puteri corespunzătoare şi obţinem

iar eroarea metodei este

Efectuând calculele prin rotunjiri şi înlocuind numărul iraţional printr-unul raţional, rezultă pentru volum valoarea . Eroarea de calcul

este .

Eroarea totală va fi: .

b) Numere aproximative. Fie mulţimea numerelor reale

5

Fig.1.1.1

Page 2: Elemente de Teoria Erorilor

R={(a1qn+ a2qn-1+ a3qn-2+...+ amqn-m+1+...)},

unde q este baza sistemului de numeraţie, ai cifrele cu 0aiq-1, qN,q2, iar m şi n sunt numere întregi, m1, care

ne arată ordinul de mărime al numerelor din R. Se ştie că orice număr real se poate scrie sub forma unui număr cu o

infinitate de cifre. Practic nu se poate lucra cu astfel de numere. Reducem numărul de cifre şi înlocuim mulţimea R

prin: ={(a1qn+ a2qn-1+ a3qn-2+...+ amqn-m+1)}, R, unde este o mulţime de numere reale cu m cifre.

Practic pentru a trece de la R la , procedăm astfel:

Fie mulţimea de numere reale

={[a1qn+ a2qn-1+ a3qn-2+...+ amqn-m+1]+qn-m(am+1+am+2q-1+am+3q-2+...)]},

Notăm = am+1+am+2q-1+am+3q-2+...

Dacă , mulţimii R facem să-i corespundă mulţimea

şi spunem că numerele din aproximează numerele corespunzătoare din prin lipsă (numere obţinute prin înlăturarea

cifrelor care urmează după am). Dacă mulţimii R facem să-i corespundă mulţimea

şi spunem că numerele din aproximează numerele din prin adaos (exces). Se observă că numerele din s-au

obţinut din numerele corespunzătoare din prin înlăturarea cifrelor care urmează după cifra am şi adunarea unităţii la

cifra am .

În cazul = q se foloseşte aproximarea prin lipsă sau exces, în funcţie de caracterul constructiv al

calculatorului.

În concluzie, mulţimea pe care se efectuează calculele este . Numerele din prin care înlocuim numerele

corespunzătoare din R se numesc numere aproximative. Prin înlocuirea lui R cu se face o eroare de calcul ( de

rotunjire ). Din cele spuse mai înainte rezultă că eroarea de rotunjire pentru obţinerea unui număr aproximativ nu

depăşeşte din unitatea ordinului ultimei cifre reţinute am (am 0).

Un număr aproximativ real şi pozitiv scris în baza 10 are forma

a=a110n+ a210n-1+ a310n-2+...+am10n-m+1,

unde 0 ai 9, , n Z, ai fiind cifrele numărului a în sistemul zecimal. Cifrele (semnele) unui număr

aproximativ a care se iau în considerare în calcul, se numesc cifre semnificative.

Definiţia 1.1.1. Se numesc cifre semnificative ale unui număr aproximativ a, cifrele nenule precum şi cele

nule care sunt situate între cifre nenule sau care indică ordinele păstrate în calcule (ordinele semnificative).

Exemplul1.1.2. La numărul a=5 -3+0 -4+0 -5+3 -6+0 -7=0,005003000 cifrele semnificative

sunt subliniate. Primele trei cifre sunt nesemnificative, ele indicând doar poziţia virgulei; ultimele două cifre nu sunt

semnificative deoarece a conţine cifre până la ordinul de mărime 10-7.

c) Erori.

Notăm cu a numărul aproximativ din a numărului exact A din R. Se mai spune că a reprezintă valoarea

aproximativă a numărului exact A şi se notează a A.

6

Page 3: Elemente de Teoria Erorilor

Dacă : a<A spunem că a este o aproximare prin lipsă a lui A;

a>A spunem că a este o aproximare prin exces a lui A.

Definiţia 1.1.2. Dacă a este o valoare aproximativă a numărului A , diferenţa =A-a se numeşte eroare a

numărului aproximativ a, iar se numeşte eroare absolută a numărului aproximativ a.

Cu definiţia modulului, putem scrie A-a= , A=a .

În calculul erorii întâlnim unul din cazurile :

Cazul 1. Se cunosc valorile a şi A. Eroarea se calculează direct, =A-a.

Cazul 2. Se cunoaşte numai aproximaţia a (nu se cunoaşte valoarea exactă A). În acest caz, eroarea nu se

poate calcula. Presupunem că pentru A cunoaştem un minorant m şi un majorant M, m<A<M. Cum a A, putem admite

că m<a<M. Prin adunarea membru cu membru a dublelor inegalităţi m<A<M, -M<-a<-m, se obţine

–(M-m)<A-a<M-m, adică <M-m sau <M-m.

Definiţia 1.1.3. Orice număr din intervalul G=[ , M-m) se numeşte eroare absolută la limită a

numărului aproximativ a.

Din această definiţie rezultă că , de fapt =inf G. Se observă de asemenea că orice număr

mai mare decât o eroare absolută la limită a unui număr a se poate numi şi el eroare absolută la limită a acelui număr.

Ne interesează practic o valoare cât mai mică pentru .

Exemplul 1.1.3. Fie A=e. Ştim că 2,71<e<2,72, m=2,71, M=2,72. Dacă a e=A, avem <M-m=2,72-2,71=0,01. Cu cele spuse

mai înainte, putem lua =0,01. Dacă ţinem seama de 2,718<e<2,72, avem <2,72-2,718=0,002, obţinem o estimare mai bună

=0,002 a erorii absolute la limită .

Definiţia 1.1.4. Se numeşte eroare relativă a numărului aproximativ a, raportul dintre eroarea absolută

şi modulul numărului A

= sau = ( când nu se cunoaşte A) (1.1.1)

de unde = sau = ( când nu se cunoaşte A) . (1.1.1I)

Definiţia 1.1.5. Se numeşte eroare relativă la limită a unui număr aproximativ a, un număr oarecare cu

proprietatea

. (1.1.2)

Din (1.1.1) şi (1.1.2) rezultă adică , de unde

= , = . (1.1.3)

Cum a A, în locul egalităţilor (1.1.3), când A este necunoscut, putem folosi egalităţile aproximative

= , = . (1.1.3I)

În anumite condiţii putem înlocui (1.1.3I) cu evaluări mai bune. Astfel, dacă a A, A>0 (evident şi a>0), iar <a,

atunci

7

Page 4: Elemente de Teoria Erorilor

= = ,

unde majorarea s-a realizat prin şi a - A a+ , astfel că

= . (1.1.4)

În aceleaşi ipoteze (A>0, a>0, <a) putem scrie

= A (a+ ) = ,

iar din egalitatea (a+ ) = rezultă prin calcul

= . (1.1.5)

Formula (1.1.4) ne dă eroarea relativă la limită a numărului aproximativ a când cunoaştem eroarea absolută la

limită şi numărul a. Formula (1.1.5) ne dă eroarea absolută la limită în funcţie de a şi .

Observaţia 1.1.1 Eroarea relativă caracterizează mai bine precizia calculelor, măsurătorilor, după cum se vede

din exemplul următor.

Exemplul 1.1.4 La măsurarea a două tije s-au obţinut lungimile aproximative t 1=142,7 cm, t2=14,5 cm cu erorile absolute la limită

= =0,1 cm.

Valorile exacte ale lungimilor celor două tije sunt T1=142,7 cm 0,1cm , T2=14,5 cm 0,1cm. Cu formula (1.1.4)obţinem :

; ;

Din faptul că cele două erori absolute la limită sunt egale s-ar putea trage concluzia că cele două măsurători au aceeaşi precizie. Falsitatea

acestei afirmaţii rezultă din faptul că erorile relative la limită sunt diferite; de aici se vede că măsurarea primei tije are precizie mai

mare.

Definiţia 1.1.6 Se spune că numărul aproximativ

a = a110n+ a210n-1+ a310n-2+...+am10n-m+1

are m cifre exacte (a1, a2, ..., am) dacă

10n-m+1 . (1.1.6)

Exemplul 1.1.5. Fie numărul exact A=257,876. Numărul aproximativ A=257,88 are cinci cifre exacte deoarece = =

= 0,004 < 102-5+1 = =0,005, unde n=2, m=5 şi s-a aplicat (1.1.6). Se observă că numărul a, care este o

aproximare a lui A la cinci cifre exacte, s-a obţinut prin rotunjire, a cincea cifră reţinută a fost mărită cu unitatea întrucât cifra înlăturată a fost 6>

10 ( rotunjire prin adaos).

Observaţia 1.1.2. Numărul a care verifică (1.1.6) se spune că aproximează numărul exact A cu m cifre

exacte în sens strict . În continuare vom folosi această definiţie pentru numărul de cifre exacte, subînţelegând “ în sens

strict ” fără a menţiona acest lucru.

În unele cazuri este comod să spunem că numărul a este o aproximare cu m cifre exacte în sens larg (slab) a

numărului A dacă se verifică

110n-m+1 (1.1.7)

8

Page 5: Elemente de Teoria Erorilor

(eroarea absolută nu depăşeşte unitatea rangului celei de a “m”-a cifră semnificativă).

Astfel, numărul a=52,67 este o aproximare la patru cifre exacte în sens larg a numărului A=52,6775,

deoarece = = =0,0075<101-4+1=10-2=0,01. Se vede că a=52,67 s-a obţinut din A prin

înlăturarea cifrelor care urmează după cifra a patra.

Numărul a=52,68 este o aproximare cu patru cifre exacte în sens srict a numărului A=52,6775, deoarece

= =0,0025< 101-4+1=0,005, (n=1, m=4). În acest caz a este o rotunjire (la 4 cifre exacte)

prin adaos a lui A.

Teorema 1.1.1. Dacă un număr aproximativ pozitiv

a=a110n + a210n-1+ a310n-2+...+ am10n-m+1

are m cifre exacte (în sens strict) atunci eroarea sa relativă verifică

101-m, (1.1.8)

unde 1 a1 9, a1 fiind prima cifră semnificativă a numărului a, m 1.

Demonstraţie Întrucât, prin ipoteză, a are m cifre exacte, se verifică (1.1.6) adică

= 10n-m+1

de unde, - 10n-m+1 A-a 10n-m+1 sau A-a - 10n-m+1

Aşadar,

A a- 10n-m+1 >a110n - 10n-m+1 = (2a1-101-m) = (a1+a1-101-m) a1, deoarece a1-101-m 0.

Ţinând seama de (1.1.1) precum şi de 10n-m+1, =A a1,

se obţine

= = 101-m

şi (1.1.8) este dovedită.

Observaţia 1.1.3 Cum membrul drept al inegalităţii (1.1.8) este mai mare decât a putem lua:

(1.1.9)

se arată că

(1.1.10)

Într-adevăr, pentru , termenul 101-m din inegalitatea obţinută în demonstraţia teoremei 1.1.1

se poate neglija astfel că

9

Page 6: Elemente de Teoria Erorilor

Apoi,

Observaţia 1.1.4 În inegalităţile (1.1.6) şi (1.1.8) şi pot fi înlocuite cu şi dacă numai

acestea sunt cunoscute.

Exemplul 1.1.6 Să aflăm numărul de cifre exacte pentru a=32,3856 ştiind că .

Cu definiţia 1.1.6, avem n=1, m necunoscut

( )

întrucât numărul a are 4 cifre exacte putem scrie a=32,39.

Ştiind că m=4 arătăm că se verifică (1.1.8). Într-adevăr luând

(A necunoscut)

obţinem inegalitatea

.

§1.2 Operaţii cu numere aproximative

a) Eroarea unei sume

Teorema1.2.1 Eroarea absolută a unei sume algebrice de numere aproximative nu depăşeşte suma erorilor

absolute ale acestor numere.

Demonstraţie Fie Xi valorile exacte, xi numerele aproximative şi erorile absolute, . Notând

şi ţinând seama că avem

care se scrie sub forma

sau

de unde cu , rezultă

(1.2.1)

şi teorema este demonstrată.

Cum din (1.2.1) rezultă

(1.2.2)

Egalitatea (1.2.2) ne spune că eroarea absolută la limită a sumei algebrice de numere aproximative este egală

cu suma erorilor absolute la limită ale acestor numere.

10

Page 7: Elemente de Teoria Erorilor

Teorema 1.2.2 Dacă termenii unei sume de numere aproximative au acelaşi semn, atunci eroarea relativă la

limită a sumei nu depăşeşte eroarea relativă la limită maximă a termenilor ei.

Demonstraţie Cu notaţiile din teorema 1.2.1, ţinând seama că Xi, , au acelaşi semn, iar

putem scrie succesiv, folosind şi (1.2.2),

astfel că

. (1.2.3)

Observaţia 1.2.1 Din (1.2.2) rezultă că eroarea absolută la limită a sumei nu poate fi sensibil inferioară erorii

absolute la limită a termenului cu eroarea absolută la limită maximă. Din această cauză ceilalţi termeni (cu eroarea

absolută la limită mai mică) nu aduc o contribuţie esenţială la precizia sumei. De aceea, pentru adunarea unor numere

cu erori absolute diferite: se lasă neschimbate numerele cu eroarea absolută maximă; se rotunjesc celelalte numere

păstrând 1-2 cifre în plus faţă de primele; se însumează numerele obţinute şi apoi se rotunjeşte rezultatul la numărul de

cifre ale numărului cu eroarea absolută la limită maximă.

Exemplul 1.2.1 Să se afle suma numerelor aproximative 256,4; 32,34; 0,084; - 0,0676; 0,00068 ştiind că cifrele lor sunt exacte.

Soluţie. Numerele date având cifrele exacte s-au obţinut prin rotunjire astfel că eroarea absolută la limită este, pentru fiecare număr, jumătate din

unitatea ordinului de mărime a ultimei cifre păstrate. Astfel, pentru primul număr, ultima cifră 4 are ordinul de mărime 10 -1.

Eroarea absolută la limită a primului număr este . Procedând analog erorile absolute la limită ale numerelor

de însumat sunt: 0,05; 0,005; 0,0005;0,00005;0,000005. Suma erorilor absolute la limită este 1=0,055555=0,0556. Rotunjim numerele păstrând o

cifră zecimală în plus în raport cu primul număr, care are eroarea absolută la limită maximă şi însumăm numerele obţinute

s=256,4+32,34+0,08-0,07+0,00=288,75.

Putem scrie s=288,7, eliminând cifra reţinută în plus faţă de primul termen. La rotunjirea termenilor sumei s-a făcut eroarea absolută

. Eroarea de rotunjire la rezultatul final este 3=0,05. Astfel

eroarea totală este =1+2+3=0,0556+0,0023+0,05=0,1079 sau . Atunci , adică . De fapt

S=288,75708, adică . (În calculul lui 2 s-a scris cu plus ce s-a omis de la termenul respectiv - rotunjire prin lipsă - şi cu

minus ce s-a adăugat la termenul respectiv - rotunjire prin adaos -).

b) Eroarea unui produs

Teorema1.2.3 Eroarea relativă a unui produs de numere aproximative diferite de zero nu depăşeşte suma

erorilor relative ale acestor numere.

Demonstraţie Fie p=x1 x2… xn produsul a n numere aproximative nenule. Cum produsul poate avea semnul +

sau - putem presupune că toate numerele sunt pozitive. Atunci lnp=lnx1+lnx2+…+lnxn. Eroarea

. Aşadar putem scrie egalitatea aproximativă:

(1.2.4)

Cu (1.2.4) avem , iar . Cu acestea, egalitatea dlnp=dlnx1+ dlnx2+…+dlnxn se scrie

11

Page 8: Elemente de Teoria Erorilor

şi cu relaţia (1.1.1) obţinem în final

. (1.2.5)

Dar astfel că . Aşadar,

(1.2.6)

Egalitatea (1.2.6) ne spune că eroarea relativă la limită a produsului este egală cu suma erorilor relative la limită a

factorilor.

Observaţia 1.2.2 Dacă se înmulţeşte un număr aproximativ x cu un număr exact k, eroarea relativă la limită a

produsului nu se schimbă, în timp ce eroarea absolută la limită este de ori mai mare.

Într-adevăr, dacă p=kx atunci cu (1.2.6) avem deoarece , k fiind număr exact. Pe

de altă parte

Observaţia 1.2.3 Din (1.2.6) se vede că eroarea relativă la limită a produsului nu poate fi semnificativ mai

mică decât eroarea relativă la limită a factorului cu cea mai mică precizie. De aceea, la efectuarea unui produs:

rotunjim factorii, lăsându-le 1-2 cifre în plus faţă de factorul cu cea mai mică precizie; după înmulţire păstrăm atâtea

cifre semnificative câte cifre exacte are factorul mai puţin precis.

Pentru a calcula numărul de cifre exacte al unui produs de factori (10) p=x1x2…x, vom presupune că

fiecare factor are cel puţin m cifre exacte (m2). Fie a1,a2,…a primele cifre semnificative a scrierii în baza 10 a celor

factori.

Factorii produsului se scriu:

(1.2.7)Cu formula (1.1.10) putem scrie

.

(1.2.8)

Scriem produsul sub forma Folosind (1.2.6) şi apoi (1.2.8) obţinem:

(1.2.9)

unde am făcut majorarea

(valoarea lui 10 este atinsă numai pentru ai=1, =10 ). (1.2.10)

Apoi,

unde am făcut majorările 1+110 şi (1.2.9). Aşadar,

(1.2.11)Ultima inegalitate şi (1.1.6) ne spune că, chiar în cazul extrem (ai=1, , =10), produsul are m-2 cifre exacte.

12

Page 9: Elemente de Teoria Erorilor

c) Eroarea unui cât

Teorema 1.2.4 Eroarea relativă a unui cât nu depăşeşte suma erorilor relative ale dempărţitului şi

împărţitorului.

Demonstraţie Fie câtul , x1, x2 numere aproximative pozitive. Avem lnq=lnx1-lnx2 , respectiv

dlnq= dlnx1-dlnx2. Folosind (1.2.4) avem egalitatea de unde, cu proprietatea modulului obţinem

inegalitatea

,

care cu (1.1.1) se scrie

(1.2.12)

Cum din (1.2.12) deducem şi astfel obţinem egalitatea

(1.2.13)

care ne arată că eroarea relativă la limită a câtului este egală cu suma erorilor relative la limită ale dempărţitului şi

împărţitorului.

Fie x1 şi x2 două numere aproximative pozitive, având fiecare cel puţin m cifre exacte (m2) şi fie a1 şi a2

primele lor cifre semnificative, iar . Cu (1.1.10) şi (1.2.13) putem scrie

unde . (1.2.14)

Apoi (1.1.3I) pentru ne dă, cu majorările şi (1.2.14),

(1.2.15)

Dacă:

10 a1 2 şi a 2 2 , avem , de unde cu (1.1.6) rezultă că q are cel puţin m-1

cifre exacte;

20 a1=1 sau a2=1, avem şi majorarea , iar , de unde cu (1.1.6) rezultă

că q are cel puţin m-2 cifre exacte.

d) Eroarea unei puteri

Fie u=xm, x>0, x1, mN, m2. Logaritmând lnu=mlnx şi luând dlnu=mdlnx cu relaţia (1.2.4) avem

13

Page 10: Elemente de Teoria Erorilor

, respectiv (1.2.16)

Aşadar, eroarea relativă (la limită) a puterii a "m"-a a unui număr aproximativ x este de m ori mai mare decât

eroarea relativă (la limită) a numărului x.

e) Eroarea unui radical.

Fie , x>0, x1, mN, m2ko. Prin ridicarea la puterea m obţinem x=um şi cu (1.2.16) avem de unde

, respectiv, . (1.2.17)

Aşadar, eroarea relativă (la limită) a rădăcinii de ordinul m a numărului aproximativ x este de m ori mai mică decât

eroarea relativă (la limită) a numărului x.

Observaţia 1.2.4 Proprietăţile privind erorile produsului, câtului, puterii, radicalului s-au dedus pentru erorile

relative. Evident, erorile absolute corespunzătoare se calculează cu formulele (1.1.1 I), (1.1.3I) sau (1.1.5).

Exemplul 1.2.2 Se dau numerele aproximative x1=12,3 şi x2=1,2564 având toate cifrele exacte. Să se calculeze p= x1 x2 şi şi să

se precizeze numărul de cifre exacte ale numerelor p şi q.

Soluţie. Numerele date având toate cifrele exacte, erorile absolute la limită au valori egale cu jumătate din unitatea ordinului de mărime al

ultimei cifre:

. Erorile relative la limită ale

celor două numere sunt (vezi(1.1.3I)):

Cu formulele (1.2.6) şi (1.2.13) avem .

Pentru calculul produsului p ţinem seama de observaţia 1.2.3 în baza căreia x 1 rămâne neschimbat (are precizia mai mică), iar pe x2 îl

rotunjim la 4 cifre, x2=1,256. Avem p=15,448 şi reţinem numărul de cifre al lui x1, p=15,4. Avem

. Utilizăm acum definiţia 1.1.6:

Aşadar, p are m=2 cifre exacte, cifrele 1 şi 5. A treia cifră este nesigură. În baza aceleiaşi observaţii produsul trebuie să aibă cel puţin o

cifră exactă ( cu două cifre mai puţin decât numărul x1).

Pentru q obţinem valoarea 9,79 dar este posibil ca q să aibă doar o cifră exactă deoarece ambele numere x 1 şi x2 au prima cifră

semnificativă egală cu 1 (vezi finele părţii c) a acestui paragraf). Pentru a preciza numărul de cifre exacte pentru q, calculăm şi aplicăm

definiţia 1.1.6:

.

Deci q are primele două cifre exacte, a treia cifră nefiind sigură.

§1.3 Formula generală a erorii

a) Problema directă a teoriei erorilor.

Calculul erorii se poate îngloba în problema determinării erorii unei funcţii diferenţiabile reale de mai multe

variabile reale, când se cunosc erorile acestor variabile.

Fie funcţia f :DR , DR n , diferenţiabilă pe D. Se presupun cunoscute erorile ale variabilelor

independente xi deci şi erorile absolute ale acestora.

14

Page 11: Elemente de Teoria Erorilor

Eroarea absolută a funcţiei f în punctul (x1,x2, …xn) este

adică

.

(1.3.1)

Membrul drept al inegalităţii (1.3.1) se majorează prin înlocuirea lui prin şi notând această

majorare cu obţinem eroarea absolută la limită a funcţiei f

.

(1.3.2) Împărţind ambii membri ai inegalităţii (1.3.1) prin avem

,

astfel că eroarea relativă a funcţiei f verifică

,

(1.3.3)

iar eroarea relativă la limită a funcţiei f este dată de egalitatea

.

(1.3.4)

Relaţiile (1.3.1)-(1.3.4) rezolvă problema directă a teoriei erorilor care constă în "determinarea erorii absolute

sau relative (la limită) într-un punct al unei funcţii reale de mai multe variabile reale când se cunosc erorile absolute (la

limită) ale variabilelor independente ale funcţiei în acel punct".

Observaţia 1.3.1. Cu formula erorii (1.3.1) putem evalua eroarea unei sume algebrice de numere

aproximative. Notând f=f(x1,x2,…,xn)= x1+x2+…+xn obţinem

ceea ce ştiam din §1.2 că eroarea absolută a unei sume de numere aproximative este majorată de suma erorilor absolute

ale acelor numere.

Observaţia 1.3.2 Pentru a deduce proprietatea privind eroarea relativă a unui produs de n numere

aproximative strict pozitive xi, , considerăm funcţia f=f(x1,x2,…,xn)= x1x2… xi…xn, logaritmăm

lnf= lnx1+lnx2+…+lnxI+… lnxn şi aplicăm (1.3.3)

.

Am obţinut astfel rezultatul cunoscut "eroarea relativă a produsului este majorată de suma erorilor relative ale factorilor

produsului".

b) Problema inversă a teoriei erorilor.

15

Page 12: Elemente de Teoria Erorilor

Această problemă se enunţă astfel: să se afle erorile absolute admisibile (adică la limită) ale variabilelor

independente ale unei funcţii reale astfel încât eroarea absolută la limită a acelei funcţii într-un punct dat să nu

depăşească o valoare dată.

Din punct de vedere matematic această problemă este nedeterminată deoarece eroarea absolută la limită dată

se poate obţine folosind o mulţime de combinaţii privind valorile erorilor absolute la limită . Rezultă că

problema pusă se poate rezolva numai dacă se impun condiţii suplimentare. Mai des întâlnite sunt următoarele condiţii

suplimentare:

10 Principiul egalităţii efectelor, prin care se presupune că cele n diferenţiale parţiale au aceeaşi

contribuţie la formarea erorii absolute la limită, adică

.

(1.3.5)

Folosind (1.3.2) avem

,

de unde

f=f(x1,x2,…,xn).

(1.3.6)

20 Erorile absolute la limită ale variabilelor xi au aceeaşi valoare

. (1.3.7)

Cu aceste condiţii suplimentare egalitatea (1.3.2) ne dă

unde

.

(1.3.8)

3 0 Erorile relative la limită ale variabilelor xi au toate aceeaşi valoare

(1.3.9)

Folosind formula de calcul (1.1.3I) a erorii relative la limită, egalităţile (1.3.9) se scriu

.

Cu formula (1.3.2) obţinem:

16

Page 13: Elemente de Teoria Erorilor

şi

.

(1.3.10)

Exemplul 1.3.1 Să se afle eroarea absolută la limită şi eroarea relativă la limită precum şi numărul de cifre exacte pentru

, ştiind că x1=1,04; x2=0,9807; x3=1,324056 au toate cifrele exacte.

Soluţie Întrucât xi, au toate cifrele exacte, erorile absolute la limită sunt:

Aplicăm formula (1.3.2) pentru calculul lui :

(1.3.11)

Găsim

Înlocuind în (1.3.11) şi efectuând calculele găsim . Cu formula (1.1.3I) avem

Pentru a preciza numărul de cifre exacte pentru f(x1,x2,x3)=1,893750609 folosim definiţia 1.1.6:

Exemplul 1.3.2 Să se afle erorile absolute la limită ale variabilelor independente care fac posibil calculul valorii funcţiei

cu patru cifre semnificative exacte, dacă x1=0,9542; x2=1,24506; x3=-0,873049 sunt numere

aproximative cu toate cifrele exacte.

Soluţie Avem

Apoi,

.

17

Page 14: Elemente de Teoria Erorilor

Pentru calculul valorilor folosim principiul egalităţii efectelor. Atunci, cu formula (1.3.6) pentru n=3 şi

, găsim

.

Determinăm acum numărul de cifre (semnificative) exacte pentru x i astfel încât să obţinem f cu 4 cifre exacte. Folosim în acest scop

definiţia 1.1.6:

.

Aşadar, luăm x1=0,954; x2=1,245; x1= -0,8730 şi obţinem :

cu 4 cifre exacte.

18