Implementace semihladké Newtonovy metody

Abstract

Cílem mé bakalářské práce je implementace semihladké Newtonovy metody. Newtonova metoda je numerickou metodou, která se využívá při řešení soustavy rovnic. Tato metoda si ovšem neporadí s řešením rovnic, které obsahují nediferencovatelné funkce. Právě z tohoto důvodu si zavedeme semihladkou Newtonovu metodu, která je založena na využití Clarkeova kalkulu. V této práci postupně rozebírám Clarkeův kalkul, pod který spadá zobecněná derivace se zobecněným gradientem a zobecněný Jacobián. Dále popisuju "klasickou" Newtonovu metodu a hlavně semihladkou Newtonovu metodu. Poté popisuju samotnou implementaci semihladké Newtonovy metody v programu Matlab, dále srovnávám řešení Newtonovou metodou pomocí analyticky spočteného zobecněného Jacobiánu a numericky spočteného zobecněného Jacobiánu. Na závěr řeším praktickou úlohu, v které zkoumám deformaci struny při kontaktu s překážkou.

The goal of this work is an implementation of semismooth Newton method. Newton method is a numerical method, which is proposed to solve the system of equations. This method is not able to solve the system of equations with nodifferentable functions. Because of it we will define semismooth Newton metod which is based on Clarke calculus. I need some basic definitions from Clarke calculus like generalized derivative, generalized gradient and generalized Jacobian. Then I describe "classical" Newton method and mainly semismooth Newton method. The next part of my work is devoted to implementation of semismooth Newthon method in program Matlab and comparison of solutions based on generalized Jacobian and numerical approximation of generalized Jacobian. In the end of my work I show you some numerical experiments with deformation of string in contact with rigid obstacle.