Lineare Algebra für Mechanik

Aus testwiki
Zur Navigation springen Zur Suche springen

Vorlage:Regal

Zusammenfassung des Projekts

Vorlage:StatusBuch

  • Zielgruppe:

Physik-Studenten in der Experimental-Physik Vorlesung Mechanik in Wikiversity

  • Lernziele:

Die Leser werden eine Anfänger-Vorlesung in Mechanik mathematisch verstehen können.

  • Buchpatenschaft / Ansprechperson:

Ansprechperson: astroboi Das Buch ist Grundmaterial zu dem Wikiversity-Kurs Mathematik für Mechnanik

  • Sind Co-Autoren gegenwärtig erwünscht?

Es wird zunächst ein Kern-Text mit den Kern-Themen gemacht. Später sind Co-Autoren erwünscht. Zunächst bitte alle Vorschläge und Textbausteine auf die Diskussionsseite.

  • Richtlinien für Co-Autoren:

Bitte sich abstimmen mit dem Wikiversity-Projekt Mathematische Methoden der Physik

  • Projektumfang und Abgrenzung zu anderen Wikibooks:

Das Buch wird nur die Themen aufgreifen, die zum Verständnis einer Anfänger-Vorlesung nötig sind. Es sollen viele Beispiele und Aufgaben aufgenommen werden. Ich werde versuchen, die Begriffsbildungen physikalsch zu motivieren.

  • Themenbeschreibung:

Die Themen behandeln die Mathematik so wie es die Mechanik eines Massenpunktes braucht.

  • Aufbau des Buches:

Erklärungen der Grundbegriffe

In dem Teilgebiet der Mathematik, das »Lineare Algebra« heißt, haben Vektoren eine etwas andere Bedeutung als in der Physik; wir beschäftigen uns hier jedoch nur mit physikalischen Vektoren und ihrer mathematischen Behandlung. Die mathematische Fassung der Begriffe ist natürlich richtig und die physikalische Auffassung schwammig. Dafür ist die physikalische Auffassung schnell zu begreifen und eigentlich wollen wir uns ja mit Physik beschäftigen. Wir werden deswegen nur so viel mathematisches Reisegepäck mitnehmen, wie wir für das Verständnis der Physik gebrauchen können.

Die Charakterisierungen in einem mathematischen Nachschlagewerk wird ein Anfänger nicht verstehen. Tatsächlich gibt es mehrere Charakterisierungen desselben Begriffes. Wir verwenden für den Anfang naive Charakterisierungen der Grundbegriffe. Zu einem anderen Zeitpunkt werden wir den Begriffsapparat nachschärfen.

Dimension

Physikalische Größen haben Dimensionen Länge, Zeit, Kraft, Arbeit... Die SI-Einheit der physikalischen Größe Länge ist Meter mit dem Symbol m. Eckige Klammer bilden die Dimension auf die Einheit ab. Man schreibt das so: [Länge] = m. Allgemein also [Physikalische Größe]= Einheit.

Skalar

Ein Skalar ist eine Zahl wie 5, -10, i, -i, 2. Skalare Größenwerte haben Einheiten angehängt 7 m, 7 s, 2V. Die Arbeit W ist zum Beispiel ein Skalar. Aus der Schule kennt man

W = Kraft x Weg = F·s

Die Formel kann man leider nur anwenden, wenn Kraft und Weg auf einer Linie liegen. Wir lernen später ein besonderes Produkt kennen, das beschäftigt sich mit solchen Fällen: das Skalarprodukt.

Skalare nennt man auch Tensoren 0.Ordnung.

Vektor

Ein Vektor hat einen Angriffspunkt, einen Betrag und eine Richtung. Er kann durch einen Pfeil dargestellt werden, er heißt auch Vektor. Den Begriff Vektor kann man noch viel präziser fassen. Für den Anfang ist das nicht notwendig. Eine vektorielle physikalische Größe ist z.B. die Kraft mit der SI-Einheit N. Bei vektoriellen physikalischen Größen sagen wir nicht Vektorbetrag, sondern besser Größenwert. Die Kraftvektor F hat z.B. den Betrag 10 und den Größenwert 10 Newton.

Für ein tieferes Verständnis von Vektoren ist es unerlässlich,charakterisieren wir ein genauer

  • ein Geometrisches Objekt,
  • das einem Vektorraum angehört,
  • das durch Koordinaten bezüglich einer Basis dargestellt werden kann,
  • das aber nicht von einer bestimmten Basis abhängt, sondern unter Basiswechsel (= Koordinatentransformation) invariant bleibt.

Tensoren

Vektoren nennt man auch Tensoren 1. Ordnung. Zwischen Tensoren verschiedener Stufen kann man Gleichungen aufstellen, wenn man Proportionalitäten einbaut. Ansonsten sind Tensoren verschiedener Stufen nicht vergleichbar. Gleichungen zwischen Tensoren gleicher Stufe können proportional sein. Die Proportionalität ist ein Tensor niedrigerer Stufe.

Beispiel: Dynamische Grundgleichung

F Kraftvektor, Tensor 1.Stufe a Beschleunigungsvektor, Tensor 1.Stufe m Masse, skalar , Tensor 0.Stufe

F= m*a

Die Kraft F ist der Beschleunigung a proportional, aber nicht gleich. Beide sind Tensoren derselben Stufe. Die Gleichheit kann hergestellt werden wenn man die Beschleunigung a mit einem Tensor 0. Stufe und zwar der skalaren Masse multipliziert. Die skalare Masse ist hier die Proportionalität.

Vorlage:W

(!Hinweis: Es handelt sich um älteren Inhalt der Wikipedia, der in der aktuellen Artikelversion anders enthalten ist!)

Exkurs:

Vorlage:W

(!Hinweis: Es handelt sich um älteren Inhalt der Wikipedia, der in der aktuellen Artikelversion anders enthalten ist!)

Klassifizierung physikalischer Vektoren

Freie Vektoren

Freie Vektoren können parallel zu sich selbst im Raum verschoben werden.

