Masurarea unui segment de drum

Stiefel - lectii tinute la ETH

La măsurarea unui segment de drum, presupunem că am efectuat 5 măsurători

AD=89m, AC=67m, BD=53m, AB=35m şi CD=20m,

Stiefellsq.png

Initializare

A=[1 1 1

1 1 0

0 1 1

1 0 0

0 0 1];

b=[89

67

53

35

20];

x=A\b

x = 3x1 double
35.1250 32.5000 20.6250

Rezolvare

x=A\b

x = 3x1 double
35.1250 32.5000 20.6250

r=b-A*x

r = 5x1 double
0.7500 -0.6250 -0.1250 -0.1250 -0.6250

norm(r)

ans = 1.1726

Matricea ecuatiilor normale este

X=A'*A

X = 3x3 double
3 2 1 2 3 2 1 2 3

y=A'*b

y = 3x1 double
191 209 162

xn=X\y

xn = 3x1 double
35.1250 32.5000 20.6250

Conditionare

[cond(A),cond(X)]

ans = 1x2 double
3.1861 10.1515