<p>This thesis consists of four papers and a summary. The common topic of the included papers are the pricing equations of path-dependent options. Various properties of barrier options and American options are studied, such as convexity of option prices, the size of the continuation region in American option pricing and pricing formulas for turbo warrants. In Paper I we study the effect of model misspecification on barrier option pricing. It turns out that, as in the case of ordinary European and American options, this is closely related to convexity properties of the option prices. We show that barrier option prices are convex under certain conditions on the contract function and on the relation between the risk-free rate of return and the dividend rate. In Paper II a new condition is given to ensure that the early exercise feature in American option pricing has a positive value. We give necessary and sufficient conditions for the American option price to coincide with the corresponding European option price in at least one diffusion model. In Paper III we study parabolic obstacle problems related to American option pricing and in particular the size of the non-coincidence set. The main result is that if the boundary of the set of points where the obstacle is a strict subsolution to the differential equation is C¹-Dini in space and Lipschitz in time, there is a positive distance, which is uniform in space, between the boundary of this set and the boundary of the non-coincidence set. In Paper IV we derive explicit pricing formulas for turbo warrants under the classical Black-Scholes assumptions.</p>

This thesis consists of four papers and a summary. The common topic of the included papers are the pricing equations of path-dependent options. Various properties of barrier options and American options are studied, such as convexity of option prices, the size of the continuation region in American option pricing and pricing formulas for turbo warrants. In Paper I we study the effect of model misspecification on barrier option pricing. It turns out that, as in the case of ordinary European and American options, this is closely related to convexity properties of the option prices. We show that barrier option prices are convex under certain conditions on the contract function and on the relation between the risk-free rate of return and the dividend rate. In Paper II a new condition is given to ensure that the early exercise feature in American option pricing has a positive value. We give necessary and sufficient conditions for the American option price to coincide with the corresponding European option price in at least one diffusion model. In Paper III we study parabolic obstacle problems related to American option pricing and in particular the size of the non-coincidence set. The main result is that if the boundary of the set of points where the obstacle is a strict subsolution to the differential equation is C1-Dini in space and Lipschitz in time, there is a positive distance, which is uniform in space, between the boundary of this set and the boundary of the non-coincidence set. In Paper IV we derive explicit pricing formulas for turbo warrants under the classical Black-Scholes assumptions.

Dans cette thèse, nous proposons des algorithmes et présentons des résultats numériques pour la résolution de jeux répétés stochastiques, à deux joueurs et somme nulle dont l'espace d'état est de grande taille. En particulier, nous considérons la classe de jeux en information complète et en horizon infini. Dans cette classe, nous distinguons d'une part le cas des jeux avec gain actualisé et d'autre part le cas des jeux avec gain moyen. Nos algorithmes, implémentés en C, sont principalement basés sur des algorithmes de type itérations sur les politiques et des méthodes multigrilles. Ces algorithmes sont appliqués soit à des équations de la programmation dynamique provenant de problèmes de jeux à deux joueurs à espace d'états fini, soit à des discrétisations d'équations de type Isaacs associées à des jeux stochastiques différentiels. Dans la première partie de cette thèse, nous proposons un algorithme qui combine l'algorithme des itérations sur les politiques pour les jeux avec gain actualisé à des méthodes de multigrilles algébriques utilisées pour la résolution des systèmes linéaires. Nous présentons des résultats numériques pour des équations d'Isaacs et des inéquations variationnelles. Nous présentons également un algorithme d'itérations sur les politiques avec raffinement de grilles dans le style de la méthode FMG. Des exemples sur des inéquations variationnelles montrent que cet algorithme améliore de façon non négligeable le temps de résolution de ces inéquations. Pour le cas des jeux avec gain moyen, nous proposons un algorithme d'itération sur les politiques pour les jeux à deux joueurs avec espaces d'états et d'actions finis, dans le cas général multichaine (c'est-à-dire sans hypothèse d'irréductibilité sur les chaînes de Markov associées aux stratégies des deux joueurs). Cet algorithme utilise une idée développée dans Cochet-Terrasson et Gaubert (2006). Cet algorithme est basé sur la notion de projecteur spectral non-linéaire d'opérateurs de la programmation dynamique de jeux à un joueur (lequel est monotone et convexe). Nous montrons que la suite des valeurs et valeurs relatives satisfont une propriété de monotonie lexicographique qui implique que l'algorithme termine en temps fini. Nous présentons des résultats numériques pour des jeux discrets provenant d'une variante des jeux de Richman et sur des problèmes de jeux de poursuite. Finalement, nous présentons de nouveaux algorithmes de multigrilles algébriques pour la résolution de systèmes linéaires singuliers particuliers. Ceux-ci apparaissent, par exemple, dans l'algorithme d'itérations sur les politiques pour les jeux stochastiques à deux joueurs et somme nulle avec gain moyen, décrit ci-dessus. Nous introduisons également une nouvelle méthode pour la recherche de mesures invariantes de chaînes de Markov irréductibles basée sur une approche de contrôle stochastique. Nous présentons un algorithme qui combine les itérations sur les politiques d'Howard et des itérations de multigrilles algébriques pour les systèmes linéaires singuliers.

