La commande à transmissions événementielles est une approche dans laquelle les instants de transmission sont définis selon un critère dépendant de l'état du système et non plus d'une horloge à l'instar des implantations périodiques. Dans cette thèse, nous nous concentrons sur la synthèse de telles lois de commande par retour de sortie. Les contributions sont les suivantes : (i) nous proposons une méthode de synthèse dite par émulation pour des systèmes non linéaires; (ii) nous présentons une méthode de synthèse jointe de la loi de commande et de la condition de déclenchement pour les systèmes linéaires; (iii) nous nous intéressons au cas de systèmes non linéaires singulièrement perturbés et nous construisons le contrôleur à partir d’approximation de la dynamique lente uniquement. / Event-triggered control is a sampling paradigm in which the sequence of transmission instants is determined based on the violation of a state-dependent criterion and not a time-driven clock. In this thesis, we deal with event-triggered output-based controllers to stabilize classes of nonlinear systems. The contributions of the presented material are threefold: (i) we stabilize a class of nonlinear systems by using an emulation-based approach; (ii) we develop a co-design procedure to simultaneously design the output feedback law and the event-triggering condition for linear systems; (iii) we propose stabilizing event-triggered controllers for nonlinear systems whose dynamics have two-time scales (in particular, we only rely on the knowledge of an approximate model of the slow dynamics)

Cette thèse traite le problème de l'estimation d'état, du diagnostic et de commande tolérante aux défauts des systèmes non linéaires représentés par un modèle de Takagi-Sugeno (T-S) à variables de prémisse non mesurables. De nombreux algorithmes pour la synthèse d'observateurs robustes vis-à-vis des perturbations, des imperfections de modélisation et des entrées inconnues sont présentés en se basant sur quatre types d'observateurs : les observateursproportionnels, les observateurs à entrées inconnues, les observateurs proportionnel intégral (PI) et multi-intégral (PMI). Par la suite, ces derniers sont utilisés pour le diagnostic de fautes affectant des systèmes non linéaires. Ceci est réalisé au moyen de trois stratégies. La première utilise l'observateur à entrée inconnue par découplage afin de rendre l'observateur insensible à certains défauts et permettre de détecter et d'isoler les défauts en construisant des bancs d'observateurs. En raison des conditions structurelles souvent insatisfaites, le découplage total des défauts de l'erreur d'estimation d'état n'est pas réalisable. Afin de s'affranchir de ces contraintes, la seconde stratégie utilise les observateurs PI et PMI pour estimer simultanément l'état et les défauts du système. La troisième stratégie qui utilise le formalisme H8 vise à concevoir un générateur de résidus minimisant l'influence des perturbations et maximisant l'influence des défauts. Un choix adéquat des paramètres du générateur de résidus permet la détection, la localisation et l'estimation des défauts. Enfin, une loi de commande tolérante aux défauts par poursuite de trajectoire d'un modèle de référence estproposée en exploitant les observateurs PI et PMI / This thesis deals with state estimation, fault diagnosis and fault tolerant control of nonlinear systems represented by a Takagi-Sugeno model with unmeasurable premise variables. The problem of state estimation of nonlinear systems with T-S model with unmeasurable premise variable is explored. Algorithms for robust observers synthesis with respect to perturbations, modeling uncertainties and unknown inputs are afterward presented. These algorithms are based on four kinds of observers called proportional, unknown input observers (UIOs), proportional-integral (PI) and multiple-integral (PMI) . The application on model-based diagnosis is studied based on three strategies. The first one uses unknown input observer to decouple some faults and makes the observers insensitive to certain faults. This allows to detect and isolate faults by constructing observers banks. Due to strong structural conditions on designing UIOs decoupling the faults on the state estimation error is not possible. To avoid this problem, the second strategy uses PI and PMI observers in order to estimate simultaneously the state and the faults of the system. The third strategy uses the H8 formalism. This aims to minimize the influence of perturbations and to maximize the effects of faults on the residual signal. An adequate choice of the residual generator parameters allows to detect, to isolate and to estimate the faults affecting the system. Lastly, a fault tolerant control law is proposed by reference trajectory tracking based on the use of PI and PMI observers

