[Users-fr] Suite de fonctions définies par récurrence

babigeon babigeon at lal.in2p3.fr
Mar 27 Aou 08:14:31 CEST 2013


Bonjour,

Quelques remarques :

- Si l'idée est effectivement de résoudre l'équation différentielle, je 
n'ai pas entièrement vérifié le script mais
j'ai tendance à faire confiance au moins à la méthode employée par l'auteur
- Si je lis de plus près, "l'idée est de générer une suite de 
fonctions", il me semble que l'objectif ne soit pas entièrement ceci. 
Est-ce vrai ?
- point mineur, mais peut-être pas puisqu'il intervient dans une e.d., 
La fonction S n'apparaît pas très bien "conditionnée". J'aurai plutôt 
écrit S comme la limite de s(x,y,t)
quand y tend vers 1 par exemple,
- Sinon, je suis beaucoup plus sensible à l'esthétique de phi et H (pour 
moi ça compte dans un modèle)
- Je suppose que vous avez exploité la symétrie en r1+r2 ? Si on le 
fait, on trouve une équation en log (phi) en fonction de S et H. Mais je 
n'ai pas trouvé d'autre combinaison simple de r1 et r2 permettant de 
transformer le problème.
- Pour moi, générer une suite de fonctions, peut consister à développer 
r1 et r2 en série de Taylor par rapport à la variable t, et injecter les 
coefficients dans les équations différentielles. On obtient alors une 
suite de fonctions dont les termes pour r1(n) et r2(n), sont les 
approximations à l'ordre n de chaque série. Par contre, il faut faire au 
préalable un développement de S à l'ordre 1 ou 2 par rapport à la 
variable t, en injectant les développements de Taylor de r1 et r2 dans S.

JLuc



Le 26/08/2013 21:35, Serge Steer a écrit :
>
> Je ne suis pas sur de comprendre ce que vous voulez faire, mais s'il 
> s'agit d'integrer le systeme d'eqations differentielles
> vous pouver le faire avec le code Scilab suivant
>
> function v=f(y,r,D,Beta,x)
>   v=exp(-(x-y)^2/D^2)*r(1)/(r(1)^Beta+r(2)^Beta)
> endfunction
>
> function v=dr(y,r,alpha,D,Beta,x)
>   v= alpha*intg(0,1,list(f,r,D,Beta,x))
> endfunction
>
> r0=[1;1];
> y0=1
> Beta=2;
> x=0;
> alpha=[1;2];
> D=1.5;
> r1=ode(r0,0,0:0.01:1,list(dr,alpha,D,Beta,x))
>
>
> Serge Steer
> INRIA
> Le 24/08/2013 23:26, pro.olive at ouvaton.org a écrit :
>> Bonjour a tous,
>>
>> Je mets en fichier joint le problème que je rencontre au sujet de 
>> Scilab : l'idée est de pouvoir générer une suite de fonctions. Je ne 
>> sais donc pas si Scilab est le meilleur logiciel qui soit...
>>
>> Merci.
>> Marc Olive
>>
>>
>> _______________________________________________
>> users-fr mailing list
>> users-fr at lists.scilab.org
>> http://lists.scilab.org/mailman/listinfo/users-fr
>
>
>
> _______________________________________________
> users-fr mailing list
> users-fr at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users-fr

-- 
repons



Jean-Luc BABIGEON
CNRS/Lal/Dacc
Laboratoire de L'accélérateur Linéaire
Centre scientifique d'Orsay - Bat 200 - BP 34
91898 - Orsay Cedex
tel 0164468247 ou 8427 -babigeon at lal.in2p3.fr 
<mailto:babigeon at lal.in2p3.fr>

-------------- section suivante --------------
Une pièce jointe HTML a été nettoyée...
URL: <https://lists.scilab.org/pipermail/users-fr/attachments/20130827/2194a4a6/attachment.htm>


Plus d'informations sur la liste de diffusion users-fr