Eine (Vektor)funktion
f
:
R
n
→
R
{\displaystyle f:\ \mathbb {R} ^{n}\rightarrow \mathbb {R} }
heißt linear,
wenn es Zahlen
a
1
,
a
2
,
…
,
a
n
∈
R
{\displaystyle a_{1},a_{2},\ldots ,a_{n}\in \mathbb {R} }
gibt, sodass
f
(
x
)
=
a
1
x
1
+
a
2
x
2
+
…
+
a
n
x
n
{\displaystyle f(\mathbf {x} )=a_{1}\mathbf {x} _{1}+a_{2}\mathbf {x} _{2}+\ldots +a_{n}\mathbf {x} _{n}}
mit
∀
x
i
∈
R
n
{\displaystyle \forall \mathbf {x} _{i}\in \mathbb {R} ^{n}}
gilt.
Eine (Vektor)funktion
f
:
R
n
→
R
{\displaystyle f:\ \mathbb {R} ^{n}\rightarrow \mathbb {R} }
Funktion ist genau dann linear, wenn
f
(
λ
x
+
μ
y
)
=
λ
f
(
x
)
+
μ
f
(
y
)
{\displaystyle f(\lambda \mathbf {x} +\mu \mathbf {y} )=\lambda f(\mathbf {x} )+\mu f(\mathbf {y} )}
mit
λ
,
μ
∈
R
;
x
,
y
∈
R
n
{\displaystyle \lambda ,\mu \in \mathbb {R} ;\;\mathbf {x} ,\mathbf {y} \in \mathbb {R} ^{n}}
gilt.
Wie bereits bekannt, kann jeder beliebige Vektor als Linearkombination
x
=
x
1
e
1
+
x
2
e
2
+
…
+
x
n
e
n
{\displaystyle \mathbf {x} =x_{1}\mathbf {e} _{1}+x_{2}\mathbf {e} _{2}+\ldots +x_{n}\mathbf {e} _{n}}
dargestellt werden.
Somit ist
y
=
f
(
x
)
=
x
1
f
(
e
1
)
+
x
2
f
(
e
2
)
+
…
+
x
n
f
(
e
n
)
{\displaystyle \mathbf {y} =f(\mathbf {x} )=x_{1}f(\mathbf {e} _{1})+x_{2}f(\mathbf {e} _{2})+\ldots +x_{n}f(\mathbf {e} _{n})}
.
Für die Bilder der Einheitsvektoren gilt
f
(
e
j
)
=
a
1
j
f
(
e
1
)
+
a
2
j
f
(
e
2
)
+
…
+
a
n
j
f
(
e
n
)
{\displaystyle f(\mathbf {e} _{j})=a_{1j}f(\mathbf {e} _{1})+a_{2j}f(\mathbf {e} _{2})+\ldots +a_{nj}f(\mathbf {e} _{n})}
.
Und somit gilt weiter
y
=
(
a
11
x
1
+
a
12
x
2
+
…
+
a
1
n
x
n
)
e
1
+
…
+
(
a
n
1
x
1
+
a
n
2
x
2
+
…
+
a
n
n
x
n
)
e
n
{\displaystyle \mathbf {y} =\left(a_{11}x_{1}+a_{12}x_{2}+\ldots +a_{1n}x_{n}\right)\mathbf {e} _{1}+\ldots +\left(a_{n1}x_{1}+a_{n2}x_{2}+\ldots +a_{nn}x_{n}\right)\mathbf {e} _{n}}
.
Natürlich kann man dies auch übersichtlicher in Matrixschreibweise ausdrücken
y
=
A
x
;
A
=
(
a
i
j
)
{\displaystyle \mathbf {y} =A\mathbf {x} ;\quad A=(a_{ij})}
.
A ist dabei eine quadratische Matrix.
Gegeben sei die lineare Vektorfunktion
y
=
A
x
{\displaystyle \mathbf {y} =A\mathbf {x} }
. Die Frage nach
der Parallelität der beiden Vektoren x und y führt zum speziellen Eigenwertproblem.
x ist dann parallel zu y , wenn
y
=
λ
x
=
λ
E
x
;
λ
∈
R
{\displaystyle \mathbf {y} =\lambda \mathbf {x} =\lambda E\mathbf {x} ;\ \lambda \in \mathbb {R} }
gilt.
(
A
−
λ
E
)
x
=
0
{\displaystyle \left(A-\lambda E\right)\mathbf {x} =0}
Dies ist ein homogenes Gleichungssystem mit der
(
n
×
n
)
{\displaystyle (n\times n)}
-Matrix
A
−
λ
E
{\displaystyle A-\lambda E}
. Dieses
Gleichungssystem hat nur dann nichttriviale Lösungen
x
≠
o
{\displaystyle \mathbf {x} \neq \mathbf {o} }
, wenn
r
g
(
A
−
λ
E
)
<
n
{\displaystyle {\rm {{rg}(A-\lambda E)<n}}}
. Es muss also gelten
P
(
λ
)
=
det
(
A
−
λ
E
)
=
0
{\displaystyle P(\lambda )=\det \left(A-\lambda E\right)=0}
Diese Gleichung wird auch als charakteristische Gleichung der Matrix A bezeichnet.
P
(
λ
)
=
det
(
A
−
λ
E
)
=
(
−
1
)
n
λ
n
+
b
n
−
1
λ
n
−
1
+
…
+
b
0
{\displaystyle P(\lambda )=\det \left(A-\lambda E\right)=(-1)^{n}\lambda ^{n}+b_{n-1}\lambda ^{n-1}+\ldots +b_{0}}
heißt charakteristisches Polynom
der Matrix A .
Die Wurzeln
λ
1
,
…
,
λ
n
{\displaystyle \lambda _{1},\ldots ,\lambda _{n}}
der charakteristischen Gleichung nennt man
die Eigenwerte der Matrix A .
Die Lösungen
x
i
≠
o
{\displaystyle \mathbf {x} _{i}\neq \mathbf {o} }
des linearen Gleichungssystems
(
A
−
λ
E
)
x
=
0
{\displaystyle \left(A-\lambda E\right)\mathbf {x} =0}
heißen Eigenvektoren der Matrix A .
Beispiel:
A
=
(
2
0
4
3
)
{\displaystyle A={\begin{pmatrix}2&0\\4&3\end{pmatrix}}}
P
(
λ
)
=
det
(
2
−
λ
0
4
3
−
λ
)
=
(
2
−
λ
)
(
3
−
λ
)
=
λ
2
−
5
λ
+
6
=
0
{\displaystyle P(\lambda )=\det {\begin{pmatrix}2-\lambda &0\\4&3-\lambda \end{pmatrix}}=(2-\lambda )(3-\lambda )=\lambda ^{2}-5\lambda +6=0}
λ
1
,
2
=
5
2
±
(
5
2
)
2
−
6
{\displaystyle \lambda _{1,2}={\frac {5}{2}}\pm {\sqrt {\left({\frac {5}{2}}\right)^{2}-6}}}
λ
1
=
3
:
{\displaystyle \lambda _{1}=3:}
(
2
−
λ
1
0
4
3
−
λ
1
)
(
x
1
x
2
)
=
(
0
0
)
{\displaystyle {\begin{pmatrix}2-\lambda _{1}&0\\4&3-\lambda _{1}\end{pmatrix}}{\begin{pmatrix}x_{1}\\x_{2}\end{pmatrix}}={\begin{pmatrix}0\\0\end{pmatrix}}}
eine mögliche Lösung:
x
(
1
)
=
(
0
1
)
{\displaystyle \mathbf {x} _{(1)}={\begin{pmatrix}0\\1\end{pmatrix}}}
Übung 1: Berechnen sie einen Eigenvektor
x
(
2
)
{\displaystyle \mathbf {x} _{(2)}}
für den
zweiten Eigenwert
λ
2
=
2
{\displaystyle \lambda _{2}=2}
des obigen Beispiels.
Übung 2: Berechnen sie Eigenwerte und Eigenvektoren für die Matrix
A
=
(
4
0
3
1
1
4
2
−
1
0
)
{\displaystyle A={\begin{pmatrix}4&0&3\\1&1&4\\2&-1&0\end{pmatrix}}}
.
Ist die
n
×
n
{\displaystyle n\times n}
-Matrix
A
=
(
a
i
j
)
{\displaystyle A=(a_{ij})}
symmetrisch
(
A
=
A
T
)
{\displaystyle (A=A^{T})}
und reell, so sind auch alle n Eigenwerte dieser Matrix reell.
Die zu verschiedenen Eigenwerten gehörenden Eigenvektoren einer reellen, symmetrischen Matrix sind orthogonal:
x
(
i
)
x
(
j
)
=
0
{\displaystyle \mathbf {x} _{(i)}\mathbf {x} _{(j)}=0}
.
Hauptachsentransformation: Die Matrix A kann mittels einer orthogonalen Matrix T auf Diagonalgestalt gebracht werden. Die Spalten i der orthogonalen Matrix T sind die auf Eins normierten Eigenvektoren
x
(
i
)
{\displaystyle \mathbf {x} _{(i)}}
.
T
−
1
A
T
=
T
T
A
T
=
(
λ
1
0
⋱
0
λ
n
)
{\displaystyle T^{-1}AT=T^{T}AT={\begin{pmatrix}\lambda _{1}&&0\\&\ddots \\0&&\lambda _{n}\end{pmatrix}}}
Übung : Gegeben ist die Matrix
A
=
(
1
0
3
0
2
−
1
3
−
1
3
)
{\displaystyle A={\begin{pmatrix}1&0&3\\0&2&-1\\3&-1&3\end{pmatrix}}}
.
Berechnen sie die Eigenwerte, Eigenvektoren und führen sie die Hauptachsentransformation durch
(Hinweis: Gram-Schmidtsches-Orthogonalisierungsverfahren).
Eine Bewegung (in der Ebene oder im Raum) ist definiert durch die Gleichung
x
′
=
A
x
+
t
{\displaystyle \mathbf {x} ^{'}=A\mathbf {x} +\mathbf {t} }
Die Matrix A sei dabei orthogonal. Dies bedeutet, dass
det
A
=
±
1
{\displaystyle \det A=\pm 1}
und
A
A
T
=
E
{\displaystyle AA^{T}\ =E}
.
det
A
=
1
{\displaystyle \det A=\ 1}
: orientierungstreue Bewegung (Schraubung = Drehung + Translation)
det
A
=
−
1
{\displaystyle \det A=\ -1}
: orientierungsumkehrende Bewegung (Umlegung)
A
=
E
{\displaystyle A\ =\ E}
: Translation
det
A
=
1
∧
t
=
o
{\displaystyle \det A=1\ \wedge \ \mathbf {t} =\mathbf {o} }
: Drehung um eine durch den Ursprung verlaufende Gerade
x
1
′
=
a
11
x
1
+
a
12
x
2
+
t
1
{\displaystyle x_{1}'=\ a_{11}x_{1}+a_{12}x_{2}+t_{1}}
x
2
′
=
a
21
x
1
+
a
22
x
2
+
t
2
{\displaystyle x_{2}'=\ a_{21}x_{1}+a_{22}x_{2}+t_{2}}
Wegen der Orthogonalität von A gilt:
a
11
2
+
a
21
2
=
1
{\displaystyle a_{11}^{2}+a_{21}^{2}=1}
a
12
2
+
a
22
2
=
1
{\displaystyle a_{12}^{2}+a_{22}^{2}=1}
a
11
a
12
+
a
21
a
22
=
0
{\displaystyle a_{11}a_{12}+a_{21}a_{22}=\ 0}
Diese Gleichungen werden durch folgende trigonometrischen Gleichungen erfüllt:
cos
2
φ
+
sin
2
φ
=
1
{\displaystyle \cos ^{2}\varphi +\sin ^{2}\varphi =1}
(
−
sin
φ
)
2
+
cos
2
φ
=
1
{\displaystyle (-\sin \varphi )^{2}+\cos ^{2}\varphi =1}
−
cos
φ
sin
φ
+
cos
φ
sin
φ
=
0
{\displaystyle -\cos \varphi \sin \varphi +\cos \varphi \sin \varphi =0}
Das bedeutet, daß die Matrix
A
=
(
cos
φ
−
sin
φ
sin
φ
cos
φ
)
{\displaystyle A={\begin{pmatrix}\cos \varphi &-\sin \varphi \\\sin \varphi &\cos \varphi \\\end{pmatrix}}}
eine Drehung der Ebene mit dem Winkel
φ
{\displaystyle \varphi }
um einen festen Punkt bewirkt.
Drehung des Koordinatensystems um den Winkel
φ
{\displaystyle \varphi }
bei fest gehaltenem Vektor:
(
x
1
x
2
)
=
(
cos
φ
−
sin
φ
sin
φ
cos
φ
)
(
x
1
′
x
2
′
)
{\displaystyle {\begin{pmatrix}x_{1}\\x_{2}\end{pmatrix}}={\begin{pmatrix}\cos \varphi &-\sin \varphi \\\sin \varphi &\cos \varphi \\\end{pmatrix}}{\begin{pmatrix}x_{1}'\\x_{2}'\end{pmatrix}}}
Drehung des Vektors um den Winkel
φ
{\displaystyle \varphi }
bei fest gehaltenem Koordinatensystem:
(
x
1
′
x
2
′
)
=
(
cos
φ
−
sin
φ
sin
φ
cos
φ
)
(
x
1
x
2
)
{\displaystyle {\begin{pmatrix}x_{1}'\\x_{2}'\end{pmatrix}}={\begin{pmatrix}\cos \varphi &-\sin \varphi \\\sin \varphi &\cos \varphi \\\end{pmatrix}}{\begin{pmatrix}x_{1}\\x_{2}\end{pmatrix}}}