Home - Rasfoiesc.com
Educatie Sanatate Inginerie Business Familie Hobby Legal
Doar rabdarea si perseverenta in invatare aduce rezultate bune.stiinta, numere naturale, teoreme, multimi, calcule, ecuatii, sisteme




Biologie Chimie Didactica Fizica Geografie Informatica
Istorie Literatura Matematica Psihologie

Matematica


Index » educatie » Matematica
» Metode iterative de solutionare a sistemelor de ecuatii liniare


Metode iterative de solutionare a sistemelor de ecuatii liniare


Metode iterative de solutionare a sistemelor de ecuatii liniare

1. Aspecte generale

Determinarea solutiei exacte a sistemului liniar de forma [A] [X] = [b] cu ajutorul unor metode de tip iterativ este posibila numai dupa efectuarea unui numar, teoretic infinit, de iteratii. Deoarece nici o metoda practica nu poate cicla la infinit, rezulta ca metodele iterative determina doar o solutie aproximativa, aproximata prin trunchiere, care se abate mai mult sau mai putin fata de solutia exacta X , in functie de precizia dorita.

Metodele iterative apeleaza la construirea unui sir de aproximatii succesive X0, X1, , Xk care, in anumite conditii, tinde catre solutia exacta X . In cazul in care pentru sirul aproximatiilor succesive nu este posibila definirea unei limite, se spune ca metoda este divergenta. In cazul in care sirul aproximatiilor succesive are o limita, se spune ca metoda este convergenta. Convergenta sirului aproximatiilor succesive este echivalenta cu posibilitatea definirii unei relatii de recurenta intre doua aproximatii succesive Xk si Xk+1.



Definirea relatiei de recurenta intre aproximatiile succesive se face pornind de la desfacerea matricei [A] in doua matrice:

(1)

Utilizand aceasta desfacere in expresia sistemului initial de ecuatii se obtine relatia:

    (2)

care permite ca, pornind de la o aproximare initiala X0, sa se construiasca un sir de aproximatii succesive pe baza relatiei de recurenta:

    (3)

In general, toate metodele iterative de rezolvare a sistemelor de ecuatii liniare, definesc matricile [N] si [P] pornind de la desfacerea standard a matricei [A]:

    (4)

unde:

[L] este o matrice strict inferior triunghiulara ale carei elemente subdiagonale sunt elementele matricei [A];

[R] este o matrice strict superior triunghiulara ale carei elemente supradiagonale sunt elementele matricei [A];

[D] este o matrice diagonala ale carei elemente nenule sunt tocmai elementele diagonale din matricea [A].

Ca in cazul metodelor directe, toate metodele iterative folosesc impartirea la un element diagonal, numit pivot. Din acest motiv, desfacerea standard definita de relatia (4) trebuie sa se caracterizeze prin elemente nenule pe diagonala matricei [A]. Daca exista cel putin un asemenea element nul este necesara efectuarea pivotarii partiale.

Fiecare metoda iterativa poate sa aiba una dintre urmtoarele trei conditii de oprire a procesului iterativ:

  1. trecerea abaterii maxime intre doua iteratii succesive sub o valoare prag impusa Emax. Abaterea maxima se defineste ca cea mai mare abatere, in valoare absoluta, intre necunoscute in doua iteratii succesive:

    (5)

  1. trecerea reziduului maxim la o iteratie sub o valoare prag impusa Rmax. Reziduul maxim este definit ca cea mai mare abatere, in valoare absoluta, intre vectorul termenilor liberi [b] si produsul matricei [A] cu vectorul aproximatiei curente [Xk]:

    (6)

  1. trecerea valorii unei functii obiectiv care depinde de toate necunoscutele Xi (i = 1, , n) calculata la fiecare iteratie, sub o valoare prag impusa Fmax:

    (7)

2. Metoda iterativa Seidel - Gauss

Metoda iterativa Seidel - Gauss foloseste o desfacere de forma:

    (8)

Folosind aceasta desfacere, relatia de recurenta (3) capata urmatoarea forma:

    (9)

Relatia (9) particularizata pentru una dintre necunsocutele Xi, conduce la formula de iterare a metodei Seidel - Gauss:

    (10)

Aparitia termenilor aii la numitorul expresiilor (7) reflecta necesitatea ca toate elementele diagonale ale matricei [A] sa fie nenule. Conditia de convergenta a metodei Seidel - Gauss o reprezinta forma diagonal dominanta a matricei coeficientilor din matricea [A].

Algoritmul asociat metodei Seidel - Gauss este prezentat in Fig. 1.

Definirea sistemului de ecuatii liniare: rangul n, matricea [A] si vectorul [b];

Definirea parametrilor de iterare: abaterea relativa maxima admisa Emax si numarului maxim de iteratii Nmax;

Calculul iterativ:

3.1. Stabilirea aproximatiei initiale, identica cu termenii liberi:

[X] [b]

3.2. Initializarea procesului iterativ:

It

3.3. Daca s-a atins numarul maxim de iteratii (it = Nmax) sau abaterea DX a

trecut sub valoarea admisa (DX £ Emax), se incheie procesul iterativ si se

trece la pasul 4;

3.4. Trecerea la o noua iteratie:

it it +1

3.5. Calculul noii aproximatii si a abaterii relative maxime, folosind o variabila

tampon z:

z Xi;

3.6. Se revine la pasul 3.3.

Stabilirea conditiilor de iesire din bucla iterativa;

- Daca DX £ Emax (metoda converge), se afiseaza solutia aproximativa si

numarul de iteratii efectuate;

Daca DX > Emax, dar it = Nmax (metoda nu converge), se afiseaza mesajul

"Depasire numar maxim iteratii".

Mersul lucrarii

Se studiaza textul lucrarii si metoda de rezolvare numerica a sistemului de ecuatii liniare.

Se elaboreaza programul de calcul.

Se testeaza programul de calcul.

Se solutioneaza cateva aplicatii numerice semnificative, prezentandu-se rezultatele.





Politica de confidentialitate





Copyright © 2024 - Toate drepturile rezervate