Discussion :

[aymenvictoire]

Bonjour à toutes et à tous,
J'ai créé une fonction pour résoudre un système en fonction de temps avec la fonction ode45, elle fonctionne mais j'aimerai que un paramétré varie lors de exécution.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
function dy=famorti(t,y)
k=2
m=1
w0=sqrt(k/m)
x0=1
dy(1)=y(2);
dy(2)=x0*w0*cos(w*t)-(w0^2)*y(1);

code d'exécution

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
[t,y] = ode45('famorti',[0 5],[0 0]);
y1=y(:,1);
y2=y(:,2);

bon, le code fonction si je le fonctionne avec w=constante, mais j'aimerai que w soit variable de 1 à 10 par exemple et que je peut présente mes résultats en fonction de w.
merci d’avance pour votre aide.

[Invité]

Bonsoir,

En utilisant le symbole @ :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
function dy=famorti(t,y,w)
k=2;
m=1;
w0=sqrt(k/m);
x0=1;
dy(1,1)=y(2);
dy(2,1)=x0*w0*cos(w*t)-(w0^2)*y(1);

(attention la fonction famorti doit retourner un vecteur colonne)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
w = ...
[t,y] = ode45(@(t,y) famorti(t,y,w),[0 5],[0 0]);
y1=y(:,1);
y2=y(:,2);

[aymenvictoire]

merci Winjerome,
j’essayai avec votre proposition, mais pas de résultats ;

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
??? Subscripted assignment dimension mismatch.
 
Error in ==> famorti at 10
dy(2,1)=X0*w0*cos(w*t)-(w0^2)*y(1);
 
Error in ==> @(t,y)famorti(t,y,w)
 
 
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
 
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

[Invité]

Le w rentré doit être un scalaire. S'il doit y avoir plusieurs valeurs, elles doivent être rentrées une par une dans une boucle.