Le travail de thèse est composé de trois thèmes principaux : le premier étudie l'existence et l'unicité pour des équations différentielles stochastiques (EDS) progressives-rétrogrades fortement couplées avec des opérateurs sous-différentiels dans les deux équations, dans l’équation progressive ainsi que l’équation rétrograde, et il discute également un nouveau type des inégalités variationnelles partielles paraboliques associées, avec deux opérateurs sous-différentiels, l’un agissant sur le domaine de l’état, l’autre sur le co-domaine. Le second thème est celui des EDS rétrogrades sans ainsi qu’avec opérateurs sous-différentiels, régies par un mouvement brownien fractionnaire avec paramètre de Hurst H> ½. Il étend de manière rigoureuse les résultats de Hu et Peng (SICON, 2009) aux inégalités variationnelles stochastiques rétrogrades. Enfin, le troisième thème met l’accent sur la caractérisation déterministe de la viabilité pour les EDS régies par un mouvement brownien fractionnaire. Ces trois thèmes de recherche mentionnés ci-dessus ont en commun d’étudier des EDS avec contraintes sur le processus d’état. Chacun des trois sujets est basé sur une publication et des manuscrits soumis pour publication, respectivement. / This PhD thesis is composed of three main topics: The first one studies the existence and the uniqueness for fully coupled forward-backward stochastic differential equations (SDEs) with subdifferential operators in both the forward and the backward equations, and it discusses also a new type of associated parabolic partial variational inequalities with two subdifferential operators, one acting over the state domain and the other over the co-domain. The second topic concerns the investigation of backward SDEs without as well as with subdifferential operator, both driven by a fractional Brownian motion with Hurst parameter H> 1/2. It extends in a rigorous manner the results of Hu and Peng (SICON, 2009) to backward stochastic variational inequalities. Finally, the third topic focuses on a deterministic characterisation of the viability for SDEs driven by a fractional Brownian motion. The three research topics mentioned above have in common to study SDEs with state constraints. The discussion of each of the three topics is based on a publication and on submitted manuscripts, respectively.

