Bonjour,
je suis dans le cas suivant:
Je veux bien que les imprécisions numériques conduisent Matlab à cela. Ce qui m'inquiète plus, c'est que x*x'/(T)-mu*mu' a une valeur propre négative très forte, alors que je m'attendais à une v.p très petite (celle qui aurait du être nulle mais qui ne l'est pas à cause des erreurs numériques)
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8 x= rand(10,15); [T,P]= size(x); mu= mean(x); z= bsxfun(@minus,x,mu); rank(z) % ->9 rank(x'*x/(T)-mu'*mu) % -> 9 rank(x*x'/(T)-mu*mu') % -> 10 !!!
Cela vous semble normal ? Qu'ais-je raté ?
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2 eig(x*x'/(T)-mu*mu')
Partager