Spelling suggestions: "subject:"opérateur"" "subject:"opérateurs""
171 |
Composition flexible par planification automatiqueMartin, Cyrille 04 October 2012 (has links) (PDF)
Nous nous positionnons dans un contexte d'informatique ambiante dans lequel il arrive que les besoins de l'utilisateur n'aient pas été prévus, notamment en situation exceptionnelle. Dans ce cas, il peut ne pas exister de système préconçu qui réponde exactement à ces besoins. Pour les satisfaire, il faut alors pouvoir composer les systèmes disponibles dans l'environnement, et le système composé doit permettre à l'utilisateur de faire des choix à l'exécution. Ainsi, l'utilisateur a la possibilité d'adapter l'exécution de la composition à son contexte. Cela signifie que la composition intègre des structures de contrôle de l'exécution, destinées à l'utilisateur : la composition est dite flexible. Dans cette thèse, nous proposons de répondre au problème de la composition flexible en contexte d'intelligence ambiante avec un planificateur produisant des plans flexibles. Dans un premier temps, nous proposons une modélisation de la planification flexible. Pour cela, nous définissons les opérateurs de séquence et d'alternative, utilisés pour caractériser les plans flexibles. Nous définissons deux autres opérateurs au moyen de la séquence et de l'alternative : l'entrelacement et l'itération. Nous nous référons à ce cadre théorique pour délimiter la flexibilité traitée par notre planificateur Lambda-Graphplan. L'originalité de Lambda-Graphplan est de produire des itérations en s'appuyant sur une approche par graphe de planification. Nous montrons notamment que Lambda-Graphplan est très performant avec les domaines se prêtant à la construction de structures itératives.
|
172 |
Contrôle du partage de l'autorité dans un système d'agents hétérogènesMercier, Stéphane 18 February 2011 (has links) (PDF)
Dans le cadre de missions réalisées conjointement par un agent artificiel (robot, drone) et un agent humain, nous présentons un contrôleur de la dynamique de l'autorité, fondé sur un graphe de dépendances entre ressources génériques contrôlables par les deux agents, dont l'objectif est d'adapter le comportement de l'agent artificiel ou de l'agent humain en cas de conflit d'autorité sur ces ressources. Nous définissons l'autorité relative de deux agents par rapport au contrôle d'une ressource, ainsi que la notion de conflit d'autorité sur une ressource. La détection d'un conflit déclenche une révision du plan de l'agent artificiel, fondée sur l'état courant des ressources et des relations d'autorité. Le contrôleur de la dynamique de l'autorité est mis en œuvre par l'outil MNMS (Multiple (Petri) Nets Management System), développé au cours de la thèse, qui permet la manipulation des ressources génériques, dépendances entre ressources et relations d'autorité, ainsi que l'intégration au sein d'une architecture robotique.
|
173 |
Evolution de la conception radio : de la radio logicielle à la radio intelligenteMoy, Christophe 08 October 2008 (has links) (PDF)
Le manuscrit de cette Habilitation à Diriger des Recherches reprend quelques thèmes de mon parcours de recherche depuis 10 ans, et que l'on peut résumer ainsi. De la conception radio classique pendant ma thèse (à base d'ASIC) j'ai vécu la transition vers la conception radio logicielle (dans un laboratoire de recherche d'un grand groupe industriel) et enfin j'étends désormais en parallèle mes travaux vers la conception radio intelligente (dans un laboratoire de recherche académique). Cinq travaux sont repris dans ce mémoire, portant sur la radio logicielle et la radio intelligente. Ils correspondent à cinq des thèses que j'ai encadrées. Ils portent tout d'abord sur la contribution à une approche de conception pragmatique pour la radio logicielle orientée composants. La radio logicielle en effet implique l'utilisation de ressources de traitements multiples et hétérogènes, générant ainsi des difficultés de conception non encore résolues dans leur ensemble par la communauté de l'architecture. L'un des principes repose sur l'indifférence qu'il doit y avoir entre traitements logiciels et matériels, quel que soit leur niveau de reconfigurabilité. Un effort de formalisation de la conception radio logicielle est proposé dans un deuxième temps, dans une approche orientée opérateurs communs. L'idée repose sur la description à différents niveaux de granularité des opérateurs de traitement du signal afin de bénéficier de la réutilisation de motifs de calcul. Le travail d'optimisation consiste à sélectionner le niveau de granularité adéquate suivant les besoins du concepteur. La gestion de reconfiguration, qui est un point central si ce n'est le point distinctif de la radio logicielle est le cœur du troisième travail présenté. Sa généralisation à la radio intelligente permet d'aboutir dans un quatrième cas d'étude à la proposition d'une architecture de gestion radio intelligente : HDCRAM. Une modélisation de cette architecture est proposée sous la forme d'un méta-modèle exécutable. Cela permet notamment de disposer d'un simulateur de scénarios du comportement de l'architecture de gestion d'une radio intelligente et de pouvoir la dimensionner pour chaque scénario. Enfin la dernière étude porte sur la définition du concept de bulle sensorielle radio intelligente qui porte particulièrement sur l'ensemble des capteurs d'un équipement de radio intelligente. Les perspectives de ces travaux sont notamment d'étudier le troisième et dernier élément du cercle cognitif simplifié (après la reconfiguration et les capteurs) qu'est la prise de décision.
|
174 |
Modélisation mathématique du poumon humainChristine, Vannier 09 July 2009 (has links) (PDF)
Nous nous intéressons à certains problèmes théoriques posés par la modélisation du poumon humain comme arbre bronchique plongé dans le parenchyme pulmonaire. L'arbre bronchique est représenté par un arbre dyadique résistif à 23 générations dans lequel un écoulement de Stokes a lieu. La loi de Poiseuille relie ainsi le débit dans chaque bronche au saut de pression à ses extrémités. Cet arbre est ensuite plongé dans un milieu visco-élastique modélisant le parenchyme. Le processus de ventilation est alors assuré par des pressions négatives, dues à une contraction du diaphragme, au niveau des alvéoles permettant l'inspiration. La première partie est consacrée à l'introduction d'un modèle d'arbre infini obtenu en faisant tendre le nombre de générations vers l'infini. Des théorèmes de trace permettent alors de modéliser le processus de ventilation comme un opérateur Dirichlet-Neumann, qui associe au champ de pression sur l'ensemble des bouts de l'arbre infini le continuum de débit sortant. La seconde partie est dédiée à l'étude de modèles du parenchyme pulmonaire. La complexité du parenchyme, milieu visco-élastique, provient de la présence de l'arbre qui relie toutes les alvéoles entre elles. Des phénomènes de dissipation non locaux sont ainsi observés dus aux couplage de toutes les sorties. Nous étudions tout d'abord un modèle monodimensionnel du parenchyme mettant en jeu une équation de type onde avec des effets non locaux. En particulier nous détaillons l'étude du comportement en temps long. Enfin, nous proposon l'ébauche d'un modèle du parenchyme en dimension supérieure prenant en compte à la fois le caractère élastique du tissu ainsi que la présence de l'arbre résistif.
|
175 |
Contrôle non destructif via des sondes courants de Foucault : nouvelles approchesJiang, Zixian 23 January 2014 (has links) (PDF)
L'objectif principal de cette thèse est de proposer et de tester quelques méthodes de l'optimisation de forme afin d'identifier et de reconstruire des dépôts qui couvrent la paroi extérieure d'un tube conducteur dans un générateur de vapeurs en utilisant des signaux courant de Foucault. Ce problème est motivé par des applications industrielles en contrôle non-destructive dans le secteur de l'énergie nucléaire. En fait, des dépôts peuvent obstruer le passage de circuit de refroidissement entre les tubes et les plaques entretoises qui les soutiennent, ce qui entraînerait une baisse de productivité et mettrait la structure en danger. On considère dans un premier temps un cas axisymétrique dans le cadre duquel on construit un modèle 2-D par des équations aux dérivées partielles pour le courant de Foucault, ce qui nous permet ensuite de reproduire des mesures d'impédances qui correspondent en réalité les signaux courants de Foucault. Pour réduire le coût de calculs de la simulation numérique, on tronque le domaine du problème en posant des conditions aux bords artificielles basées sur des études sur le comportement de la solution, notamment sur un calcul semi-analytique de l'opérateur Dirichlet-to-Neumann dans la direction axiale. Pour la partie identification et reconstruction, on distingue deux sortes de dépôts et établit pour chacun une méthode d'inversion spécifique. Le premier cas concernent des dépôts dont la conductivité est relativement faible (d'environs 1.e4 S/m). On utilise la méthode d'optimisation de forme qui consiste à exprimer explicitement la dérivée des mesures d'impédance par rapport au domaine du dépôt en fonction d'une déformation et à représenter le gradient d'un fonctionnel de coût à minimiser par l'intermédiaire d'un état adjoint proprement défini. Motivée par la présence des dépôts et des plaques de maintient non-axisymétriques, on fait aussi une extension en 3-D des méthodes précédentes. Pour le deuxième cas, des dépôts sont fortement conducteurs et sous forme de couche mince d'une épaisseur de l'ordre de micron. La méthode adaptée à la première sorte de dépôts devient ici trop coûteuse à cause du degré de liberté très élevé des éléments finis sur un maillage extrêmement raffiné à l'échelle de la couche mince. Pour relever cette difficulté, les études sont portées sur plusieurs modèles asymptotiques qui remplace la couche mince par des conditions d'interface sur la surface du tube portante du dépôt. Le nom de modèle asymptotique vient du fait que les conditions d'interface effectives sont obtenues par le développement asymptotique de la solution en fonction d'un paramètre caractérisant la conductivité et l'épaisseur de la couche mince. La validation numérique a permis de retenir un modèle asymptotique qui est à la fois suffisamment précis et simple à inverser. L'inversion (recherche de l'épaisseur du dépôt) consiste alors à rechercher des paramètres dans les conditions d'interface (non standard). Pour les deux cas, la validation et des exemples numériques sont proposés pour le modèle direct et l'inversion.
|
176 |
Méthodes stochastiques en dynamique moléculaire / Stochastic methods in molecular dynamicPerrin, Nicolas 20 March 2013 (has links)
Cette thèse présente deux sujets de recherche indépendants concernant l'application de méthodes stochastiques à des problèmes issus de la dynamique moléculaire. Dans la première partie, nous présentons des travaux liés à l'interprétation probabiliste de l'équation de Poisson-Boltzmann qui intervient dans la description du potentiel électrostatique d'un système moléculaire. Après avoir introduit l'équation de Poisson-Boltzmann et les principaux outils mathématiques utilisés, nous nous intéressons à l'équation linéaire parabolique de Poisson-Boltzmann. Avant d’énoncer le résultat principal de la thèse, nous étendons des résultats d'existence et unicité des équations différentielles stochastiques rétrogrades. Nous donnons ensuite une interprétation probabiliste de l'équation non-linéaire de Poisson-Boltzmann sous la forme de la solution d'une équation différentielle stochastique rétrograde. Enfin, dans une seconde partie prospective, nous commençons l'étude d'une méthode proposée par Paul Malliavin de détection des variables lentes et rapides d'une dynamique moléculaire. / This thesis presents two independent research topics. Both are related to the application of stochastic problems to molecular dynamics. In the first part, we present a work related to the probabilistic interpretation of the Poisson-Boltzmann equation. This equation describes the electrostatic potential of a molecular system. After an introduction to the Poisson-Boltzmann equation, we focus on the parabolic and linear equation. After extending an existence and uniqueness result for backward stochastic differential equations, we establish a probabilistic interpretation of the nonlinear Poisson-Boltzmann equation with backward stochastic differential equations. Finally, in a more prospective second part, we initiate a study of a slow and fast variables detection method due to Paul Malliavin.
|
177 |
Solveur parallèle pour l’équation de Poisson sur mailles superposées et hiérarchiques, dans le cadre du langage Python / Parallel solver for the Poisson equation on a hierarchy of superimposed meshes, under a Python frameworkTesser, Federico 11 September 2018 (has links)
Les discrétisations adaptatives sont importantes dans les problèmes de fluxcompressible/incompressible puisqu'il est souvent nécessaire de résoudre desdétails sur plusieurs niveaux, en permettant de modéliser de grandes régionsd'espace en utilisant un nombre réduit de degrés de liberté (et en réduisant letemps de calcul).Il existe une grande variété de méthodes de discrétisation adaptative, maisles grilles cartésiennes sont les plus efficaces, grâce à leurs stencilsnumériques simples et précis et à leurs performances parallèles supérieures.Et telles performance et simplicité sont généralement obtenues en appliquant unschéma de différences finies pour la résolution des problèmes, mais cetteapproche de discrétisation ne présente pas, au contraire, un chemin faciled'adaptation.Dans un schéma de volumes finis, en revanche, nous pouvons incorporer différentstypes de maillages, plus appropriées aux raffinements adaptatifs, en augmentantla complexité sur les stencils et en obtenant une plus grande flexibilité.L'opérateur de Laplace est un élément essentiel des équations de Navier-Stokes,un modèle qui gouverne les écoulements de fluides, mais il se produit égalementdans des équations différentielles qui décrivent de nombreux autres phénomènesphysiques, tels que les potentiels électriques et gravitationnels. Il s'agitdonc d'un opérateur différentiel très important, et toutes les études qui ontété effectuées sur celui-ci, prouvent sa pertinence.Dans ce travail seront présentés des approches de différences finies et devolumes finis 2D pour résoudre l'opérateur laplacien, en appliquant des patchsde grilles superposées où un niveau plus fin est nécessaire, en laissant desmaillages plus grossiers dans le reste du domaine de calcul.Ces grilles superposées auront des formes quadrilatérales génériques.Plus précisément, les sujets abordés seront les suivants:1) introduction à la méthode des différences finies, méthode des volumes finis,partitionnement des domaines, approximation de la solution;2) récapitulatif des différents types de maillages pour représenter de façondiscrète la géométrie impliquée dans un problème, avec un focussur la structure de données octree, présentant PABLO et PABLitO. Le premier estune bibliothèque externe utilisée pour gérer la création de chaque grille,l'équilibrage de charge et les communications internes, tandis que la secondeest l'API Python de cette bibliothèque, écrite ad hoc pour le projet en cours;3) la présentation de l'algorithme utilisé pour communiquer les données entreles maillages (en ignorant chacune l'existence de l'autre) en utilisant lesintercommunicateurs MPI et la clarification de l'approche monolithique appliquéeà la construction finale de la matrice pour résoudre le système, en tenantcompte des blocs diagonaux, de restriction et de prolongement;4) la présentation de certains résultats; conclusions, références.Il est important de souligner que tout est fait sous Python comme framework deprogrammation, en utilisant Cython pour l'écriture de PABLitO, MPI4Py pour lescommunications entre grilles, PETSc4py pour les parties assemblage et résolutiondu système d'inconnues, NumPy pour les objets à mémoire continue.Le choix de ce langage de programmation a été fait car Python, facile àapprendre et à comprendre, est aujourd'hui un concurrent significatif pourl'informatique numérique et l'écosystème HPC, grâce à son style épuré, sespackages, ses compilateurs et pourquoi pas ses versions optimisées pour desarchitectures spécifiques. / Adaptive discretizations are important in compressible/incompressible flow problems since it is often necessary to resolve details on multiple levels,allowing large regions of space to be modeled using a reduced number of degrees of freedom (reducing the computational time).There are a wide variety of methods for adaptively discretizing space, but Cartesian grids have often outperformed them even at high resolutions due totheir simple and accurate numerical stencils and their superior parallel performances.Such performance and simplicity are in general obtained applying afinite-difference scheme for the resolution of the problems involved, but this discretization approach does not present, by contrast, an easy adapting path.In a finite-volume scheme, instead, we can incorporate different types of grids,more suitable for adaptive refinements, increasing the complexity on thestencils and getting a greater flexibility.The Laplace operator is an essential building block of the Navier-Stokes equations, a model that governs fluid flows, but it occurs also in differential equations that describe many other physical phenomena, such as electric and gravitational potentials, and quantum mechanics. So, it is a very importantdifferential operator, and all the studies carried out on it, prove itsrelevance.In this work will be presented 2D finite-difference and finite-volume approaches to solve the Laplacian operator, applying patches of overlapping grids where amore fined level is needed, leaving coarser meshes in the rest of the computational domain.These overlapping grids will have generic quadrilateral shapes.Specifically, the topics covered will be:1) introduction to the finite difference method, finite volume method, domainpartitioning, solution approximation;2) overview of different types of meshes to represent in a discrete way thegeometry involved in a problem, with a focuson the octree data structure, presenting PABLO and PABLitO. The first one is anexternal library used to manage each single grid’s creation, load balancing and internal communications, while the second one is the Python API ofthat library written ad hoc for the current project;3) presentation of the algorithm used to communicate data between meshes (beingall of them unaware of each other’s existence) using MPI inter-communicators and clarification of the monolithic approach applied building the finalmatrix for the system to solve, taking into account diagonal, restriction and prolongation blocks;4) presentation of some results; conclusions, references.It is important to underline that everything is done under Python as programmingframework, using Cython for the writing of PABLitO, MPI4Py for the communications between grids, PETSc4py for the assembling and resolution partsof the system of unknowns, NumPy for contiguous memory buffer objects.The choice of this programming language has been made because Python, easy to learn and understand, is today a significant contender for the numerical computing and HPC ecosystem, thanks to its clean style, its packages, its compilers and, why not, its specific architecture optimized versions.
|
178 |
Résultats de convergence pour les inéquations variationnelles et applications en mécanique du contact / Convergence results for variational inequalities and applications in contact mechanicsBenraouda, Ahlem 06 June 2018 (has links)
Le sujet de cette thèse porte sur quelques résultats de convergence pour les inéquations variationnelles avec applications dans l'étude des problèmes aux limites décrivant le contact entre un corps déformable et une fondation. La thèse est composée de deux parties. Dans la première partie, nous nous intéressons à l'analyse des inéquations quasivariationnelles, avec ou sans opérateurs de mémoire, dans un espace de Hilbert. Nous prouvons plusieurs résultats de convergence liés à la perturbation de l'ensemble des contraintes ainsi qu'à une méthode de pénalisation. Aussi, pour une classe d'inéquations quasivariationnelles avec opérateurs de mémoire nous étudions une formulation duale pour laquelle nous présentons des résultats d'existence, d'unicité et d'équivalence. La deuxième partie est consacrée à l'application de ces résultats abstraits dans l'étude de six problèmes de contact pour des matériaux élastiques, viscoélastiques et viscoplastiques, dans le cas statique ou quasistatique. Les lois de contact considérées sont la loi de Signorini, la loi de contact avec compliance normale et contrainte unilatérale et la loi de contact avec contrainte unilatérale et seuil critique. Enfin, nous étudions un nombre de problèmes de contrôle optimal associés aux certains modèles de contact. Pour ces problèmes nous obtenons des résultats d'existence et de convergence. / The topic of this thesis concerns some convergence results for variational inequalities with applications in the study of boundary value problems which describe the contact between a deformable body and a foundation. The thesis is divided into two parts. In the first part, we are interested in the analysis of quasivariational inequalities, with or without history-dependent operators, in Hilbert spaces. We prove some convergence results related to a perturbation of the set of constraints and a penalty method, as well. Moreover, for a class of history-dependent quasivariational inequalities we study a dual formulation for which we present existence, uniqueness and equivalence results. The second part is devoted to applications of these abstract results in the study of six contact problems with elastic, viscoelastic and viscoplastic materials, both in the static or quasistatic case. The contact conditions we consider are the Signorini condition, the normal compliance condition with unilateral constraint, the unilateral constraint condition with yield limit. Finally, we study a number of optimal control problems associated to some contact models. For these problems we provide existence and convergence results.
|
179 |
Apport de l’analyse temps-fréquence combinée à l’analyse de formes pour le traitement ISARCorretja, Vincent 30 January 2013 (has links)
Dans le cadre de la surveillance maritime, les opérationnels ont de plus en plus recours à l'imagerie radar pour classifier à grande distance un objet marin. Le traitement ISAR (Inverse Synthetic Aperture Radar) répond à ce besoin. Il repose en particulier sur l'analyse des mouvements propres de l'objet marin. Une fois l'objet détecté, il s'agit d'afficher sur la console tactique la représentation de la fréquence Doppler en fonction de la distance, aussi appelée image range-Doppler. Le travail présenté dans ce mémoire s'inscrit dans une perspective d'évolution opérationnelle de la chaîne de traitement existante. Il vise à produire de manière automatique la « meilleure » image range-Doppler. Dans cette thèse, nos contributions s'appuient sur l'idée de reconsidérer la chaîne de traitement en tenant compte de l'a priori que l'objet marin est un objet rigide dont la géométrie structure l'évolution du signal radar. Ainsi, dans une première contribution, nous proposons une nouvelle méthode d'analyse temps-fréquence du signal radar afin d'obtenir une image instantanée où l'opérationnel peut distinguer « au mieux » les superstructures de l'objet marin. Cette dernière est fondée sur la fusion de plusieurs représentations temps-fréquence issues de la classe de Cohen en faisant l'hypothèse que les composantes temps-fréquence sont des trajectoires structurées 2D dans le plan temps-fréquence, contrairement aux termes d'interférences induits par la propriété de bilinéarité des membres de cette classe. Une étude comparative sur données synthétiques et ISAR est menée pour confirmer la pertinence de notre approche, notamment du point de vue de la résolution temps-fréquence et de la suppression des termes d'interférences.Dans une seconde contribution, nous établissons une nouvelle procédure pour qualifier chaque image range-Doppler, obtenue à l'issue de l'analyse temps-fréquence, avec des mesures d'irrégularité de formes que nous fusionnons à l'aide d'un opérateur d'agrégation. Des simulations sur données réelles sont réalisées. Les résultats concordent avec une analyse subjective menée par des opérationnels, ce qui confirme l'efficacité de notre méthode. / In maritime surveillance, radar imaging plays a key role to classify a maritime object. ISAR processing is one of the solutions, which takes advantage of the object rotational motion to provide a range-Doppler image.The work, presented in this report, is an evolution of the existing ISAR processing chain. Therefore, our contributions are based on the processing chain reconsideration by taking into account the fact that the maritime object is a rigid object, the geometry of which influences the radar signal evolution.In a first contribution, we propose a new time-frequency analysis method based on the aggregation of some time-frequency representations obtained with Cohen class members. It consists in differentiating the signal, assumed to be characterized by 2-D near-linear stable trajectories in the time-frequency plane, and the cross-terms, assumed to be geometrically unstructured. A comparative study is then carried out on ISAR synthetic data to confirm the efficiency of our approach.In a second contribution, we present a new procedure to characterize each range-Doppler image, obtained from a time-frequency analysis, by means of shape irregularity measures that are combined with a fuzzy logic operator. To validate our approach, simulations on real data are done. The results are compared to a subjective analysis carried out with practionners.
|
180 |
Etude de l'asymptotique du phénomène d'augmentation de diffusivité dans des flots à grande vitesse / The asymptotic of the phenomenon of enhancement of diffusivity in high speed flowNguyen, Thi-Hien 29 September 2017 (has links)
En application, on souhaite générer des nombres aléatoires avec une loi précise (méthode de Monte Carlo par chaines de Markov - MCMC (Markov Chaine Monte Carlo)). La méthode consiste à trouver une diffusion qui a la loi invariante souhaitée et à montrer la convergence de cette diffusion vers son équilibre avec une vitesse exponentielle. L’exposant de cette convergence est le trou spectral du générateur. Il a été montré par Chii-Ruey Hwang, Shu-Yin Hwang-Ma, et Shuenn-Jyi Sheu qu’on peut agrandir le trou spectral, en rajoutant un terme non-symétrique au générateur auto-adjoint (souvent utilisé en MCMC). Ceci correspond à passer d’une diffusion réversible (en detailed balance) à une diffusion non réversible. Un moyen de construire une diffusion non-réversible avec la même mesure invariante est de rajouter un flot incompressible à la dynamique de la diffusion réversible.Dans cette thèse, nous étudions le comportement de la diffusion lorsqu’on accélère le flot sous-jacent en multipliant le champ des vecteurs qui le décrit par une grande constante. P. Constantin, A.Kisekev, L.Ryzhik et A.Zlatoš (2008) ont montré que si le flot était faiblement mélangeant alors l’accélération du flot suffisait pour faire converger la diffusion vers son équilibre en un temps fini. Dans ce travail, on explicite la vitesse de ce phénomène sous une condition de corrélation du flot. L’article de B. Franke, C.-R.Hwang, H.-M. Pai et S.-J. Sheu (2010) donne l’expression asymptotique du trou spectral lorsque le flot sous-jacent est accéléré vers l’infini. Ici aussi, on s’intéresse à la vitesse avec laquelle le phénomène se manifeste. Dans un premier temps, nous étudions le cas particulier d’une diffusion du type Ornstein-Uhlenbeck qui est perturbée par un flot préservant la mesure gaussienne. Dans ce cas, grâce à un résultat de G. Metafune, D. Pallara et E. Priola (2002), nous pouvons réduire l’étude du spectre du générateur à des valeurs propres d’une famille de matrices. Nous étudions ce problème avec des méthodes de développement limité des valeurs propres. Ce problème est résolu explicitement dans cette thèse et nous donnons aussi une borne pour le rayon de convergence du développement. Nous généralisons ensuite cette méthode dans le cas d’une diffusion générale de façon formelle. Ces résultats peuvent être utiles pour avoir une première idée sur les vitesses de convergence du trou spectral décrites dans l’article de Franke et al. (2010). / In application, we would like to generate random numbers with a precise law MCMC (Markov Chaine Monte Carlo). The method consists in finding a diffusion which has the desired invariant law and in showing the convergence of this diffusion towards its equilibrium with an exponential rate. The exponent of this convergence is the spectral gap of the generator. It was shown by C.-R. Hwang, S.-Y. Hwang-Ma and S.-J. Sheu that the spectral gap can grow up by adding a non-symmetric term to the self-adjoint generator.This corresponds to passing from a reversible diffusion to a non-reversible diffusion. A means of constructing a non-reversible diffusion with the same invariant measure is to add an incompressible flow to the dynamics of the reversible diffusion.In this thesis, we study the behavior of diffusion when the flow is accelerated by multiplying the field of the vectors which describes it by a large constant. In 2008, P. Constantin, A. Kisekev, L. Ryzhik and A. Zlatoˇs have shown that if the flow was weakly mixing then the acceleration of the flow was sufficient to converge the diffusion towards its equilibrium after finite time. In this work, the speed of this phenomenon is explained under a condition of correlation of the flow. The article by B. Franke, C.-R.Hwang, H.-M. Pai and S.-J.Sheu (2010) gives the asymptotic expression of the spectral gap when the large constant goes to infinity. Here we are also interested in the speed with which the phenomenon manifests itself. First, we study the special case of an Ornstein-Uhlenbeck diffusion which is perturbed by a flow preserving the Gaussian measure. In this case, thanks to a result of G. Metafune, D. Pallara and E. Priola (2002), we can reduce the study of the generator spectrum to eigenvalues of a family of matrices. We study this problem with methods of limited development of eigenvalues. This problem is solved explicitly in this thesis and we also give a boundary for the convergence radius of the development. We then generalize this method in the case of a general diffusion in a formal way. These results may be useful to have a first idea on the speeds of convergence of the spectral gap described in the article by Franke et al. (2010).
|
Page generated in 0.077 seconds