1 / 25

Metody rozwi ązywania układów równań nie liniowych

Metody rozwi ązywania układów równań nie liniowych. Ogólne sformułowanie procesu iteracyjnego. Twierdzenie o punkcie stałym dla układów równań nieliniowych. Niech D Í D f będzie obszarem domkniętym i ograniczonym w R n i niech f ( x ) Î D dla wszystkich x Î D odzworowuje D na siebie.

eara
Download Presentation

Metody rozwi ązywania układów równań nie liniowych

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Metody rozwiązywania układów równań nieliniowych Ogólne sformułowanie procesu iteracyjnego

  2. Twierdzenie o punkcie stałym dla układów równań nieliniowych • Niech DÍDf będzie obszarem domkniętym i ograniczonym w Rn i niech f(x)ÎD dla wszystkich xÎD odzworowuje D na siebie. • Jeżeli istnieje taka stała 0 £L <1 i norma wektorowa ||·|| taka, że dla każdego x, x’ÎD jest spełniony warunek Lipschitza ||f(x)-f(x’)|| £L||x-x’||, wtedy: • Istnieje dokładnie jeden punkt stały x* w D. • Proces iteracyjny jest zbieżny do punktu x* dla każdego przybliżenia początkowego x(0). • Słuszne są następujące oszacowania błędów: oszacowanie aprioryczne oszacowanie aposterioryczne

  3. Normy wektorów n-norma 2-norma (norma euklidesowa) 1-norma norma maksimum

  4. Macierz Jacobiego w szacowaniu zbieżności

  5. Rząd zbieżności Proces iteracyjny ma rząd zbieżności r jeżeli istnieje taka stała 0 £ M £ ¥, że

  6. Metoda Newtona - wyprowadzenie

  7. Formuła iteracyjna wielowymiarowej metody Newtona

  8. Metoda Newtona • Przyjąć przybliżenie początkowe x(0) • W p-tej iteracji obliczyć macierz Jacobiego Jf(x(p)). • Rozwiązać układ równań na Dx(p). • Jf(x(p))Dx(p)=-f(x(p)). • 4. Obliczyć x(p+1)=x(p)+Dx(p). • 5. Jeżeli ||Dx(p+1)||<e lub ||f(x(p+1))||<d lub przekroczono największą dopuszczalną liczbę iteracji proces iteracyjny kończy się, w przeciwnym wypadku zwiększyć licznik iteracji o 1 i przejść do punktu 2. • Rząd zbieżności metody Newtona wynosi 2 jeżeli macierz Jacobiego w punkcie odpowiadającym rozwiązaniu układu jest nieosobliwa: • ||x(p+1)-x*||<M||x(p)-x*||2

  9. Tłumiona metoda Newtona W punkcie 4 podstawiamy x(p+1)=x(p)+2-kDx(p) gdzie k jest najmniejszą liczbą całkowitą taką, że ||f(x(p+1)||<||f(x(p))|| • Inne modyfikacje metody Newtona: • Dla oszczędności czasu można macierz Jacobiego liczyć nie w każdym kroku ale co parę kroków. • Pochodne można przybliżać symetrycznymi ilorazami różnicowymi

  10. Metoda siecznych dla układów równań nieliniowych Macierz Jacobiego przybliżamy macierzą ilorazów różnicowych liczonych pomiędzy iteracją p a p-1.

  11. Metoda największego spadku

  12. Metoda sprzężonych gradientów • Wybieramy punkt startowy x0. Obliczamy d0=-g0=-JTf(0). Jeżeli g0=0 procedura jest zakończona. • W kroku k=0,1,…,n-1 obliczamy Wstawiamy xk+1=xk+akdk Wstawiamy gk+1=gk+akJdk. Jeżeli ||gk+|||<e kończymy proces, w przeciwnym wypadku wstawiamy a następnie wstawiamy dk+1=-gk+1+bkdk

  13. Zastosowanie metody Newtona do obliczania stężeń równowagowych w układach wieloskładnikowych Przykład: równowagi ustalające się w amoniakalnym roztworze soli srebra 6 niewiadomych (stężeń równowagowych) 6 równań

  14. Algebraiczny zapis równowag chemicznych i bilansu masy.

  15. Obliczanie macierzy Q na podstawie macierzy A.

  16. Sposób pierwszy: rozwiązujemy układ równań na stężenia (Alcock, R.M.,Hartley, F.R. & Rogers, D.E. J. Chem. Soc. Dalton Trans.115, 1978)

  17. a jest dobierane tak aby nowe stężenia były dodatnie.

  18. Sposób 2: wyrażamy zmiany stężeń w wyniku reakcji przez postęp reakcji. Zależności po lewej można uzyskać w czysto formalny sposób. Ponieważ macierz stechiometryczna A jest ortogonalna do macierzy bilansu Q, wartość wyrażenia QC nie zmieni się, jeżeli dodamy do stężeń niej dowolną liniową kombinację wierszy macierzy A.

  19. jest dobierane tak aby nowe stężenia były dodatnie. Przybliżeniem początkowym jest jakikolwiek zestaw dodatnich stężeń spełniających równania bilansu masy.

  20. Przykład: równowagi w roztworze nad osadem siarczku żelaza (II).

  21. log (stezenie) H+ OH- Fe2+ FeOH+ Numer iteracji log(stezenie) S2- HS- H2S Numer iteracji

  22. Sposób 3: zmiany logarytmów stężeń wyrażamy poprzez parametry związane z równaniami bilansu masy Jeżeli macierz Q można podzielić na (n-m)x(n-m)-wymiarową część jednostkową odpowiadającą komponentom i pozostałą odpowiadającą kompleksom to parametry t mają sens zmian potencjałów chemicznych komponentów podzielonych przez RT.

  23. jest dobierane tak aby zmiany stężeń nie były niefizycznie duże. Przybliżeniem początkowym jest jakikolwiek zestaw stężeń spełniających równania na stałe równowag.

More Related