Discussion :

[miniclems]

Bonjour,

Je souhaiterai générer une courbe de weibull mais inversée.

Je m'explique: j'ai fait un tirage aléatoire grâce à la fonction wblrnd qui me génère les valeurs aléatoires suivant la distribution de weibull et l'histogramme mon montre bien la forme de la courbe comme suit:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
a=1;
b=2;
x = wblrnd(a,b,500,1);
hist (x1,20)
hold on 
 
xx = 0:.1:4 ;
yy = 10*(b1.*(a1.^(-b1)).*xx.^(b1-1).*exp(-(xx./a1).^b1));    %pour tracer %la courbe réélle
plot(xx,yy,'g-','LineWidth',2)
hold off
%cf. weibull1.JPG

ce que je voudrai maintenant, c'est d'avoir un histogramme la même allure que la courbe qui ressemble à l'image suivante (courbe rouge dans weibull2.JPG)

Je ne réussis que quand on fait :

Code :

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

x= -wblrnd(a,b,500,1)

mais les valeurs sont négatives, ce qui ne m'arrange pas.. (cf hist.jpg)

Pourriez vous me dire comment faire pour avoir un histogramme de ce type mais à valeurs positives?

Merci d'avance pour vos réponses.

[svagrim]

il faut que tu fasses un changement de variable.

0 devienne 4, 4 devienne 0 etc ...
du coup si une fois ta valeur x tirée tu dis
y = 4 - x ou un truc du genre ça devrait faire ce que tu veux nan ?

ou alors tu réalises toi même le tirage aléatoire avec un petit code monte carlo