Zum Inhalt springen

Benutzerin:Gabriele Hornsteiner/Test-Baustelle

Aus Wikibooks

D:\Phili\FH\Lehre\Mathe\Lineare Algebra\Vorlesung

Berechnung einer Inversen

[Bearbeiten]

Inverse einer Matrix

Es gibt zum reziproken Wert einer Zahl eine Entsprechung, nämlich die Inverse einer Matrix :

Die Entsprechung ist so zu verstehen:

bzw.

mit als nichtsingulärer ()-Matrix und als Einheitsmatrix der Ordnung .

Die Inverse einer Matrix wird vor allem beim Berechnen linearer Modelle verwendet.

Wie berechnet sich ? Nicht, indem man etwa für alle Elemente einfach berechnet. Sondern es muß bei der Multiplikation die Einheitsmatrix resultieren.

Meistens ist die Ermittlung der Inverse sehr aufwendig und wird in der Regel nur auf Computern durchgeführt. Es gibt mehrere Verfahren zur Berechnung der Inversen. Hier wird die manuelle Berechnung mittels eines Tableaus und der vollständigen Elimination gezeigt.

Beispiel:

Die Matrix soll invertiert werden. Man beginnt mit einem Tableau, das senkrecht in zwei Hälften geteilt ist: Links steht die aktuelle Matrix , rechts dagegen die Einheitsmatrix. Ziel ist es, die linke Matrix unter Einbeziehung des gesamten Tableaus in eine Einheitsmatrix umzuwandeln. Es steht dann auf der rechten Seite die Inverse von .

Tableau:

Mit der vollständigen Elimination wird das Tableau umgeformt. Wenn links die Einheitsmatrix steht, ist rechts die Inverse gegeben.

Wir arbeiten nun die Spalten des Tableaus nun von links nach rechts ab:

Die 1.Spalte soll sein .

Ein Vielfaches der 1. Zeile wird von einem Vielfachen der 2. Zeile subtrahiert:

1. Zeile ist ok.

2. Zeile

3. Zeile ist ok.

Die erste Spalte ist fertig und wir wollen nun ein neues Tableau als Zwischenergebnis festhalten.

Tableau:

Die 2.Spalte soll sein .

2. Zeile

1. Zeile:

3. Zeile:

Neues Tableau:

Die 3.Spalte soll sein .

3. Zeile

1. Zeile:

2. Zeile:

Ergebnistableau:


Also ist .

Die Probe ergibt

.

Da eine Matrix nur eine Inverse haben kann, haben wir tatsächlich die Inverse von gefunden.

Inverse einer 2x2-Matrix

[Bearbeiten]

Für den Spezialfall einer ()-Matrix mit den Elementen ergibt sich für die Inverse

mit als Determinante von .

Beispiel:

Gegeben ist :.

Wir erhalten

.

Umformung von Matrizengleichungen

[Bearbeiten]

Die vereinfachte Darstellung ganzer Zahlensysteme ermöglicht die knappe Wiedergabe und elegante Berechnung komplexer Ausdrücke wie z.B.

oder . Ebenso kann man auch Matrizengleichungen umformen und bezüglich einer unbekannten Matrix oder eines Vektors auflösen.

Regeln für das Umformen von Gleichungen:

a) Man kann eine Matrizengleichung mit einer Matrix additiv (also mit + oder -) von links oder rechts erweitern, wobei natürlich die Ordnung der Matrizen übereinstimmen muss.

Beispiel:

Gegeben sind , und .



b) Man kann eine Matrizengleichung mit einer Matrix multiplikativ erweitern, wobei natürlich die Ordnung der Matrizen übereinstimmen muss:

Von links:

Beispiel:

Gegeben sind , und .

.

Ausklammern:

.

Von rechts:

Beispiel:

Gegeben sind , und .

.

Ausklammern:

.

Umformungen wie sind im allgemeinen nicht zulässig und oft auch gar nicht definiert. c) Es gilt beispielsweise mit , mit , beide invertierbar, und der Nullmatrix :

.

.

. .)

Skalare werden wie Zahlen behandelt. Insbesondere können die Seiten der Gleichung wahlweise von links und von rechts mit einem Skalar multipliziert werden.

Beispiel:

Beispiele für Umformungen:

Gegeben ist das lineare Gleichungssystem . muss invertierbar sein. Wir wollen den Lösungsvektor ermitteln:

Wir multiplizieren das Gleichungssystem von links mit :

.

