Aufgabe5.tex 2.6 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768
  1. \section*{Aufgabe 5}
  2. \subsection*{Teilaufgabe a}
  3. Eine Quadraturformel $(b_i, c_i)_{i=1, \dots, s}$ hat die Ordnung
  4. $p$, falls sie exakte Lösungen für alle Polynome vom Grad $\leq p -1$
  5. liefert.
  6. \subsection*{Teilaufgabe b}
  7. Die Ordnungsbedingungen, mit denen man zeigen kann, dass eine Quadraturformel
  8. mindestens Ordnung $p$ hat, lautet:
  9. \[\forall p \in \Set{1, \dots, p}: \sum_{i=1}^s b_i c_i^{q-1} = \frac{1}{q}\]
  10. \subsection*{Teilaufgabe c}
  11. \paragraph{Aufgabe} Bestimmen Sie zu den Knoten $c_1 = 0$ und $c_2 = \frac{2}{3}$ Gewichte, um eine Quadraturformel
  12. maximaler Ordnung zu erhalten. Wie hoch ist die Ordnung?
  13. \paragraph{Lösung}
  14. Nach VL kann bei Vorgabe von $s$ Knoten auch die Ordnung $s$ durch
  15. geschickte Wahl der Gewichte erreicht werden. Nach Satz 27 ist diese
  16. Wahl eindeutig.
  17. Also berechnen wir die Gewichte, um die Ordnung $p=2$ zu sichern.
  18. Dazu stellen wir zuerst die Lagrange-Polynome auf:
  19. \begin{align}
  20. L_1(x) &= \frac{x-x_2}{x_1 - x_2} = \frac{x-c_2}{c_1-c_2} = \frac{x-\nicefrac{2}{3}}{-\nicefrac{2}{3}} = -\frac{3}{2} x + 1\\
  21. L_2(x) &= \frac{x-x_1}{x_2 - x_1} = \frac{x-c_1}{c_2-c_1} = \frac{x}{\nicefrac{2}{3}} = \frac{3}{2} x
  22. \end{align}
  23. Nun gilt für die Gewichte:
  24. \begin{align}
  25. b_i &= \int_0^1 L_i(x) \mathrm{d}x\\
  26. b_1 &= \int_0^1 -\frac{3}{2} x + 1 \mathrm{d}x = \left [ -\frac{3}{4}x^2 + x \right ]_0^1 = \frac{1}{4}\\
  27. b_2 &= \frac{3}{4}
  28. \end{align}
  29. Nun sind die Ordnungsbedingungen zu überprüfen:
  30. \begin{align}
  31. \nicefrac{1}{1} &\stackrel{?}{=} b_1 c_1^0 + b_2 c_2^0 = \nicefrac{1}{4} + \nicefrac{3}{4} \text{\;\;\cmark}\\
  32. \nicefrac{1}{2} &\stackrel{?}{=} b_1 c_1^1 + b_2 c_2^1 = \frac{1}{4} \cdot 0 + \frac{3}{4} \cdot \frac{2}{3} \text{\;\;\cmark}\\
  33. \nicefrac{1}{3} &\stackrel{?}{=} b_1 c_1^2 + b_2 c_2^2 = \frac{1}{4} \cdot 0 + \frac{3}{4} \cdot \frac{4}{9} \text{\;\;\cmark}\\
  34. \nicefrac{1}{4} &\stackrel{?}{=} b_1 c_1^3 + b_2 c_2^3 = \frac{1}{4} \cdot 0 + \frac{3}{4} \cdot \frac{8}{27} \text{\;\;\xmark}\\
  35. \end{align}
  36. Die Quadraturformel mit den Knoten $c_1 = 0$, $c_2 = \nicefrac{2}{3}$ sowie
  37. den Gewichten $b_1 = \nicefrac{1}{4}$, $b_2 = \nicefrac{3}{4}$ erfüllt
  38. also die 1., 2. und 3. Ordnungsbedingung, nicht jedoch die 4.
  39. Ordnungsbedingung. Ihre maximale Ordnung ist also $p=3$.
  40. \textbf{Anmerkungen:} Da $c_1 = 0$ kann es sich nicht um die Gauß-QF handeln.
  41. Somit können wir nicht Ordnung $p=4$ erreichen.
  42. Bei der Suche nach den Gewichten hätte man alternativ auch das folgende
  43. LGS lösen können:
  44. \begin{align}
  45. \begin{pmatrix}
  46. c_1^0 & c_2^0\\
  47. c_1^1 & c_2^1
  48. \end{pmatrix}
  49. \cdot x
  50. =
  51. \begin{pmatrix}
  52. 1\\
  53. \nicefrac{1}{2}
  54. \end{pmatrix}
  55. \end{align}