Bonjour à tous,

Je suis débutant en Matlab, que j'utilise pour le compte de mes études.
J'ai un soucis que je n'arrive pas à résoudre malgré de multiples recherches.
J'ai simplifié mon programme de base pour que ce soit davantage compréhensible.

Je souhaite résoudre une équation différentielle avec un paramètre qui varie périodiquement, et ça j'y arrive. Le soucis c'est que j'aimerai pouvoir extraire le paramètre qui varie périodiquement de la fonction dans laquelle il se trouve.

Voici le programme principale :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
duree=[0 5];
a=0.5;
[t,y] = ode45(@(t,y)testfonction(t,y,a), duree, [0 0.01]);
plot(t, y(:,1), t, y(:,2))


et la fonction :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
function dydt = testfonction(t,y,a)
b=a*cos(2*pi*(1/0.5)*t)+a;
dydt = zeros(2,1);
dydt(1) = y(2);
dydt(2) = b.*t.*y(1);
end

Finalement, ce que je souhaite faire, c'est tracer b en fonction de t. Je ne vais pas étaler tout ce que j'ai essayé de faire, ce serait trop long et puis de toute façon ça n'a pas marché...

Donc voilà, si quelqu'un a une idée de comment s'y prendre je suis preneur !
Merci à vous