Linienflüchtige Vektoren

Linienflüchtige Vektoren (zum Beispiel Kräfte) sind an ihre Wirkungslinie gebunden und nur längs dieser verschiebbar.

Ortsgebundene Vektoren

Ortsgebundene Vektoren können überhaupt nicht verschoben werden. Dazu gehören die Feldvektoren, die einem bestimmten Punkt zugeordnet sind, und die Ortsvektoren, die immer im Ursprung des Basissystems beginnen.

Klassifizierung der Vektoren nach Transformationsverhalten

Vektoren kann man auch durch ihr Verhalten gegen Rauminversion klassifizieren. Eine Rauminversion ist die Spiegelung aller Raumpunkte an einem bestimmten, vorgegebenen Punkt (z.B. Koordinatenursprung).

Polare Vektoren werden negativ bei Inversion. Axiale Vektoren erhalten sich.

Frage: Wie teilst Du diese vektoriellen Größen ein ?

Kraft, Drehmoment, Geschwindigkeit, Beschleunigung

Frage: Was ist ein Pseudoskalar ?

Ein Nullvektor an einen Vektor angeheftet verändert ihn nicht. Frage: Wie klassifiziert man einen Nullvektor ?

Das Transformationsverhalten von Vektoren beim Basiswechsel wird durch Matrizen vermittelt. Dabei wird eine Matrix von links an den Vektor ranmultipliziert.

Exkurs:

Transformationsverhalten bei Tensoren

Nullvektor

Ein Vektor und sein entgegengesetzter Vektor angeheftet erzeugt den Nullvektor.

Fragen:

  • Gibt es Mengen von Vektoren, die man so einander heften kann, dass sie denn Nullvektor ergeben ?
  • Kann man einen Vektor und einen entgegengesetzten Vektor zu einem Nullvektor aneinanderheften, wenn sie auf verschiedenen Ebenen liegen.
  • Kann man eine Gruppe von Vektoren durch Strecken und Stauchen einzelner Vektoren immer zu einem Nullvektor zusammensetzen ?
  • Was geschieht wenn man den Nullvektor streckt oder staucht ?

Lineare Unabhängigkeit

Wir haben eine Menge von Vektoren gegeben auf einer Ebene. Jeder Vektor wird eindeutig mit einem Index i gekennzeichnet. Jeder der Vektoren läßt sich mit einem Streckfaktor ai strecken oder stauchen, je nachdem ob es ein echter Bruch oder ein unechter Bruch ist.

Wenn es unmöglich ist die Vektoren durch aneinanderheften zum Nullvektor zu machen trotz Strecken und Stauchen der Vektoren, dann sind die Vektoren linear unabhängig. Natürlich ist es verboten, dass alle Streckfaktoren ai gleichzeitig verschwinden.

Koordinatensysteme

Das bekannteste Koordinatensystem sollte das kartesische sein. Hier stehen die Hauptachsen senkrecht zueinander. Daneben gibt es schiefwinklige Koordinatensysteme. Besonders interessant sind die Koordinatensysteme mit den senkrechten Achsen, ihre Basis heißt Orthonormal-Basis (ONB). Mit der Vektoraddition können wir erst anfangen zu rechnen, wenn wir die Komponenten eines Vektors erklärt haben. In der Schule hat bestimmt jeder schon einmal das kartesische Koordinatensystem kennengelernt. Das sind einfach zwei Achsen senkrecht zueinander. Der Schnittpunkt heißt Koordinatenursprung. Wenn man einen Vektor einzeichnet, dann findet man die kartesischen Komponenten des Vektors durch das Lot auf die Abzisse und die Ordinatenachse. Die Abzisse und die Ordinatenachse sind selber Vielfache von Vektoren und zwar ganz besondere mit den Namen Einheits- und Basisvektoren. Außer den kartesischen Basisvektoren gibt es noch andere Sorten.

Kartesisches Koordinatensystem Koordinaten Polarkoordinaten Schiefwinklige Koordinaten

Übung:

Stelle den Punkt (1,1) im Kartesischen Koordinatensystem in Polarkoordinaten dar.

Einheitsvektoren

Wir können eine Schar von Vektoren auswählen. Jeder Vektor aus einer Menge soll sich durch aneinanderheften aus dieser Schar entstehen können. Diese Vektoren aus der Schar nennen wir Einheitsvektoren. Einheitsvektoren haben den Betrag 1. Orthogonale Einheitsvektoren stehen aufeinander senkrecht.

Übung: Was sind die Einheitsvektoren bei den Polarkoordinaten ?

Basis

Die Einheitsvektoren bilden eine Menge genannt Basis.

Komponentendarstellung, Koordinaten

Die Projektion eines Vektors auf einen Einheitsvektor bedeutet das Lot vom Anfangs- und Endpunkt auf den Einheitsvektor finden. Es schneidet einen Abschnitt auf dem Einheitsvektor aus. Wenn der Einheitsvektor einen Maßstab hat, kann man die Endpunkte des Abschnittes mit Skalaren bezeichnen. Ein Anfangs- und ein Endpunkt ergibt in der Ebene zwei Skalare oder ein Paar. Im Raum wäre es ein 3-Tupel mit 3 Skalaren. Die Skalare aus der Projektion des Anfangs- und Endpunktes tragen wir in ( , ,...) ein mit dem Komma als Trennzeichen.

Prototyp: Kartesisches Koordinatensystem

2D: Ebene

(0,0) , (-1,5),...

3D: Raum

(0,0,0) , (-1,4, 5),...

Die Skalare nennen wir Komponenten oder Koordinaten

Satz: (Invarianz gegenüber Basiswechsel)

Ein Vektor ist unabhängig von seiner Darstellung durch die Basis.

Koordinatensysteme auf einer gekrümmten Kurve

Allgemeine Koordinaten Koordinatenlinien Koordinatenraster Geradlinige- & Krummlinige Koordinatenlinien