Les broncho-pneumopathies chroniques obstructives (BPCO) constituent un groupe de maladies des poumons caractérisées par le blocage du passage de l’air, rendant la respiration de plus en plus difficile. L’emphysème et la bronchite chronique sont les deux principales affections parmi les BPCO, mais les BPCO peuvent également être provoquées par les dégâts causés par des bronchites chroniques asthmatiques. L’emphysème pulmonaire est une maladie pulmonaire caractérisée par l’élargissement des espaces aériens distaux en amont des bronchioles terminales non respiratoires, accompagné de la destruction des parois alvéolaires. Ces modifications du parenchyme pulmonaire sont pathognomoniques de l’emphysème. La bronchite chronique est une forme de bronchite caractérisée par une production excessive d’expectoration, entraînant l’apparition d’une toux chronique et l’obstruction des voies respiratoires. Dans toutes ces affections, les dégâts causés aux voies respiratoires finissent par affecter les échanges gazeux dans les poumons. L’emphysème est généralement diagnostiqué de façon indirecte, sur la base d’un examen clinique, d’explorations de la fonction respiratoire (EFR), et d’une évaluation visuelle subjective des scanners des tomodensitogrammes. Ces tests présentent une valeur limitée dans les cas d’emphysème léger ou modéré. La présente étude aborde la possibilité d’appliquer une démarche d’analyse non linéaire à la répartition de la densité de l’air au sein de l’arbre des voies respiratoires des poumons à un quelconque niveau des ramifications. Les images sources de tomodensitométrie (TDM) du poumon sont traitées en deux phases, afin de produire un coefficient fractal de répartition de la densité de l’air. Au cours de la première phase, les valeurs brutes de pixel des images sources correspondant à toutes les densités d’air possibles sont traitées par un outil logiciel, mis au point pour construire une image cible. On y parvient par suppression en cascade des éléments indésirables (SCEI) : une étape de prétraitement dans l’analyse de l’image source. Celle-ci permet d’identifier les valeurs de densité d’air au sein de l’arbre des voies respiratoires, tout en éliminant toutes les valeurs non relatives à la densité de l’air. La seconde phase consiste en une réduction itérative de la résolution (RIR). Chaque réduction de la résolution produit un nouvel histogramme. Chaque histogramme ainsi produit comporte un certain nombre de pics, chacun d’entre eux correspondant à un ensemble de densités d’air. La courbe mettant en relation chaque réduction de la résolution avec le nombre de pics correspondant, obtenus à la résolution concernée, est tracée. Ceci permet de calculer la dimension fractale par une régression linéaire sur un graphique log – log. / Chronic Obstructive Pulmonary Disease (COPD) refers to a group of lung diseases that block airflow and make it increasingly difficult for you to breathe. Emphysema and chronic bronchitis are the two main conditions that make up COPD, but COPD can also refer to damage caused by chronic asthmatic bronchitis. Pulmonary emphysema is defined as a lung disease characterized by “abnormal enlargement of the air spaces distal to the terminal, non-respiratory bronchiole, accompanied by destructive changes of the alveolar walls”. These lung parenchymal changes are pathognomonic for emphysema. Chronic bronchitis is a form of bronchitis characterized by excess production of sputum leading to a chronic cough and obstruction of air flow. In all cases, damage to your airways eventually interferes with the exchange of oxygen and carbon dioxide in your lungs. Habitual techniques of emphysema’s diagnosis are based on indirect features, such as clinical examination; Pulmonary Function Tests (PFT) and subjective visual evaluation of CT scans. These tests are of limited value in assessing mild to moderate emphysema. The presented work discusses the possibility of applying a nonlinear analysis approach on air density distribution within lung airways tree at any level of branching. Computed Tomography (CT) source images of the lung are subjected to two phases of treatment in order to produce a fractal coefficient of the air density distribution. In the first phase, raw pixel values from source images, corresponding to all possible air densities, are processed by a software tool, developed in order to, construct a product image. This is done through Cascading Elimination of Unwanted Elements (CEUE): a preprocessing analysis step of the source image. It identifies values of air density within the airways tree, while eliminating all non-air-density values. Then, during the second phase, in an iterative manner, a process of Resolution Diminution Iterations (RDI) takes place. Every resolution reduction produces a new resultant histogram. A resultant histogram is composed of a number of peaks, each of which corresponding to a cluster of air densities. A curve is plotted for each resolution reduction versus the number of peaks counted at this particular resolution. It permits the calculation of the fractal dimension from the regression slope of log-log power law plot.

