Numerické metody pro řídké matice
V této diplomové práci se zabýváme numerickými metodami řešení parciálních diferenciálních rovnic eliptického typu. Metody konečných diferencí a konečných prvků vytváří soustavy lineárních rovnic s velkou a řídkou maticí soustavy. Tyto soustavy jsou řešeny iteračními metodami, nejdříve se aplikují J...
Uloženo v:
| Hlavní autor: | |
|---|---|
| Další autoři: | |
| Typ dokumentu: | VŠ práce nebo rukopis |
| Jazyk: | Čeština |
| Vydáno: |
2014
|
| Témata: | |
| On-line přístup: | http://is.muni.cz/th/357617/prif_m/ |
| Shrnutí: | V této diplomové práci se zabýváme numerickými metodami řešení parciálních diferenciálních rovnic eliptického typu. Metody konečných diferencí a konečných prvků vytváří soustavy lineárních rovnic s velkou a řídkou maticí soustavy. Tyto soustavy jsou řešeny iteračními metodami, nejdříve se aplikují Jacobiho a Gaussova - Seidelova metoda, což je spojeno s určením spektrálních poloměrů iteračních matic obou metod pro zjištění, zda-li posloupnost aproximací řešení soustav bude konvergovat ke skutečnému řešení. V případě konvergence posloupnosti aproximací se zjistí, že tato konvergence je pomalá. Proto jsou použity Krylovovské metody a metoda sdružených gradientů, které generují posloupnost aproximací řešení, jejíž konvergence ke skutečnému řešení je rychlejší, a které jsou nenáročné pro paměť počítače. Ke všem metodám jsou popsány teoretické základy. In this diploma thesis we are concerned with numerical methods for solutions to partial elliptic differential equations. The finite difference method and the finite element method create systems of linear equations with a matrix, that is large and sparse. These systems are solved by iterative methods. At the beginning the Jacobi method and the Gauss - Seidel method are used, but their usage is associated with determination spectral radius of both iterative matrices. This is necessary because of convergence of sequence approximations of solution to true solution of systems. If sequence approximations is convergent, this convergence is too slow. That's why there are used the conjugate gradients method and Krylovov's methods, which generate sequence approximations of solution with higher speed of convergence. These methods are not expensive for computer memory. Theory of these methods is described in this thesis. |
|---|---|
| Popis jednotky: | Vedoucí práce: Jiří Zelinka |
| Fyzický popis: | 104 l. |