In jedem Punkt einer gekrümmten Kurve kann man eine ONB konstruieren. Auf diese ONB kann man dann Ortsvektoren in dem Punkt beziehen und in Komponenten zerlegen. Natürlich ändert sich die ONB dann von Punkt zu Punkt und bewegt sich mit der Kurve mit. Das soll hier nur erwähnt werden.

Vektoren

Die Kraft F ist eine Größe mit einer Richtung und einem Betrag, solche Größen nennen wir Vektoren. Ein Skalar ist einfach eine reelle Zahl. In der Tensorsprache heißt ein Skalar auch Tensor 0.ter Ordnung und ein Vektor Tensor 1.Ordnung nur um das einmal gehört zu haben.

Ein Vektor in der Physik bedeutet etwas anderes als ein Vektor in der Mathematik. Ein physikalischer Vektor hat einen Angriffspunkt, einen Betrag und eine Richtung, in der Mathematik wird der Vektor-Begriff axiomatisch gefaßt. Wir bleiben aber im Rahmen des Buches bei dem anschaulicheren Begriff.

Damit wir mit den Vektoren etwas anfangen können, müssen wir damit Zahlen errechnen können. Anders als in der Mathematik geht es uns nicht nur um mathematische Zusammenhänge, sondern wir haben in der Physik auch Zahlwerte aus Messungen.

Elementare Rechenregeln

Gruppe der Vektoren

Wir erklären dem Leser am Anfang den Gruppenbegriff. Das Kapitel hätte auch Rechenregeln für Vektoren heißen können, aber wir erläutern den Gruppenbegriff hier, weil der Leser schon mal die Begriffe in diesem Kurs gehört haben sollte, auf die er ständig in den umfangreicheren Kurs Lineare Algebra treffen wird. Wir empfehlen dem Leser sich gründlich mit der Lineare Algebra zu beschäftigen, auch wenn sie trocken sind. Denn in der Linearen Algebra finden Begriffsbildungen statt, auf die man später immer wieder in der Mathematik zur Quantenmechanik trifft. Es erleichtert das Verständnis, wenn man da Analogie-Brücken von der Begrifflichkeit in der Linearen Algebra aus bauen kann. Übrigens ist Quantenmechanik auch eine ziemlich trockene Angelegenheit.

Die erste Gruppe der Rechenregeln, die wir für die Menge der Vektoren haben, heißt auch Gruppe. Eine Gruppe nennt man eine Menge mit gewissen Verknüpfungsregeln darauf.

Wir bezeichnen einmal die Menge der Vektoren mit V und mit der Verknüpfung + haben wir die Gruppe (V,+) der Vektoren. Das Pluszeichen steht hier für die Gruppen-Operation, im Falle der Vektoren ist das gleichbedeutend mit der Addition. Es hätte auch z.B. * stehen können.

  • (G1) Abgeschlossenheit

Zwei Vektoren a und b ergeben wieder einen Vektor a+b. Es ist nicht möglich bei der Addition zweier Vektoren ein Skalar z.B. 1 zu bekommen.

  • (G2) Assoziativität

Klammern zählen in der Vektoraddition nicht.

(a+b)+c = a+(b+c)
  • (G3) Neutrales Element

Das neutrale Element ist der Nullvektor. Nicht verwechseln mit der Zahl 0.

Nullvektor O.

Diese drei Strukturmerkmale bei dem Paar (Menge, Operation) bilden schon den Gruppenbegriff. Wenn noch weitere Merkmale, wird die Gruppe z.B. eine kommutative Gruppe.

  • (G4*)Kommutativität
Die Vektoraddition kann man vertauschen.
a+b = b+a

Die Vertauschbarkeit der Operation gegen die Gruppenelemente kann man nicht aus den drei Gruppen-Geetzen herleiten.

Zur Übung kann der Leser den ersten Satz aus den Gruppen-Gesetzen ableiten.

Satz über das Inverses Element:

Das inverse Element zu a ist erklärt a + inverses Element = neutrales Element (hier O). Das inverse Element bekommt hier das Symbol -a .

a + (-a) = O

Beweis: -> Leser

Aufgabe: Warum sind die nätürlichen Zahlen und die Operation - keine Gruppe ?

Wenn man die inversen Elemente erklärt hat, dann kann man eine abkürzende Schreibweise für die Anordnung

a+(-b) zu dem kürzeren a-b erklären. Das bekommt den neuen Namen Vektor-Subtraktion. Das ist eigentlich keine neue Operation, sondern einfach eine Abkürzung. Die Vektor-Subtraktion ist wie die Vektor-Addition erklärt.

Der Leser hat jetzt begriffen, was eine Gruppe ist. Wir werden in diesem Buch später noch weitere Gruppen kennelernen, nämlich die Gruppe der orthogonalen Drehungen; wenn wir das Koordinatensystem drehen werden.

Um den Begriff Gruppen gibt es die Gruppentheorie. Davon braucht der Leser noch nichts zu wissen.

Die nächste Gruppe der Rechenregeln erklärt die Verknüpfung mit Elementen aus einer anderen Menge als die ursprüngliche Gruppe. Wir nennen sie Multiplikation und die jetzt zusätzlich auftretende Menge nennen wir Skalarkörper. Denn es wird die Menge der reelen Zahlen sein und ihre Elemente heißen Skalare in der Vektorsprache für die richtungslosen Größen.

Die Rechenregeln dazu nennt man Körperaxiome.

  • (K1) Skalarmultiplikation

Skalarprodukt

Das skalare Produkt zweier Vektoren soll ein Skalar erzeugen. Das Skalarprodukt mißt die Ähnlichkeit zweier Vektoren

Charakterisierung:

𝐕𝐖=VWcos(𝐕,𝐖).

Der Cosinus-Faktor bringt einen Winkel im Bogenmaß hinein. Je kolinearer die Vektoren sind desto näher ist der Cosinus an 1 dran, so ist das Ähnlichkeitsmaß zu verstehen.

Ein Vektorraum mit einem Skalarprodukt heißt unitärer Raum.

