Semestr zimowy 2008/09


Konsultacje:

Patrz: Plan. (pokój 5010 - IV piętro - wieża północna przy schodach). w czasie sesji czy ferii tylko po uprzednim kontakcie.

Numeryczna algebra liniowa

(wyklad czw 1215-1345) - sala 3130.
Egzamin II termin zakończony!
Wyniki:
USOS

Lista tematow na egzamin ustny.
Należy wybrać jeden temat i zgłosić się na egzamin. Poza wybranym tematem należy się ogólnie orientować w pozostałych zagadnieniach.

W praktycznych obliczeniach naukowych czy rozwiązujemy zadanie różniczkowe, czy optymalizacji, czy nieliniowe - praktycznie zawsze natkniemy sie na problem rozwiązania układu równań liniowych czy zagadnienia własnego często jako część algorytmu właściwego np wielowymiarowa metoda Newtona czy schematy zamknięte dla RRZw czy równania siatkowe dla RRcz - jeśli chcesz się dowiedzieć o algorytmach rozwiążywania takich zadań - ich własnościach - zaletach i wadach - ten wykład jest dla ciebie. Postaramy się w przystępny sposób opisać wyprowadzenie algorytmów i podać ich podstawowe własności - część wykładu będzie przeznaczona na omówienie problemów implementacyjnych np format macierzy rzadkich (sparse) oraz bibliotek zawierających zaimplementowane algorytmu omówione na wykładzie.

Program wykladu


Z uwagi na brak ćwiczeń prawdopodobnie ostatnie kilka wykładów się nie zmieści.
  1. Wstęp czym zajmuje sie numeryczna algebra liniowa - podstawowe zadania NAL, co będzie na wykładzie, literatura itp. Przypomnienie podstawowych pojęć: normy wektorowe i macierzowe.
  2. Rozkład P1AP2=LU czyli eliminacja Gaussa, wyprowadzenie algorytmu
  3. Dokończenie rozkładu LU arytmetyka fl głowne własności (przypomnienie), współczynnik uwarunkowania zadania, algorytm numerycznie poprawny, uwaga o własnościach numerycznych rozkładu LU z różnymi typami wyboru elementu głównego (same wyniki bez szczegółowych dowodów)
  4. Dokończenie własności fl dla rozkładu LU, algorytm Choleskiego jako szczeg. przypadek rozkładu LU

  5. LU dla macierzy pasmowych. Macierze rzadkie: różne formaty ich przechowywania, najprostsze operacje na macierzach rzadkich np mnożenie wektora przez macierz w odp. formacie (0,5-1 wykład).
  6. Metody iteracyjne rozwiązywania układów równań liniowych: klasyczne (stacjonarne) Gaussa-Seidla, Jacobi, Richardson, SOR.
  7. Zbieżność stacjonarnych metod iteracyjnych: Richardsan, Gaussa-Seidla w przypadku macierzy dodatnio określonej i Jakobiego i G-S w normie max.
  8. Metody wielomianowe - w szczególności metoda cykliczna Czebyszewa (2poziomowa) i 3-poziomowa metoda Czebyszewa. Zbieżność met. Czebyszewa dla macierzy dodatnio określonych.
  9. Dokończenie implementacji 3-poziomowej met. Czebyszewa. Metody gradientowe np najszybszego spadku, ich zbieżnośc,
  10. Metoda CG - definicja i wyprowadzenie wzorów.
  11. Zbieżność CG. PCG tzn CG z prekonditionerem lewostronnym i prawostronym. Implementacja i zbieżność.
  12. Dokończenie o prekonditioningu prawostronnym dla CG. GSMR - zbieżność i efektywna implementacja. (o ile starczy czasu mniej lub więcej dokładnie) - 1 wykład.
  13. Dokończnie implementacji GMRES i krótka informacja o zbieżności Numeryczna zagadnienie własne - sprowadzanie macierzy do macierzy podobnej postaci Hessenberga (dla mac. symetrycznej do trójdiagonalnej)
  14. Numeryczna zagadnienie własne Metoda potęgowa, odwrotna potęgowa, metoda QR klasyczna wersja i wersja praktyczna z przesunięciami. informacja o metodzie z ilorazem Rayleigha.
  15. Dokończenie metody QR - szkic dowodu zbieżności (klasycznej wersji bez shiftów), SVD + zastosowanie SVD do rozwiązania LZNK w przypadku ogólnym.

Literatura

Metody numeryczne algebry liniowej

Podręczniki - podstawy metod numerycznych w których podstawowe metody algebry liniowej też są omówione

Geometria z algebrą liniową

Dotyczy osób dopuszczonych do egz ustnego w II terminie: egz ustny w pokoju 5010 we wtorek 17 marca 2009 - 10-12; można też się zgłosić do dr Bechlera w tej sprawie.

ćwiczenia 2 grupy - wtorek 10-12 i 12-14 - sala 3130.
Do uzyskania 100punktów (2x40p za dwa kolokwia i 40p za zadania domowe - kilka serii).
Zaliczenie od 50p + obecności.
Zadania domowe będą zadawane seriami - termin oddania (zawsze w czasie ćwiczeń) zazwyczaj 3 tygodnie. (terminy są nieprzekraczalne)

Serie zadań domowych. (pliki pdf)

  1. Pierwsza seria - termin oddania : wtorek 4-11-2008
  2. Druga seria - termin oddania : wtorek 25-11-2008
    Doprecyzowałem treść niektórych zadań i dodałem dwa dodatkowe dla chętnych
  3. Trzecia seria - termin oddania : wtorek 16-12-2008
  4. Czwarta seria - termin oddania : wtorek 13-01-2009
  5. Piąta seria - dla chętnych
Wyniki Kolokiwum 2, ZD 4ta seria i zaliczenia ćwiczeń: plik pdf
Zadania z egzaminu + szkic rozwiązań: plik pdf
Strona do zajęć wykładowcy prof Leszka Plaskoty

Powrót do mojej strony domowej


Ostatnia aktualizacja: 16 marca 2009

Dziś jest