Nous étudions la convergence de systèmes dynamiques vers des équilibres. En particulier, nous nous intéressons à deux types d'équilibres. D'une part, les solutions d'inéquations variationnelles sous contraintes qui interviennent aussi dans la résolution de problèmes d'optimisation hiérarchique. D'autre part l'état stable d'un système dynamique, c'est à dire l'état où l'énergie du système est nulle. Cette thèse est divisée en deux parties principales, chacune focalisée sur la recherche d'un de ces équilibres. Dans la première partie nous étudions une classe d'algorithmes explicite-implicites pour résoudre certaines inéquations variationnelles sous contraintes. Nous introduisons un algorithme proximal-gradient pénalisé, "splitting forward-backward penalty scheme". Ensuite, nous prouvons sa convergence ergodique faible vers un équilibre dans le cas général d'un opérateur maximal monotone, et sa convergence forte vers l'unique équilibre si l'opérateur est de plus fortement monotone. Nous appliquons aussi notre algorithme pour résoudre des problèmes d'optimisation sous contrainte ou hiérarchique dont les fonctions objectif et de pénalisation sont formées d'une partie lisse et d'une autre non lisse. En effet, nous démontrons la convergence faible de l'algorithme vers un optimum hiérarchique lorsque l'opérateur est le sous-différentiel d'une fonction convexe semi-continue inférieurement et propre. Nous généralisons ainsi plusieurs algorithmes connus et nous retrouvons leurs résultats de convergence en affaiblissant les hypothèses utilisées dans nombre d'entre eux.Dans la deuxième partie, nous étudions l'action d'un contrôle interne local sur la stabilisation indirecte d'un système dynamique couplé formé de trois équations d'ondes, le système de Bresse. Sous la condition d'égalité des vitesses de propagation des ondes, nous montrons la stabilité exponentielle du système. En revanche, quand les vitesses sont différentes, nous prouvons sa stabilité polynomiale et nous établissons un nouveau taux de décroissance polynomial de l'énergie. Ceci étend des résultats présents dans la littérature au sens où le contrôle est localement distribué (et non pas appliqué à tout le domaine) et nous améliorons le taux de décroissance polynomial de l'énergie pour des conditions au bord de type Dirichlet et Dirichlet-Neumann. / We study the convergence of dynamical systems towards equilibria. In particular, we are interested in two types of equilibria. On one hand solutions of constrained variational inequations that are also involved in the resolution of hierarchical optimization problems. On the other hand the stable state of a dynamical system, i.e. the state when the energy of the system is zero. The thesis is divided into two parts, each focused on one of these equilibria. In the first part, we study a class of forward-backward algorithms for solving constrained variational inequalities. We consider a splitting forward-backward penalty scheme. We prove the weak ergodic convergence of the algorithm to an equilibrium for a general maximal monotone operator, and the strong convergence to the unique equilibrium if the operator is an addition strongly monotone. We also apply our algorithm for solving constrained or hierarchical optimization problems whose objective and penalization functions are formed of a smooth and a non-smooth part. In fact, we show the weak convergence to a hierarchical optimum when the operator is the subdifferential of a closed convex proper function. We then generalize several known algorithms and we find their convergence results by weakening assumptions used in a number of them. In the second part, we study the action of a locally internal dissipation law in the stabilization of a linear dynamical system coupling three wave equations, the Bresse system. Under the equal speed wave propagation condition we show that the system is exponentially stable. Otherwise, when the speeds are different, we prove the polynomial stability and establish a new polynomial energy decay rate. This extends results presented in the literature in the sense that the dissipation law is locally distributed (and not applied in the whole domain) and we improve the polynomial energy decay rate with both types of boundary conditions, Dirichlet and Dirichlet-Neumann.

