Re: [Debian]:DGL's
Andreas Mueller <andrmuel@uos.de> writes:
> Mein Problem sieht so aus:
>
> u' = Au + b
>
> Ich dachte da an ein Tool, dem ich A und b übergebe und welches mir
> u ausgibt.
>
> Ich brauche kein Grafisches Frontend oder ähnliches ist alles sekundär.
Also ein Anfangswertproblem.
Tja, wenn u, A und b skalar sind, würde ich folgende Lösung
vorschlagen:
mc:axeld@rts-pc7> maxima
GCL (GNU Common Lisp) Version(2.2.1) Wed Apr 7 23:15:26 CEST 1999
Licensed under GNU Public Library License
Contains Enhancements by W. Schelter
Maxima 5.2 Wed Apr 7 22:27:56 CEST 1999 (with enhancements by W. Schelter).
(C1) 'diff(u,t)=A*u+b;
dU
(D1) -- = A U + B
dT
(C2) ode2(c1,u,t);
- A T
B %E A T
(D2) U = (%C - ---------) %E
A
d.h. u(t) = ( C - b/A*exp(-A*t) )*exp(A*t)
Nun brauchen wir nur noch den Anfangswert u0 = u(t=0), um die
unbekannte Konstante "C" zu bestimmen:
u(t=0) = u0 = (c-b/A*exp(-A*0))*exp(A*0)
= (c-b/A*1)*1
= c-b/A
=> c = u0 + b/a
Das ganze jetzt noch in C zu programmieren,
dürfte kein großes Problem sein.
MfG,
AxelD
--
Axel Dürrbaum / Universität Gh Kassel / FB 15 - RTS Regelungstechnik
Mönchebergstraße 7 / 34109 Kassel / Germany / Technik I/II / Raum 2510
phone:+49 561 804 3261 Email:axeld@rts.maschinenbau.uni-kassel.de
------------------------------------------------
Um sich aus der Liste auszutragen schicken Sie
bitte eine E-Mail an majordomo@jfl.de die im Body
"unsubscribe debian-user-de <deine emailadresse>"
enthaelt.
Bei Problemen bitte eine Mail an: Jan.Otto@jfl.de
------------------------------------------------
Anzahl der eingetragenen Mitglieder: 773
Reply to: