Discussion :

[annibel]

Bonjour,


Qu'elle est la fonction "random" sous matlab, qui permet de générer des valeurs différentes qui ont une moyenne ( par exemple moyenne=200)
Par exemple, elle peut générer : 60, 100, 220, 500, 120
Ce sont 5 nombres qui ont une moyenne égale à 200

Merci pour vous tous

[FLB]

La fonction rand te permet de faire un tirage uniforme entre 0 et 1:

Code :

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

doc rand

Pour un tirage d'entiers entre 100 et 300, de moyenne 200 :

Code :

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

r =  100+round((300-100).*rand)

Tu préféreras peux être un tirage gaussien

Code :

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

doc randn

dans ce cas

Code :

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

r =  200+round(sigma*randn)

Permet de tirer des entiers, centrés sur 200, avec une répartition gaussienne sigma.

[shaiHulud]

J'ai également le souvenir de fonctions dans la financial toolbox qui permettent de générer un échantillon dont le moyenne EMPIRIQUE est égal à l'espérance ; si toutefois c'est ce que tu cherches à faire. Un tel échantillon ne saurait être indépendant !

[ArnaudFu13]

Bonjour,

cette fonction pourrait t'aider (sauf si tu as trouvé avant)

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
function [val] = randCentre (valCent, nbVal)
    for i=1:nbVal
        val(i) = abs(randi(1000));
    end
 
    rajout = nbVal - mod (sum(val),nbVal);
    val(1) = val(1) + rajout;
 
    moy = mean(val);  
    delta = valCent - moy;
    for i=1:nbVal
        val(i) = val(i) + delta;
    end
 
    moy = mean (val);
    minimum = min (val);
    while minimum < 0
        indice = find (val == minimum);
        rajout = abs(randi(50));
        minimum = abs(minimum) + rajout;  
        val(indice)  = val(indice)+ minimum;
 
        maximum = max (val);
        indiceMax = find (val == maximum);
        val(indiceMax) = val(indiceMax) - minimum;
 
        minimum = min(val);
    end
 
    val
    moy = mean(val)