Gegeben ist das Gleichungssystem . Es soll nach aufgelöst werden. Also

Alle Ausdrücke mit auf die linke Seite bringen:
nach rechts ausklammern:
Von links mit multiplizieren:
Wegen

erhalten wir , also .

Beispiel 3:

Gegeben ist die Matrix , wobei ist. Gesucht ist :

.

Man kann nun zusammenfassen:

, was ergibt.

Man nennt übrigens eine Matrix mit der Eigenschaft idempotent, was bedeutet, dass diese Matrix beim Potenzieren sich selbst hervorbringt.

Falls Sie meinen, das sei nur so zum Spaß ausgedacht: In der multiplen Regression wird diese Matrix für das Berechnen der Regressionskoeffizienten benötigt.

=== anker ===\ Der Rang gibt die "Information" an, die in einer Matrix steckt.

Wir betrachten noch einmal das zweite Beispiel von Lisa und Familie Meier:

Im Juli kauft Lisa einen Schokoriegel und zwei Tüten Chips und zahlt 8 Euro. Familie Meier kauft vier Schokoriegel und acht Tüten Chips und zahlt 32 Euro.

Die Koeffizientenmatrix ist . Man könnte auch als Zusammenfassung von zwei Spaltenvektoren auffassen: mit und .

Hier sind zwei Spalten (und auch zwei Zeilen ) linear abhängig voneinander.

Lineare Abhängigkeit bedeutet hier bzw. , allgemein ( const.), was bedeutet, dass im Grund schon eine Spalte die nötige Information enthält.

Nun wollen wir ein Beispiel mit drei Spaltenvektoren ansehen.

Gegeben sind drei Spaltenvektoren , und . Falls gilt

,

wobei und Skalare (= Konstanten) sind und mindestens ein Skalar ist, nennt man eine Linearkombination aus und . , und sind dann linear abhängig.

Beispiel: Gegeben sind die beiden Spaltenvektoren 2. Ordnung (also ) und . und sind linear unabhängig, denn es gibt kein , das die Beziehung ermöglicht.

Aus Linearkombinationen mit und können nun alle restlichen Vektoren der Ordnung 2 gebildet werden.

Beispielsweise soll der Vektor mit Hilfe von und dargestellt werden. Wir erreichen das mit

also .

Übrigens kriegt man die Werte raus, indem man ein Gleichungssystem

bzw.

löst.

Man nennt und Basisvektoren. Sie spannen den zweidimensionalen Raum auf.

Wie ist das gemeint? Stellen wir uns ein Koordinatensystem aus und vor. Mit Hilfe der --Paare können wir jeden Punkt in diesem Koordinatensystem darstellen. und spannen also den zweidimensionalen Raum auf. Das Gleiche haben wir mit und gemacht.

Basisvektoren im engeren Sinn sind die Einheitsvektoren und . Sie bilden das kartesische Koordinatensysten. In der Regel werden mit dem Begriff Basisvektoren diese Einheitsvektoren gemeint.

Werden Spaltenvektoren der Ordnung betrachtet, können maximal Vektoren linear unabhängig sein.

Sind beispielsweise die Vektoren , , und der Ordnung 3 gegeben, müssen diese Vektoren insgesamt linear abhängig sein, es gilt also hier

oder auch   (, , , const.) usw.

Betrachten wir die drei Spaltenvektoren 3. Ordnung , und . , und sind linear unabhängig. Jeder Vektor der Ordnung 3 kann als Linearkombination dieser drei Vektoren dargestellt werden, z.B.

.

Entsprechendes gilt übrigens auch für die Zeilenvektoren.

etappe

[Bearbeiten]

Besteht eine -Matrix aus genau unabhängigen Vektoren , hat die Matrix den Rang , also .

Im obigen Beispiel von Lisa und den Meiers ist , weil - wie wir im Kapitel über Lösbarkeit von Gleichungssystemen nur 1 aussagefähige Gleichung existiert. In der ursprünglichen Konstellation des Beispiels ist

,

denn hier sind beide Spaltenvektoren linear unabhängig, kann nicht als Linearkombination von dargestellt werden, beide Zeilen enthalten eigenständige Informationen.

Wie berechnet man den Rang einer Matrix?

Beispielsweise mit Hilfe des Gaußschen Algorithmus. Man bringt die Matrix so weit wie möglich auf Dreiecksgestalt (also auf Trapezgestalt). Die Zahl der Zeilen, die keine Nullvektoren sind, geben den Rang an.

