DGL Lösung < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Aufgabe | Deine Aufgabe ist es, die Bewegung [mm] \theta(t) [/mm] für gegebene Werte L,A, w zu programmieren, mit [mm] \theta(0) [/mm] und [mm] \theta'(0) [/mm] zu initialisieren. Nutze dazu die Matlab ODE45.
Die Funktion lautet: [mm] \theta''=3/2L*(g-A*w^2*\sin(wt))*\sin(\theta) [/mm] |
Wenn ich es richtig gesehen habe, arbeitet ode45 mit DGLs 1. Ordnung, hier habe ich aber eien DGL 2. Ordnung. Wie wandle ich diese um ? Und wie nutze ich diese ode45, bitte gebt mir einen Ansatz, danke.
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 17:10 Mi 11.07.2007 | | Autor: | leduart |
Hallo
eine Dgl 2ter Ordnung kann man immer in ein System von Dgl. erster Ordnung umschreiben:
y''=f(x,y)
y1=y
y2=y'
dann hast du das System:
y1'=y2
y2'=f(x,y1)
mit y1(0) und y2(0) als Anfangswerte.
in Ode45 musst dus dann selbst einsetzen.nutz doch den Kommentar der routine !
Gruss leduart
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 18:28 Mi 11.07.2007 | | Autor: | Phlipper |
Also das mit der Ode45 Funktion, das bekomme ich dann schon hin. Mein Problem ist wirklich diese Fkt. korrekt umzuformen.
Habe es jetzt so:
function dteta = teta(t,teta)
A=0.5;
w=5.3;
dteta = zeros(2,1);
dteta (2) = 3/(2*10)*(386.09 - [mm] A*w^2*sin(w*t))*sin(teta(1));
[/mm]
dteta(1) = teta(2);
Ich glaube aber, dass da noch ein kleiner Fehler drin ist.
Würde mich riesig über einen Hinweis freuen. Philipp
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 12:20 Sa 14.07.2007 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
