[Users-fr] Utilisation de la fonction impl

[Calixte Denizet]

On 16/10/2012 15:55, LE BELGUET Alix 226276 wrote:
>
> Bonjour,
>
> Je souhaite utiliser la fonction 'impl' de Scilab afin de résoude un 
> système d'équations différentielles linéaire du type : A(t,y) 
> dy/dt=g(t,y), y(t0)=y0.
>
> Cependant, je n'arrive pas à comprendre ce que signifie l'argument 
>  'adda' définie par : r=A(t,y)+p.
>
> Comment choisir p ?
>

Bonjour Alix,

Par exemple, vous cherchez a resoudre en dimension 2:
Id * y' = y avec y(0) = [1;2]

function r=res(t,y,ydot)
r = y-eye(2,2)*ydot;
endfunction

function r=adda(t,y,p)
r=eye(2,2)+p
endfunction

puis:
t=0:0.1:1;
y=impl([1;2],[0;0],0,t,res,adda)

On verifie que y vaut a peu pres [1;2]*exp(t)
max(abs(y-[1;2]*exp(t)))

retourne chez moi 1.3e-6.

Dans le code ci-dessus, il suffit de remplacer eye(2,2) par A(t,y) et 
dans res de remplacer y par g(t,y). Vous n'avez pas a vous souciez de ce 
que vaut p...

Calixte

> Je ne parviens pas non plus à trouver d'exemple parlant de cette 
> fonction, c'est pourquoi je me tourne vers vous...
>
> En vous remerciant par avance pour votre réponse,
>
> Cordialement,
>
> Alix LB
>
>
>
> _______________________________________________
> users-fr mailing list
> users-fr at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users-fr


-- 
Calixte Denizet
Software Development Engineer
-----------------------------------------------------------
Scilab Enterprises
143bis rue Yves Le Coz - 78000 Versailles, France
http://www.scilab-enterprises.com

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