Es gilt:

Bei einer ()-Matrix ist der Rang höchstens der kleinere Wert von und .

Ist der Rang einer quadratischen Matrix (also ) gleich , hat vollen Rang, ist regulär, nichtsingulär, invertierbar. Alle Bezeichnungen sind äquivalent.

Beispiele für Ergebnistableaus von oben. Gesucht ist der Rang von links im Tableau:

I

II

III

IV

V

VI

VII

VIII

Abgearbeitet

[Bearbeiten]

Nehmen wir uns noch einmal das Beispiel mit dem Tante-Emma-Laden in der Einführung vor:

Es sollen beide Gleichungen zugleich gelten, also schreibt man eigentlich korrekt

.

Sind Missverständnisse ausgeschlossen, können die „und“ weggelassen werden.

In Matrizenschreibweise kann man das System darstellen als

,

allgemein

mit der Koeffizientenmatrix , dem Variablenvektor und dem Lösungsvektor .


Lösung linearer Gleichungssysteme

[Bearbeiten]

Gaußscher Algorithmus

[Bearbeiten]

In der Matrizenrechnung wird aus praktischen Gründen das Eliminationsverfahren verwendet.

Wie löst man idealerweise ein lineares Gleichungssystem (kurz: LGS)?

Betrachten wir ein Beispiel:

,


in Matrizenschreibweise

, bzw.
,

Wie lösen wir dieses Gleichungssystem? Am besten rekursiv von unten nach oben. Rekursiv heißt, dass man die Informationen der vorherigen Schritte für den nächsten Lösungsschritt verwendet.

,

Wir beginnen mit Zeile 4:

Wir erhalten

.

Wir nehmen uns nun die dritte Gleichung vor und berechnen :

.

Für die zweite Gleichung ergibt sich

.

Schließlich bestimmen wir :

Unser Ergebnis ist also

.

Dieses Gleichungssystem war einfach zu lösen, denn die Gleichungen bauten aufeinander auf und ermöglichten ein rekursives Lösungsverfahren. Wir können die Rekursivität unmittelbar der Koeffizientenmatrix entnehmen, denn sie ist eine Dreiecksmatrix.

Wie sieht es nun beispielsweise aus mit der Lösung des Gleichungssystems

in Matrizenschreibweise

bzw.

.

Nett wäre es, wenn man auch hier auf Dreiecksgestalt wie oben in der Form bringen könnte. Das kann man mit dem Eliminationsverfahren erreichen. Man wendet nun so genannte Zeilenoperationen an. Zulässige Zeilenoperationen sind:

  1. Multiplikation einer Zeile mit einer Konstanten.
  2. Addition zweier Zeilen.
  3. Vertauschen zweier Zeilen.

Für die manuelle Durchführung der Zeilenoperationen fasst man alle Zahlen zweckmäßigerweise in einem so genannten Tableau zusammen, d.h. man trägt die Koeffizienten und die Lösungen in einer Tabelle zusammen.

Bemerkung: Wenn man nun Klammern darumsetzte, hätte man eine erweiterte Koeffizientenmatrix, was uns inhaltlich aber nicht weiter bringt. Das hervorgehobene "manuell" deutet an, dass es hier vor allem darum geht, die Funktionsweise des Lösungsverfahren zu erfassen. Normale Menschen verwenden spezielle Programme für die Lösungen. Studierende werden aber gequält. Sie müssen sich ja das Examen verdienen.

Tableau:

Die Spalten des Tableaus werden nun von links nach rechts abgearbeitet:

Die 1.Spalte soll sein . Ein Vielfaches der 1. Zeile wird von einem Vielfachen der 2. Zeile subtrahiert:

Ein Vielfaches der 1. Zeile wird von einem Vielfachen der 3. Zeile subtrahiert:

Die erste Spalte ist fertig und wir wollen nun ein neues Tableau als Zwischenergebnis festhalten.

Nun wird die zweite Spalte transformiert. Die 2.Spalte soll sein .

Ein Vielfaches der 2. Zeile wird von einem Vielfachen der 3. Zeile subtrahiert:


Wir speichern Berechnungen als Tableau ab und sehen, dass wir mit den Zeilentransformationen fertig sind, denn unsere Koeffizientenmatrix hat Dreiecksgestalt.

Das entspricht dem Gleichungssystem

.

das nun rekursiv von unten nach oben gelöst werden kann:

Wir erhalten den Lösungsvektor bzw. etwas platzsparender .

