Discussion :

[bibed]

Bonjour,
je suis assez débutant en Matlab et je cherche à résoudre une équation différentielle du 2eme ordre qui est la suivante:

(1-x^2)y'' - x*y' + 49y =0

avec y'=dy/dx
y(x=0)=0
y'(x=0)=-7

et je dois tracer la courbe entre 0 et 1

Après avoir fait des recherches sur internet je tombe souvent sur la fonction ode23 mais j'arrive même pas à l'utiliser pour un premier ordre . Ma question est la suivante si quelqu'un peux m'expliquer avec un exemple concret comment fonctionne la fonction ode23. Ou bien existe il une autre fonction pour résoudre une équation differentielle du deuxième ordre.

Merci infiniment d'avoir prit le temps de me lire j'attends vos réponses avec impatience.

[Dam2227]

tu peux essayer ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
function bibed
 
[X,Y] = ode45(@bibed2,[0:0.01:1],[0 -7]);
 
plot(X,Y)
 
    function [dy]=bibed2(x,y)
        dy=zeros(2,1);
        dy(1)=y(2);                                     %y(1)=y     y(2)=y'     dy(1)=y'    dy(2)=y''
        dy(2)=(x*y(2)-49*y(1))/(1-x*x);
    end
end

Par contre, en 1, ta fonction diverge et devient infinie. Donc il faudra que tu t'arrêtes avant (0.99 dans ce cas)

++