METODE NUMERICE DE REZOLVARE A ECUATIILOR DIFERENTIALE

1. CURS 6 Introducere METODE NUMERICE DE REZOLVARE A ECUATIILOR DIFERENTIALE O ecuatie diferentiala poate sa nu posede solutie sau, chiar daca are solutie, nu totdeauna aceasta se poate explicita. In multe situatii, mai ales in cazul ecuatiilor diferentiale neliniare, trebuie sa ne consideram multumiti daca obtinem o aproximatie a solutiei ecuatiei diferentiale. In cele ce urmeaza, utilizand metode numerice, se vor obtine seturi de puncte care, atunci cand se utilizeaza discretizari suficient de fine se poate aproxima solutia ecuatiei diferentiale considerate, asa cum se prezinta in fig. 1.

2. d2x dw d2w dp dx dt2 dt dx2 EI = - M(x) si  = dy = f(x,y) dx ECUATII DIFERENTIALE DE ORDINUL 1 Forma generala a acestei ecuatii este: dat fiind faptul ca ecuatiile diferentiale de ordin superior pot fi reduse la sisteme de ecuatii diferentiale de ordinul 1. EXEMPLE 1) Ecuatia de miscare a unei particule de masa m într-un câmp de forte F(x): m = F(x)  = F(x) 2) Ecuatia diferentiala a fibrei medii deformate: Intr-adevar, considerand ecuatia diferentiala: si introducand variabilele:

3. rezulta ca: In acest fel ecuatia diferentiala de ordinul n scrisa sub forma (1) se poate echivala cu urmatorul sistem de n ecuatii diferentiale de ordinul 1: ,

4. a) METODAEULER (METODA LINIILOR POLIGONALE) Aceasta este una dintre cele mai simple tehnici de solu\ionare numerica a ecuatiilor diferentiale, fiind cunoscuta si sub denumirea de metoda tangentelor. Se considera ecuatia diferentiala (inclusiv conditia initiala aferenta): Daca se considera pe axa Ox o diviziune echidistanta de pas h se poate gasi un punct (x1, y1) = (x0+h, y1) pe tangenta la curba ce reprezinta solutia ecuatiei diferentiale in punctul (x0, y0), asa cum se prezinta in figura alaturata. Se poate scrie: unde:

5. . Daca se noteaza x0+h = x1, atunci punctul de coordonate (x1, y1) situat pe tangenta considerata reprezinta a aproximatie a punctului de coordonate (x1 , y(x1)) situat pe curba ce reprezinta solutia ecuatiei diferentiale. Este evident faptul ca eroarea metodei este cu atat mai redusa cu cat valoarea pasului diviziunii considerate - h este mai mic. Se poate construi astfel un proces iterativ: a carui reprezentare grafica se prezinta in figura alaturata.

6. Aproximarea solutiei ecuatiei diferentale

7. APLICATIA 1 Considerand ecuatia diferentiala: si conditia initiala y(1)=1, sa se obtina o aproximatie pentru a gasi valoarea lui y in punctul de abscisa x = 1.5, utilizand un pas h cu valorile 0.1, 0.05 si 0.01, apoi sa se calculeze si erorile relative raportate la valoarea exacta, cu patru zecimale. Solutie Solutia analitica a ecuatiei diferentiale considerate este:

8. Tab. 2 Valorile obtinute pentru h = 0.05 Tab. 1 Valorile obtinute pentru h = 0.10 xn yn Exact Eroare Eroare % xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2155 1.5044 1.8955 2.4311 3.1733 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0182 0.0483 0.0982 0.1806 0.3171 0.00 1.47 3.11 4.93 7.98 9.08 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2000 1.4640 1.8154 2.2874 2.9278 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0337 0.0887 0.1784 0.3244 0.5625 0.00 2.73 5.71 8.95 12.42 16.12

9. Tab. 3 Valorile obtinute pentru h = 0.01 xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2298 1.5423 1.9723 2.5719 3.4197 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0039 0.0104 0.0214 0.0398 0.0707 0.00 0.31 0.67 1.07 1.52 2.03

10. b) METODA EULER-HEUN Aceasta metoda este cunoscuta si sub denumirea de metoda Euler imbunatatita, situatia fiind prezentata in figura de mai jos pentru primul pas de integrare. Procesul iterativ pentru solutionarea numerica este descris de relatia urmatoare:

11. Analizand figura se poate constata ca valoarea lui y1 este mai apropiata de solutie decat valoarea obtinuta pentru y1*, ceea ce poate incadra aceasta metoda in clasa metodelor predictor-corector. Intr-adevar, aceasta situatie se poate interpreta in felul urmator: Pasul predictor: Pasul corector: In continuare, cu referire la ecuatia diferentiala considerata anterior, se prezinta rezultatele obtinute prin aplicarea acestei metode pentru diferiti pasi de integrare si programul folosit pentru solutionarea numerica.