Orientador: Samuel da Silva / Resumo: Muitas aplicações da engenharia envolvem estruturas essencialmente não-lineares onde várias técnicas têm sido recentemente estudadas e investigadas por muitos pesquisadores. Dentre as várias abordagems, as que usam séries de Volterra têm apresentado propriedades úteis para fornecer um melhor entendimento para identificação e análise. Neste contexto, a presente tese propõem novas contribuições em como usar as séries de Volterra para caracterização, identificação e análise dinâmica de sistemas não-lineares usando sinais de entrada e saída e sinais somente de saída. Inicialmente, apresenta-se uma metodologia para análise de sistemas mecânicos não-lineares através das funções de resposta em frequência de alta-ordem (HOFRFs) e o conceito de HOFRFs estendidas com dados apenas de saída é introduzido e descrito em detalhes. Após isso, uma abordagem para identificação de sistemas não-lineares com base nas séries de Volterra através da expansão na base ortonormal de Kautz é proposta. Essa técnica permite identificar os seus núcleos mais facilmente e permite separar as contribuições dos termos lineares e não-lineares usando somente sinais de saída. Além disso, uma metodologia para análise modal de sistemas fracamente não-lineares sujeito a excitações com vários níveis de amplitude é também apresentada. A contribuição desse novo método reside no fato de que as HOFRFs são simplesmente estimadas como função das FRFs lineares. Basicamente, essa metodologia estende o conceito ... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Most recent engineering applications involve structures essentially nonlinear where several techniques have been recently studied and investigated by many researchers. Among them, the methods based on Volterra series expansion have presented powerful properties to provide a better understanding for identification and analysis. In this context, the present thesis proposes new contributions in how to use Volterra series for characterization, identification and dynamical analysis of nonlinear systems based on input and output signals and output-only signals. Initially, a methodology for analysis of nonlinear mechanical systems through higher-order frequency response functions (HOFRFs) is presented and the concept of extended HOFRFs based on output-only is introduced and described in detail. Afterwards, an approach for identification of nonlinear systems based on Volterra series through the expansion onto orthonormal Kautz basis is proposed. This technique allows to identify the Volterra kernels easily and enable to split the contribution of the linear and nonlinear terms using input-output as well as output-only signals. Furthermore, a methodology for modal analysis of weakly nonlinear systems under multilevel excitation is also proposed. The contribution of this new approach lies in the fact that HOFRFs are simply computed as functions of the linear FRFs. Basically, it extends the conventional experimental modal analysis methods in order to characterize and treat no... (Complete abstract click electronic access below) / Doutor

