Bonjour à tous !

J'aimerais savoir s'il existe un moyen simple (du compliqué, je peux aussi en faire ^^) de tirer un nombre aléatoire appartenant non pas à un intervalle, mais à un parmi deux intervalles.

Par exemple, je voudrais tirer un nombre appartenant à [-1;-0.5[U]0.5;1] : est-ce possible ?

Merci d'avance

-
Edité par Turambar 2 février 2014 à 15:32:08

ça ne me dit rien mais en attendant on peut tricher ;-)

from random import uniform

intervalle = [(-1, -0.5), (0.5, 1)]

def tirage(intervalle):
    while 1:
        a = uniform(intervalle[0][0], intervalle[1][1])
        if a > intervalle[0][1] and a < intervalle[1][0]:
            tirage(intervalle)
        else: return a
    
print(tirage(intervalle))

Tu peux définir deux intervalles tout simplement en les nommant de deux façons différentes

Oui, ça je sais faire, mais je voudrais une façon de le faire en une seule ligne, s'il existe une fonction permettant de le faire... mais merci quand même pour votre aide

def tirage(*args):
    return uniform(*choice(args))

tirage((-1, -0.5), (0.5, 1))

Non testé

Ne pas oublier d'importer uniform et choice du module random...

-
Edité par Anonyme 2 février 2014 à 16:55:55

Ça marche parfaitement, merci énormément

J'ai remodifié je pense qu'on peut faire plus léger

PS : par contre, pourrais-tu, stp, m'expliquer le fonctionnement du *choice(args) ?

choice est une fonction qui dans ce cas, va retourner un intervalle de type tuple

Le problème c'est que la fonction uniform prend deux arguments de type float, int et pas un tuple

En mettant l'étoile devant c'est comme si je faisais

x, y = choice(args)
n = uniform(x, y)

Merci, je ne connaissais pas ça

Bonjour.

Choisir un intervalle puis tirer aléatoirement un nombre dedans n'aura pas la même distribution que de tirer aléatoirement un nombre dans l'union des 2 intervalles !!!

Par exemple pour des listes d'entiers (problème similaire), pour {1} U [[0;100]], tu auras 1 chance sur 2 de tirer 1 (dans l'ensemble {1}) alors que tu auras 1 chance sur 200 (=2*100) de tirer 99 (dans l'ensemble [[0;100]]).

La solution est de joindre les 2 intervalles :

import random

I1 = (-1, -0.5)
I2 = (0.5, 1)

I3 = (I1[0], I1[1]+I2[1]-I2[0])
res = random.uniform(I3[0], I3[1])
if I1[1] < res:
	res += I2[0]
print(res)

Il n'y a ici aucun contrôle sur les intervalles. Si tu veux faire quelque chose de propre il faudra vérifier qu'ils se superposent et qu'ils sont corrects.

Effectivement ; toutefois, dans le cas présent, mes deux intervalles sont de même taille, on devrait donc bien avoir une (pseudo-)équiprobabilité, non ?

Salut, à la limite, je comprends pas pourquoi tu te prends la tête comme ça. Un petit

nb=random.uniform(-1;0)
nb+=(nb>=-0.5)

suffit, non ?

-
Edité par Anonyme 2 février 2014 à 17:52:44

Salut,

Pour aller plus loin, on peut généraliser à plusieurs intervalles de largeur différente:

import random
import matplotlib.pyplot as plt

def tirage(n, *intervalles):
    largeurs = [ max(i) - min(i) for i in intervalles ]
    s = sum(largeurs)
    def eggs(x):
        for idx, interv in enumerate(intervalles):
            if x < sum(largeurs[:idx+1]):
                return min(interv) + x - sum(largeurs[:idx])
    return [ eggs(random.uniform(0, s)) for _ in range(n) ]

def main():
    intervalles = ( (-1., -0.6), (-0.4, -0.1), (0.4, 1.), (1.1, 2.) )
    res = tirage( int(1e6), *intervalles)

    plt.hist(res, bins=60)
    plt.show()

if __name__ == '__main__':
    main()

-
Edité par azertymen 2 février 2014 à 19:24:11

Turambar a écrit:

Effectivement ; toutefois, dans le cas présent, mes deux intervalles sont de même taille, on devrait donc bien avoir une (pseudo-)équiprobabilité, non ?

@dri1 a écrit:

Salut,
à la limite, je comprends pas pourquoi tu te prends la tête comme ça. Un petit

nb=random.uniform(-1;0)
nb+=(nb>=-0.5)

suffit, non ?

Comme tu l'as dit l'union d'intervalles que tu as cité n'était qu'un exemple, donc la solution que l'on te propose doit correspondre au problème et non à l'exemple