Répondre à la discussion
Affichage des résultats 1 à 3 sur 3

Matlab equadiff



  1. #1
    maki00

    Matlab equadiff


    ------

    Bonjour, je suis embêté par matlab qui me renvoit une erreur bizarre:
    je cherche à résoudre:
    y1''=-2*y1'/x+2*y1/x²+sin(ln(x))/x² pour x entre 1 et 2
    avec les conditions: y1(1)=1 ; y1'(1)=0

    voilà les codes:
    clear all;
    xspan=[1 2];
    y1ini=[1 0]';
    y2ini=[0 1]';
    [x,y1]=ode45(@func1,xspan,y1ini);


    function ydot1=func1(x,y1)
    ydot(1)=y1(2);
    ydot(2)=-2*y1(2)/x+2*y1(1)/x^2+sin(log(x))/x^2;
    ydot=ydot(: );



    et ça me renvoie:
    Error using ==> C:\matlab\toolbox\matlab\funfu n\private\odearguments
    Solving FUNC1 requires an initial condition vector of length 0.



    -----

  2. #2
    ChromoMaxwell

    Re : Matlab equadiff

    Je vois pas l'erreur. Peut-être que le nom de ta fonction est ydot et non ydot1. J'ai des doutes sur le fait que xspan soit en ligne et non en colonne. Bref, c'est de la syntaxe.

    Je te propose le code suivant, qui fonctionne apparemment.




    f=inline('[x(2) ; -2.*x(2)./(t+2)+2.*x(1)./t^2+sin(log(t))/t^2 ]','t','x');

    [tsol,ysol]=ode45(f,[1 2]',[1 0]');

    plot(tsol,ysol)
    figure
    a=ysol(:,1);
    b=ysol(:,2);
    plot(a,b)

  3. #3
    maki00

    Re : Matlab equadiff

    ok merci pour l'info; problème, jamais vu inline en cours, ca risque de ne pas passer dans un contrôle...
    enfin je vais faire avec ça parce qu'au moins là ca marche

Discussions similaires

  1. Equadiff et écoulement
    Par alias_sg1 dans le forum Physique
    Réponses: 6
    Dernier message: 23/07/2009, 00h17
  2. Résoudre rigoureusement une équadiff
    Par Bleyblue dans le forum Mathématiques du supérieur
    Réponses: 23
    Dernier message: 26/02/2007, 19h43
  3. EquaDiff, comment faire ?
    Par TuTu57 dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 26/11/2006, 09h34
  4. equadiff
    Par bastien8312 dans le forum Mathématiques du collège et du lycée
    Réponses: 8
    Dernier message: 28/10/2006, 11h11
  5. Resolution equadiff
    Par juanita dans le forum Mathématiques du supérieur
    Réponses: 5
    Dernier message: 24/09/2005, 20h23