Bonjour,

je suis dans le cas suivant:
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 !!!
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
 
eig(x*x'/(T)-mu*mu')
Cela vous semble normal ? Qu'ais-je raté ?