Entsprechend würde ein Sinus-Faktor ausdrücken wie "verschieden" die Vektoren von der kolinearen Lage wären. Das nutzen wir bei dem vektorwertigen oder Vektorprodukt aus.

abab

Übung: Beweise die Ungleichung, indem Du die Charakterisierung des Skalarproduktes mit dem Cos-Faktor ausnutzt !

Anwendung

Für eine ortsunabhängige Kraft F, die entlang eines geraden Weges s wirkt, ist die Arbeit definiert als das Skalarprodukt

W=Fs=Fscosα,

wobei α der Winkel zwischen der Richtung der Kraft und der Richtung des Weges ist und F und s die Beträge der entsprechenden Vektoren sind.

Wenn die wirkende Kraft in Richtung des zurückgelegten Weges angreift und konstant ist, dann vereinfacht sich dieser Ausdruck zu

W=Fs.

Übung: Zerlege die Kraft F in einen parallelen Anteil zur Wirkungslinie also kolinear zu s.

Der Ausdruck sollte sein

F=1F2(Fs)F

Kosinussatz

c=ab
c2=(ab)2=a22ab+b2
c2=a2+b22abcosα

Anwendungen des Skalarproduktes

Norm

Das Skalarprodukt erzeugt eine Norm ||V|| charakterisiert als

||V|| = V*V

Mit der Norm können wir Vektoren normieren, indem wir den Vektor V durch seine Norm ||V|| teilen. Einheitsvektoren sind immer normiert.

Ein Vektorraum mit einer Norm heißt normierter Raum.

Mit der Norm können wir Vektoren normieren, indem wir den Vektor V durch seine Norm ||V|| teilen. Einheitsvektoren sind immer normiert.

Die Norm hat immer 3 Eigenschaften

  1. x=0x=0 (Definitheit);
  2. αx=|α|x (Homogenität);
  3. x+yx+y (die w: Dreiecksungleichung).

Ein Vektorraum mit einer Norm heißt normierter Vektorraum oder normierter Raum.

Anwendung

Arbeitsbegriff

Für eine ortsunabhängige Kraft F, die entlang eines geraden Weges s wirkt, ist die Arbeit definiert als das Skalarprodukt

W=Fs=Fscosα,

wobei α der Winkel zwischen der Richtung der Kraft und der Richtung des Weges ist und F und s die Beträge der entsprechenden Vektoren sind.

Wenn die wirkende Kraft in Richtung des zurückgelegten Weges angreift und konstant ist, dann vereinfacht sich dieser Ausdruck zu

W=Fs.

Projektion

Projektion

Eine Projektion eines Vektors x auf die Richtung eines anderen Vektors n. Diese Projektion meint das man die Anfangspunkte von a und n aneinanderheften soll. Das Lot von a auf n ergibt einen Schnittpunkt.

Die Projektion von a auf n ist ein Vektor mit den gemeinsamen Anfangspunkt von a,n und den Schnittpunkt als Endpunkt. Die Projektion gibt sozusagen den Anteil von n am Vektor a wieder.

Übung: Zerlege die Kraft F in einen parallelen Anteil zur Wirkungslinie also kolinear zu s.

Der Ausdruck sollte sein

F=1F2(Fs)F

Kosinussatz

Kosinussatz

c=ab
c2=(ab)2=a22ab+b2
c2=a2+b22abcosα

Projektion, Norm

Satz von Schwarz

Cauchy-Schwarzscher Satz

abab
0(a+αb)(a+αb)
a2+α2b2+αba+αab
a2+α2b2+2αab
α
α=abb2
0a2b2(ab)2

Satz von Thales

Übung: Zeige den Satz von Thales. Hinweis: Identifiziere dabei die Terme in (a+b)(ba) mit den geometrischen Größen im Satz.

Dreiecks-Ungleichung

Dreiecks-Ungleichung

Es muss heißen:

|ab|a±ba+b

aba+ba+b
a2+b22aba2+b2+2aba2+b2+2ab
(ab)2(a+b)2(a+b)2 normieren
ababab
aba+ba+b

Mehr in der Wikipedia zum Normierten Raum Skalarmultiplikation Skalarprodukt Vektorprodukt Mehrfache Produkte Normierter Vektorraum Komponentenschreibweise

Vektorprodukt

Vektorprodukt

Vorlage:Merke

Es heißt auch äußeres Produkt, Kreuzprodukt oder Vektorprodukt. Die Charakterisierung unterscheidet sich von dem Skalarprodukt durch die Verwendung eines Sinusfaktors. Übrigens nennt man das Skalarprodukt auch inneres Produkt.

W=|𝐔×𝐕|=||U||||V||sinφ=UVsin(𝐔,𝐕).

Sinussatz

Anwendung: das Drehmoment M

M=r×F

Übung: Projiziere die Kraft F in ihren senkrechten Anteil zu dem Weg s.

Lösung sollte

F=1F2F×(s×F)F

Kreuzprodukt

Skalarmultiplikation Skalarprodukt Vektorprodukt Mehrfache Produkte Normierter Vektorraum Komponentenschreibweise

Mehrfache Vektorprodukte

Basisvektoren

Orthonormalität und Vollständigkeit des Basissystems ONB

Dass man einen Vektor so schreiben kann oder nach einer Basis entwickeln kann ist ein immer wieder auftauchendes Verfahren. Wir wenden das immer an, um bequem rechnen zu können. Am komfortabelsten ist natürlich die Rechnung in einer Basis zu machen, die möglichst einfache Koordinaten liefert. Unter den einfachen Basissystemen ist am prominentesten die orthonormalen Basissystemen, hier stehen die Basisvektoren senkrecht aufeinander, was man mit dem Skalarprodukt feststellen kann, denn das Skalarprodukt orthogonaler Basisvektoren verschwindet.

Die wichtigsten Zusammenhänge sind die Vollständigkeit- , die Orthonormalitäts-Relationen und ein Entwicklungssatz.

Orthonormalität

Das Skalarprodukt der Basisvektoren ergibt paarweise Null.

