Spelling suggestions: "subject:"toeplitz 2matrices"" "subject:"toeplitz cicatrices""
41 |
Tests non paramétriques minimax pour de grandes matrices de covariance / Non parametric minimax tests for high dimensional covariance matricesZgheib, Rania 23 May 2016 (has links)
Ces travaux contribuent à la théorie des tests non paramétriques minimax dans le modèle de grandes matrices de covariance. Plus précisément, nous observons $n$ vecteurs indépendants, de dimension $p$, $X_1,ldots, X_n$, ayant la même loi gaussienne $mathcal {N}_p(0, Sigma)$, où $Sigma$ est la matrice de covariance inconnue. Nous testons l'hypothèse nulle $H_0:Sigma = I$, où $I$ est la matrice identité. L'hypothèse alternative est constituée d'un ellipsoïde avec une boule de rayon $varphi$ autour de $I$ enlevée. Asymptotiquement, $n$ et $p$ tendent vers l'infini. La théorie minimax des tests, les autres approches considérées pour le modèle de matrice de covariance, ainsi que le résumé de nos résultats font l'objet de l'introduction.Le deuxième chapitre est consacré aux matrices de covariance $Sigma$ de Toeplitz. Le lien avec le modèle de densité spectrale est discuté. Nous considérons deux types d'ellipsoïdes, décrits par des pondérations polynomiales (dits de type Sobolev) et exponentielles, respectivement.Dans les deux cas, nous trouvons les vitesses de séparation minimax. Nous établissons également des équivalents asymptotiques exacts de l'erreur minimax de deuxième espèce et de l'erreur minimax totale. La procédure de test asymptotiquement minimax exacte est basée sur une U-statistique d'ordre 2 pondérée de façon optimale.Le troisième chapitre considère une hypothèse alternative de matrices de covariance pas nécessairement de Toeplitz, appartenant à un ellipsoïde de type Sobolev de paramètre $alpha$. Nous donnons des équivalents asymptotiques exacts des erreurs minimax de 2ème espèce et totale. Nous proposons une procédure de test adaptative, c-à-d libre de $alpha$, quand $alpha$ appartient à un compact de $(1/2, + infty)$.L'implémentation numérique des procédures introduites dans les deux premiers chapitres montrent qu'elles se comportent très bien pour de grandes valeurs de $p$, en particulier elles gagnent beaucoup sur les méthodes existantes quand $p$ est grand et $n$ petit.Le quatrième chapitre se consacre aux tests adaptatifs dans un modèle de covariance où les observations sont incomplètes. En effet, chaque coordonnée du vecteur est manquante de manière indépendante avec probabilité $1-a$, $ ain (0,1)$, où $a$ peut tendre vers 0. Nous traitons ce problème comme un problème inverse. Nous établissons ici les vitesses minimax de séparation et introduisons de nouvelles procédures adaptatives de test. Les statistiques de test définies ici ont des poids constants. Nous considérons les deux cas: matrices de Toeplitz ou pas, appartenant aux ellipsoïdes de type Sobolev / Our work contributes to the theory of non-parametric minimax tests for high dimensional covariance matrices. More precisely, we observe $n$ independent, identically distributed vectors of dimension $p$, $X_1,ldots, X_n$ having Gaussian distribution $mathcal{N}_p(0,Sigma)$, where $Sigma$ is the unknown covariance matrix. We test the null hypothesis $H_0 : Sigma =I$, where $I$ is the identity matrix. The alternative hypothesis is given by an ellipsoid from which a ball of radius $varphi$ centered in $I$ is removed. Asymptotically, $n$ and $p$ tend to infinity. The minimax test theory, other approaches considered for testing covariance matrices and a summary of our results are given in the introduction.The second chapter is devoted to the case of Toeplitz covariance matrices $Sigma$. The connection with the spectral density model is discussed. We consider two types of ellipsoids, describe by polynomial weights and exponential weights, respectively. We find the minimax separation rate in both cases. We establish the sharp asymptotic equivalents of the minimax type II error probability and the minimax total error probability. The asymptotically minimax test procedure is a U-statistic of order 2 weighted by an optimal way.The third chapter considers alternative hypothesis containing covariance matrices not necessarily Toeplitz, that belong to an ellipsoid of parameter $alpha$. We obtain the minimax separation rate and give sharp asymptotic equivalents of the minimax type II error probability and the minimax total error probability. We propose an adaptive test procedure free of $alpha$, for $alpha$ belonging to a compact of $(1/2, + infty)$.We implement the tests procedures given in the previous two chapters. The results show their good behavior for large values of $p$ and that, in particular, they gain significantly over existing methods for large $p$ and small $n$.The fourth chapter is dedicated to adaptive tests in the model of covariance matrices where the observations are incomplete. That is, each value of the observed vector is missing with probability $1-a$, $a in (0,1)$ and $a$ may tend to 0. We treat this problem as an inverse problem. We establish the minimax separation rates and introduce new adaptive test procedures. Here, the tests statistics are weighted by constant weights. We consider ellipsoids of Sobolev type, for both cases : Toeplitz and non Toeplitz matrices
|
42 |
Sur les méthodes rapides de résolution de systèmes de Toeplitz bandes / Fast methods for solving banded Toeplitz systemsDridi, Marwa 13 May 2016 (has links)
Cette thèse vise à la conception de nouveaux algorithmes rapides en calcul numérique via les matrices de Toeplitz. Tout d'abord, nous avons introduit un algorithme rapide sur le calcul de l'inverse d'une matrice triangulaire de Toeplitz en se basant sur des notions d'interpolation polynomiale. Cet algorithme nécessitant uniquement deux FFT(2n) est manifestement efficace par rapport à ses prédécésseurs. ensuite, nous avons introduit un algorithme rapide pour la résolution d'un système linéaire de Toeplitz bande. Cette approche est basée sur l'extension de la matrice donnée par plusieurs lignes en dessus, de plusieurs colonnes à droite et d'attribuer des zéros et des constantes non nulles dans chacune de ces lignes et de ces colonnes de telle façon que la matrice augmentée à la structure d'une matrice triangulaire inférieure de Toeplitz. La stabilité de l'algorithme a été discutée et son efficacité a été aussi justifiée. Finalement, nous avons abordé la résolution d'un système de Toeplitz bandes par blocs bandes de Toeplitz. Ceci étant primordial pour établir la connexion de nos algorithmes à des applications en restauration d'images, un domaine phare en mathématiques appliquées. / This thesis aims to design new fast algorithms for numerical computation via the Toeplitz matrices. First, we introduced a fast algorithm to compute the inverse of a triangular Toeplitz matrix with real and/or complex numbers based on polynomial interpolation techniques. This algorithm requires only two FFT (2n) is clearly effective compared to predecessors. A numerical accuracy and error analysis is also considered. Numerical examples are given to illustrate the effectiveness of our method. In addition, we introduced a fast algorithm for solving a linear banded Toeplitz system. This new approach is based on extending the given matrix with several rows on the top and several columns on the right and to assign zeros and some nonzero constants in each of these rows and columns in such a way that the augmented matrix has a lower triangular Toeplitz structure. Stability of the algorithm is discussed and its performance is showed by numerical experiments. This is essential to connect our algorithms to applications such as image restoration applications, a key area in applied mathematics.
|
43 |
Inégalités de von Neumann sous contraintes, image numérique de rang supérieur et applications à l’analyse harmonique / Constrained von Neumann inequalities, higher rank numarical range and applications to harmonic analysisGaaya, Haykel 05 December 2011 (has links)
Cette thèse s’inscrit dans le domaine de la théorie des opérateurs. L’un des opérateurs qui m’a particulièrement intéressé est l’opérateur modèle noté S(Φ) qui désigne la compression du shift unilatéral S sur l’espace modèle H(Φ) où Φ est une fonction intérieure. L’étude du rayon numérique de S(Φ) semble être importante comme l’illustre bien un résultat dû à C. Badea et G. Cassier qui ont montré qu’il existe un lien entre le rayon numérique de tels opérateurs et l’estimation des coefficients des fractions rationnelles positives sur le tore. Nous fournissons une extension de leur résultat et nous trouvons une expression explicite du rayon numérique de S(Φ) dans le cas particulier où Φ est un produit de Blaschke fini avec un unique zéro. Dans le cas général où Φ est un produit de Blaschke fini quelconque, une estimation du rayon numérique de S(Φ) est aussi donnée. Dans la deuxième partie de cette thèse on s’est intéressé à l’image numérique de rang supérieur Λk(T) qui est l’ensemble de tous les nombres complexes λ vérifiant PTP = λP pour une certaine projection orthogonale P de rang k . Cette notion a été introduite récemment par M.-D. Choi, D. W. Kribs, et K. Zyczkowski et elle est utilisée pour certains problèmes en physique. On montre que l’image numérique de rang supérieur du shift n-dimensionnel coïncide avec un disque de rayon bien déterminé / This thesis joins in the field of operator theory. We are specially interested by the extremal operator S(Φ) defined by the compression of the unilateral shift S to the model subspace H(Φ) where Φ is an inner function on the unit disc. The numerical radius of S(Φ) seems to be important and have many applications to harmonic analysis. C. Badea and G. Cassier showed that there is a relationship between the numerical radius of such operators and the Taylor coefficients of positive rational functions. We give an extension of C. Badea and G. Cassier result and an explicit formula of the numerical radius of S(Φ) in the particular case where Φ is a finite Blaschke product with unique zero. An estimate in the general case is also established. The second part is devoted to the study of the higher rank-k numerical range denoted by Λk(T) which is the set of all complex number λ satisfying PTP = λP for some rank-k orthogonal projection P. This notion was introduced by M.-D. Choi, D. W. Kribs, et K. Zyczkowski motivated by a problem in Physics. We show that if Sn is the n-dimensional shift then its rank-k numerical range is the circular discentered in zero and with a precise radius
|
Page generated in 0.0557 seconds