Submitted by Adalberto Zanatta Neder Lazarini null (adalberto.lazarini@yahoo.com) on 2017-02-02T20:05:36Z No. of bitstreams: 1 Dissertação Adalberto - Final.pdf: 1297717 bytes, checksum: 4cea64840d9e8f7f2841989ef7ea349a (MD5) / Approved for entry into archive by LUIZA DE MENEZES ROMANETTO (luizamenezes@reitoria.unesp.br) on 2017-02-06T17:12:01Z (GMT) No. of bitstreams: 1 lazarini_azn_me_ilha.pdf: 1297717 bytes, checksum: 4cea64840d9e8f7f2841989ef7ea349a (MD5) / Made available in DSpace on 2017-02-06T17:12:01Z (GMT). No. of bitstreams: 1 lazarini_azn_me_ilha.pdf: 1297717 bytes, checksum: 4cea64840d9e8f7f2841989ef7ea349a (MD5) Previous issue date: 2017-01-19 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Neste trabalho são feitas análises sobre quando as condições de existência dos teoremas apresentados por (GUEDES, 2015), que propõem condições necessárias e suficientes para a estabilidade de sistemas não lineares de tempo contínuo descritos através de modelos fuzzy Takagi-Sugeno (TS), que transformam sistemas não lineares em um conjunto convexo de sistemas lineares a partir de regras se-então, baseadas em Funções de Lyapunov Fuzzy (FLF), são satisfeitas. Primeiramente, são analisados os casos particulares de 3 e 4 modelos locais. Logo após, são tratados casos genéricos, com quantidades de regras tanto pares quanto ímpares, identificando quando as condições impostas pelos teoremas citados são satisfeitas. Para os casos nos quais as condições não são satisfeitas, é proposto também um algoritmo para obtenção do “pior caso” possível, varrendo todas as possibilidades que o sistema apresenta. Este algoritmo oferece condições necessárias e suficientes para o problema e também pode ser utilizado quando as condições exigidas pelos dois teoremas apresentados em (GUEDES, 2015) são satisfeitas. Por fim, são analisadas as contribuições do uso do algoritmo, onde são considerados parâmetros como número de LMIs, número de variáveis matriciais simétricas nxn simétricas utilizadas na resolução das LMIs e tempo computacionalmente necessário. / In this work,the required conditions for the theorems presented by(GUEDES,2015), which propose necessary and sufficient conditions for the stability of continuous-time nonlinear systems described by Takagi-Sugeno (TS) fuzzy models based on Lyapunov fuzzy functions (LFF), are analysed. First of all, the cases with 3 and 4 local models are considered. After that, generic cases are studied, with odd or even number os rules, identifying when the conditions imposed by these theorem shold. If the conditions are not satisfied, an algorithm is proposed as way of finding the “worst case scenario” for the system, considering all possible options. This algorithm offers necessary and sufficient conditions for solving this problem and also can be used when the conditions for the application of the theorems presented in(GUEDES,2015) hold. Finally,the algorithm’s contributions are analyzed, considering parameters such as number of LMIs,number of symmetric matricial variables used on solving the LMIs and computational time.

Estimar a região de estabilidade de um ponto de equilíbrio assintoticamente estável é importante em aplicações tais como sistemas de potência, economia e ecologia. A compreensão da estrutura qualitativa da fronteira da região de estabilidade é fundamental para estimar com eficiência a região de estabilidade. Caracterizações topológicas e dinâmicas da fronteira da região de estabilidade foram desenvolvidas ao longo das últimas décadas. Estas caracterizações foram desenvolvidas sob hipóteses de hiperbolicidade dos pontos de equilíbrio na fronteira e transversalidade. Para sistemas que dependem de parâmetros, a condição de hiperbolicidade pode ser violada em pontos de bifurcações. Estaremos interessados em estimar a região de estabilidade, para sistemas sujeitos a variações de parâmetros, onde ocorre a violação da condição de hiperbolicidade dos pontos de equilíbrio na fronteira da região de estabilidade devido ao aparecimento de uma bifurcação sela-nó do tipo zero nesta fronteira. Apresentaremos neste trabalho uma caracterização completa da fronteira da região de estabilidade na presença de um ponto de equilíbrio não hiperbólico sela-nó do tipo zero. Motivados também em oferecer um algoritmo conceitual para obter estimativas da região de estabilidade perturbada via conjunto de nível de uma dada função energia na vizinhança de um parâmetro de bifurcação sela-nó do tipo zero, buscaremos exibir resultados que permitam compreender o comportamento da região de estabilidade e de sua fronteira sob a influência das variações do parâmetro, incluindo variações do parâmetro próximo a um parâmetro de bifurcação sela-nó do tipo zero. / Estimating the stability region of an asymptotically stable equilibrium point is fundamental in applications such as power systems, economy and ecology. The knowledge of the qualitative structure of the stability boundary is essential to estimate with efficiency the stability region. Topological and dynamical characterizations of the stability boundary have been developed over the past decades. These characterizations were developed under assumptions of hyperbolicity of equilibrium points on the stability boundary and transversality. For systems that depend on parameters, the condition of hyperbolicity can be violated at points of bifurcations. We will be primarily interested in estimating the stability region, for systems subjected to parameter variations, when the condition of hyperbolicity of equilibrium points on the stability boundary is violated due to the appearance of a type-zero saddle-node bifurcation on the stability boundary. We will develop in this work, a complete characterization of the stability boundary in the presence of a type-zero saddle-node non-hyperbolic equilibrium point. Also, motivated to providing a conceptual algorithm to obtain estimates of the perturbed stability region via level sets of a given energy function in the neighborhood of a type-zero saddle-node bifurcation parameter, we offer results that explain the behavior of the stability region and its boundary under the influence of parameter variations, including variations of the parameter close to a type-zero saddle-node bifurcation parameter.

