Wenige Zeilen VBA-Code genügen, um mit Excel lineare Gleichungssysteme zu lösen:

Tabellenblatt

Wir schauen uns das Tabellenblatt zum obigen Makro an, mit Eingabe und Lösung des linearen Gleichungssystems

Die Koeffizienten geben wir also ab Zelle B2 in den grauen Bereich und die Konstanten in den grünen Bereich ein. Das Pluszeichen muss man nicht eintippen, und Leerstellen werden automatisch durch Nullen ersetzt (siehe Programm-Code). Die Lösung   x₁ = 3,   x₂ = 4,   x₃ = -2 in Spalte O ist gelb unterlegt. Der graue Bereich bietet Platz für maximal 10 Gleichungen mit 10 Unbekannten - Makro und Tabellenblatt kann man aber leicht modifizieren, so dass dann auch wesentlich größere Systeme gelöst werden können

Anmerkungen zum Programm-Code

'Gleichungen und Unbekannte zählen - also n und m berechnen :
Für jede Zeile und Spalte des Eingabebereichs B2:K11 wird überprüft, ob Einträge existieren, ggf. erhöhen sich die Werte von n bzw. m um jeweils 1. Dabei verwenden wir die Signum-Funktion "Sgn" und ersparen uns somit IF-Prozeduren. Die Standard-Methoden zum Zählen von Zeilen und Spalten sind hier ungeeignet, weil die eingegebenen Koeffizienten mitunter keinen zusammenhängenden Bereich bilden.

Hinweis: Auf Excel-Tabellenfunktionen kann man in VBA mit dem Zusatz "Application.WorksheetFunction" zugreifen. Wann dieser Code tatsächlich erforderlich ist, hängt von der Objekt-Bibliothek ab, in der sich die Funktion befindet. In den unterschiedlichen Excelversionen ist das nicht einheitlich geregelt. Wir haben im obigen Makro nur "Application" vorangestellt - in den aktuellen Excelversionen funktioniert das einwandfrei.

'Namen festlegen:
Mit den Zahlen n und m lässt sich nun der genaue Koeffizienten-Bereich (im obigen Beispiel B2:D4) bestimmen und benennen ("matrix"), ebenso die Zellen in Spalte L ("vektor_b").

'Leerstellen im Gleichungssystem durch Nullen ersetzen :
Leerstellen in den Bereichen "matrix" und "vektor_b" werden nun durch Nullen ersetzt, damit kann man beide Bereiche mathematisch als Matrix bzw. Vektor verwenden.

'Vorherige Arrayformeln löschen :
Teile von Arrayformeln lassen sich nicht überschreiben - wir löschen einfach die Inhalte der reservierten Bereiche für die inverse Matrix und die Lösung.

'Eingabe auf Fehler/Lösbarkeit prüfen:
Im Fall "If Not IsNumeric ..." stimmt die Anzahl n der Gleichungen nicht mit der Anzahl m der Unbekannten überein, oder es liegen nichtnumerische Einträge vor. Im Fall "Else-IF" hat die Determinante der Koeffizientenmatrix den Wert Null, und dann gibt es keine oder unendlich viele Lösungen.

'Lösung berechnen:
Die drei letzten Programmzeilen sind der wesentliche Teil des Makros. Wenn MDeterm eine Zahl ungleich Null liefert, hat das LGS eine eindeutige Lösung, die sich durch Multiplikation der inversen Matrix mit den Zahlen aus Spalte L ergibt. Die entsprechende Array-Formel in der Spalte O liefert dann die gesuchte Lösung.

Mathematische Erläuterungen

Die Koeffizienten eines linearen Gleichungssystems mit n Gleichungen und n Unbekannten kann man zu einer quadratischen Matrix A anordnen und die Determinante det(A) berechnen. Im Fall det(A) ≠ 0 existiert die zu A inverse Matrix A^-1, und das Gleichungssystem hat genau eine Lösung.

Hierzu schauen wir uns das Gleichungssystem aus dem Tabellenblatt an:

Weil die Determinante von A ungleich Null ist, existiert die inverse Matrix A^-1. Somit können wir A·x = b von links mit A^-1 multiplizieren. Die Lösung des linearen Gleichungssystem lautet also: x = A^-1·b - in VBA übersetzt: "=MMULT(inverse,vektor_b)".

Bei sehr umfangreichen Berechnungen ist es wegen der unvermeidlichen Rechenungenauigkeit von Computersystemen ratsam, die Lösung mit einer Fehlerabschätzung zu überprüfen und ggf. automatisch zu korrigieren.

Wenn es weniger Gleichungen als Unbekannte gibt oder im Fall det(A) = 0, existieren entweder unendlich viele oder überhaupt keine Lösungen. Die Matrix A ist dann nicht invertierbar, und der Programmieraufwand steigt erheblich. Zu den unendlich vielen Lösungen bleibt anzumerken, dass dies natürlich nicht auf Systeme zutrifft, die sich auf endliche Zahlenstrukturen beziehen (etwa bei der Nachrichtenverschlüsselung).

Falls es mehr Gleichungen als Unbekannte gibt, ist das LGS im Allgemeinen nicht lösbar; solche überbestimmten Gleichungssysteme resultieren etwa aus Messfehlern. Zwar existiert keine exakte Lösung, aber unter allen "falschen Lösungen" ist diejenige gesucht, die den theoretischen Vorgaben am besten entspricht. Diese sogenannten Ausgleichsprobleme lassen sich mit Excel und VBA grundsätzlich gut lösen - auch in nichtlinearen Fällen.