Methods for enforcing non-negativity of solution in Krylov regularization
Metody vynucení nonnegativity řešení v krylovovské regularizaci
diploma thesis (DEFENDED)

View/ Open
Permanent link
http://hdl.handle.net/20.500.11956/127505Identifiers
Study Information System: 226592
Collections
- Kvalifikační práce [10594]
Author
Advisor
Referee
Pozza, Stefano
Faculty / Institute
Faculty of Mathematics and Physics
Discipline
Numerical and computational mathematics
Department
Department of Numerical Mathematics
Date of defense
23. 6. 2021
Publisher
Univerzita Karlova, Matematicko-fyzikální fakultaLanguage
English
Grade
Excellent
Keywords (Czech)
lineární inverzní problém, regularizace, iterační metody, krylovovský prostor, nonnegativitaKeywords (English)
linear inverse problem, regularization, iterative methods, Krylov subspace, non-negativityCílem této práce je studovat nezáporné inverzní úlohy a komplikace, které nastávají při jejich řešení pomocí standardních metod Krylovových podprostorů. Nejdříve shrneme teorii vázající se k nezáporným inverzním problémům. Poté se zaměříme na vybrané modifikace metod Krylovových podprostorů, které vedou k značným vylepšením řešení uvažovaných úloh. Popíšeme jejich vlastnosti, nastíníme jejich implementaci a navrhneme vylepšení pro jednu z metod. Dále provedeme numerické experimenty pro srovnání jednotlivých algoritmů, kde se zaměříme speciálně na analýzu vlivu výběru parametrů na chování řešičů. V práci je názorně ukázáno, že metody vynucující podmínku nezápornosti konstruují obecně lepší aproximaci neznámého přesného řešení. Navíc nově navržená metoda vede v některých případech k úspoře celkového výpočetního času při zachování dobré kvality aproximace.
The purpose of this thesis is to study how to overcome difficulties one typically encounters when solving non-negative inverse problems by standard Krylov subspace methods. We first give a theoretical background to the non-negative inverse problems. Then we concentrate on selected modifications of Krylov subspace methods known to improve the solution significantly. We describe their properties, provide their implementation and propose an improvement for one of them. After that, numerical experiments are presented giving a comparison of the methods and analyzing the influence of the present parameters on the behavior of the solvers. It is clearly demonstrated, that the methods imposing nonnegativity perform better than the unconstrained methods. Moreover, our improvement leads in some cases to a certain reduction of the number of iterations and consequently to savings of the computational time while preserving a good quality of the approximation.