www.matheraum.de
Das Matheforum.
Das Matheforum des MatheRaum.

Für Schüler, Studenten, Lehrer, Mathematik-Interessierte.
Hallo Gast!einloggen | registrieren ]
Startseite · Forum · Wissen · Kurse · Mitglieder · Team · Impressum
Forenbaum
^ Forenbaum
Status Mathe
  Status Schulmathe
    Status Primarstufe
    Status Mathe Klassen 5-7
    Status Mathe Klassen 8-10
    Status Oberstufenmathe
    Status Mathe-Wettbewerbe
    Status Sonstiges
  Status Hochschulmathe
    Status Uni-Analysis
    Status Uni-Lin. Algebra
    Status Algebra+Zahlentheo.
    Status Diskrete Mathematik
    Status Fachdidaktik
    Status Finanz+Versicherung
    Status Logik+Mengenlehre
    Status Numerik
    Status Uni-Stochastik
    Status Topologie+Geometrie
    Status Uni-Sonstiges
  Status Mathe-Vorkurse
    Status Organisatorisches
    Status Schule
    Status Universität
  Status Mathe-Software
    Status Derive
    Status DynaGeo
    Status FunkyPlot
    Status GeoGebra
    Status LaTeX
    Status Maple
    Status MathCad
    Status Mathematica
    Status Matlab
    Status Maxima
    Status MuPad
    Status Taschenrechner

Gezeigt werden alle Foren bis zur Tiefe 2

Navigation
 Startseite...
 Neuerdings beta neu
 Forum...
 vorwissen...
 vorkurse...
 Werkzeuge...
 Nachhilfevermittlung beta...
 Online-Spiele beta
 Suchen
 Verein...
 Impressum
Das Projekt
Server und Internetanbindung werden durch Spenden finanziert.
Organisiert wird das Projekt von unserem Koordinatorenteam.
Hunderte Mitglieder helfen ehrenamtlich in unseren moderierten Foren.
Anbieter der Seite ist der gemeinnützige Verein "Vorhilfe.de e.V.".
Partnerseiten
Mathe-Seiten:Weitere Fächer:

Open Source FunktionenplotterFunkyPlot: Kostenloser und quelloffener Funktionenplotter für Linux und andere Betriebssysteme
StartseiteMatheForenGewöhnliche DifferentialgleichungenLösung eines lin. DGL-Systems
Foren für weitere Schulfächer findest Du auf www.vorhilfe.de z.B. Philosophie • Religion • Kunst • Musik • Sport • Pädagogik
Forum "Gewöhnliche Differentialgleichungen" - Lösung eines lin. DGL-Systems
Lösung eines lin. DGL-Systems < gewöhnliche < Differentialgl. < Analysis < Hochschule < Mathe < Vorhilfe
Ansicht: [ geschachtelt ] | ^ Forum "Gewöhnliche Differentialgleichungen"  | ^^ Alle Foren  | ^ Forenbaum  | Materialien

Lösung eines lin. DGL-Systems: Aufgabenhilfe
Status: (Frage) beantwortet Status 
Datum: 21:16 Do 07.05.2009
Autor: Ultio

Aufgabe
Bestimmen Sie die Lösung des Systems von linearen Differentialgleichungen
x'=2x+y
y'=2y
zu den Anfangswerten x(0)=1, y(0)=3 indem sie exp(t [mm] \pmat{ 2 & 1 \\ 0 & 2 }) [/mm] berechnen.

Hi,
Ich wollte mal fragen, ob jemand mir bitte weiterhelfen kann und soweit wie ich es fertig habe kontrollieren könnte.
Vielen Dank im Voraus.

Die Matrix folgt aus den Gleichungen. Dann folgt folgende Überlegung:
[mm] \vektor{x'(t) \\ y('t)}= \pmat{ 2 & 1 \\ 0 & 2 } [/mm] * [mm] \vektor{x(0) \\ y(0)} [/mm] = [mm] \pmat{ 2 & 1 \\ 0 & 2 } [/mm] * [mm] \vektor{1 \\ 3} [/mm]
daraus folgt:
[mm] \vektor{x \\ y} [/mm] (t) = exp(t [mm] \pmat{ 2 & 1 \\ 0 & 2 }) \vektor{x(0) \\ y(0)} [/mm] =exp(t [mm] \pmat{ 2 & 1 \\ 0 & 2 }) \vektor{1 \\ 3} [/mm]
= [mm] \pmat{ e^{2t} & e^{t} \\ 0 & e^{2t} } [/mm]
(doch was sagt mir diese Matrix?Denn damit sollten wir ja rechnen laut Aufgabe)

