Discussion :

[khalil.ajmi]

bonjour a tous ;
je suis entrain d'étudier quelque aspects sur le comportements des valeurs propres des matrices aléatoires alors j'ai écris ce petit code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
%------------------------------------
clear all;
close all;
n=100;
m=100;
for i=1:10000
 
h=rand(n,m);
r=(1/n).*(h*h');
m=max(eig(r));
l(i)=trace(r)/m;
 
end 
figure;
plot(l);
axis([0.0001,5,0.0001,1.5]);
 
%-------------------------------------

j'ai simulé le code mais j'ai bien eu une figure qui décrit un comportement sur la somme des valeur propres deviseés par la valeur propre maximale mais j'ai pas pu identifié la loi de probabilité .est ce que quelqu'un pourrait m'aider a l'identifier ?
prière de me répondre .
merci

[Invité]

Bonsoir khalil.ajmi,

Attention, dans ta boucle for, tu modifies la valeur de m en prenant la valeur maximale de tes valeurs propres.
Essaies plutôt:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
for i=1:10000
    h=rand(n,m);
    r=(1/n).*(h*h');
    M=max(eig(r));
    l(i)=trace(r)/M;
end

Tu as surement du avoir le message suivant pour cela car m n'est plus entier:

Warning: Size vector should be a row vector with integer elements

[khalil.ajmi]

merci bien pour votre réponse .
en fait je n'ai pas eu le message que vous m'aviez indiqué plutôt j'ai eu

Warning: Divide by zero


a part ca le script tourne et me lance la figure .