Matlab - Runge-Kutta de ordinul 5

Ultima postare: apr 02 2014 13:58, Iniţiat de lupulalbastru1 , mar 31 2014 16:27

#1
lupulalbastru1

Publicat: 31 martie 2014 - 16:27

Junior Member

Grup: Members
Mesaje: 110
Înscris: 24.05.2013

Buna!
Am o intrebare, la o problema in Matlab, daca poate cineva sa ma ajute, va rog.
Să se rezolve următoarea problemă Cauchy, folosind funcţia Matlab care se bazează şi pe metoda Runge-Kutta de ordinul 5:
2⋅y "'+3⋅y = t2⋅et,
y(2)=y '(2)=y "(2)=-1, t∈[2,4]

Codul l-am scris asa:

function dy=ecdif(t,y)
dy=(t.^2.*exp(t)-3*y)./2;
%conditia initiala
y0=-1;
%domeniul
dom=[2,4];
%rezolvarea ecuatiei diferentiale
[xval,yval]=ode45('ecdif',dom,y0)
%reprezentarea grafica a solutiei
plot(xval,yval)

Dar cand ii dau RUN imi da eroarea:

??? Input argument "t" is undefined.
Error in ==> ecdif at 2
dy=(t.^2.*exp(t)-3*y)./2;

Nu stiu de ce da eroarea asta si cum sa o rezolv.

#2
flegma

Publicat: 31 martie 2014 - 17:27

Senior Member

Grup: Senior Members
Mesaje: 2.294
Înscris: 09.08.2009

Pai nu ai definit t-ul. Cel putin, asa suna mesajul de eroare. Si nici nu il vad in codul tau.

#3
dan-s

Publicat: 31 martie 2014 - 19:33

Active Member

Grup: Banned
Mesaje: 1.699
Înscris: 02.07.2013

Nu a definit nici y. Din câte știu nu e nevoie să dai run la o funcție. Ci o apelezi direct (eventual te asiguri că fișierul este adăugat în calea în care caută MATLAB funcțiile).

#4
soarce

Publicat: 01 aprilie 2014 - 21:44

Active Member

Grup: Members
Mesaje: 1.467
Înscris: 02.10.2008

Functia o definesti intr-un fisier separat, apoi o chemi co ode45. Uitate pe exemplele din matlab, in fereastra de comezi ii dai "help ode45".

#5
lupulalbastru1

Publicat: 02 aprilie 2014 - 13:58

Junior Member

Grup: Members
Mesaje: 110
Înscris: 24.05.2013

Da, asta era. Multumesc @soarce

Anunturi

Bun venit pe Forumul Softpedia!

Matlab - Runge-Kutta de ordinul 5

#1 lupulalbastru1 Publicat: 31 martie 2014 - 16:27

#2 flegma Publicat: 31 martie 2014 - 17:27

#3 dan-s Publicat: 31 martie 2014 - 19:33

#4 soarce Publicat: 01 aprilie 2014 - 21:44

#5 lupulalbastru1 Publicat: 02 aprilie 2014 - 13:58