Problème Résolution système équation différentiel en utilisant ODE45

Bonsoir à tous,
je souhaiterais résoudre sur Matlab ce système d'équation :

Pièce jointe 528783

J'ai essayé le code ci dessous :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30

function ypoint = ab(t,y) global mumax1 mumax2 ks1 ks2 ki2 k1 k2 k3 k4 ypoint(1) = (mumax1*y(3)*y(1)/(y(3)+ ks1)); ypoint(3) = -k1*(mumax1*y(3)*y(1)/(y(3)+ks1)); ypoint(2) = (mumax2*y(4)*y(2)/(y(4)+ks2+(y(4)^2/ki2))) ypoint(4) = k2*(mummax1*y(3)*y(1)/(y(3)+ks1))-k3*(mumax2*y(4)*y(2)/(y(4)+ks2+(y(4)^2/ki2))); ypoint(5) = k4*(mumax2*y(4)*Y(2)/(y(4)+ks2+(y(4)^2/ki2))); ypoint = ypoint(:); end   global mumax1 mumax2 ks1 ks2 ki2 k1 k2 k3 k4 q mumax1=0.4; mumax2= 0.4; ks1= 35; ks2=4; ki2= 170; k1= 50; k2= 50; k3= 15 k4= 75;% paramètres tfinal= 100; % temps final y01= 0.4; y02= 0.01; y03= 10; y04= 2; y05= 0;% conditions initiales [t,y]= ode45('ab',[0 tfinal],[y01 y02 y03 y04 y05])% résolution y1 = y(:,1);%extractions de Y1 Y2 Y3 Y4 et y5 y2 = y(:,2); y3 = y(:,3); y4 = y(:,4); y5 = y(:,5); q = k4*(mummax2*y4*y2/(y4+ks2+(y4^2/ki2))); plot(t,y1,'r') % y1 fonction de temps plot(t,y2,'g') % y2 fonction de temps plot(t,y3,'o') % y3 fonction de temps plot(t,y4,'a') % y4 fonction de temps plot(t,y1,t,y2,t,y3,t,y4)% y4 fonvtion de temps plot(t,y5)% y4 fonvtion de temps plot(t,q,'z'); plot(t,y5,'r',t,q,'p');

---

Mais je recois cette erreur :

>>Untitled2
Error : File : Untitled2.m Line : 10 column : 1
This statement is not inside any function.
(It follows the END that terminates the definition of the function "ab".)

Pouvez vous m'aider svp

Bonjour,

La définition de la fonction ab doit se trouver à la fin du script, et non au début.
Si tu n'es pas familier de la fonction ode45, je te conseille de commencer par tenter de résoudre une EDO plus simple, voire de te baser sur un example publié par Mathworks. e.g. https://fr.mathworks.com/help/matlab/ref/ode45.html
Note par ailleurs qu'il serait préférable d'éviter les variables globales. Mathworks propose diverses manière de procéder : https://fr.mathworks.com/help/matlab...5.html#bu3uhuk

Citation:

---

Envoyé par Antoane

Bonjour,

La définition de la fonction ab doit se trouver à la fin du script, et non au début.
Si tu n'es pas familier de la fonction ode45, je te conseille de commencer par tenter de résoudre une EDO plus simple, voire de te baser sur un example publié par Mathworks. e.g. https://fr.mathworks.com/help/matlab/ref/ode45.html
Note par ailleurs qu'il serait préférable d'éviter les variables globales. Mathworks propose diverses manière de procéder : https://fr.mathworks.com/help/matlab...5.html#bu3uhuk

---

bonsoir,
je voudrais que vous m'aider; résoudre cette équation différentielle à l'aide de runde kutta 4
voici le code de mon programme de la fonction de l'autre coté de l'égalité

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

function [g]= foc(t,y) gama=10; f=1500; alpha=0.01; Mo=14; V=1.4e-6; tg=alpha/gama;te=alpha/Mo*V^2; b=10e2; d=10e8; w=2*3.14*f; A=[0.5-3*d*V^2/16 0 -0.25+d*V^2/8 0 -d*V^2/32]; B=[0 -0.25*b*V 0 b*V/12 0]; % t= linspace(0,20,1000) ; s=0; e=V*sin(w*t); W=Mo*(0.5-e*b/3-d/4*e.^2).*e.^2; figure(1) plot(t,W) for n=2:5; s=s+A(1,n)*cos(n*w.*t)+B(1,n).*sin(n*w.*t); end r=0.5*A(1,1)+s; % g=W/(te*Mo*V^2)-y/tg; g=1/te*r-y/tg; end

---

voici mon programme runge kutta 4

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

function [y,t] = RK4(f,tmin,tmax,Nint,y0) % Méthode de Runge Kutta d’ordre 4 % Nint - nombre de sous intervalles % tmin - temps t0 % tmax - temps t0 + T % f est une fonction avec comme arguments t et y : f(t,y(t)) % y0 contient les valeurs des conditions limites h = (tmax-tmin)/Nint ; % valeur du pas (G) t = linspace(tmin,tmax,Nint+1) ; % vecteur de t discrétisé t=[tmin,tmax] n = length(t); y = zeros(1,n); % on fixe la taille du tableau y y(1,1) = y0 ; t(1,1)=0;% y contient les solutions de y(tn)n = 1, ..., Nint + 1 for n = 2:Nint-1 t(n+1)=t(n)+h; k1=h*f(t(n),y(n)); k2=h*f(t(n)+(h/2),y(n)+(k1/2)); k3=h*f(t(n)+(h/2),y(n)+(k2/2)); k4=h*f(t(n)+h,y(n)+k3); y(n+1)=y(n)+(k1/6)+(k2/3)+(k3/3)+(k4/6); end end

---

Pièce jointe 536226

Bonjour,

Et alors ?
Comment ce code est-il exécuté ? Qu'est ce que ça donne ? Quel éventuel message d'erreur ?
etc.