[mm] det(\pmat{ 2 & 1 \\ 0 & 2} [/mm] - [mm] \lambda \pmat{ 1 & 0 \\ 0 & 1 }) [/mm] = [mm] \pmat{ 2-\lambda & 1 \\ 0 & 2-\lambda} [/mm] = [mm] (2-\lambda)^{2} [/mm] = [mm] \lambda^{2} [/mm] - [mm] 4\lambda [/mm] +4
[mm] \Rightarrow \lambda_{1,2} [/mm] = 2  [mm] \pm \sqrt{4-4} [/mm]
[mm] \Rightarrow \lambda_{1} [/mm] = 2 (Eigenwert)

so jetzt bin ich mir absolut nicht mehr sicher:
[mm] exp(\pmat{ 2 & 0 \\ 0 & 2}) [/mm] = [mm] \pmat{ e^{2t} & 0 \\ 0 & e^{2t}}ist [/mm] dann was? Ist das meine Lösung? Oder fehlt da noch was?



Dankeschön.

        
Bezug
Lösung eines lin. DGL-Systems: Antwort
Status: (Antwort) fertig Status 
Datum: 22:31 Do 07.05.2009
Autor: MathePower

Hallo Ultio,

> Bestimmen Sie die Lösung des Systems von linearen
> Differentialgleichungen
> x'=2x+y
>  y'=2y
>  zu den Anfangswerten x(0)=1, y(0)=3 indem sie exp(t [mm]\pmat{ 2 & 1 \\ 0 & 2 })[/mm]
> berechnen.
>  Hi,
> Ich wollte mal fragen, ob jemand mir bitte weiterhelfen
> kann und soweit wie ich es fertig habe kontrollieren
> könnte.
>  Vielen Dank im Voraus.
>  
> Die Matrix folgt aus den Gleichungen. Dann folgt folgende
> Überlegung:
>  [mm]\vektor{x'(t) \\ y('t)}= \pmat{ 2 & 1 \\ 0 & 2 }[/mm] *
> [mm]\vektor{x(0) \\ y(0)}[/mm] = [mm]\pmat{ 2 & 1 \\ 0 & 2 }[/mm] * [mm]\vektor{1 \\ 3}[/mm]
> daraus folgt:
>  [mm]\vektor{x \\ y}[/mm] (t) = exp(t [mm]\pmat{ 2 & 1 \\ 0 & 2 }) \vektor{x(0) \\ y(0)}[/mm]
> =exp(t [mm]\pmat{ 2 & 1 \\ 0 & 2 }) \vektor{1 \\ 3}[/mm]
>  = [mm]\pmat{ e^{2t} & e^{t} \\ 0 & e^{2t} }[/mm]


Das stimmt nicht ganz.

Berechne hier sämtliche Potenzen von A

[mm]e^{t*A}:=\summe_{k=0}^{\infty}\bruch{t^{k}}{k!}*A^{k}[/mm]

In der Regel reichen hier die ersten Glieder um
eine Systematik dahinter zu erkennen.


>  
> (doch was sagt mir diese Matrix?Denn damit sollten wir ja
> rechnen laut Aufgabe)


Die Matrix [mm]e^{t*A}[/mm] ist dann die Fundamentallösung dieses DGL-Systems.


>  
> [mm]det(\pmat{ 2 & 1 \\ 0 & 2}[/mm] - [mm]\lambda \pmat{ 1 & 0 \\ 0 & 1 })[/mm]
> = [mm]\pmat{ 2-\lambda & 1 \\ 0 & 2-\lambda}[/mm] = [mm](2-\lambda)^{2}[/mm]
> = [mm]\lambda^{2}[/mm] - [mm]4\lambda[/mm] +4
>  [mm]\Rightarrow \lambda_{1,2}[/mm] = 2  [mm]\pm \sqrt{4-4}[/mm]
>  [mm]\Rightarrow \lambda_{1}[/mm]
> = 2 (Eigenwert)


Zunächst benötigst Du einen Eigenvektor zum Eigenwert 2.

Dieser ergibt sich aus dem Gleichungssystem

[mm]\pmat{0 & 1 \\ 0 & 0}*ev_{1}=\pmat{0 \\ 0}[/mm]


