maximum sur maple
bonjour,
j'ai un petit probleme,
j'utilise maple, j'ai une fonction un peu complique, et je voudrais savoir son maximum approximatif (enivron 0.0001pres)
j'ai essayer maximize mais la fonction est trop complexe meme apres 1heure... cependant j'arrive sans peobleme a avoir le graphique..
y aurait-il un methode (numerique peut etre, pour la rapidite) pour trouver la valeur de ce maximum..
merci
heu... c'est quoi cette fonction ?
Shokin
Pardon, humilité, humour, hasard, tolérance, partage, curiosité et diversité => liberté et sérénité.
je ne sais pas si ca nous avance beaucoup de connaitre la fonction, c'est la procedure, je pense, qui manque...Envoyé par shokin
la fonction represente un comportement normal centre sur le lever du jour que l'on observerait sans tenir compte des changement d'heure du lever du soleil au cours de l'annee...
voici la fonction:
g:=x->1/n*sum(1/(sd*sqrt(2*Pi))*exp(-0.5*((x-((12-((arccos((-0.01454-sin((arcsin(0.3978*sin((280*ra d+(1.914*rad*sin((357*rad+0.98 56*rad*j))+0.02*rad*sin(2*(357 *rad+0.9856*rad*j)))+0.9856*ra d*j)))))*sin(la*rad))/(cos((arcsin(0.3978*sin((280*r ad+(1.914*rad*sin((357*rad+0.9 856*rad*j))+0.02*rad*sin(2*(35 7*rad+0.9856*rad*j)))+0.9856*r ad*j)))))*cos(la*rad))))/rad)/15+(4*((1.914*rad*sin((357*rad +0.9856*rad*j))+0.02*rad*sin(2 *(357*rad+0.9856*rad*j)))+(-2.466*rad*sin(2*(280*rad+(1.91 4*rad*sin((357*rad+0.9856*rad* j))+0.02*rad*sin(2*(357*rad+0. 9856*rad*j)))+0.9856*rad*j))+0 .053*sin(4*(280*rad+(1.914*rad *sin((357*rad+0.9856*rad*j))+0 .02*rad*sin(2*(357*rad+0.9856* rad*j)))+0.9856*rad*j))))))))/sd)^2),j=1..n):
avec n,rad,la,sd connu..
merci.
peut etre que tu pourrais faire ton programme toi meme ???
Arg... gloup... ya pas moyen de simplifier cette fonction ?
d'habitude, pour trouver le maximum d'une fonction, je dérive celle-ci. Je calcule le x tel que f'(x)=0. (x;f(x)) représente alors un extrêmum de la fonction. Pour m'assurer de cet extrêmum, je calcule f''(x), si f''(x)>0, j'ai affaire un minimum, un creux, si f''(x)<0, j'ai affaire à un maximum, à une bosse.
Heu... pour dériver ta fonction... ça doit pas être compliqué, juste un peu long (les fautes de calculs sont probables)...
Si tu sais dériver sin, cos, exp, ax^n, le produit de deux fonctions, le quotient de deux fonctions, la fonction de fonction, pas de problème !
Ya juste ton j qui semble pouvoir prendre plusieurs valeurs...
Shokin
Pardon, humilité, humour, hasard, tolérance, partage, curiosité et diversité => liberté et sérénité.
j'avais pense a ca, f(x)'=0
"rien de plus simple que de resoudre cette equation pour maple.."
cependant, c'est pareil.. l'ordinateur me prend des lustres... et j'aurais besoin de repeter le calcul...
le probleme de faire ce calcul maluellement est qu'il comprend une somme.., et cette somme va differer selon les different parametre.. le plus simple donc serait une approche numerique de la chose je pense..
mais maple a t il une telle fonction?
Désolé, mais je ne connais pas du tout maple.
Shokin
Pardon, humilité, humour, hasard, tolérance, partage, curiosité et diversité => liberté et sérénité.
bon c'est pas grave
je crois que j'ai la solution, c'est barbare mais je vais programmer une procedure de maximum a differente echelle de precision pour resserrer l'interval et ameliorer la precision..
merci encore..