Man nennt dieses Verfahren, nämlich das Gleichungssystem in Dreiecksgestalt umzuwandeln,

  • Gaußscher Algorithmus,
  • Gaußsches Iterationsverfahren,
  • Gaußsches (teilweises) Eliminationsverfahren.

Bemerkung:

Oftmals sind in der Koeffizientenmatrix Nullen enthalten, die bei geschickter Platzierung das Eliminationsverfahren erleichtern. Hier können bei Bedarf die Spalten von vertauscht werden. Allerdings muß auf die Reihenfolge der Variablen geachtet werden, z.B.


Vollständiges Eliminationsverfahren

[Bearbeiten]

Das Gaußsche Eliminationsverfahren ermöglicht eine rekursive Berechnung der Variablen von unten nach oben. Das lineare Gleichungssystem kann aber weiter umgeformt werden, so dass z.B. ein Tableau der Form

entsteht. Hier kann die Lösung direkt abgelesen werden:

, , .

Man nennt dieses Verfahren vollständiges Gaußsches Eliminationsverfahren. Das Verfahren ist zwar aufwendiger, wird aber aus verschiedenen praktischen Gründen häufig angewendet, vor allem im EDV-Einsatz. Beispiel: Gegeben ist das Tableau

das mit Hilfe der vollständigen Elimination umgeformt werden soll. Man geht spaltenweise von links nach rechts vor, zeilenweise von oben nach unten.

Die 1.Spalte soll sein .

Man nennt die Zeile, in der schließlich 1 stehen soll, Pivotzeile.

1. Zeile (Pivotzeile):

2. Zeile:


3. Zeile: Ist schon o.k.

Neues Tableau als Zwischenergebnis:

Die 2.Spalte soll sein . Pivotzeile ist also nun Zeile 2.

1. Zeile ist schon o.k.

3. Zeile:


Neues Tableau als Zwischenergebnis:

Die 3.Spalte soll sein . Pivotzeile ist nun Zeile 3.

1. Zeile:

2. Zeile:

Neues Tableau als Ergebnis:

Wir erhalten die Lösungen

, , bzw. den Lösungsvektor . Für unsere Zwecke dieses Mathekurses genügt i.a. das teilweise Gaußsche Eliminationsverfahren.

Lösbarkeit linearer Gleichungssysteme

[Bearbeiten]

Aus den bisherigen Ausführungen wurde ersichtlich, daß z.B. für die Bestimmung von 3 Unbekannten i.a. drei Gleichungen benötigt werden, d.h. für die Bestimmung von Unbekannten braucht man möglicherweise Gleichungen.

Fragen: Genügen immer Gleichungen? Was ist, wenn mehr als Gleichungen vorhanden sind? Ist das System dann noch lösbar?

Im Juni haben wir

als Gleichungssystem

Wir haben also zwei Gleichungen und zwei Unbekannte und lösen nun das LGS wie gewohnt. Wir gehen aus vom Tableau

und wollen die Koeffizientenmatrix zur Dreiecksmatrix umwandeln:

Tableau:

Ein Schokoriegel kostet also 2€ und eine Tüte Chips 3€. Es gibt nur eine Lösung für und . Dieses Gleichungssystem ist eindeutig lösbar.

Im Juli kauft Lisa wie gewohnt einen Schokoriegel und zwei Tüten Chips und zahlt 8 Euro. Familie Meier kauft vier Schokoriegel und acht Tüten Chips und zahlt 32 Euro. Wir interessieren uns wieder für die Preise der Artikel.

Wir gehen aus vom Tableau

und wollen die Koeffizientenmatrix umwandeln:

Hmmm...

Das neue Tableau wird nicht besser:

Wenn wir unser Gleichungssystem genauer betrachten, sehen wir, dass in diesem Fall die zweite Gleichung eigentlich nur das Vierfache der ersten Gleichung ist. In II ist keine neue Information enthalten, so dass wir im Grunde nur eine Gleichung haben. Genau das verrät uns das neue Tableau.

Wir erhalten bzw. . Mehr Infos sind nicht zu bekommen. Wir haben hier letztlich zwei Unbekannte, aber nur eine Gleichung. Wir können allerdings für einen Wert festlegen und dann entsprechend berechnen: oder usw. Es gibt theoretisch unendlich viele Lösungen, je nach Wert von . Dieses System ist also mehrdeutig lösbar. Eine spezielle Lösung ist die so genannte Basislösung, also die Lösung für :.

