E-machten van matrices
-
- Berichten: 101
E-machten van matrices
In deze opgave willen we e-machten van matrices berekenen. Om dit te doen moeten we eerst definieren wat een e-macht van een matrix is. Laat nu A een n × n-matrix zijn, in analogie met de Taylorreeks definieren we
e^A = (som van k=0 tot oneindig) 1/k! * A^k
a) Bereken voor een n × n diagonaalmatrix A zijn e-macht en schrijf deze als 1 matrix
b) Toon aan dat e^(CAC^-1) = C*(e^A)*C^-1
c) Bereken e^A (Hint: Breng A op diagonaalgedaante en gebruik onderdeel b)
met A:
-1 4 -2
-3 4 0
-3 1 3
Iemand een idee hoe je dit moet aanpakken?
e^A = (som van k=0 tot oneindig) 1/k! * A^k
a) Bereken voor een n × n diagonaalmatrix A zijn e-macht en schrijf deze als 1 matrix
b) Toon aan dat e^(CAC^-1) = C*(e^A)*C^-1
c) Bereken e^A (Hint: Breng A op diagonaalgedaante en gebruik onderdeel b)
met A:
-1 4 -2
-3 4 0
-3 1 3
Iemand een idee hoe je dit moet aanpakken?
- Berichten: 7.556
Re: E-machten van matrices
In de definitie van de e-macht van een matrix komt de macht van die matrix voor. Wat weet je over de macht van een diagonaalmatrix?a) Bereken voor een n × n diagonaalmatrix A zijn e-macht en schrijf deze als 1 matrix
Never express yourself more clearly than you think.
- Niels Bohr -
- Niels Bohr -
-
- Berichten: 101
Re: E-machten van matrices
A^n bij diagonaalmatrix wordt
a11^n 0 0
0 a22^n 0
0 0 a33^n
a11^n 0 0
0 a22^n 0
0 0 a33^n
- Berichten: 24.578
Re: E-machten van matrices
Dus wat geeft dat voor de e-macht van een diagonaalmatrix...? Probeer een stap verder te denken.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
e^a11 1 1
1 e^a22 1
1 1 e^a33
Zoiets? . ik heb geen idee eigenlijk
1 e^a22 1
1 1 e^a33
Zoiets? . ik heb geen idee eigenlijk
- Berichten: 24.578
Re: E-machten van matrices
Geen idee? Omdat je gokt en niet de tijd neemt om zelf met pen en papier eens wat te proberen...?
\(e^A = \sum\limits_{k = 0}^{ + \infty } {\frac{{A^k }}{{k!}}} = \sum\limits_{k = 0}^{ + \infty } {\left( {\frac{1}{{k!}}\left( {\begin{array}{*{20}c} a \hfill & 0 \hfill & 0 \hfill \\ 0 \hfill & b \hfill & 0 \hfill \\ 0 \hfill & 0 \hfill & \ddots \hfill \\\end{array}} \right)^k } \right)} = \sum\limits_{k = 0}^{ + \infty } {\left( {\begin{array}{*{20}c} {\frac{{a^k }}{{k!}}} \hfill & 0 \hfill & 0 \hfill \\ 0 \hfill & {\frac{{b^k }}{{k!}}} \hfill & 0 \hfill \\ 0 \hfill & 0 \hfill & \ddots \hfill \\\end{array}} \right)} = \ldots\)
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
Ja, ik volg het zo wel allemaal. Maar wat ik vervolgens zou moeten doen..
Je kan hem nu uitschrijven maar dat is niet echt de bedoeling denk ik.
Voor die term linksboven in de matrix krijg je dan 1, a, a^2/2!, a^3/3! enz Zo voor de hele diagonaal. Maar daar zit geen e-macht in:P
Je kan hem nu uitschrijven maar dat is niet echt de bedoeling denk ik.
Voor die term linksboven in de matrix krijg je dan 1, a, a^2/2!, a^3/3! enz Zo voor de hele diagonaal. Maar daar zit geen e-macht in:P
- Berichten: 24.578
Re: E-machten van matrices
Nu staat er een som van matrices, maar die kan je samennemen tot een matrix (als je weet hoe de optelling voor matrices werkt). Begin dat toch eens uit te schrijven, maar kijk goed wat er gebeurt. Concentreer je bijvoorbeeld op het eerste diagonaalelement, wat komt daar te staan als je de matrices allemaal optelt?
In plaats van uitschrijven, kan je ook symbolisch verder werken. Opnieuw, als je weet hoe de optelling voor matrices werkt, wat gebeurt er dan met de sommatie die nu buiten de matrix staat? Als je dit 'inziet', doe het dan zo. Indien niet, schrijf het dan toch eens (een stukje) uit zodat je ziet wat er gebeurt.
In plaats van uitschrijven, kan je ook symbolisch verder werken. Opnieuw, als je weet hoe de optelling voor matrices werkt, wat gebeurt er dan met de sommatie die nu buiten de matrix staat? Als je dit 'inziet', doe het dan zo. Indien niet, schrijf het dan toch eens (een stukje) uit zodat je ziet wat er gebeurt.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
Die 1e term is gelijk aan e^a als je hem uitschrijft en optelt
- Berichten: 24.578
Re: E-machten van matrices
Inderdaad, en dus de tweede diagonaalterm? En algemeen? Dus hoe ziet de hele matrix eruit?
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
2e is dan e^b 3e e^c
ik snap niet hoe je dit dan kunt generaliseren
ik snap niet hoe je dit dan kunt generaliseren
- Berichten: 24.578
Re: E-machten van matrices
Het is toch duidelijk dat je gewoon de e-macht van elk diagonaalelement moet nemen...? Verder symbolisch:
\(\sum\limits_{k = 0}^{ + \infty } {\left( {\begin{array}{*{20}c} {\frac{{a^k }}{{k!}}} \hfill & 0 \hfill & 0 \hfill \\ 0 \hfill & {\frac{{b^k }}{{k!}}} \hfill & 0 \hfill \\ 0 \hfill & 0 \hfill & \ddots \hfill \\\end{array}} \right)} = \left( {\begin{array}{*{20}c} {\sum\limits_{k = 0}^{ + \infty } {\frac{{a^k }}{{k!}}} } \hfill & 0 \hfill & 0 \hfill \\ 0 \hfill & {\sum\limits_{k = 0}^{ + \infty } {\frac{{b^k }}{{k!}}} } \hfill & 0 \hfill \\ 0 \hfill & 0 \hfill & \ddots \hfill \\\end{array}} \right) = \left( {\begin{array}{*{20}c} {e^a } \hfill & 0 \hfill & 0 \hfill \\ 0 \hfill & {e^b } \hfill & 0 \hfill \\ 0 \hfill & 0 \hfill & \ddots \hfill \\\end{array}} \right)\)
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
Ja dat zei ik toch:P. Maar dat is dan dus al de matrix?
- Berichten: 24.578
Re: E-machten van matrices
Ja toch.... Kijk eens naar het begin (linkerlid) van m'n uitwerking, dat was e^A.
Let wel, dit geldt alleen als A een diagonaalmatrix is! Dit is nog steeds deel a) dus...
Let wel, dit geldt alleen als A een diagonaalmatrix is! Dit is nog steeds deel a) dus...
"Malgré moi, l'infini me tourmente." (Alfred de Musset)
-
- Berichten: 101
Re: E-machten van matrices
Oke, dat was a dan. Snap ik wel.
Bij b dan.
Voor die c te krijgen moet je toch de eigenwaarden en eigenvectoren uitrekenen. Maar hoe zou dat hier moeten?
Bij b dan.
Voor die c te krijgen moet je toch de eigenwaarden en eigenvectoren uitrekenen. Maar hoe zou dat hier moeten?