Wir betrachten Optimalsteuerungsprobleme, die von partiellen Differentialgleichungen beziehungsweise Variationsungleichungen mit Differentialoperatoren erster Ordnung abhängen. Wir führen die Reformulierung eines Tagebauplanungsproblems, das auf stetigen Funktionen beruht, ein. Das Resultat ist ein Optimalsteuerungsproblem für Viskositätslösungen einer Eikonalgleichung. Die Existenz von Lösungen dieses und bestimmter Hilfsprobleme, die von semilinearen PDG‘s mit künstlicher Viskosität abhängen, wird bewiesen, Stationaritätsbedingungen hergeleitet und ein schwaches Konsistenzresultat für stationäre Punkte präsentiert. Des Weiteren betrachten wir Optimalsteuerungsprobleme, die von stationären Variationsungleichungen erster Art mit linearen Differentialoperatoren erster Ordnung abhängen. Wir diskutieren Lösbarkeit und Stationaritätskonzepte für diese Probleme. Für letzteres vergleichen wir Ergebnisse, die entweder durch die Anwendung von Penalisierungs- und Regularisierungsansätzen direkt auf Ebene von Differentialoperatoren erster Ordnung oder als Grenzwertprozess von Stationaritätssystemen für viskositätsregularisierte Optimalsteuerungsprobleme unter passenden Annahmen erhalten werden. Um die Konsistenz von ursprünglichem und regularisierten Problemen zu sichern, wird ein bekanntes Ergebnis für Lösungen von VU’s mit degeneriertem Differentialoperator erweitert. In beiden Fällen ist die erhaltene Stationarität schwächer als W-stationarität. Die theoretischen Ergebnisse werden anhand numerischer Beispiele verifiziert. Wir erweitern diese Ergebnisse auf Optimalsteuerungsprobleme bezüglich zeitabhängiger VU’s mit Differentialoperatoren erster Ordnung. Hierfür wird die Existenz von Lösungen bewiesen und erneut ein Stationaritätssystem mit Hilfe verschwindender Viskositäten unter bestimmten Beschränktheitsannahmen hergeleitet. Die erhaltenen Ergebnisse werden anhand von numerischen Beispielen verifiziert. / We consider problems of optimal control subject to partial differential equations and variational inequality problems with first order differential operators. We introduce a reformulation of an open pit mine planning problem that is based on continuous functions. The resulting formulation is a problem of optimal control subject to viscosity solutions of a partial differential equation of Eikonal Type. The existence of solutions to this problem and auxiliary problems of optimal control subject to regularized, semilinear PDE’s with artificial viscosity is proven. For the latter a first order optimality condition is established and a mild consistency result for the stationary points is proven. Further we study certain problems of optimal control subject to time-independent variational inequalities of the first kind with linear first order differential operators. We discuss solvability and stationarity concepts for such problems. In the latter case, we compare the results obtained by either utilizing penalization-regularization strategies directly on the first order level or considering the limit of systems for viscosity-regularized problems under suitable assumptions. To guarantee the consistency of the original and viscosity-regularized problems of optimal control, we extend known results for solutions to variational inequalities with degenerated differential operators. In both cases, the resulting stationarity concepts are weaker than W-stationarity. We validate the theoretical findings by numerical experiments for several examples. Finally, we extend the results from the time-independent to the case of problems of optimal control subject to VI’s with linear first order differential operators that are time-dependent. After establishing the existence of solutions to the problem of optimal control, a stationarity system is derived by a vanishing viscosity approach under certain boundedness assumptions and the theoretical findings are validated by numerical experiments.

Tikhonov-type regularization of linear and nonlinear ill-posed problems in abstract spaces under sparsity constraints gained relevant attention in the past years. Since under some weak assumptions all regularized solutions are sparse if the l1-norm is used as penalty term, the l1-regularization was studied by numerous authors although the non-reflexivity of the Banach space l1 and the fact that such penalty functional is not strictly convex lead to serious difficulties. We consider the case that the sparsity assumption is narrowly missed. This means that the solutions may have an infinite number of nonzero but fast decaying components. For that case we formulate and prove convergence rates results for the l1-regularization of nonlinear operator equations. In this context, we outline the situations of Hölder rates and of an exponential decay of the solution components.

