Technische Universität München

Lineare Algebra
für Informatiker [MA 0901]

Übungsblatt 9

Tutorium
R1
T9.1 Gegeben sei der Polynomraum R[x]2 mit dem Skalarprodukt ⟨p, q⟩ = 0 p(x)q(x) dx und der
Unterraum W = ⟨1 + x2 ⟩.

(a) Bestimmen Sie eine Basis von W ⊥ .

(b) Bestimmen Sie mit dem Gram-Schmidt-Verfahren aus der Basis

p1 (x) = 1, p2 (x) = x, p3 (x) = x2

von R[x]2 eine Orthonormalbasis von R[x]2 .

1
Lösung T9.1:

(a) p ∈ W ⊥ ⇔ p(x) = ax2 + bx + c ∈ R[X]2 und ⟨p, q⟩ = 0 für alle q ∈ W , d.h. ⟨p, 1 + x2 ⟩ = 0
D E D E
0 = p, 1 + x2 = ax2 + bx + c, 1 + x2
Z1 Z1
2 2
= (ax + bx + c)(1 + x ) dx = ax4 + bx3 + (a + c)x2 + bx + c dx
0 0
1
a b a+c 3 b 2 a b a+c b
= x5 + x4 + x + x + cx = + + + +c
5 4 3 2 0
5 4 3 2
8 3 4
= a+ b+ c
15 4 3
⇒ a = − 45
32
b − 52 c      
⇒ p(x) = ax2 + bx + c = − 45
32
b − 5
2
c x 2
+ bx + c = − 45 2
32
x + x b + − 5 2
2
x + 1 c
n o
Damit ist − 45
32
x2 + x, − 52 x2 + 1 ein linear unabhängiges Erzeugendensystem, also Basis von W ⊥ .

(b)
1
p1 Z
q1 := = p1 , da ||p1 ||2 = 1 dx = 1
||p1 ||
0
v
Z1
u 1
1 uZ
1 1
q2′ := p2 − ⟨q1 , p2 ⟩ q1 = x − x dx = x − mit ||q2′ || = t x2 −x+ dx = √
u
2 4 12
0 0
q2′ √1
  √
q2 := = 12 x −
= 3(2x − 1)
||q2′ || 2
D E D√ E √
q3′ := x2 − 1, x2 · 1 − 3(2x − 1), x2 · 3(2x − 1)
Z1 Z1
2 2 1 1 1
=x − x dx − (6x − 3) 2x3 − x2 dx = x2 − − (6x − 3) = x2 − x +
3 6 6
0 0
q3′ √ 
2 1
 √
q3 := ′ = 180 x − x + = 5(6x2 − 6x + 1)
||q3 || 6
√ √
Es ist also {1, 3(2x − 1), 5(6x2 − 6x + 1)} eine Orthonormalbasis von R[X]2 .

T9.2

(a) Zu den Messwerten

i 1 2 3 4
xi −1 0 1 2
yi 1 0 1 2
soll eine Gerade der Form y(x) = α + βx so gelegt werden, dass die Abweichung
4
(y(xi ) − yi )2
X

i=1

minimal wird. Bestimmen Sie die optimalen Parameter α und β unter Verwendung der Normalglei-
chung.
2
(b) Lösen Sie das lineare Ausgleichsproblem ∥b − Ax∥ = min für
   
1 1 2
A = 1 1.0001 und b = 0.0001 .
   

1 1.0001 4.0001

Lösung T9.2: (a) Wir verwendenunser Rezept:

1 −1
⊤
1 0 
1 1 .
(1) Setze b = (1, 0, 1, 2) und A =  

1 2
4 2
 
⊤ ⊤ ⊤
(2) Wir lösen die Normalengleichung A A x = A b, wobei A A = 2 6 :

 1
 
−1
4 2 1 1 1 1   1 −2
0= 6 4 4/5
     
