[scipy] Séries mathématiques - fonction comme argument d'une autre

Version imprimable

03/03/2007, 11h45
sebsainte

[scipy] Séries mathématiques - fonction comme argument d'une autre

Bonjour,

Je suis amené à calculer des séries mathématiques pour des applications
scientifiques. (ce sont des sommes de termes qui dépendent d'un paramètre)

l'algo serait le suivant :

------------------------------------------------

somme = 0

alpha =150
beta=12
...

pour n de 0 à 100
somme = somme + f(n, alpha, beta, ...)
fin pour

Resultat = somme

------------------------------------------------

avec f une fonction qui dépend d'un entier n (entier sur lequel on
effectue la sommation) mais aussi d'autres paramètes qui eux sont fixés
avant la somme (alpha, beta, ...) .

je cherche à créer une fonction qui permet de réaliser directement la somme

ex :

Resultat = MaSuperSomme(f (n,alpha, beta, .... ), ParametreDeSomme = n, n=0, n=100)

Je suis en galère :aie: :

- pour passer la fonction f en paramètre de MaSuperSomme
- pour déterminer sur lequel des arguments de f on effecue la somme.

Merci pour vos idées !

Seb

Salut,

Voici une solution.
c'est certainement pas la meilleure, mais ca te permettra peut être de te débloquer...

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
def toto(n,alpha,beta) :
    return n+alpha+beta
 
def MaSuperSomme(f,args,parametre,debut,fin):
    somme = 0
    for i in range(debut,fin) :
        args[parametre]=i
        fct = f.__name__+"("
        for cle,valeur in args.items(): fct+=cle+"="+str(valeur)+","
        fct=fct[:-1] # pour virer la virgule
        fct += ")"
        exec "somme+="+fct
    return somme
 
 
print MaSuperSomme(toto,{'n':0,'alpha':3,'beta':4},parametre='n',debut=0,fin=100)

Si tu as des questions sur le code, je suis là :yaisse2:

See++

03/03/2007, 13h23
sebsainte

Bonjour see++

La solution que tu proposes est vraiment intéressante. Elle ouvre
un certain nombre d'applicatons sur les manipulations génériques
de fonctions.

Merci pour l'effort d'avoir construit un exemple fonctionnel

Le script est très clair

Merci et à bientôt sur le forum !

++

Seb
03/03/2007, 14h43
see++

De rien, ca fait toujours plaisir de pouvoir aider quelqu'un

A+ sur le forum

See++
05/03/2007, 11h28
Matthieu Brucher
On devrait pouvoir s'arranger sans passer par le exec.
Si tu passes ta fonction en argument, tu peux directement l'appeler, sans pb. Le réel problème reste les arguments.
Pour cela, tu peux faire qqch du style :
Code:

1 2 3 4 5 def MaSuperFonction(f, debut, fin, **parameters): somme = 0 # for i in range(debut, fin): somme += f(parameters) return somme
l'intérêt du **parameters est que c'est un dictionnaire contenant toutes les paramètres passés à la fonction qui ne font pas partie des attributs définis explicitements. Ensuite, il te reste à faire un :

Code:

print MaSuperFonction(f = maFonction, debut = 0, fin = 10, alpha = 30, beta = 45, gamma = 56, ...)

On peut aussi se servir des fonctions sum et map:
(y'a peut-être moyen d'être plus générique)

exemple:

Code:

1
2
3
4
5
6
7
8
9
alpha =150
beta=12
 
def somme(f, debut, fin, a, b):
    return sum(map(f, range(debut, fin+1), (fin-debut+1) * [a], (fin-debut+1) * [b]))
 
def f(n, a, b): return n*a+b
 
print somme(f, 0, 10, a = alpha, b = beta)

Edit: finalement avec les listes compréhensions, ce sera plus générique:

Code:

1
2
3
4
5
6
7
8
9
alpha =150
beta=12
 
def somme(f, debut, fin, **param):
    return sum([f(i, **param) for i in range(debut, fin+1)])
 
def f(n, a, b): return n*a+b
 
print somme(f, 0, 10, a = alpha, b = beta)

20/09/2007, 00h59
SRT__

Bonjour, (j'ai changé de paramètres de connexion... anciennement sebsainte)
Je me suis finalement rapproché de fonctions présentes dans la librairie scipy,
qui procède globalement comme les deux dernières réponses (notamment pour les fonctions d'intégration)

le code est très proche de ces deux dernières

Ca marche très bien !

Merci pour vos réponses