Chaque matrice définie semi-positive correspond-elle à une matrice de covariance?

12

Il est bien connu qu'une matrice de covariance doit être définie semi-positive, mais l'inverse est-il vrai?

Autrement dit, chaque matrice définie semi-positive correspond-elle à une matrice de covariance?

Jingjings
la source

Réponses:

12

En suivant les définitions de PD et PSD ici , oui, je pense que oui, car nous pouvons le faire par construction. Je suppose pour un argument un peu plus simple que vous voulez dire pour les matrices avec des éléments réels, mais avec des changements appropriés, cela s'étendrait aux matrices complexes.

Soit une vraie matrice PSD; d'après la définition à laquelle je suis lié, il sera symétrique. Toute matrice réelle symétrique définie positive A peut être écrit sous la forme A = L L T . Cela peut être fait par L = Q UNEUNEUNE=LLTsiA=QDQTavec orthogonaleQet diagonaleDetL=QQTUNE=QQTQ en tant que matrice de racines carrées sages composantes deD. Ainsi, il n'a pas besoin d'être complet.

Soit une variable aléatoire vectorielle, de dimension appropriée, avec la matrice de covariance I (qui est facile à créer).Zje

Puis a covariance matrice A .LZUNE

[Au moins c'est en théorie. Dans la pratique, il y aurait divers problèmes numériques à résoudre si vous vouliez de bons résultats et - en raison des problèmes habituels avec le calcul en virgule flottante - vous n'obtiendrez qu'environ ce dont vous avez besoin; qui est, la variance de la population d'une calculée ne serait habituellement pas exactement A . Mais ce genre de chose est toujours un problème quand on arrive à calculer les choses]LZ UNE

Glen_b -Reinstate Monica
la source
S'il est vrai qu'une décomposition est possible sans rang complet, l'algorithme de Cholesky ne fonctionne qu'avec A normal . Donc, sans rang complet, il ne peut pas s'agir d'une décomposition de Cholesky. Calculativement, on pourrait faire cette décomposition au cas singulier par diagonalisation. (Bien que ce soit beaucoup plus cher)UNE=LLUNE
Horst Grünbusch
@Horst: Pourquoi être triangulaire inférieur? L=QQT
amibe dit Réintégrer Monica
1
@amoeba Bien que l'on puisse l'organiser tel quel, il n'est pas nécessaire qu'il soit triangulaire inférieur pour que l'argument fonctionne - c'est une caractéristique du Cholesky mais ce n'est pas nécessaire pour que le résultat fonctionne.
Glen_b -Reinstate Monica
1
@Glen Est-ce que la symétrie est une condition nécessaire pour être PSD ou cette définition est-elle une des nombreuses?
114
1
@ 114 pour la relation entre symétrique et PSD voir math.stackexchange.com/questions/516533/…
Frank