Hier brauchst Du noch einen sogenanten Eigenvektor 2. Stufe.

Dieser ist Lösung von

[mm]\pmat{0 & 1 \\ 0 & 0}*ev_{2}=ev_{1}[/mm]

Diese Gleichung ergibt sich unmittelbar aus

[mm]\pmat{0 & 1 \\ 0 & 0}^{2}*ev_{1}=\pmat{0 \\ 0}[/mm]


>  
> so jetzt bin ich mir absolut nicht mehr sicher:
>  [mm]exp(\pmat{ 2 & 0 \\ 0 & 2})[/mm] = [mm]\pmat{ e^{2t} & 0 \\ 0 & e^{2t}}ist[/mm]
> dann was? Ist das meine Lösung? Oder fehlt da noch was?
>  


Wenn, wie hier. ein doppelter Eigenwert auftritt,
dann ergeben sich die Lösungen zu

[mm]\pmat{x\left(t\right) \\ y\left(t\right) }=C_{1}*ev_{1}*e^{2t}+C_{2}*ev_{2}*t*e^{2t}[/mm]


>
>
> Dankeschön.


Gruß
MathePower

Bezug
                
Bezug
Lösung eines lin. DGL-Systems: weitere Hilfe Bitte
Status: (Frage) beantwortet Status 
Datum: 14:44 Fr 08.05.2009
Autor: Ultio


> Das stimmt nicht ganz.
>  
> Berechne hier sämtliche Potenzen von A
>  
> [mm]e^{t*A}:=\summe_{k=0}^{\infty}\bruch{t^{k}}{k!}*A^{k}[/mm]
>  
> In der Regel reichen hier die ersten Glieder um
> eine Systematik dahinter zu erkennen.

Ist das so richtig?:
[mm] e^{t*A}:=\summe_{k=0}^{\infty}\bruch{t^{k}}{k!}*A^{k}=1 [/mm] + [mm] t*\pmat{ 2 & 1 \\ 0 & 2} [/mm] + t/2 [mm] \pmat{ 2 & 1 \\ 0 & 2}^{2} [/mm] + t/6 [mm] \pmat{ 2 & 1 \\ 0 & 2}^{3}+t/24 \pmat{ 2 & 1 \\ 0 & 2}^{4}... [/mm]
=1 + t * [mm] \pmat{ 2 & 1 \\ 0 & 2} [/mm] + t/2 [mm] \pmat{ 4 & 4 \\ 0 & 4} [/mm] + t/6 [mm] \pmat{ 8 & 12 \\ 0 & 8} [/mm] + t/24 [mm] \pmat{ 16 & 32 \\ 0 & 16} [/mm]
weiter als zu verallgemeiner  mit [mm] \pmat{ 2^{k} & ??? \\ 0 & 2^{k}} [/mm] komme ich aber gerade nicht oder was meintest du mit Systematik?


> Die Matrix [mm]e^{t*A}[/mm] ist dann die Fundamentallösung dieses
> DGL-Systems.
>  

stimmt ich erinnere mich, das hatten wir kurz angesprochen.




>
> Zunächst benötigst Du einen Eigenvektor zum Eigenwert 2.
>  
> Dieser ergibt sich aus dem Gleichungssystem
>  
> [mm]\pmat{0 & 1 \\ 0 & 0}*ev_{1}=\pmat{0 \\ 0}[/mm]
>  

du musst entschuldigen, Eigenvektoren und Einheitswerte hatten wir in der linearen Algebra nicht und in Analysis 2 ebenso wenig, daher bin ich hier ein bisschen fragebedürftig bzw. hilfebedürftig.
ist das denn richtig?:
[mm] \pmat{ 0 & 1 \\ 0 & 0} [/mm] * [mm] \pmat{ x \\ y} [/mm] = [mm] \pmat{ 0\\ 0} [/mm]
dann ist
0 = 0*x + 1*y
0 = 0*x + 0*y
aber damit kann ich irgendwie nicht anfangen für alle x ist es 0 da multiplikation mit 0 und y = 0.


> Hier brauchst Du noch einen sogenanten Eigenvektor 2.
> Stufe.
>  
> Dieser ist Lösung von
>  
> [mm]\pmat{0 & 1 \\ 0 & 0}*ev_{2}=ev_{1}[/mm]
>  
> Diese Gleichung ergibt sich unmittelbar aus
>  
> [mm]\pmat{0 & 1 \\ 0 & 0}^{2}*ev_{1}=\pmat{0 \\ 0}[/mm]
>  

