Bonjour,
Je travaille sur de l'optimisation de modèle numérique à partir de Scilab. L'un des problèmes concerne la résolution d'équations différentielles (ED).
Pour certaines ED le vecteur y solution n'a pas la même dimension que l'intervalle de temps sur lequel l'ED doit être résolue.
Ex :
clear
t=0:0.1:%pi;
y0=0;t0=0;
function ydot=f(t,y), ydot=(2*y)^2-y*sin(t)+cos(-0.5*t), endfunction
y=ode(y0,t0,t,f);
Y=y'
pour cette ED, dim(Y)=9 alors que dim(t)=32 !
Or si je donne cette équation
function ydot=f(t,y), ydot=(1*y)^2-y*sin(t)+cos(1*t), endfunction
Y et t ont la même dimension.
Je crois que l'explication vient du fait que la résolution de la 1ère ED est un problème raide... et c'est pourquoi je n'obtiens pas toutes les valeurs de Y correspondant à t. Pourriez vous m'aidez à récupérer toutes ses valeurs ?
Autre question : Est-ce que les premières valeurs de Y correspondent aux 1ères valeurs de t ?
Merci beaucoup de votre aide.
-----