Industrial robot manipulators are general-purpose machines used for industrial automation in order to increase productivity, flexibility, and quality. Other reasons for using industrial robots are cost saving, and elimination of heavy and health-hazardous work. Robot motion control is a key competence for robot manufacturers, and the current development is focused on increasing the robot performance, reducing the robot cost, improving safety, and introducing new functionalities. Therefore, there is a need to continuously improve the models and control methods in order to fulfil all conflicting requirements, such as increased performance for a robot with lower weight, and thus lower mechanical stiffness and more complicated vibration modes. One reason for this development of the robot mechanical structure is of course cost-reduction, but other benefits are lower power consumption, improved dexterity, safety issues, and low environmental impact. This thesis deals with three different aspects of modeling and control of flexible, i.e., elastic, manipulators. For an accurate description of a modern industrial manipulator, the traditional flexible joint model, described in literature, is not sufficient. An improved model where the elasticity is described by a number of localized multidimensional spring-damper pairs is therefore proposed. This model is called the extended flexible joint model. This work describes identification, feedforward control, and feedback control, using this model. The proposed identification method is a frequency-domain non-linear gray-box method, which is evaluated by the identification of a modern six-axes robot manipulator. The identified model gives a good description of the global behavior of this robot. The inverse dynamics control problem is discussed, and a solution methodology is proposed. This methodology is based on a differential algebraic equation (DAE) formulation of the problem. Feedforward control of a two-axes manipulator is then studied using this DAE approach. Finally, a benchmark problem for robust feedback control of a single-axis extended flexible joint model is presented and some proposed solutions are analyzed.