𝐞1𝐞2=𝐞2𝐞3=𝐞3𝐞1=0.

Nach der Charakterisierung des Skalarproduktes mit dem cos-Faktor stehen die Einheitsvektoren also senkrecht aufeinander. Wenn wir die Einheitsvektoren verlängern dann begrenzen sie einen rechtwinkliges Raster.

Normierung

Normierte Basisvektoren heißen Einheitsvektoren, weil die Norm 1 ist.

𝐞1𝐞1=𝐞2𝐞2=𝐞3𝐞3=1

Vollständigkeit

Die Basisvektoren reichen aus, um den gesamten Vektorraum aufzuspannen. Die Vollständigkeitsrelation erklären wir mit dem Kroneckerdelta

δ ist das sogenannte Kronecker-Symbol. Es verschwindet, wenn die Indizes i = j , und es ist 1 bei i = j.

δ={1,wenn i=j0,wenn i=j

Der Zusammenhang mit unserem Thema steckt in der Gleichung:

k=13ek,i ek,j=δij,

wobei der erste Index - wie schon oben - den "Namen" des Basisvektors und der zweite Index die Komponente (z.B. x-, y-, oder z-Koordinate im kartesischen Koordinatensystem) bezeichnet: (𝐞k)i=ek,i

Entwicklungssatz

Jeder Vektor kann in der orthonormierten Basis dargestellt werden. Orthonormiert ist ein Kunstwort aus normiert und orthogonal.

𝐕=(V1,V2,V3)=𝐕1+𝐕2+𝐕3=V1𝐞1+V2𝐞2+V3𝐞3.

Die Darstellung eines Vektors in einem ONB ist eindeutig. ONS steht Orthonormiertes Basis-System.

Die Komponenten/Koordinaten eines Vektors 𝐚 bezeichnen wir mit einem angehängten Index a1, a2 oder a3. Dabei bedeuten sie das Produkt aus dem Skalarprodukt Vektor 𝐚 mit dem Basisvektor 𝐞i, i= 1,2,3.

ai=𝐚𝐞i

Die Darstellung des Vektors a in der Basis wird dann zu einer Summe.

𝐚=i=13ai 𝐞i=a1 𝐞1+a2 𝐞2+a3 𝐞3

Zusammenhang mit Vollständigkeit:

Setzt man nun in den Audruck für die k-te Komponente des Vektors 𝐚 den obigen Ausdruck für die Koordinaten ai ein, so erhält man:

(𝐚)k=ak=i=13ai (𝐞i)k=i=13(𝐚𝐞i) (𝐞i)k=i=13j=13aj ei,j ei,k=j=13aji=13ei,j ei,k=j=13aj δjk=ak

Beim vorletzten Gleichzeichen wurde die Vollständigkeitsrelation verwendet. Das heißt, dass ein Vektor nur dann vollständig nach einer Basis entwickelt werden kann, wenn die Basisvektoren die Vollständigkeitsrelation erfüllen.

Komponentendarstellungen

Rechnen mit Vektoren in Komponentendarstellung

Wir zeigen, wie man mit die Vektor-Operationen rechnet in der Komponentendarstellung. Wir tun das für das kartesische (kanonische) Basissystem.

Summe und Differenz zweier Vektoren

Man addiert zwei Vektoren, indem man bei gleicher Basis ihre skalaren Komponenten addiert.

ei(a+b)=ai+bi=ci für i = 1,2,3
c=(a1+b1,a2+b2,a3+b3)

Addition und Subtraktion

Die Summe der beiden Vektoren

a=a1i+a2j+a3k und b=b1i+b2j+b3k

berechnet sich als:

a+b=(a1+b1)i+(a2+b2)j+(a3+b3)k=(a1+b1a2+b2a3+b3).

Skalarprodukt zweier Vektoren

Skalarprodukt

Das Skalarprodukt ist die Summe der Komponentenprodukte

(ab)
= i=13(aiei) j=13(ajej)
= i,j=13(aibj)(eiej)
= i=13(aibj)δij
(ab) = i=13(aibi)
δ ist das sogenannte Kronecker-Symbol. Es verschwindet, wenn die Indizes i := j , und es ist 1 bei i = j.
δ={1,wenn i=j0,wenn i=j

Im Kartesischen Koordinatensystem berechnet sich das Skalarprodukt als inneres Produkt

ab=(a1a2a3)(b1b2b3)=a1b1+a2b2+a3b3

insbesondere gilt für das Quadrat eines Vektors

aa=(a1a2a3)(a1a2a3)=a12+a22+a32

Die Wurzel daraus ist die Norm des Vektors a.

Vektorprodukt zweier Vektoren

Rechtssystem

Diese Formeln charakterisieren ein Rechtssstem.

𝐞1×𝐞2=𝐞3,𝐞2×𝐞3=𝐞1,𝐞3×𝐞1=𝐞2,𝐞2×𝐞1=𝐞3,𝐞3×𝐞2=𝐞1,𝐞1×𝐞3=𝐞2,𝐞1×𝐞1=𝐞2×𝐞2=𝐞3×𝐞3=𝟎.

Diese Formeln für die Basisvektoren in der ONB haben die angenehme Eigenschaft, dass man das Kreuzprodukt recht einfach auf eine komponentenweise Berechnung der Koordinaten zurückführen kann.

Das Prinzip ist wie beim Skalarprodukt.

Vergleich:

a =i=13aei = ae1+ae2+ae3

Hier hatten wir die Projektionen auf die Richtung der ei aufsummiert. Nun müssen wir für das Kreuzprodukt a×b den entstehenden Vektor auf die Richtungen von ei×ej projizieren und dann die Projektionen aufsummieren. Es sind die Richtungen zu berechnen auf das wir den Produktvektor a×b projizieren wollen.

Die Kreuzprodukte der Faktoren sind einfach und bilden auf die Koordinaten des Produktvektors ab, wie man an der rechten Seite erkennt.

e1×e2=e3 ,
e2×e3=e1 ,
e3×e1=e2

Es ist nur noch die Projektion zu leisten, dass wir zur Vorbereitung mit einem Basisvektor vornehmen.

ei(ej×ek) =: εijk

Dieses Produkt ist das bekannte Spatprodukt aus den Basisvektoren und ist eine Prominente. Epsilon-Tensor - total antisymmetrischer Tensor. Ein Hinweis darauf, dass die Basis ein Rechtssystem ist.

εijk={+1falls (i,j,k,) eine gerade Permutation von (1,2,3,) ist,1falls (i,j,k,) eine ungerade Permutation von (1,2,3,) ist,0wenn mindestens zwei Indizes gleich sind.

Kreuzprodukt

Damit fassen wir die Kreuzprodukte der Basis-Vektoren übersichtlich zusammen

(ei×ej) = k=13εijkek

Mit dem Kreuzprodukt der Basis-Vektoren können wir nun auch komponentenweise das Kreuzprodut anschreiben.

a×b=εijkajbkei=εijkaibjek
(a×b)i=j=13k=13εijkajbk.

Das ist im Grunde eine Kurzdarstellung der drei Komponenten-Gleichungen

c1 = a2b3a3b2
c2 = a3b1a1b3
c3 = a1b2a2b1

Das soll der Leser selber mit den obigen Formeln nachrechnen, um den Umgang mit dem εijk zu lernen. Später werden wir diese Fertigkeit bei Beweisen gut gebrauchen können.

Epsilon-Tensor

Das Spatprodukt

Das Spatprodukt kann man einfach formulieren mit den früheren Formeln für das Kreuzprodukt

a(b×c)
=i,j,kaibjckei(ej×ek)
=i,j,kεi,j,kaibjck

Eine wichtige Eigenschaft des Spatproduktes ist, das man × und vertauschen kann.

Nachrechnen am Epsilon-Tensor !

Das Spatprodukt ist positiv, wenn die Permutation gerade ist und negativ falls ungerade.

Die Reihenfolge der Indizes ist {1,2,3}. Die Reihenfolge kann auch vertauscht sein {2,1,3}, hier ist genau eine Vertauschung der Indizes mit seinem Nachbarn geschehen. Die Permutation ist deswegen ungerade. Wenn ich bei {2,1,3} 1 -> 3 vertausche entsteht die gerade Permutation {2,3,1}.

Zyklische Vertauschung meint im Uhrzeigersinn vertauschen. Antizyklische Vertauschung gegen den Uhrzeigersinn. Dabei schreib man zuerst den Index 1 an im Uhrzeigersinn den Index 2 und dann den Index

Im Gegensatz zur zyklischen Vertauschung tritt bei der antizyklischen Vertauschung ein Vorzeichenwechsel auf:

(a,b,c)=(b,a,c)

Zwei gleiche Vektoren lassen a×a=0 verschwinden:

(a,a,b)=0

Sonst gelten die vertrauten Eigenschaften von Produkten.

Skalare Multiplikation α ist assoziativ:

(αa,b,c)=α(a,b,c)

Distributiv:

(a,b,c+d)=(a,b,c)+(a,b,d)

Spatprodukt

Vektorprodukt dreier Vektoren (»Entwicklungssatz«)

Wir brauchen die k-te Komponente des doppelten Vektorproduktes a x (b x c). Für uns sind das zwei Vektorprodukte

a×(b×c)

Jedes läßt sich mit dem ϵijk Tensor in eine Komponenten-Gleichung schreiben

[a×(b×c)]k
= i,jai(b×c)j
= i,jl,mεijkεlmjaiblcm
= i,jl,mεijkεjlmaiblcm

Mit der Formel

εijkεjlm = δilδkmδimδkl

Die Herkunft dieser Formel erleben wir in den Übungen. Setze Sie an dieser Stelle einfach ein ohne sich Gedanken zu machen.

[a×(b×c)]k
= i,l,maiblcm(δi,mδk,lδi,lδk,m)
= i(aibkciaibick)
= bk(ac)ck(ab)
= [b(ac)c(ab)]k

Die Merkregel dafür ist einprägsam und hat einen Namen Bac-cab-Regel

Vektoriell ist dies auch als die Graßmann-Identität bekannt:

a×(b×c)=(ac)b(ab)c
a×(b×c)+b×(c×a)+c×(a×b)=0

Weitere Produkte mit vektoriellen Faktoren

Lagrange-Identität

(a×b)(c×d)=(ac)(bd)(bc)(ad)

Matrizen

Rechenregeln für Matrizen

Koordinatentransformation

e1,e2,e3
r = (x1,x2,x3) in
r = (x1,x2,x3) in
j=13xjej = j=13xjej
ej = j=13djkek
ej = kdjkek
djm = ejem=cosϕjm
D=(d11d12d13d21d22d23d31d32d33)
eiej=δij=k,mdikdjm(ekem=mdimdjm)
mdimdjm=mcos(ϕim)cos(ϕjm)=δij
xi=j=13xj(ejei)=j=13cosϕijxj; i= 1,2,3
(x1x2x3) = (d11d12d13d21d22d23d31d32d33) (x1x2x3) r () = Dr()
D1D=DD1=E
D1r () = D1Dr () = Er ()
D1 (x1x2x3) = D1D (x1x2x3) = E (x1x2x3) = (x1x2x3)
xi=j=1xj(ejei) = j=13cosϕjixj i = 1,2,3
(x1x2x3) = (d11d12d13d21d22d23d31d32d33) (x1x2x3) r() = D1r ()
D1=DT=((dij1=dji)
δij = mdim(d1)mj=dimdjm
δij = m(d1)imdmj=dmidmj

Drehungen in der Ebene

Mehrfache Drehungen in der Ebene

Drehung im Raum der 3-Achse

Determinanten

Rechenregeln für Determinanten

Für eine nur aus einem Koeffizienten bestehende 1×1-Matrix A ist

det(A)=det(a11)=a11

Geometrische Deutun: Parallelogramm-Fläche

Ist A eine 2×2-Matrix, dann ist

det(A)=det(a11a12a21a22)=a11a22a21a12

Geometrische Deutung: Spatvolumen

Für eine 3×3-Matrix A gilt die folgende Formel:

det(A)=det(a11a12a13a21a22a23a31a32a33)=a11a22a33+a12a23a31+a13a21a32a13a22a31a11a23a32a12a21a33

Will man diese Determinante von Hand berechnen, so stellt die w: Regel von Sarrus dafür ein einfaches Schema zur Verfügung.

Epsilontik

Zwischen dem Epsilon-Tensor und dem Kronecker-Delta gilt die Beziehung

εijkεlmn=|δilδimδinδjlδjmδjnδklδkmδkn|=δilδjmδkn+δimδjnδkl+δinδjlδkmδimδjlδknδilδjnδkmδinδjmδkl

Übung: Zeige dass gilt

εijkεimn=δjmδknδjnδkm
εijkεijn=2δkn
εijkεijk=6

(wiederum mit Summenkonvention). Sie ist nützlich bei der Vektorrechnung im 3.

Die Determinante einer n×n-Matrix A=(Aij) kann mit dem Levi-Civita-Symbol und der Summenkonvention wie folgt geschrieben werden:

detA=εi1i2inA1i1A2i2Anin.
  • Berechne:
Die Determinante dieser Matrix
(100010001)
  • Berechne das Kreuzprodukt
(100)×(010)
  • Berechne das Spatprodukt
(100)(010)×(001)

Inverse Matrix

Als inverse Matrix A1 wird die Matrix bezeichnet, die die Beziehung

A1A=AA1=E

erfüllt. In diesem Fall wird A als invertierbar oder auch regulär bezeichnet. Aufgrund der Kommutativität der obigen Beziehung ist A in jedem Fall eine quadratische Matrix. Die inverse Matrix einer quadratischen Matrix A existiert genau dann, wenn die Determinante dieser Matrix nicht verschwindet, wenn also det(A)0 gilt.

Die Berechnung der inversen Matrix kann mittels des Verfahrens von Gauß-Jordan erfolgen. Hierzu wird die Matrix durch Addition und Subtraktion einzelner Zeilen der Matrix mit anderen Zeilen der Matrix sowie Multiplikation und Division einzelner Zeilen mit passenden reellen Zahlen in die Einheitsmatrix umgeformt (dies ist im Fall einer invertierbaren Matrix stets möglich). Indem nun die gleichen Rechenschritte (Addition und Subtraktion einzelner Zeilen sowie Multiplikation und Division einzelner Zeilen mit reellen Zahlen), die zur Umformung der Matrix in die Einheitsmatrix nötig waren auf eine Einheitsmatrix angewandt werden, ergibt sich die inverse Matrix aus der auf diese Weise umgeformten Einheitsmatrix:

(a11......a1n....................................an1......ann)(10......001..............................100......01)(10......001..............................100......01)(b11......b1n....................................bn1......bnn)

Es gilt:

(a11......a1n....................................an1......ann)(b11......b1n....................................bn1......bnn)=(10......001..............................100......01)


Beispiel

Gegeben sei die quadratische Matrix (2101). Offensichtlich lässt sich diese Matrix durch Subtraktion der zweiten Zeile von der ersten Zeile und anschließendes Dividieren der ersten Zeile durch 2 in die Einheitsmatrix umformen. Wendet man die selben Rechenschritte (Subtraktion der zweiten Zeile von der ersten Zeile und anschließendes Dividieren der ersten Zeile durch 2) auf die Einheitsmatrix an, ergibt sich die inverse Matrix: (1001)(1101)(121201).

Wie man sofort nachrechnen kann ist (2101)(121201)=(1001).

Rotation

Drehmatrix

Lineare Gleichungssysteme

Koordinatensysteme

Koordinatenlinien

Alle Koordinatensysteme haben etwas gemeinsam. Sie "überziehen" die Zeichenebene mit einem Muster aus Linien, den Koordinatenlinien. Eine Koordinatenlinie fixen Koordinate und eine veränderliche Koordinate. Die Koordinatenlinien bilden Muster "Koordinaten-Raster". Es erinnert an Längen- und Breitengrade.

Im Fall kartesischer und schiefwinkeliger Koordinaten sind alle Koordinatenlinien gerade. Der Koordinaten-Raster besteht aus zwei Scharen von zueinander parallelen Geraden. Diese Koordinaten werden daher geradlinig genannt.

Polarkoordinaten sind ein Beispiel für krummlinige Koordinaten. Ihre Koordinatenlinien sind

einerseits alle Strahlen ( = Halbgeraden) durch den Ursprung (jene Linien, entlang derer sich der Wert der Koordinate ϕ nicht ändert),
andererseits alle Kreise um den Ursprung (jene Linien, entlang derer sich der Wert der Koordinate r nicht ändert).

Für diese Koordinaten macht der Begriff der Achsen überhaupt keinen Sinn!

Probleme können eine typische Symmetrie haben und in einem angepassten Koordinatensystem einfacher berechnet werden.

Die Koordinaten eines Punktes in Bezug auf ein Koordinatensystem können in jene bezüglich eines anderen Koordinatensystems umgerechnet werden. Man spricht dann von einer Koordinatentransformation.

Funktionaldeterminante

Jacobideterminante

Die Determinante der Jacobi-Matrix spielt z. B. bei Transformationen von Integralen eine wichtige Rolle und wird Funktionaldeterminante genannt.

Für den Fall m=n ist f eine n×n-Abbildung und die Jacobi-Matrix ist quadratisch. Hierfür kann man dann die Jacobi-Determinante berechnen.

Diese Determinante gibt zu einem gegebenen Punkt wichtige Informationen über das Verhalten der Funktion f in der Nähe dieses Punktes. Wenn beispielsweise die Jacobideterminante einer stetig differenzierbaren Funktion in einem Punkt p nicht null ist, so ist die Funktion in einer Umgebung von p invertierbar. Weiterhin gilt, dass bei positiver Determinante in p die Funktion ihre Orientierung beibehält und bei negativer Jacobideterminante die Orientierung umkehrt. Der absolute Wert der Determinante im Punkt p gibt den Wert an, mit dem die Funktion in der Nähe von p expandiert oder schrumpft.

Die Jacobi-Matrix (nach Carl Gustav Jacob Jacobi; auch Funktionalmatrix genannt) einer differenzierbaren Funktion f:nm ist die m×n-Matrix sämtlicher erster partieller Ableitungen. Genutzt wird sie z. B. in der näherungsweisen Berechnung/Approximation oder Minimierung mehrdimensionaler Funktionen in der Mathematik.

Die Jacobi-Matrix bildet die Matrix-Darstellung der ersten Ableitung der Funktion f. Sie wird mit Jf, Df oder fx bezeichnet.

Bezeichnet man die Koordinaten im Urbildraum n mit x=(x1,,xn) und die Komponentenfunktionen von f mit f1,,fm, so lautet die Jacobi-Matrix

Jf=(fixj)i=1,,m; j=1,,n,

beziehungsweise ausführlich

Jf=fx=(f1x1f1x2f1xnfmx1fmx2fmxn).

Sie kann, wenn man sie für einen Punkt p=(p1,,pn) ausrechnet, zur Näherung der Funktionswerte von f in der Nähe von p verwendet werden:

f(x1,,xn)f(p1,,pn)+Jf(p1,,pn)(x1p1xnpn).

Diese affine Abbildung entspricht der Taylor-Approximation erster Ordnung (Linearisierung).

Für m=1 entspricht die Jacobi-Matrix dem Gradienten von f. Je nach Definition des Gradienten, der manchmal als Zeilenvektor und manchmal als Spaltenvektor definiert wird, unterscheidet sich jedoch in diesem Fall die Jacobi-Matrix als Zeilenvektor vom Gradienten.

Ein Beispiel für eine Rechnung mit der Jacobi-Matrix ist die Transformation in Polarkoordinaten.

Beispiel

Zylinderkoordinaten

Die Umrechnungsformeln von Zylinderkoordinaten (r, φ, h) in kartesische Koordinaten lauten:

x=rcosφ
y=rsinφ
z=h

Die Funktionaldeterminante lautet also:

det(x,y,z)(r,φ,h)=|cosφrsinφ0sinφrcosφ0001|=r.

Folglich ergibt sich für das Volumenelement dV:

dV=|(x,y,z)(r,φ,h)|drdφdh=rdrdφdh.

Genauso gut hätte man eine andere Reihenfolge der Polarkoordinaten wählen können. Die Funktionaldeterminante lautet dann beispielsweise:

det(x,y,z)(r,h,φ)=|cosφ0rsinφsinφ0rcosφ010|=r.

In das Transformationsgesetz geht jedoch immer nur der Betrag der Determinante ein, also ist das Ergebnis dann unabhängig von der gewählten Reihenfolge der Variablen, nach denen abgeleitet wird.

Kugelkoordinaten

Die Umrechnungsformeln von Kugelkoordinaten (r,θ,φ) in kartesische Koordinaten lauten:

x=rsinθcosφ,
y=rsinθsinφ und
z=rcosθ.

Die Funktionaldeterminante lautet also:

det(x,y,z)(r,θ,φ)=|sinθcosφrcosθcosφrsinθsinφsinθsinφrcosθsinφrsinθcosφcosθrsinθ0|=r2sinθ.

Folglich ergibt sich für das Volumenelement dV:

dV=|(x,y,z)(r,θ,φ)|drdθdφ=r2sinθdrdθdφ.

Krummlinige Koordinaten

Zylinderkoordinaten

Polarkoordinaten

Für die Umrechnung von ebenen Polarkoordinaten in kartesische Koordinaten gilt demnach:

x=rcosφ
y=rsinφ
r=rcosφex+rsinφey

Winkels im Intervall (−π, π]

Mit Hilfe des Arkuskosinus kommt man mit nur zwei Fallunterscheidungen aus:

φ={+arccosxrfu¨r y0arccosxrfu¨r y<0

Winkels im Intervall [0, 2π)

Die Formel mit dem Arkuskosinus kommt auch in diesem Fall mit nur zwei Fallunterscheidungen aus:

φ={arccosxrfu¨r y02πarccosxrfu¨r y<0

Wenn man ein kartesisches Koordinatensystem so ausrichtet, dass die z-Achsen zusammenfallen und die x-Achse in Richtung φ=0 zeigt, dann ergeben sich die folgenden Umrechnungsformeln:

x=ρcosφ
y=ρsinφ
z=z

Für die Umrechnung von kartesischen Koordinaten in Zylinderkoordinaten ergeben sich für ρ und φ die gleichen Formeln wie bei den Polarkoordinaten.

Polarkoordinaten

Für die Umrechnung von ebenen Polarkoordinaten in kartesische Koordinaten gilt demnach:

x=rcosφ
y=rsinφ
r=rcosφex+rsinφey

Winkels im Intervall (−π, π]

Mit Hilfe des Arkuskosinus kommt man mit nur zwei Fallunterscheidungen aus:

φ={+arccosxrfu¨r y0arccosxrfu¨r y<0

Winkels im Intervall [0, 2π)

Die Formel mit dem Arkuskosinus kommt auch in diesem Fall mit nur zwei Fallunterscheidungen aus:

φ={arccosxrfu¨r y02πarccosxrfu¨r y<0

Kugelkoordinaten

Die Transformationsgleichungen von kartesischen in Kugelkoordinaten lauten

r=x2+y2+z2;
φ={arccosxx2+y2fu¨r y0,[,5em]2πarccosxx2+y2fu¨r y<0;
θ=arccotzx2+y2 = π2arctanzx2+y2.

Die Rücktransformation erfolgt nach den Gleichungen

x=rsinθcosφ
y=rsinθsinφ
z=rcosθ.