Metody równań
Metoda Gaussa
- Przygotowanie macierzy: System równań jest przedstawiany w postaci macierzy, gdzie każdy wiersz odpowiada
równaniu, a kolumny reprezentują współczynniki przy zmiennych.
a11 |
a12 |
a1n |
a21 |
a22 |
a2n |
an1 |
an2 |
ann |
- Wybór wiodącego elementu: W pierwszym kroku należy wybrać wiodący element, czyli element wiodący w pierwszej
kolumnie. Może to być największa wartość bezwzględna lub inny wybrany kryterium.
- Eliminacja: W tej fazie dokonuje się eliminacji współczynników w kolumnach poniżej wiodącego elementu.
Pozostałe równania są redukowane poprzez odejmowanie odpowiednich wielokrotności pierwszego równania od
pozostałych równań, w celu wyeliminowania współczynników przy zmiennej.
- Powtarzanie kroków 2 i 3: Proces wyboru wiodącego elementu i eliminacji jest powtarzany dla pozostałych
kolumn, aż do uzyskania macierzy trójkątnej górnej.
a11 |
a12 |
a1n |
0 |
a22 |
a2n |
0 |
0 |
ann |
- Rozwiązanie systemu: Gdy macierz jest w postaci trójkątnej górnej, można łatwo obliczyć wartości nieznanych
zmiennych poprzez proces zwany "odwrotną substytucją". Zaczynając od ostatniego równania, wartości są
obliczane wstecz.
Metoda Gaussa (kroki)
Newtona, która polega na iteracyjnym przybliżaniu rozwiązania. (kroki):
- Ustalenie warunku początkowego iteracji x(0),
- Wyznaczenie wektora F(x(i)) w i-tym kroku,
- Macierz Jacobiego J(x(i)) w i-tym kroku,
- Rozwiązanie układu równań liniowych 𝐽(𝑥(i))∆𝑥(i) = −𝐹(𝑥(i)) w każdym
kroku ze względu na wektor ∆𝑥(i),
- Rozwiązanie wzorem 𝑥(i+1) = 𝑥(i) + ∆𝑥(i),
- Sprawdzenie kryterium stopu np. ‖𝐹(𝑥(i))‖ = max ‖𝑓k(𝑥(i))‖ < 𝜀 ,
- Jeżeli 6. Spełniona to otrzymamy wektor x(i) jest rozwiązaniem po i-iteracji dla zadanej
dokładności 𝜀,
jeśli nie to 𝑖 ∶= 𝑖 + 1 oraz przejęcie do punktu 2.
Ekstrapolacja
Jest techniką numeryczną wykorzystywaną do oszacowania wartości funkcji w punkcie, który jest poza zakresem
dostępnych danych. Pozwala na przewidywanie wartości funkcji na podstawie dostępnych danych numerycznych.
Kroki ekstrapolacji Richardsona
To metoda polegająca na ulepszaniu przybliżeń numerycznych poprzez zastosowanie różnych poziomów dokładności
obliczeń.
Oto ogólny opis kroków Ekstrapolacji Richardsona:
- Przygotowanie siatki wartości: Pierwszym krokiem jest utworzenie siatki wartości funkcji w danym obszarze,
które są dostępne do obliczeń. Im gęstsza siatka, tym lepsza jakość ekstrapolacji.
- Wybór odstępów: Następnie wybieramy odstępy między wartościami funkcji w siatce. Odstępy powinny być dobrze
zrównoważone, aby zapewnić równomierne rozłożenie punktów.
- Ekstrapolacja: Korzystając z dostępnych wartości funkcji, stosujemy formułę ekstrapolacyjną, która bierze
pod uwagę różne poziomy dokładności obliczeń. Możemy użyć wielomianu interpolacyjnego, np. wielomianu
Lagrange'a lub wielomianu Newtona, aby ekstrapolować wartość funkcji w punkcie poza dostępnymi danymi.
- Ulepszanie przybliżeń: Następnie można zastosować kroki ulepszania, które zwiększają dokładność przybliżeń.
Jest to możliwe poprzez obliczenie wartości funkcji przy wykorzystaniu większej liczby punktów lub
zastosowanie bardziej zaawansowanych technik interpolacji.
- Iteracyjne powtarzanie: Jeśli uzyskane przybliżenie jest wystarczająco dokładne, proces można zakończyć. W
przeciwnym razie można powtórzyć kroki Ekstrapolacji Richardsona, korzystając z bardziej precyzyjnych danych
lub innych metod numerycznych, aby dalej poprawiać przybliżenia.
Ekstrapolacja Richardsona jest użyteczna w przypadkach, gdy potrzebujemy estymacji wartości funkcji w punktach,
dla których brakuje dokładnych danych. Może być również stosowana do ulepszania przybliżeń numerycznych, co
prowadzi do bardziej precyzyjnych wyników obliczeń.
Całkowanie numeryczne
znane również jako metody numeryczne całkowania, to zbiór
technik wykorzystywanych do obliczania przybliżonych wartości całek oznaczonych funkcji.
Pozwala to na rozwiązanie problemów całkowych, gdy nie jest możliwe wykonanie dokładnych
obliczeń analitycznych.
-
Metoda prostokątów
Metoda ta polega na przybliżaniu obszaru pod krzywą za
pomocą prostokątów. W zależności od wyboru punktu, na którym jest wartość funkcji
(lewy, prawy lub środkowy punkt przedziału), metody prostokątów są nazywane
odpowiednio metodą prostokątów lewych, prawych lub środkowych. Całka jest
przybliżana poprzez sumowanie pól prostokątów.
-
Metoda trapezów
Metoda ta polega na przybliżeniu obszaru pod krzywą za pomocą
trapezów. Całka jest przybliżana przez sumowanie pól trapezów, które mają
wierzchołki na krzywej. Metoda trapezów może być również ulepszona poprzez
zastosowanie reguły trapezów złożonych, gdzie przedział całkowania jest podzielony
na wiele mniejszych podprzedziałów, a następnie stosuje się regułę trapezów na
każdym z podprzedziałów.
-
Metoda Simpsona
, znana również jako reguła Simpsona, jest bardziej zaawansowaną
metodą całkowania numerycznego, która oferuje wyższą dokładność przybliżenia całek
w porównaniu do prostszych metod, takich jak metoda prostokątów czy metoda
trapezów. Metoda ta polega na przybliżeniu obszaru pod krzywą za pomocą
wielomianów drugiego stopnia (parabol).
Ważne jest, aby pamiętać, że metody całkowania numerycznego są przybliżone, a wynik zależy
od wyboru kroku dyskretyzacji (np. podziału przedziału na mniejsze podprzedziały) i liczby
punktów użytych do przybliżania wartości funkcji. Im większa liczba podprzedziałów lub
punktów, tym dokładniejsze przybliżenie całki.
Obliczanie szeregów Fouriera
jest techniką matematyczną wykorzystywaną do
rozbicia okresowej funkcji na składowe sinusoidalne o różnych częstotliwościach. Szeregi
Fouriera są szeroko stosowane w dziedzinach takich jak analiza harmoniczna, przetwarzanie
sygnałów, teoria sterowania i wiele innych. Oto krótki opis kroków obliczania szeregów
Fouriera:
- Określenie okresu funkcji: Pierwszym krokiem jest ustalenie okresu funkcji, dla której
chcemy obliczyć szereg Fouriera. Okres to najmniejsza dodatnia wartość T, dla której funkcja powtarza się.
- Rozbicie funkcji na składowe: Następnie, na podstawie okresu funkcji, możemy rozłożyć ją na
składowe harmoniczne. Szereg Fouriera jest sumą nieskończonej liczby sinusów i kosinusów o różnych
częstotliwościach, których amplitudy i fazy zależą od funkcji pierwotnej.
- Obliczenie współczynników Fouriera: Kolejnym krokiem jest obliczenie współczynników
Fouriera, które określają amplitudy i fazy składowych harmonicznych. Istnieją różne wzory i metody
obliczania tych współczynników, takie jak wzory Eulera, wzory trigonometryczne lub wykorzystanie
transformaty Fouriera.
- Wykorzystanie składowych harmonicznych: Na podstawie obliczonych współczynników Fouriera,
można odtworzyć funkcję pierwotną jako sumę sinusów i kosinusów o odpowiednich częstotliwościach,
amplitudach i fazach. Im więcej składowych harmonicznych uwzględniamy w szeregu Fouriera, tym dokładniejsze
przybliżenie funkcji.
- Analiza wyników: Po obliczeniu szeregu Fouriera, można dokonać analizy wyników, takiej jak
określenie najważniejszych składowych harmonicznych, częstotliwości dominujących w funkcji, określenie
energii w poszczególnych składowych itp.
- Szybkie przekształcenie Fouriera według algorytmu Cooleya-Tukeya.
Układ równań różniczkowych
to zbiór równań, które opisują zależności między
pochodnymi wielu funkcji. Ogólnie, rozwiązanie układu równań różniczkowych polega na
znalezieniu funkcji, które spełniają wszystkie równania jednocześnie.
-
Metoda Eulera: : jest jedną z najprostszych metod numerycznych do rozwiązywania
układów równań różniczkowych. Polega ona na przybliżeniu rozwiązania przez
skończony krok czasowy. W tej metodzie, równania różniczkowe są zamieniane na
równania różnicowe, a pochodne są przybliżane za pomocą różnic skończonych.
Następnie, dla danego punktu początkowego, kolejne wartości funkcji są obliczane
iteracyjnie na podstawie wartości poprzednich.
-
Metoda Rungego-Kutty: jest bardziej zaawansowaną metodą numeryczną do
rozwiązywania układów równań różniczkowych. Jest ona oparta na rozwinięciu Taylora
funkcji w punkcie początkowym, aby obliczyć wartości funkcji w kolejnych krokach.
Istnieje wiele wariantów metody Rungego-Kutty, takich jak metoda RK2, RK4 itd., które
różnią się stopniem dokładności przybliżenia. Im wyższy stopień metody Rungego-Kutty, tym bardziej dokładne
przybliżenie można uzyskać, ale zwykle kosztem większej
złożoności obliczeniowej.
Podstawą wszystkich metod Rungego-Kutty jest wyrażenie zależności całkowej:
Należy pamiętać, że metody numeryczne do rozwiązywania układów równań różniczkowych
mają swoje ograniczenia i wymagają odpowiednich warunków początkowych oraz
określonego przedziału czasowego. Wybór odpowiedniej metody zależy od charakterystyki
układu równań i oczekiwanej dokładności rozwiązania.