Vandermonde-Matrix

Unter einer Vandermonde-Matrix (nach A.-T. Vandermonde) versteht man in der Mathematik eine Matrix, die eine im Folgenden beschriebene spezielle Form hat.

Für ein <math>n</math>-Tupel <math> (x_1, x_2, \ldots , x_n) </math> reeller Zahlen oder allgemeiner von Elementen in einem Körper ist die Vandermonde-Matrix definiert durch:

<math>

V (x_1, x_2, \ldots , x_n) = \begin{pmatrix}

 1      & x_1    & x_1^2  & \cdots & x_1^{n-1} \\
 1      & x_2    & x_2^2  & \cdots & x_2^{n-1} \\
 1      & x_3    & x_3^2  & \cdots & x_3^{n-1} \\
 \vdots & \vdots & \vdots & \ddots & \vdots    \\
 1      & x_n    & x_n^2  & \cdots & x_n^{n-1}

\end{pmatrix} </math>

Die Determinante wird auch Vandermonde-Determinante genannt, sie hat den Wert

<math> \det V(x_1,x_2, \ldots, x_n) = \prod_{1 \leq i < j \leq n} (x_j - x_i) </math>.<ref>{{#ifexist:Vorlage:bibISBN/{{#invoke:URIutil|plainISBN|9783834817242}}

| {{bibISBN/{{#invoke:URIutil|plainISBN|9783834817242}}

  |record = Vorlage:bibISBN/{{#invoke:URIutil|plainISBN|9783834817242}}
  |format = Literatur
  |Autor = 
  |Titel = 
  |TitelErg = 
  |Band = 
  |Auflage = 
  |Kommentar= 
  |Kapitel = 
  |Seite = 113–116
  |Seiten = 
  |Spalten = 
  |ArtikelNr = 
  |Fundstelle = 
  |DOI = 
  |Online = 
  |URL = 
  |Linktext = 
  |Format = 
  |KBytes = 
  |Abruf = 
  |Typ = 

}}{{#ifeq: 0 | 0

   | {{#invoke:TemplatePar|check
       |all= 1=
       |opt= 2= format= Autor= Titel= TitelErg= Hrsg= Sammelwerk= WerkErg= Band= Nummer= Auflage= Datum= Sprache= NummerReihe= BandReihe= HrsgReihe= Kommentar= Kapitel= Seite= Seiten= Spalten= ArtikelNr= Fundstelle= DOI= Online= URL= Linktext= Format= KBytes= Abruf= Typ=

|template=Vorlage:bibISBN |cat=Wikipedia:Vorlagenfehler/Vorlage:BibISBN}}

     }}

| {{#if:||{{#ifeq: 0 | 0 | }}}}}}</ref>

Insbesondere ist die Vandermonde-Matrix genau dann regulär, wenn die <math> x_i </math> paarweise verschieden sind.

Anwendung: Polynominterpolation

Die Vandermonde-Matrix spielt bei der Interpolation von Funktionen eine wichtige Rolle: Wenn an den Stützstellen <math> (x_1, x_2, \ldots , x_n) </math> die Funktionswerte <math>(f_1, f_2, \ldots, f_n) </math> durch ein Polynom <math>p</math> vom Grad <math> n - 1 </math> (oder kleiner) interpoliert werden sollen, dann führt der Ansatz

<math>p(x) = a_0 + a_1x^1 + a_2x^2 + \cdots + a_{n-1} x^{n-1}</math>

auf das lineare Gleichungssystem

<math>

\begin{pmatrix}

 1      & x_1    & x_1^2  & \cdots & x_1^{n-1} \\
 1      & x_2    & x_2^2  & \cdots & x_2^{n-1} \\
 1      & x_3    & x_3^2  & \cdots & x_3^{n-1} \\
 \vdots & \vdots & \vdots & \ddots & \vdots    \\
 1      & x_n    & x_n^2  & \cdots & x_n^{n-1}

\end{pmatrix} \cdot \begin{pmatrix} a_0 \\ \vdots \\ a_k \\ \vdots \\ a_{n-1} \end{pmatrix} = \begin{pmatrix} f_1 \\ \vdots \\ f_i \\ \vdots \\ f_n \end{pmatrix} </math>

mit einer Vandermonde-Matrix als Koeffizientenmatrix. Aus der oben genannten Eigenschaft der Vandermonde-Determinante folgt daher insbesondere, dass das Interpolationsproblem genau dann eindeutig lösbar ist, wenn alle Stützstellen paarweise verschieden sind.

In der Standardbasis der Polynome ist die Matrix allerdings sehr schlecht konditioniert und die Auflösung mit Standardmethoden mit einer Laufzeit in <math>O(n^3)</math> recht teuer, weswegen man andere Darstellungen für die Polynome wählt. Näheres bei Polynominterpolation und unten.

Weitere Eigenschaften

Für paarweise verschiedene Stützstellen diagonalisiert die Vandermonde-Matrix <math>V</math> aus dem obigen Gleichungssystem die Begleitmatrix <math>C</math> des Polynoms

<math>p(x) = \prod_{i=1}^n(x-x_i) = b_0+b_1x+\ldots+b_{n-1}x^{n-1}+x^n</math>,

es gilt:

<math> V C V^{-1} = \operatorname{diag} (x_1, x_2, \dots, x_n)</math>.

Für große Anzahlen <math>n</math> kann man das Gleichungssystem oben auch über den folgenden Zusammenhang lösen, durch den die Inverse der Vandermonde-Matrix eng mit ihrer Transponierten verbunden ist. Mit den eingeführten Polynomkoeffizienten bildet man die Hankel-Matrix

<math>H:=\begin{pmatrix}

 b_1     & b_2 & \cdots & b_{n-1} & 1 \\
 b_2     & b_3 & \cdots & 1       & 0 \\
 \vdots  &     & .\cdot               \\ 
 b_{n-1} & 1   &        & 0           \\
 1       & 0   &        &         & 0

\end{pmatrix}</math> und die Diagonalmatrix <math>D:= \operatorname{diag}(p'(x_i))</math>. Wenn alle Stützstellen paarweise verschieden sind, ist <math>D</math> regulär. Damit gilt

<math> VHV^T=D,\quad V^{-1}=HV^TD^{-1}.</math>

Literatur

• Uwe Luther, Karla Rost: Matrix exponentials and inversion of confluent Vandermonde matrices. In: Electronic Transactions on Numerical Analysis. Bd. 18, 2004, {{#invoke:URIutil|{{#ifeq:1|1|linkISSN|targetISSN}}|1068-9613|0}}{{#ifeq:1|0|[!]

}}{{#ifeq:0|1

        |{{#switch:00
                  |11= (print/online)
                  |10= (print)
                  |01= (online)
          }}

}}{{#ifeq:0|0

        |{{#ifeq:0|0
              |{{#if:{{#invoke:URIutil|isISSNvalid|1=1068-9613}}
                    |
                    |{{#invoke:TemplUtl|failure|ISSN ungültig}}}}}}

}}, S. 91–100.

• Martin Hermann: Numerische Mathematik, Band 2: Analytische Probleme. 4., überarbeitete und erweiterte Auflage. Walter de Gruyter Verlag, Berlin und Boston 2020. ISBN 978-3-11-065765-4.

Einzelnachweise

<references/>

Weblinks

{{#if: Weisstein, Eric W. | Weisstein, Eric W. | Eric W. Weisstein }}: Vandermonde Matrix. In: MathWorld (englisch). {{#if: VandermondeMatrix | {{#ifeq: {{#property:P2812}} | VandermondeMatrix | | {{#if: {{#property:P2812}} | {{#ifeq: 0 | 0 | }} | {{#ifeq: 0 | 0 | }} }} }} }}