Discussion :

[Jack_dev]

Bonjour,

Je veux définir une fonction avec plusieurs éléments, par exemple:

def function((a,b,c), (d, e, f), (m, n, o, p, q), (r, s, t, u, v)):
où les item (r, s, t, u, v) correspondent respectivement à une probabilité pour chacun des items (m, n, o, p, q)

la fonction effectue différents calculs et donne un résultat "z".

Je veux créer une boucle de façon à obtenir un "z" maximum en utilisant un (m, n, o, p, q) de départ et en faisant varier les données par 0.10 à chaque itération.

La somme des items (m, n, o, p, q) doit toujours être égale à 1 et je veux incrémenter les valeurs de 0.10 uniquement, pour limiter les possibilités:

Exemples de valeurs permises: (0.1, 0.2, 0.4, 0.1, 0.2), ....(0.3, 0.0, 0.2, 0.2, 0.3)...et ainsi de suite.....

Chaque itération me donnera un nouveau "z", et je désire sortir de la boucle avec le plus grand "z".....

Comment faire varier uniquement les données (m, n, o, p, q)...

Merci pour l'aide.

Jacques

[__dardanos__]

Bonjour.

Comment faire varier uniquement les données (m, n, o, p, q)...

En utilisant un générateur et avec itertools.product, tu peux obtenir toutes les données d'entrées :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
import itertools
from numpy import arange
 
 
def get_one():
    items = arange(0.0, 1.0, 0.1)
    for x in itertools.product(items, repeat=5):
        if sum(x) == 1.0:
            yield x
 
 
if __name__ == "__main__":
    for i, x in enumerate(get_one()):
        print(i, *x)

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
0 0.0 0.0 0.0 0.1 0.9
1 0.0 0.0 0.0 0.2 0.8
2 0.0 0.0 0.0 0.3 0.7
3 0.0 0.0 0.0 0.4 0.6
4 0.0 0.0 0.0 0.5 0.5
5 0.0 0.0 0.0 0.6 0.4
6 0.0 0.0 0.0 0.7 0.3
7 0.0 0.0 0.0 0.8 0.2
8 0.0 0.0 0.0 0.9 0.1
9 0.0 0.0 0.1 0.0 0.9
10 0.0 0.0 0.1 0.1 0.8
11 0.0 0.0 0.1 0.2 0.7
12 0.0 0.0 0.1 0.3 0.6
13 0.0 0.0 0.1 0.4 0.5
14 0.0 0.0 0.1 0.5 0.4
15 0.0 0.0 0.1 0.6 0.3
16 0.0 0.0 0.1 0.7 0.2
17 0.0 0.0 0.1 0.8 0.1
18 0.0 0.0 0.1 0.9 0.0
19 0.0 0.0 0.2 0.0 0.8
20 0.0 0.0 0.2 0.1 0.7
...

[Jack_dev]

Bonjour Dardanos,

Merci pour le truc avec numpy. Ça semble si simple quand on s'y connait...

Maintenant, pour chaque ligne d'entrée, je vais calculer un "Z" qui est le résultat de certains calculs sur les autres données. Je veux m'assurer que le résultat final soit le plus grand de tous les "Z" calculés avec toutes les lignes d'entrée obtenues à partir de numpy: autrement dit, optimiser le "Z".

J'ai pensé obtenir le nombre de lignes fournies par numpy, l'utiliser pour faire tourner la boucle x fois, et toujours stocker le plus grand "Z"....sortir le résultat final.

Il y a probablement une solution plus simple...

Merci encore

[Jack_dev]

Bonjour,

J'ai testé la fonction et je reçois un message d'erreur de syntaxe sur la dernière ligne: print (i, *x).

J'utilise Python 2.7. Si j'enlève le * devant le x, ça fonctionne, mais le nombre de décimal est trop, comment fixer les décimales à 0.1?

Merci

Jacques

[__dardanos__]

En PYTHON2, print est une instruction ; en PYHTON3, c'est une fonction.
Le module __future__ permet d'importer la fonction en PYTHON2.

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
#!/usr/bin/python2
# -* coding: utf-8 -*-
from __future__ import print_function
import itertools
from numpy import arange
 
def get_one():
    items = arange(0.0, 1.0, 0.1)
    for x in itertools.product(items, repeat=5):
        if sum(x) == 1.0:
            yield x
 
if __name__ == "__main__":
    for i, x in enumerate(get_one()):
        print(i, *x)

Si tu débutes, pourquoi ne pas utiliser PYTHON3 ?

[Jack_dev]

Salut,

Pourquoi le 2 au lieu du 3, simplement parce que j'utilise python dans le cadre d'une formation en finance et que tout le reste du matériel est en 2. C'est assez compliqué pour quelqu'un qui n'est pas programmeur, et en plus, je dois jongler avec deux logiciels différents: R et python!!!

Jacques

[alexdevl]

Bonjour,

Je te propose ce code en 3 parties :

1. Calcul de listes de combinaisons pour 3 puis 5 paramètres
2. Création de listes en combinant les 3 et 5 paramètres
3. calcul du max de la fonction f avec toutes ces combinaisons

Je n'ai pris en compte que les éléments (ab,c) et (m,n,o,p,q)


Est ce ok ?

Alex

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
33
34
35
36
37
38
39
40
41
42
43
44
 
def f((a,b,c), (d, e, f), (m, n, o, p, q), (r, s, t, u, v)):
    return a+b+c+m+n+o+p+q
 
# Creation des listes de valeur ou la somme=1
# ------pour 3 parametres
valeurs_dim3=[] 
valeurs_dixieme=[0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0]
for a in valeurs_dixieme:
    for b in valeurs_dixieme:
        for c in valeurs_dixieme:
            if a+b+c==1:
                valeurs_dim3.append((a,b,c))
print "nombre de cas pour 3 parametres",len (valeurs_dim3)
# ------pour 5 parametres
valeurs_dim5=[] 
for a in valeurs_dixieme:
    for b in valeurs_dixieme:
        for c in valeurs_dixieme:
            for d in valeurs_dixieme:
                for e in valeurs_dixieme:
                    if a+b+c+d+e==1:
                        valeurs_dim5.append((a,b,c,d,e))
print "nombre de cas pour 5 parametres",len (valeurs_dim5)
 
# Creation de la liste en combinant les 3 parametre et 5 parametres
valeurs_total=[]
for p1 in valeurs_dim3:
    for p2 in valeurs_dim5:
        valeurs_total.append((p1,p2)) 
print "nombre de cas pour 3+5 parametres",len(valeurs_total)
 
# Recherche du maxi
max=0
for i,valeur in enumerate(valeurs_total):
    val_cal=f(valeur[0],(0.1,0.2,0.3),valeur[1],(0.1,0.2,0.3,0.4,0.5))
    if val_cal>max:
        max=val_cal
        item_max=i
        print max,valeurs_total[i]
print "---resultat---"
print "maxi=",max
print "Numero d'element=",item_max
print "parametre pour le max=",valeurs_total[item_max]