Le contrôle par mode glissant est une technique d'automatique qui possède une longue histoire, la littérature remontant jusqu'au année 50. Son essence est la suivante : le contrôle est définit comme étant l'image d'une fonction discontinue de la variable de glissement, contraignant le système à évolué sur une variété, le système glisse alors dessus, d'où le nom. Cette variable de glissement est elle définie à partir de l'état du système. Les développements ont mené à la constitution d'une théorie bien établie à propos de cette technique, avec de nombreuses propriétés théoriques fort intéressante. Toutefois ceci ne porte que sur la version continue, c'est à dire quand le contrôle peut changer de valeur à chaque instant. En comparaison la version discrète du ce contrôleur est définie par le fait que la valeur du contrôle ne peut changer qu'à des instants isolés discrets. On a alors une fonction en escalier, constante sur la période d'échantillonnage. Cette situation est rencontrée par exemple lorsque le contrôleur est implémenté à l'aide d'un micro-contrôleur, ce qui est le cas dans nombre d'applications industrielles. Le principal problème avec le mode glissant est l'apparition d'un phénomène largement indésirable, le chattering (ou broutement) avec la version discrète du contrôleur, où même déjà en simulation. Dans ce dernier cas, nous appelons ceci du chattering numérique que nous attribuons à une mauvaise discrétisation du contrôle. L'approche développée ici se focalise sur ce point et est largement inspirée par les travaux effectués en mécanique non régulière, où ce type de comportement a aussi été observé lors de la simulation de système avec frottements et/où impacts. L'idée principale est de discretisé le contrôle de manière implicite et non explicite. Ceci permet d'éliminer le chattering numérique dans les cas simples (systèmes linéaires par exemple) où bien de le réduire grandement. Pour mener à bien l'analyse, des outils provenant de l'analyse convexe ainsi que des inégalités variationnelles en dimension finie sont utilisés. Le contrôleur proposé possède des propriétés intéressantes et proches de celles du temps continu. Ainsi on peut montrer que la variable de glissement est régie par une dynamique stable en temps finie, avec une fonction de Lyapunov. Le contrôle discret convergence vers celui du cas continu quand la période d'échantillonnage tends vers 0. Une atténuation d'éventuelles perturbations de type "matching" peut être établie. Ces travaux ont essentiellement portés sur le contrôle par mode glissant classique. L'algorithme dit twisting a pu être discrétisé avec la même technique et sa stabilité en temps finie grâce à une fonction de Lyapunov a pu être montrée. Ces propriétés ont été vérifiée en simulation, mais aussi de manière expérimentale. Ainsi des essais ont pu être menés sur deux banc d'essai: le premier est basé sur un système electropneumatique où à la fois le contrôle par mode glissant classique ainsi que le twisting ont pu être implémentés. L'objectif étant de suivre une trajectoire de référence. Le second système est un pendule inverse où le système doit être stabilisé à la position d'équilibre instable. Ici seul le contrôleur classique a été testé. L'analyse des données expérimentales a permis de mettre en lumière les performances supérieures des contrôleurs proposés par rapport à ceux classiquement usités. Les objectifs de contrôle sont mieux atteint et le chattering est grandement diminué. / Sliding Mode Control is a control technique with a long history, with research efforts dating back to the 50's. The basic idea is to define the control input as a discontinuous function of the sliding variable, which solely depends on the state, and to constraint the system to evolve on a manifold, hence the term sliding. Over the years a strong theory was build around this technique, but only in continuous time. In our context, this means that control input value can change value at any time. The discrete-time case is when the control input can only change at isolated time instants and the dynamical system on which the control is still a continuous-time process. The control input is therefore a step function. This case appears when the controller is digitally implemented, for instance with the help of a microcontroller. This kind of setup is nowadays ubiquitous in benchmarks and industrial applications. One of the main limitation of the applicability of sliding mode control is the chattering phenomenon that is witnessed when this control technique is applied in practice, but already in simulations. In contrast to previous approaches, we single out the chattering that is already witnessed in simulation, even with no disturbance and with perfect knowledge of the dynamics. This is called the numerical chattering and one of its distinct feature is the constant chattering, or high-frequency bang-bang behavior, of the control input. This naturally induces a chattering of the sliding variable. The claim that this type of chattering is usually predominant and that it is due to a bad discretization of the signum multifunction. The approach developed in this work was inspired by the research effort in the nonsmooth mechanical to properly simulate some systems like those with dry friction and/or unilateral constraints. The main point is to discretize the signum in an implicit fashion, that is its argument is the value of the sliding variable at the end of the next sampling period. With this change, the numerical chattering can be removed in the simplest cases, largely attenuated. The research effort was focused on classical sliding mode controller, rather than the higher order ones. The frameworks used to perform the analysis are convex analysis and variational inequalities. This discrete-time controller enjoys several interesting theoretical properties. First it is finite-time Lyapunov stable: the sliding variable goes to 0 in finite-time. The discrete-time control input converges to the continuous-time one as the sampling period goes to 0. The control action also attenuates the effect of matched perturbations. Also the increase of the gain of the controller does not affect the performances when the system is sliding. The twisting controller can be discretized in the same way and is also finite-time Lyapunov stable. This good theoretical properties have been verified in simulations, but also on experimental setups. Two tests were conducted: the first one on an electropneumatic system, where both the classical first-order sliding mode controller and the twisting algorithm were tested. The objective was to track a reference trajectory. The second one was an inverted pendulum on a cart with only the classical SMC. The goal was to stabilize the system at the unstable equilibrium. The analysis from the data collected during those experiments shows that the proposed controllers perform better than the their explicitly discretized versions. The performances are better and the chattering is effectively reduced.