x = (A⊤ A)−1 A⊤ b = 2 6 −1 0 1 2 1 20 −2 4 4 = 2/5 .
2
4
(3) Die gesuchte Ausgleichsgerage lautet damit y = 5
+ 25 x.
(b) Mit δ = 0.0001 gilt

1 1
 
3 3 + 2δ
 
A = 1 1 + δ , b = (2, δ, 4 + δ)⊤ , A⊤ A = .
1 1+δ 3 + 2δ 3 + 4δ + 2δ 2

Für die Lösung ergibt sich

2
 
−1 
3 3 + 2δ 1 1 1
 
⊤ −1 ⊤
x = (A A) A b = δ 
3 + 2δ 3 + 4δ + 2δ 2

1 1+δ 1+δ
4+δ
2
 
1 3 + 4δ + 2δ 2 −3 − 2δ 1 1 1
   
= 2  δ 
2δ −3 − 2δ 3 1 1+δ 1+δ
4+δ
1 2δ + 2δ 2 −δ −δ  2 
 
1 2δ 2 1
     
= 2 δ = 2 2 = 1 .
2δ −2δ δ δ 4+δ 2δ 2δ

T9.3 Im R3 seien zwei Geraden

G = {a + λu | λ ∈ R} und H = {b + µv | µ ∈ R}

mit        
1 −2 1 2
a = 3 , b = −2 , u = 1 , v = 3
       

2 1 2 4
gegeben. Bestimmen Sie auf den Geraden G und H die beiden Punkte p = a+λ0 u ∈ G und q = b+µ0 v ∈ H
mit minimalem Abstand ||p − q||.

3
Lösung T9.3: Aus der Bedingung
!
λ !
||p − q|| = ||a + λu − b − µv|| = ||(u, −v) − (b − a)|| = min
µ

erhalten wir ein lineares Ausgleichsproblem mit Matrix

   
1 −2 −3
A = (u, −v) = 1 −3 und rechter Seite c := b − a = −5 .
   

2 −4 −1
 
λ
Die Lösung x = µ
von A⊤ Ax = A⊤ c erhalten wir mit
   
1 −2 ! ! −3 ! !
⊤ 1 1 2  6 −13 ⊤ 1 1 2   −10
A A= 1 −3 = und A c= −5 =

−2 −3 −4 −13 29 −2 −3 −4 25
2 −4 −1

mit Hilfe des Gauß-Algorithmus aus

! ! !
6 −13 −10 6 −13 −10 −1 3 5
→ →
−13 29 25 −1 3 5 5 20

zu µ = 4 und λ = 7. Dies entspricht den Punkten

     
1 1 8
p = a + λu = 3 + 7 1 = 10
     

2 2 16

und      
−2 2 6
q = b + µv = −2 + 4 3 = 10 .
     

1 4 17

4
 Zusätzliche Übungen

Z9.1 Schreiben Sie ein Programm, das bei Eingabe von n Zeitpunkten t = (t1 , . . . , tn )⊤ ∈ Rn und n
Messwerten b = (y1 , . . . , yn )⊤ ∈ Rn die Koeffizienten λ1 , λ2 der Ausgleichsgerade f (x) = λ1 + λ2 x ausgibt
und einen Plot mit der Punktwolke (t1 , y1 ), . . . , (tn , yn ) und dem Graph der Ausgleichsgeraden f zeichnet.
Testen Sie Ihr Programm mit t=(0:0.1:10)’; b=t+rand(101,1) .* sign(randn(101,1));

Lösung Z9.1: Der folgende Code erfüllt alle Wünsche:

1 function [ lambda ] = ausgleichsgerade ( t , b )
2 % bestimmt die ausgleichsgerade f und gibt einen plot mit der punktwolke
und
3 % dem graph der ausgleichsgerade aus
4 n = length ( t ) ;
5 A =[ ones (n ,1) , t ];
6 lambda =( A '* A ) \( A '* b ) ;
7 f =@( x ) lambda (1) + lambda (2) * x ;
8 hold on
9 plot (t ,b , 'o ') ;
10 ezplot (f , [ min ( t ) , max ( t ) ]) ;
11 hold off
12 end