damit ist das gleiche Siel:
[mm] \pmat{0 & 1 \\ 0 & 0}^{2}*ev_{1}=\pmat{0 \\ 0} [/mm] = pmat{0 & 1 [mm] \\ [/mm] 0 & 0}*pmat{0 & 1 [mm] \\ [/mm] 0 & 0} = pmat{0 & 0 [mm] \\ [/mm] 0 & 0}
das heißt für alle x und alle y oder wie ist das zu verstehen?

>
> Wenn, wie hier. ein doppelter Eigenwert auftritt,
> dann ergeben sich die Lösungen zu
>  
> [mm]\pmat{x\left(t\right) \\ y\left(t\right) }=C_{1}*ev_{1}*e^{2t}+C_{2}*ev_{2}*t*e^{2t}[/mm]
>  
> >

und damit komme ich ja gar nicht klar, entsprechend meinen oberen Nichtverständnisses.


Danke für die Antwort,
Gruß

Bezug
                        
Bezug
Lösung eines lin. DGL-Systems: Antwort
Status: (Antwort) fertig Status 
Datum: 18:26 Fr 08.05.2009
Autor: MathePower

Hallo Ultio,

>
> > Das stimmt nicht ganz.
>  >  
> > Berechne hier sämtliche Potenzen von A
>  >  
> > [mm]e^{t*A}:=\summe_{k=0}^{\infty}\bruch{t^{k}}{k!}*A^{k}[/mm]
>  >  
> > In der Regel reichen hier die ersten Glieder um
> > eine Systematik dahinter zu erkennen.
>  
> Ist das so richtig?:
>  [mm]e^{t*A}:=\summe_{k=0}^{\infty}\bruch{t^{k}}{k!}*A^{k}=1[/mm] +
> [mm]t*\pmat{ 2 & 1 \\ 0 & 2}[/mm] + t/2 [mm]\pmat{ 2 & 1 \\ 0 & 2}^{2}[/mm] +
> t/6 [mm]\pmat{ 2 & 1 \\ 0 & 2}^{3}+t/24 \pmat{ 2 & 1 \\ 0 & 2}^{4}...[/mm]
>  
> =1 + t * [mm]\pmat{ 2 & 1 \\ 0 & 2}[/mm] + t/2 [mm]\pmat{ 4 & 4 \\ 0 & 4}[/mm]
> + t/6 [mm]\pmat{ 8 & 12 \\ 0 & 8}[/mm] + t/24 [mm]\pmat{ 16 & 32 \\ 0 & 16}[/mm]


Das muß hier so lauten:

[mm]=\pmat{1 & 0 \\ 0 & 1} + t * \pmat{ 2 & 1 \\ 0 & 2} + t^{\red{2}}/2 \pmat{ 4 & 4 \\ 0 & 4} + t^{\red{3}}/6 \pmat{ 8 & 12 \\ 0 & 8} + t^{\red{4}}/24 \pmat{ 16 & 32 \\ 0 & 16}+ ...[/mm]


Sicher habt ihr schon die Fundamentallösung eines DGL-Systems mit doppelten Eigenwerten behandelt.

Interessant ist hier der Eintrag in der 1. Zeile, 2. Spalte:

[mm]a_{12}=0+t+\bruch{4}{2!}t^{2}+\bruch{12}{6}t^{3}+\bruch{32}{24}*t^{4}+...[/mm]


[mm]=t+\bruch{4}{2!}t^{2}+\bruch{12}{3!}t^{3}+\bruch{32}{4!}*t{4}+...[/mm]

[mm]=t*\left(1+\bruch{4}{2!}t+\bruch{12}{3!}t^{2}+\bruch{32}{4!}*t{3}+...\right)[/mm]

[mm]=t*\left(1+\bruch{2}{2}*\bruch{2t}{1!}+\bruch{3}{3}*\bruch{4}{2!}t^{2}+\bruch{4}{4}\bruch{8}{3!}*t{3}+...\right)[/mm]

[mm]=t*\left(\left(2t\right)^{0}+\bruch{\left(2t\right)^{1}}{1!}+\bruch{\left(2t\right)^{2}}{2!}+\bruch{\left(2t\right)^{3}}{3!}+...¸\right)[/mm]

Die Reihe in der Klammer dürfte Dir bekannt sein.