Effiziente und zuverlässige a posteriori Fehlerabschätzungen sind eine Hauptzutat für die effiziente numerische Berechnung von Lösungen zu Variationsungleichungen durch die Finite-Elemente-Methode. Die vorliegende Arbeit untersucht zuverlässige und effiziente Fehlerabschätzungen für beliebige Finite-Elemente-Methoden und drei Variationsungleichungen, nämlich dem Hindernisproblem, dem Signorini Problem und dem Bingham Problem in zwei Raumdimensionen. Die Fehlerabschätzungen hängen vom zum Problem gehörenden Lagrange Multiplikator ab, der eine Verbindung zwischen der Variationsungleichung und dem zugehörigen linearen Problem darstellt. Effizienz und Zuverlässigkeit werden bezüglich eines totalen Fehlers gezeigt. Die Fehleranschätzungen fordern minimale Regularität. Die Approximation der exakten Lösung erfüllt die Dirichlet Randbedingungen und die Approximation des Lagrange Multiplikators ist nicht-positiv im Falle des Hindernis- und Signoriniproblems, und hat Betrag kleiner gleich 1 für das Bingham Problem. Dieses allgemeine Vorgehen ermöglicht das Einbinden nicht-exakter diskreter Lösungen, welche im Kontext dieser Ungleichungen auftreten. Aus dem Blickwinkel der Anwendungen ist Effizienz und Zuverlässigkeit im Bezug auf den Fehler der primalen Variablen in der Energienorm von großem Interesse. Solche Abschätzungen hängen von der Wahl eines effizienten diskreten Lagrange Multiplikators ab. Im Falle des Hindernis- und Signorini Problems werden postive Beispiele für drei Finite-Elemente Methoden, der konformen Courant Methode, der nicht-konformen Crouzeix-Raviart Methode und der gemischten Raviart-Thomas Methode niedrigster Ordnung hergeleitet. Partielle Resultate liegen im Fall des Bingham Problems vor. Numerischer Experimente heben die theoretischen Ergebnisse hervor und zeigen Effizienz und Zuverlässigkeit. Die numerischen Tests legen nahe, dass der aus den Abschätzungen resultierende adaptive Algorithmus mit optimaler Konvergenzrate konvergiert. / Efficient and reliable a posteriori error estimates are a key ingredient for the efficient numerical computation of solutions for variational inequalities by the finite element method. This thesis studies such reliable and efficient error estimates for arbitrary finite element methods and three representative variational inequalities, namely the obstacle problem, the Signorini problem, and the Bingham problem in two space dimensions. The error estimates rely on a problem connected Lagrange multiplier, which presents a connection between the variational inequality and the corresponding linear problem. Reliability and efficiency are shown with respect to some total error. Reliability and efficiency are shown under minimal regularity assumptions. The approximation to the exact solution satisfies the Dirichlet boundary conditions, and an approximation of the Lagrange multiplier is non-positive in the case of the obstacle and Signorini problem and has an absolute value smaller than 1 for the Bingham flow problem. These general assumptions allow for reliable and efficient a posteriori error analysis even in the presence of inexact solve, which naturally occurs in the context of variational inequalities. From the point of view of the applications, reliability and efficiency with respect to the error of the primal variable in the energy norm is of great interest. Such estimates depend on the efficient design of a discrete Lagrange multiplier. Affirmative examples of discrete Lagrange multipliers are presented for the obstacle and Signorini problem and three different first-order finite element methods, namely the conforming Courant, the non-conforming Crouzeix-Raviart, and the mixed Raviart-Thomas FEM. Partial results exist for the Bingham flow problem. Numerical experiments highlight the theoretical results, and show efficiency and reliability. The numerical tests suggest that the resulting adaptive algorithms converge with optimal convergence rates.