Die Ausgleichungsfunktion erhalten wir durch

lambda= ausgleichsgerade( t,b );

f = @(x) lambda(1) + lambda(2) * x;

Wir erhalten mit der Eingabe t=(0:0.1:10)'; b=t+rand(101,1).* sign(randn(101,1)); beispiels-

weise (da es sich bei den Variablen t und b um zufällig erstellte Vektoren handelt) die Funktion

f = @(x) − 0.1145 + 1.0112 ∗ x

und den Plot

5
 λ(1)+λ(2) x

10

0 2 4 6 8 10
x

Z9.2 Gezeitenprognose: Messungen an einer Küste ergeben die Tabelle

t 0 2 4 6 8 10
h 1.0 1.6 1.4 0.6 0.2 0.8

für den Wasserstand h (Meter) zur Tageszeit t (Stunden). Schaetzen Sie unter der natuerlichen Annahme,
dass h(t) durch eine harmonische Schwingung
π π
h(t) = x1 + x2 cos t + x3 sin t
6 6
beschrieben wird, mittels linearer Ausgleichsrechnung ab, wie groß hmax und hmin sind.

Lösung Z9.2: Die “Modellschwingung” h(t) = x1 + x2 cos π6 t + x3 sin π6 t ist hier bestmöglich anzupassen.
Dazu lesen wir zunächst A ∈ R6×3 und b ∈ R6 aus der Datentabelle ab:
1 1 0
     
√ h1 1.0
1 3
1 2
 

 √2



 h2 


 1.6 

1 − 21 3
h3 1.4
     
   
A=
 2 
und b =   =  .
1 −1 0√ h4 0.6
     
     
 
1 − 2 − √23
1 h5 0.2
     
     
1 12 − 23 h6 0.8

Erfreulicherweise reduziert sich hierbei die Normalgleichung auf

    
6 0 0 x1 5.6
⊤ ⊤
   
0.8 
A Ax =  0 3 0   x2  =  √  = A b,
0 0 3 x3 3

6
 
14 4 1
√ 
woraus wir die Parameter x = , ,
15 15 3
3 ≈ (0.933, 0.267, 0.577) gewinnen. Für die Amplitude A der
Schwingung erhält man q
A= x22 + x23 ≈ 0.63596 .
Folglich erhalten wir bei Flut hmax = x1 + A ≈ 1.569 m sowie bei Ebbe hmin = x1 − A ≈ 0.297 m.

Z9.3 Bestimmen Sie die Ausgleichsgerade zu den 5 Messpunkten:

ti 1 2 3 4 5
yi 5,0 6,5 9,5 11,5 12,5

Fertigen Sie eine Zeichnung an.

Lösung Z9.3: Gesucht ist eine Ausgleichsgerade f (x) = λ1 + λ2 x. Analog zur Vorlesung definieren wir
   
1 1 5, 0
1 2  6, 5 
   
   
1
A :=  3 , sowie b :=  9, 5 
 
1 4 11, 5
   
1 5 12, 5

und loesen das lineare Gleichungssystem

   
1 1 5, 0
1
! 2  6, 5 
!  ! !  
λ1 1 1 1 1 1   λ1 1 1 1 1 1 
A⊤ A = A⊤ b ⇔
 
1 3 = 9, 5 
λ2 1 2 3 4 5   λ
2 1 2 3 4 5  
1 4 11, 5
 

1 5 12, 5
! ! ! ! !
5 15 λ1 45 λ1 3
⇔ = ⇔ = .
15 55 λ2 155 λ2 2

Die Ausgleichsgerade hat damit die Gleichung f (x) = 3 + 2x.

Aktuelle Informationen zu Vorlesung und Übungen finden Sie unter: [Link]

7