>  
> weiter als zu verallgemeiner  mit [mm]\pmat{ 2^{k} & ??? \\ 0 & 2^{k}}[/mm]
> komme ich aber gerade nicht oder was meintest du mit
> Systematik?
>  
>
> > Die Matrix [mm]e^{t*A}[/mm] ist dann die Fundamentallösung dieses
> > DGL-Systems.
>  >  
> stimmt ich erinnere mich, das hatten wir kurz
> angesprochen.
>  
>
>
>
> >
> > Zunächst benötigst Du einen Eigenvektor zum Eigenwert 2.
>  >  
> > Dieser ergibt sich aus dem Gleichungssystem
>  >  
> > [mm]\pmat{0 & 1 \\ 0 & 0}*ev_{1}=\pmat{0 \\ 0}[/mm]
>  >  
>
> du musst entschuldigen, Eigenvektoren und Einheitswerte
> hatten wir in der linearen Algebra nicht und in Analysis 2
> ebenso wenig, daher bin ich hier ein bisschen
> fragebedürftig bzw. hilfebedürftig.
>   ist das denn richtig?:
>  [mm]\pmat{ 0 & 1 \\ 0 & 0}[/mm] * [mm]\pmat{ x \\ y}[/mm] = [mm]\pmat{ 0\\ 0}[/mm]
>  
> dann ist
> 0 = 0*x + 1*y
>  0 = 0*x + 0*y
>  aber damit kann ich irgendwie nicht anfangen für alle x
> ist es 0 da multiplikation mit 0 und y = 0.


Nun, daraus folgt die Lösungsmenge

[mm]\pmat{x \\ y }=s*\pmat{1 \\ 0}[/mm]

, das heißt [mm]\pmat{1 \\ 0}[/mm] ist der erste Eigenvektor.


>  
>
> > Hier brauchst Du noch einen sogenanten Eigenvektor 2.
> > Stufe.
>  >  
> > Dieser ist Lösung von
>  >  
> > [mm]\pmat{0 & 1 \\ 0 & 0}*ev_{2}=ev_{1}[/mm]
>  >  
> > Diese Gleichung ergibt sich unmittelbar aus
>  >  
> > [mm]\pmat{0 & 1 \\ 0 & 0}^{2}*ev_{1}=\pmat{0 \\ 0}[/mm]
>  >  
>
> damit ist das gleiche Siel:
>  [mm]\pmat{0 & 1 \\ 0 & 0}^{2}*ev_{1}=\pmat{0 \\ 0}[/mm]Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)

= pmat{0 &

> 1 [mm]\\[/mm]Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)

Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)

0 & 0}*pmat{0 & 1 [mm]\\[/mm]Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)
Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)

0 & 0} = pmat{0 & 0 [mm]\\[/mm]Eingabefehler: "{" und "}" müssen immer paarweise auftreten, es wurde aber ein Teil ohne Entsprechung gefunden (siehe rote Markierung)

0 & 0}

>  das heißt für alle x und alle y oder wie ist das zu
> verstehen?


Der Eigenvektor 2. Stufe muß natürlich in [mm]\operatorname{ker}\left(A-2*I\right)^{2}[/mm] liegen.

Also Lösung von

[mm]\pmat{0 & 1 \\ 0 & 0}^{2}*ev_{2}=\pmat{0 \\ 0}[/mm]

ist.

Wie schon erwähnt reduziert sich diese Gleichung auf

[mm]\pmat{0 & 1 \\ 0 & 0}*ev_{2}=ev_{1}=\pmat{1 \\ 0}[/mm]


>  
> >
> > Wenn, wie hier. ein doppelter Eigenwert auftritt,
> > dann ergeben sich die Lösungen zu
>  >  
> > [mm]\pmat{x\left(t\right) \\ y\left(t\right) }=C_{1}*ev_{1}*e^{2t}+C_{2}*ev_{2}*t*e^{2t}[/mm]
>  
> >  

> > >
> und damit komme ich ja gar nicht klar, entsprechend meinen
> oberen Nichtverständnisses.
>  
>
> Danke für die Antwort,
>  Gruß


Gruß
MathePower

Bezug
Ansicht: [ geschachtelt ] | ^ Forum "Gewöhnliche Differentialgleichungen"  | ^^ Alle Foren  | ^ Forenbaum  | Materialien


^ Seitenanfang ^
www.matheforum.net
[ Startseite | Forum | Wissen | Kurse | Mitglieder | Team | Impressum ]