Zusammenfassung Diese Arbeit ist ein Versuch, verschiedene klassische und neuere Methodender glatten bzw. nichtglatten Optimierung zu verallgemeinern und in ihrem Zusammenhang darzustellen. Als Hauptinstrument erweist sich dabei die sogenannte verallgemeinerte Kojima-Funktion. Neben reichlichen Beispielen setzen wir einen besonderen Akzent auf die Betrachtung von Variationsungleichungen, Komplementaritaetsaufgaben und der Standartaufgabeder mathematischen Programmierung. Unter natuerlichen Voraussetzungen an diese Probleme kann man u.a. Barriere-, Straf- und SQP-Typ-Methoden, die auf Newton-Verfahrenbasieren, aber auch Modelle, die sogenannte NCP-Funktionen benutzen, mittelsspezieller Stoerungen der Kojima-Funktion exakt modellieren. Daneben werdendurch explizite und natuerliche Wahl der Stoerungsparameter auch neue Methoden dieser Arten vorgeschlagen. Die Vorteile solcher Modellierungsind ueberzeugend vor allem wegen der direkt moeglichen (auf Stabilitaetseigenschaften der Kojima-Gleichung beruhendenden)Loesungsabschaetzungen und weil die entsprechenden Nullstellen ziemlich einfach als Loesungen bekannter Ersatzprobleme interpretiert werden koennen. Ein weiterer Aspekt der Arbeit besteht in der genaueren Untersuchungder "nichtglatten Faelle". Hier wird die Theorie von verschiedenen verallgemeinerten Ableitungen und dadurch entstehenden verallgemeinerten Newton-Verfahren, die im Buch "Nonsmooth Equations in Optimization" von B. Kummer und D. Klatte vorgeschlagen und untersucht wurde, intensiv benutzt. Entscheidend ist dabei, dass die benutzten verallgemeinerten Ableitungen auch praktisch angewandt werden koennen, da man sie exakt ausrechnen kann. / This work attempts to generalize various classical and new methods of smooth or nonsmooth optimization and to show them in their interrelation. The main tool for doing this is the so-called generalized Kojima-function. In addition to numerous examples we specialy emphasize the consideration of variational inequalities, complementarity problems and the standard problem of mathematical programming. Under natural assumptions on these problems we can model e.g. barrier-, penalty-, and SQP-Type-methods basing on Newton methods, and also methods using the so-called NCP-function exactly by means of special perturbations of the Kojima-function. Furthermore, by the explicit and natural choice of the perturbation parameters new methods of these kinds are introduced. The benefit of such a modelling is obvious, first of all due to the direct solution estimation (basing on stability properties of the Kojima-equation) and because the corresponding zeros can easily be interpreted as solutions of known subproblems. A further aspect considered in this paper is the detailed investigation of "nonsmooth cases". The theory of various generalized derivatives and resulting generalized Newton methods, which is introduced and investigated in the book "Nonsmooth Equations in Optimization" of B. Kummer and D. Klatte, is intensely used here. The crucial point is the applicability of the used generalized derivatives in practice, since they can be calculated exactly.