26.12.2013 Aufrufe

Einführung in die Numerische Mathematik - Lehrstuhl Numerische ...

Einführung in die Numerische Mathematik - Lehrstuhl Numerische ...

Einführung in die Numerische Mathematik - Lehrstuhl Numerische ...

MEHR ANZEIGEN
WENIGER ANZEIGEN

Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.

YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.

3.5 <strong>Numerische</strong> Integration<br />

mit e<strong>in</strong>em Polynom p(x) möglichst exakt <strong>in</strong>tegrieren. Die Funktion ω(x) nennen wir Gewichtsfunktion.<br />

Diese Gewichtsfunktion fließt nun <strong>in</strong> das Konstruktionspr<strong>in</strong>zip von speziellen<br />

Gauß-Quadraturformeln e<strong>in</strong>. Zunächst formulieren wir den hierfür zentralen Satz:<br />

Satz 3.66 (Gewichtetes Skalarprodukt). Durch ω ∈ L ∞ ([−1, 1]) sei e<strong>in</strong>e positive<br />

ω(x) > 0<br />

fast überall,<br />

Gewichtsfunktion gegeben. Dann ist durch<br />

e<strong>in</strong> Skalarprodukt gegeben.<br />

Beweis: Übung!<br />

(f, g) ω :=<br />

∫ 1<br />

−1<br />

ω(x) f(x) g(x) dx,<br />

Die bisher bewiesenen Sätze 3.56, 3.57, 3.58, 3.59 sowie 3.61 gelten alle auch bei Verwendung<br />

<strong>die</strong>ses gewichteten Skalarprodukts. Das Übertragen der Beweise überlassen wir als<br />

Übung. Die Gauß-Tschebyscheff Quadratur wählt als spezielles Gewicht <strong>die</strong> Funktion<br />

ω(x) =<br />

1<br />

√<br />

1 − x 2 .<br />

Diese Funktion legt e<strong>in</strong> starkes Gewicht auf <strong>die</strong> beiden Intervallenden. Es gilt der folgende<br />

zunächst erstaunliche Satz:<br />

Satz 3.67 (Tschebyscheff-Polynome). Die Tschebyscheff-Polynome T n ∈ P n mit<br />

T n (x) := cos(n arccos x), −1 ≤ x ≤ 1,<br />

s<strong>in</strong>d <strong>die</strong> bezüglich des (·, ·) ω Skalarprodukts orthogonalisierten Monome {1, x, . . . } mit der<br />

Gewichtsfunktion<br />

1<br />

ω(x) = √ . 1 − x 2<br />

Die n paarweise verschiedenen Nullstellen des n-ten Tschebyscheff-Polynoms T n (x) s<strong>in</strong>d<br />

gegeben durch<br />

( ) 2k + 1<br />

x k = cos<br />

2n π , k = 0, . . . , n − 1.<br />

Beweis: (i) Wir müssen zunächst nachweisen, dass durch T n (x) überhaupt Polynome<br />

gegeben s<strong>in</strong>d. Wir führen den Beweis <strong>in</strong>duktiv durch Herleiten e<strong>in</strong>er Iterationsformel. Es<br />

gilt:<br />

T 0 (x) = 1, T 1 (x) = x,<br />

also <strong>in</strong>sbesondere T 0 ∈ P 0 und T 1 ∈ P 1 . Aus dem Additionstheorem für <strong>die</strong> Kos<strong>in</strong>usfunktion<br />

cos((n + 1)y) + cos((n − 1)y) = 2 cos(nt) cos(y),<br />

□<br />

91

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!