12. Tab. 4 Valorile obtinute pentru h = 0.1 Tab.5 Valorile obtinute pentru h=0.05 xn yn Exact Eroare Eroare % xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2332 1.5514 1.9909 2.6060 3.4795 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0004 0.0013 0.0029 0.0057 0.0108 0.00 0.04 0.08 0.14 0.22 0.31 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2320 1.5479 1.9832 2.5908 3.4509 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0017 0.0048 0.0106 0.0209 0.0394 0.00 0.14 0.31 0.53 0.80 1.13

13. Tab.6 Valorile obtinute pentru h = 0.01 xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2337 1.5527 1.9936 2.6115 3.4899 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0000 0.0000 0.0001 0.0002 0.0005 0.00 0.00 0.00 0.005 0.010 0.025

14. c) METODA DEZVOLTARII IN SERIE TAYLOR Aceasta metoda nu prezinta o larga aplicabilitate deoarece, in principal, rezultatele sunt comparabile cu cele obtinute prin aplicarea metodei Euler-Heun in condi\iile unor calcule mai laborioase. Asa cum se stie, dezvoltarea unei functii y(x) in jurul unui punct x = a are expresia: Pentru a apropia expresia de mai sus de problema curenta, daca se considera cazul in care a = xn si x = xn+h, se obtine: Daca se presupune ca functia y(x) este solutia ecuatiei diferentiale:

15. Oprind numai doi termeni din se poate scrie: Se observa faptul ca daca in ecuatia de mai sus se inlocuieste y(xn+h) cu yn+1 si y(xn) cu yn se obtine formula ce caracterizeaza metoda Euler: Oprind trei termeni si facand inlocuirile antementionate se obtine urmatoarea formula: Metoda poate fi considerata drept cheia de bolta a majoritatii metodelor de rezolvare numerica a ecuatiilor diferentiale.

16. d) METODE TIP RUNGE-KUTTA Aceasta clasa de metode reprezinta una dintre cele mai folosite in abordarea numerica a ecuatiilor diferentiale, imbinand numarul relativ redus de operatii elementare cu acuratetea rezultatelor. Metoda Runge-Kutta de ordinul II consta in ga[sirea constantelor a, b astfel incat expresia: cu: sa se apropie de dezvoltarea in serie Taylor pentru cat mai multi termeni posibili. Meritul principal al acestei clase de metode rezida deci in aceea ca se apropie de acuratetea unei dezvoltari in serie Taylor fara insa a fi nevoie sa se calculeze si derivatele de ordin superior.

17. Se fac urmatoarele observatii: a) Daca: atunci relatia anterioara reprezinta dezvoltarea in serie Taylor, oprind primii trei termeni, a functiei considerate. b) Daca: atunci relatia anterioara se reduce la expresia metodei Euler-Heun. c) Se poate constata ca metoda Euler este de fapt o procedura Runge-Kutta de ordinul I. Folosind aceleasi consideratii se determina astfel algoritmul de calcul si in cazul metodei Runge-Kutta de ordinul IV (se opresc primii patru termeni din dezvoltarea in serie Taylor):

18. In continuare se prezinta rezultatele numerice obtinute in cazul aplicarii acestei metode pentru rezolvarea ecuatiei diferentiale considerate anterior si programul folosit in acest scop.

19. xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2337 1.5527 1.9937 2.6117 3.4903 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.00 0.00 0.00 0.00 0.00 0.0028 xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2337 1.5527 1.9937 2.6116 3.4902 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0000 0.0000 0.0000 0.0001 0.0002 0.00 0.00 0.00 0.00 0.0038 0.0076 Tab. 7 Valorile obtinute pentru h = 0.1 Tab. 8 Valorile obtinute pentru h = 0.05

20. Tab. 9 Valorile obtinute pentru h = 0.01 xn yn Exact Eroare Eroare % 1.00 1.10 1.20 1.30 1.40 1.50 1.0000 1.2337 1.5527 1.9937 2.6117 3.4903429 1.0000 1.2337 1.5527 1.9937 2.6117 3.4904 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000571 0.00 0.00 0.00 0.00 0.00 0.0016

21. e) METODA MILNE Aceasta este o metoda[ care lucreaza in doi pasi. • Pasul predictor: • Pasul corector:

22. Se remarca faptul ca valorile pentru y0, y1, y2 si y3 trebuie cunoscute pentru demararea procesului iterativ. Acestea se pot calcula cu oricare dintre metodele anterioare, recomandandu-se insa ca acest calcul initial si se faca cu un nivel inalt de acuratete. CONCLUZII Aplicarea uneia sau alteia dintre aceste metode in vederea solutionarii unei probleme concrete implica o analiza atenta a erorilor induse, in vederea realizarii unui compromis acceptabil intre precizie si efortul de calcul cerut.