Meistens finden die Leute eine eindeutige Lösung schicker als eine mehrdeutige. Allerdings hat man bei einer mehrdeutigen mehr Spielraum, etwa bei der Produktionsplanung.

Im August haben wir die Konstellation

als Gleichungssystem

als Tableau

Wir erhalten durch Zeilentransformation

So haben wir hier die Situation, dass ist. Sogar einem Meerschweinchen würde auffallen, dass wir hier eine inkonsistente Gleichung haben, dass also keine Lösung möglich ist.

Es sind bei einem Linearen Gleichungssystem folgende Lösungsmöglichkeiten denkbar:

1. Eine (einzige) eindeutige Lösung

2. Unendlich viele Lösungen

3. Keine Lösung

Es gibt verschiedene Methoden, herauszufinden, inwieweit eine Lösung existiert. Das Gaußsche Eliminationsverfahren ist beispielsweise geeignet: Man bringt die erweiterte Koeffizientenmatrix, also zusammen mit bzw. das entsprechende Tableau auf Trapezgestalt, d.h. man formt so weit wie möglich auf Dreiecksgestalt um. Schließlich ergeben sich Konstellationen wie oben erläutert.

ergebnisse

[Bearbeiten]

Beispiele für Ergebnistableaus:

I

II

III

IV

V

VI

VII

VIII

Lösung für die Tableaus:

  • I eindeutig lösbar
  • II eindeutig lösbar
  • III eindeutig lösbar, vierte Zeile überflüssig
  • IV alle trivial
  • V mehrdeutig lösbar
  • VI mehrdeutig lösbar
  • VII mehrdeutig lösbar
  • VIII keine Lösung

halt2

[Bearbeiten]

2. Beispiel von oben:

konnte bestimmt werden. Was aber ergibt ? Wir versuchen es:

Wir haben zwei Spalten bei , aber drei Zeilen von . Hier ist eine Multiplikation nicht definiert. PQ konnte bestimmt werden. Was aber ergibt QP?

style="text-align:center"|{|

|

|} Was wäre aber P‘Q?

Also ist die Anordnung der Matrizen von Bedeutung.

Beachten:

Im allgemeinen ist

Es gilt AA = A2, AAA = A3, ... (falls die Multiplikation erlaubt ist).

und das aktuelle Produktionsprogramm der Unternehmen in der Matrix

angegeben. Wir suchen die Menge von Produktionsfaktoren , die für das Produktionsprogramm beschafft werden müssen. Es ergibt sich

als

halt1

[Bearbeiten]
  A         B    
a a a   b b b b  
a a a   b b b b  
        b b b b  
                 
2 × 3   3 × 4   Spaltenzahl von A = Zeilenzahl von B!!
         
    2 × 4       neue Matrix

class="wikitable"

| 2 ||×|| colspan="3" style="text-align:center"|3       3|| × || 4|| | 2 ||×|| 3 ||  || 3|| × | class="hintergrundfarbe8"| Y2

| colspan="2" style="text-align:center"| Zelle 2

,|| colspan="3"||| 3|| × || 4||

  A     B
2 ×
3 ... 3
× 4
m ×
r   r
× n

Beispiel

1

halt2

[Bearbeiten]

.

Beispiel mit etwas mehr Verallgemeinerung: Zwei Unternehmen U1 und U2 mit gleichem Produktangebot stellen drei Produkte x, y, z her. Dazu benötigen sie die Produktionsfaktoren a, b, c. Die Inputmenge, die für die Herstellung einer Einheit von x, y oder z benötigt werden, sind in der Produktionsmatrix

und das aktuelle Produktionsprogramm der Unternehmen in der Matrix

angegeben. Wir suchen die Menge von Produktionsfaktoren , die für das Produktionsprogramm beschafft werden müssen. Es ergibt sich

als

Wir multiplizieren:

und erhalten

Lineare Gleichungssysteme

Lisa kauft im Tante-Emma-Laden 1 Schokoriegel und 2 Tüten Chips und zahlt 8 Euro. Familie Meier kauft im selben Laden 4 Schokoriegel und 5 Tüten Chips und zahlt 23 Euro. Was kostet ein Schokoriegel, was eine Tüte Chips?

Gegeben ist also

Schokoriegel Chips Summe
Lisa 1 2 8 €
Fam. Meier 4 5 23 €


Wir bezeichnen den Preis eines Schokoriegels als und den einer Tüte Chips als und erhalten das lineare Gleichungssystem