Matlab - Runge-Kutta de ordinul 5

#1
lupulalbastru1

Posted 31 March 2014 - 16:27

Junior Member

Grup: Members
Posts: 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.

Back to top

#2
flegma

Posted 31 March 2014 - 17:27

Senior Member

Grup: Senior Members
Posts: 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.

Back to top

#3
dan-s

Posted 31 March 2014 - 19:33

Active Member

Grup: Banned
Posts: 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).

Back to top

#4
soarce

Posted 01 April 2014 - 21:44

Active Member

Grup: Members
Posts: 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".

Back to top

#5
lupulalbastru1

Posted 02 April 2014 - 13:58

Junior Member

Grup: Members
Posts: 110
Înscris: 24.05.2013

Da, asta era. Multumesc @soarce

Back to top

Anunturi

Second Opinion

Folosind serviciul second opinion ne puteți trimite RMN-uri, CT -uri, angiografii, fișiere .pdf, documente medicale.

Astfel vă vom putea da o opinie neurochirurgicală, fără ca aceasta să poată înlocui un consult de specialitate. Răspunsurile vor fi date prin e-mail în cel mai scurt timp posibil (de obicei în mai putin de 24 de ore, dar nu mai mult de 48 de ore). Second opinion – Neurohope este un serviciu gratuit.

www.neurohope.ro

Matlab - Runge-Kutta de ordinul 5

#1 lupulalbastru1 Posted 31 March 2014 - 16:27

#2 flegma Posted 31 March 2014 - 17:27

#3 dan-s Posted 31 March 2014 - 19:33

#4 soarce Posted 01 April 2014 - 21:44

#5 lupulalbastru1 Posted 02 April 2014 - 13:58