loi binomiale

[morane]

bonjour.
Normalement, il y a 95% des valeurs entre moyenne -2sigma et moyenne + 2sigma.
Le fichier montre 100 expériences d'un mouvement brownien de 60 pas, en partant de l'origine.
Or il s'avère que 100% des positions finales sont dans l'intervalle à 95%
Tous les chemins semblent partir d'un même côté.
Qui peut m'expliquer ? y a-t-il une erreur ?
merci c'est urgent.
Morane.

 

[morane]

Re : loi binomiale

le fichier que j'ai envoyé ne fonctionne pas bien.
Voici le bon

 

[herve62]

Re : loi binomiale

Bonsoir
Là il y a un truc à piger , est ce basé sur un plan d'expériences ( binomiale pure) ou une approximation de la binomiale par la loi de Poisson ?
Je ne retrouve pas la formule P(k)=C (n k) . P^k etc ... (on peut pas l'écrire ici)
C'est déjà pas évident ( que des souvenirs de stat.) à moins de tomber sur un pro , et si tu donnes pas assez de détails sur la structure de ton VBA .... ca va pas le faire , du style les col en xx ca fait CA , les plages YYY c'est pour ca ...etc , mon erreur est LA !!! explique met des commentaires dans ton fichier
C'est pour une stat d'un parc de machines qui doit pas tomber en panne ? et trouver le risque ?

 

[morane]

Re : loi binomiale

fini...
cela correspond à la dispersion des distances entre la position de départ et la position d'arrivée d'un mouvement aléatoire à 4 possibilités sur 100 pas.
L'erreur provenait d'une double boucle avec random...

 

[ROGER2327]

Re : loi binomiale

Bonsoir à tous.

Peut-être avec ça ?​

Sub Macro1()
Dim i, k, ex As Integer
  Randomize '***
  i = 1
  Columns("C:D").ClearContents
  Columns("M").ClearContents
  ex = Range("J31").Value
  Cells(1, 3).Value = 0
  Cells(1, 4).Value = 0
  For c = 1 To ex
    x = 0: y = 0
    For t = 2 To 60
      x = x + 1 + 2 * (Rnd < 0.5)
      y = y + 1 + 2 * (Rnd < 0.5)
      i = i + 1
      Cells(i, 3).Resize(1, 2).Value = Array(x, y)
    Next t
    i = i + 1
    Cells(i, 3).Resize(2, 2).Value = Array(0, Empty)
    Cells(c, 13).Value = Sqr(x * x + y * y)
  Next c
  Range("K15").Value = x
  Range("L15").Value = y
End Sub

Bonne soirée.


ℝOGER2327
#7347

Mercredi 25 Palotin 141 (Saints Potassons, acolythes - fête Suprême Quarte)
25 Floréal An CCXXII, 9,1191h - carpe
2014-W20-3T21:53:09Z


P.s. : On ne rafraîchit jamais assez ! Je n'avais pas vu le message précédent...

 

[ROGER2327]

Re : loi binomiale

Suite...

Une variation sur le thème...​

ℝOGER2327
#7348

Jeudi 26 Palotin 141 (Sainte Prétentaine, rosière - fête Suprême Quarte)
26 Floréal An CCXXII, 1,1508h - fusain
2014-W20-4T02:45:43Z

 

[Dranreb]

Re : loi binomiale

Bonjour.

Information.
Je pense que pour simuler un mouvement brownien d'un point en 2 dimensions, j'utiliserais ces fonctions en spécifiant comme Moy les valeurs X et Y du dernier dernier point trouvé.

Function AléatN(ByVal ÉTyp As Double, ParamArray Moy() As Variant) As Variant
Dim TRés() As Double, N As Long, Np As Long, AléN1 As Double: Static AléN2 As Double, DéjàDonné As Boolean
If UBound(Moy) <= 0 Then
   GoSub Donner: AléatN = AléN1 * ÉTyp
   If UBound(Moy) = 0 Then AléatN = Moy(0) + AléatN
Else
   ReDim TRés(1 To UBound(Moy) + 1) As Double
   Np = 0: N = 1: Do:
      GoSub Donner: TRés(N) = Moy(Np) + AléN1 * ÉTyp: If N = UBound(TRés) Then Exit Do
      Np = N: N = N + 1: Loop
   AléatN = TRés
   End If
Exit Function
Dim AléU1 As Double, AléU2 As Double, S As Double
Donner: If DéjàDonné Then
   AléN1 = AléN2: DéjàDonné = False
Else:                 AléU1 = Aléat * 2 - 1
   Do: AléU2 = AléU1: AléU1 = Aléat * 2 - 1: S = AléU1 * AléU1 + AléU2 * AléU2: Loop Until S <= 1
   S = Sqr(-2 * Log(S) / S): AléN1 = AléU1 * S: AléN2 = AléU2 * S: DéjàDonné = True: End If
Return
End Function
'

Function Aléat(Optional ByVal G As Double = -1) As Double
Static X As Double
If G >= 0 Then If G > 0 Then X = G Else X = Now / 2958466
X = (X + 1.35198775424545) ^ 7: X = X - Int(X)
Aléat = X
End Function

Remarque: théoriquement Aléat pourrait être remplacé par Rnd. J'avais écrit cette fonction Aléat en constatant, avec Rnd, de nombreuses coïncidences insolentes nuisant à la crédibilité d'un ciel étoilé, qui ne se produisaient pas avec mon Aléat. Mais j'ai gardé un souvenir très vague de ce que celui ci m'a fait un peu le même coup dans un autre contexte. Enfin bref, ma foncion AléatN est faite pour une génération aléatoire en distribution normale autour d'une moyenne réelle, complexe, spatiale, hyper complexe etc.