This Thesis comprises of two parts containing similar studies of Nonlinear Localization induced Base Isolation of structural systems. The present method of base isolation,like other nonlinear vibration isolation methods, enjoys certain merits like capability of absorbing broad band vibrations, attenuating heavy shocks etc. The research in this thesis is an extension of this base isolation strategy first proposed by Vakakis and co-author. The strategy involves augmenting an appendage referred to as the secondary system with the main structural unit or the primary system, which we want to isolate from disturbances at the base. The primary system is coupled to the secondary system through a stiffness element. Both the primary and secondary systems have nonlinear dynamic behavior. It is seen that for certain choice of values of the coupling element, steady state vibration of very small magnitude is induced in the primary system. This result was established by considering a general discrete nonlinear system with viscous damping. Now it is a well known fact that viscous damping, though being widely used in literature as well as in practice doesn't turn out to be accurate enough to capture structural damping behaviors. Moreover, the actual damping mechanism if governed by some nonlinear function of the system variables, may influence the physics governing the nonlinear localization phenomenon in a manner rendering the present method not suitable for structural systems at the very outset. So in the present study we focus our attention in establishing the robustness and hence utility of the method by considering technically more defensible models of structural damping. These models efficiently capture certain complex phenomena which structures are known to exhibit. The occurrence of localization induced vibration isolation in structural systems in the presence of these damping models is taken as a proof of the efficacy of the method and its applicability to a wide range of situations. The present study establishes existence of localization through relevant analytical and numerical exercises. In the first part of the thesis we take up the study of nonlinear localization induced base isolation of a three degrees of freedom system having cubic nonlinearities under sinusoidal base excitation. The damping forces in the system are hysteretic in nature. In the present setting this is captured by Bouc-Wen model of hysteresis. Bouc-Wen model is one of the most widely used phenomenological model of hysteresis to have a ready-to-use mathematical description of hysteretic patterns appearing in structural engineering systems. The nature of responses of the different degrees of freedom as excitation frequency varies is a better way of analyzing the performance of the vibration isolation system. We adopt this line of approach for the present study. Normally Harmonic Balance Method (HBM) serves this purpose very well but in the present case as the hysteretic variable is not explicitly related to the system variables, HBM cannot be straightway implemented. Moreover, the hysteretic variable is related to other state variables through a relation which contains non-smooth terms. As a result, Incremental Harmonic Balance (IHB) method is used to obtain amplitude frequency relationship of the system response. The stability analysis of the solution branches is done by using Floquet Theory. Direct numerical simulation is then made use of to support our results that are obtained from this approximate numeric-analytic estimate of the amplitudefrequency relationships of the system, which helps us to analyze the efficacy of this method of base isolation for a broad class of systems. In the next part we consider a similar system where the damping forces in the system are described by functions of fractional derivative of the instantaneous displacements. Fractional Derivative based damping model has been found to be very effective in describing structural damping. We adopt half-order fractional derivative for our study, which can capture damping behavior of polymeric material very well. Typically linear and quadratic damping is considered separately as these are the two most relevant representations of structural damping. Under the assumption of smallness of certain system parameters and nonlinear terms an approximate estimate of the response at each degree of freedom of the system is obtained using Method of Multiple Scales. We then consider a situation where the nonlinear terms and certain other system parameters are no longer small. For the case where asymptotic methods are no longer valid, the assessment of performance of the vibration isolation system is made from amplitude-frequency relations. As a result, we take recourse to the Harmonic Balance Method in conjunction with arc length based continuation technique for obtaining the frequency amplitude plot for linear damping and Incremental Harmonic Balance method for quadratic damping, each of which is validated against results obtained from direct numerical simulation of the system. It needs to be appreciated that base isolation obtained this way has no counterpart in the linear theory.

<p>Industrial robot manipulators are general-purpose machines used for industrial automation in order to increase productivity, flexibility, and quality. Other reasons for using industrial robots are cost saving, and elimination of heavy and health-hazardous work. Robot motion control is a key competence for robot manufacturers, and the current development is focused on increasing the robot performance, reducing the robot cost, improving safety, and introducing new functionalities. Therefore, there is a need to continuously improve the models and control methods in order to fulfil all conflicting requirements, such as increased performance for a robot with lower weight, and thus lower mechanical stiffness and more complicated vibration modes. One reason for this development of the robot mechanical structure is of course cost-reduction, but other benefits are lower power consumption, improved dexterity, safety issues, and low environmental impact.</p><p>This thesis deals with three different aspects of modeling and control of flexible, i.e., elastic, manipulators. For an accurate description of a modern industrial manipulator, the traditional flexible joint model, described in literature, is not sufficient. An improved model where the elasticity is described by a number of localized multidimensional spring-damper pairs is therefore proposed. This model is called the extended flexible joint model. This work describes identification, feedforward control, and feedback control, using this model.</p><p>The proposed identification method is a frequency-domain non-linear gray-box method, which is evaluated by the identification of a modern six-axes robot manipulator. The identified model gives a good description of the global behavior of this robot.</p><p>The inverse dynamics control problem is discussed, and a solution methodology is proposed. This methodology is based on a differential algebraic equation (DAE) formulation of the problem. Feedforward control of a two-axes manipulator is then studied using this DAE approach.</p><p>Finally, a benchmark problem for robust feedback control of a single-axis extended flexible joint model is presented and some proposed solutions are analyzed.</p>

Για την αντιμετώπιση του προβλήματος ρύθμισης ενός συστήματος μη ελάχιστης φάσης είναι γνωστοί δύο τρόποι από τη θεωρία των γραμμικών συστημάτων. Ο ένας αφορά στην επιλογή βέλτιστης συνθετικής εξόδου ως προς την οποία το σύστημα είναι ελάχιστης φάσης. Ο δεύτερος τρόπος περιλαμβάνει άμεση κατασκευή βέλτιστου νόμου ανάδρασης καταστάσεων ως προς ένα σύνθετο δείκτη απόδοσης. Στην παρούσα εργασία αρχικά αναπτύσσεται μέθοδος για τη σύνθεση βέλτιστου νόμου ανάδρασης καταστάσεων για μη γραμμικές διεργασίες, όπου η είσοδος υπεισέρχεται μη γραμμικά στις διαφορικές εξισώσεις, με βάση ένα σύνθετο τετραγωνικό δείκτη απόδοσης. Ο δείκτης αυτός εξαρτάται τόσο από τη ρυθμιστική απόκλιση, όσο και από την απόκλιση της μεταβλητής χειρισμού. Για την επίλυση του προβλήματος δυναμικής βελτιστοποίησης χρησιμοποιούνται οι εξισώσεις Hamilton – Jacobi μέσω των οποίων υπολογίζεται ο βέλτιστος νόμος ανάδρασης καταστάσεων. Η λύση των εξισώσεων Hamilton – Jacobi υπολογίζεται με βάση την επαναληπτική μέθοδο Newton – Kantorovich. Σε κάθε βήμα της επανάληψης επιλύεται προσεγγιστικά μια μερική διαφορική εξίσωση τύπου Zubov με τη βοήθεια αναπτύγματος σε δυναμοσειρά. Στο Νοστό βήμα της επανάληψης η μέθοδος παράγει τη Νοστής τάξης προσέγγιση του αναπτύγματος κατά Taylor του βέλτιστου νόμου ανάδρασης καταστάσεων. Η παραπάνω μέθοδος εφαρμόζεται σε προβλήμα ρύθμισης της συγκέντρωσης προϊόντος σε σύστημα δύο μη ισοθερμοκρασιακών αντιδραστήρων CSTR, όπου λαμβάνει χώρα εξώθερμη αντίδραση, στην περίπτωση που η είσοδος υπεισέρχεται μη γραμμικά στις δυναμικές εξισώσεις της διεργασίας. Επίσης μελετώνται οι ιδιότητες σύγκλισης της επαναληπτικής μεθόδου Newton – Kantorovich, όταν αυτή εφαρμόζεται για την επίλυση της εξίσωσης Hamilton – Jacobi – Bellman που αντιστοιχεί στο πρόβλημα βελτιστοποίησης ενός σύνθετου τετραγωνικού δείκτη απόδοσης υπό τους περιορισμούς μιας μη γραμμικής δυναμικής όπου η είσοδος υπεισέρχεται γραμμικά στις διαφορικές εξισώσεις. Στη συνέχεια, για τη βέλτιστη ρύθμιση μη γραμμικών συστημάτων με ασταθή δυναμική μηδενιστών (συστήματα μη ελάχιστης φάσης), χρησιμοποιείται ο συνήθης τετραγωνικός δείκτης απόδοσης ISE. Στην περίπτωση αυτή το πρόβλημα δυναμικής βελτιστοποίησης είναι ιδιόμορφο. Για την επίλυση του προβλήματος αυτού το μη γραμμικό σύστημα μετασχηματίζεται στην κανονική μορφή Byrnes-Isidori, εφαρμόζεται η θεωρία Hamilton – Jacobi και υπολογίζεται στατικά ισοδύναμη συνθετική έξοδος με ευσταθή δυναμική μηδενιστών. Η ρύθμιση της συνθετικής εξόδου στο προκαθορισμένο σημείο επιτυγχάνεται με γραμμικοποίηση εισόδου/εξόδου. Για την επίλυση των σχετικών εξισώσεων Hamilton–Jacobi αναπτύσσεται η επαναληπτική μέθοδος Newton – Kantorovich, η οποία περιλαμβάνει την επίλυση μιας μερικής διαφορικής εξίσωσης τύπου Zubov σε κάθε βήμα της επανάληψης. Η μέθοδος εφαρμόζεται σε πρόβλημα ρύθμισης της συγκέντρωσης του επιθυμητού προϊόντος σε μη ισοθερμοκρασιακό αντιδραστήρα CSTR με κινητική Van de Vusse που παρουσιάζει ασταθή δυναμική μηδενιστών. Τέλος, οι δύο μέθοδοι συγκρίνονται με βάση τους επιμέρους δείκτες απόδοσης ISE και ISC, των οποίων ο γραμμικός συνδυασμός συνιστά το σύνθετο δείκτη απόδοσης της πρώτης μεθόδου, ενώ τα αποτελέσματά τους συγκρίνονται όταν αυτές εφαρμόζονται σε πρόβλημα ρύθμισης της συγκέντρωσης του επιθυμητού προϊόντος σε μη ισοθερμοκρασιακό αντιδραστήρα CSTR με κινητική Van de Vusse. / For the control of nonlinear nonminimum – phase systems, there are two possible lines of attack, originating from linear systems theory: a) direct calculation of the optimal state feedback with respect to a quadratic performance index that represents a combination of an error measure and a control effort measure (composite index), and b) calculation of the ISE-optimal minimum-phase output and subsequent input/output linearization on that output. This work develops a numerical algorithm for the calculation of an optimal nonlinear state feedback law for nonlinear systems. A quadratic performance index is used, which contains quadratic error terms and quadratic input penalty terms. The optimization problem is solved using the Hamilton-Jacobi equations, which determine the optimal nonlinear state feedback law. A Newton-Kantorovich iteration is developed for the solution of the pertinent Hamilton-Jacobi equations, which involves solving a Zubov partial differential equation at each step of the iteration, using a power series method. At step N of the iteration, the method generates the (N+1)-th order truncation of the Taylor series expansion of the optimal state feedback function. The method is applied to the problem of controlling a system of two non-isothermal continuous stirred tank reactors (CSTR), where an exothermic reaction takes place. Convergence properties of the algorithm are also developed independently of Kantorovich’s theorem, and the results are illustrated in a numerical example. For the optimal regulation of nonminimum-phase nonlinear systems, the performance index ISE (Integral of the Square of the Error) is used. The problem of minimizing ISE subject to the dynamics of the system and closed-loop stability is singular. The problem of calculation of an ISE-optimal, statically equivalent, minimum-phase output for nonminimum-phase compensation is formulated using Hamilton-Jacobi theory and the Byrnes-Isidori normal form representation of the nonlinear system. An input/output linearizing state feedback law is applied to regulate the synthetic output to a constant set point. A Newton-Kantorovich iteration is developed for the solution of the pertinent Hamilton-Jacobi equations, which involves solving a Zubov equation at each step of the iteration. The method is applied to the problem of controlling a nonisothermal CSTR with Van de Vusse kinetics, which exhibits nonminimum-phase behaviour. Finally, the two methods are compared with respect to the constituent indexes ISE and ISC (Integral of the Square of the Control), whose linear combination forms the composite performance index. The numerical results from both methods are compared in the control of a nonisothermal CSTR with Van de Vusse kinetics.