Mathematik online lernen im Mathe-Forum. Nachhilfe online

Startseite » Forum » Ausgleichsebene für 4 gemessene xyz Punkte

Ausgleichsebene für 4 gemessene xyz Punkte

Universität / Fachhochschule

Tags: Ausgleichsebene, Lineare Regression, Punktwolke

swali aktiv_icon

18:07 Uhr, 27.12.2016

Hallo Zusammen,

ich habe einen Z-Taster an einer CNC-Anlage mit dem ich 4 "XYZ-Punkte" im Raum abtaste. Die Positionen für

X

und

Y

sind fest definiert. Position

Z

wird jeweils ausgemessen. Mir ist bewusst, dass man zur Bestimmung einer Ebene nur 3 Punkte im Raum brauch. Ich brauche jedoch Aufgrund der Genauigkeit, die Ausgleichsebene der 4 Punkte. Ich habe hier auch schon einen Lösungsansatz gefunden, jedoch ist mir noch nicht klar, wie der genaue Rechenweg dahin aussieht?

Z = a 0 + a 1 x 1 + a 2 x 2

Wie genau komme ich auf die 4 neu errechneten "Z-Werte"? Die Berechnung soll später über ein AHK-Script (AutoHotkey) automatisiert werden.
Für Hilfe wäre ich sehr dankbar.

Ausgemessene Z-Positionen:

P 1

X : - 15

Y : - 15

Z : - 5.335

P 2

X : - 15

Y : 15

Z : - 3.332

P 3

X : 15

Y : 15

Z : - 4.971

P 4 :

X : 15

Y : - 15

Z : - 6.802

Anhang:
- Grafik zur veranschaulichung (jedoch mit verdrehtem Vorzeichen bei

Z ...)

- Bisher gefundener Lösungsansatz

Für alle, die mir helfen möchten (automatisch von OnlineMathe generiert):
"Ich möchte die Lösung in Zusammenarbeit mit anderen erstellen."

pwmeyer aktiv_icon

18:39 Uhr, 27.12.2016

Hallo,

Du bildest eine

4 \times 4

-Matrix

A (4 :

Zahl der Messpunkte, 3: Zahl der zu bestimmenden Parameter. 1. Spalte: alle

1, 2

. Spalte: x-Werte, 3. Spalte: y-Werte.

Dann löst Du das Gleichungssystem

A^{T} \cdot A v = A^{T} \cdot z (z :

Z-Werte).
Die Lösung

v

enthält dann

a_{0}, a_{1}, a_{2}

Gruß pwm

anonymous

19:09 Uhr, 27.12.2016

Hallo
Oder...
Ich vermute, das hätte heissen wollen:

z = a_{0} + a_{1} \cdot x + a_{2} \cdot y

Wenn du definierst:
Abweichung

b =

z_ideal - z_gemessen

= a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i}

Fehlerquadratsumme:

F = \sum b^{2} = \sum {(a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i})}^{2}

und die Fehlerquadratsumme minimierst, dann bekommst du drei Gleichungen für deine drei Parameter

a_{0}, a_{1}, a_{2}

anonymous

19:16 Uhr, 27.12.2016

PS:
Für dein Zahlenbeispiel ergibt sich:

z = - 5.1100 - 0.051766666666 \cdot x + 0.0639 \cdot y

Roman-22

21:20 Uhr, 27.12.2016

Sei

z = a \cdot x + b \cdot y + c

die gesuchte Ebenengleichung.

Seien

X, Y

und

Z

die

4 \times 1

Vektoren mit den

x -, y -

und

z -

Koordinaten deiner vier Punkte.
Da deine

x

und

y

Koordinaten bereits zentriert sind und sich nicht ändern, kannst du die gesuchten Koeffizienten wie folgt errechnen:

Seit

D = | \begin{matrix} X X & X Y \\ X Y & Y Y \end{matrix} | = X X \cdot Y Y - {(X Y)}^{2}

Dabei soll zB

X Y

für das Skalarprodukt der beiden Vektoren

X

und

Y

stehen.

Nun gilt

a = \frac{| \begin{matrix} X Z & X Y \\ Y Z & Y Y \end{matrix} |}{D}

b = \frac{| \begin{matrix} X X & X Z \\ X Y & Y Z \end{matrix} |}{D}

c = \bar{Z}

und

\bar{Z}

ist der Mittelwert der z-Werte.

Zum Nachlesen zB
http//www.gfai.de/%7Egoetze/links/Approxim_2_5.pdf

swali aktiv_icon

17:40 Uhr, 29.12.2016

Vielen Dank für die Hilfe. Ich verstehe jetzt wie ich mit der Ebenengleichung auf meine 4 neuen Z-Werte kommen müsste.

Was ich nicht verstehe ist, wie ich jetzt auf die Koeffizienten

a, b

und

c

komme? Ich bin leider nicht so fit in Mathe. Wenn ich den Beitrag von "Roman-22" nehme und versuche "D" "a" und "b" zu errechnen, frage ich mich welche

X

und

Y

und

Z

Werte ich jetzt von all denen die mir zur Verfügung stehen nehmen soll?

Könnte mir jemand zur Verdeutlichung die Gleichung mit Zahlen vom Beispiel füllen, so dass ich den genauen Rechenweg bis zum neuen Z-Wert sehen kann?
Bis jetzt schon mal vielen Dank für alle Beiträge!!

Roman-22

18:26 Uhr, 29.12.2016

>

frage ich mich welche

X

und

Y

und

Z

Werte ich jetzt von all denen die mir zur Verfügung stehen nehmen soll?
Lies den Anfang meines Beitrags nochmals durch, dort steht es!

\to

"Seien $X, Y$ und $Z$ die 4×1 Vektoren mit den x−,y− und z− Koordinaten deiner vier Punkte."
UNd weiter unten steht auch noch, das zB XY für das Skalarprodukt aus den beiden Vektoren

X

und

Y

steht. Analog natürlich XX, XZ, etc.

Aber da sich bei dir die

x -

und y-Werte (also die

\pm 15)

ja nicht ändern, kann ich es dir noch ein weniger einfacher liefern.
Seien

z_{1}

bis

z_{4}

die z-Koordinaten deiner vier Punkte und zwar genau in der Reihenfolge, in der du ursprünglich die Punkte hier angegeben hast.
Also

z_{1}

is die z-Koordinate des Punktes

(- 15 / - 15 / z_{1}),

etc.

Dann ist

a = \frac{1}{60} \cdot (z_{3} + z_{4} - z_{1} - z_{2})

b = \frac{1}{60} \cdot (z_{2} + z_{3} - z_{1} - z_{4})

c = \frac{1}{4} \cdot (z_{1} + z_{2} + z_{3} + z_{4})

Das lässt sich sicher leicht in deinem Skript umsetzen und wenn du die Werte, die du hier gepostet hast, einsetzt, solltest auf das Ergebnis kommen, welches kreador schon genannt hatte, nämlich

a = - 0,0517 \bar{6}

b = 0,0639

c = - 5,11

Zur Erinnerung, es geht um die Gleichung

z = a \cdot x + b \cdot y + c

anonymous

18:29 Uhr, 29.12.2016

Ich hatte oben schon den Ansatz hingeschrieben:
Fehlerquadratsumme:

F = \sum {(a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i})}^{2}

Minimieren, das heisst die partiellen Ableitungen sind Null:
dF/da_0

= 0 = \sum 2 \cdot (a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i}) \cdot 1

\to

0 = n \cdot a_{0} + a_{1} \cdot \sum x + a_{2} \cdot \sum y - \sum z_{i}

(Hierbei ist

n

die Anzahl an Stützpunkten. In deinem Beispielfall also

n = 4)

dF/da_1

= 0 = \sum 2 \cdot (a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i}) \cdot x

\to

0 = a_{0} \cdot \sum x + a_{1} \cdot \sum x^{2} + a_{2} \cdot \sum x \cdot y - \sum x \cdot z_{i}

dF/da_2

= 0 = \sum 2 \cdot (a_{0} + a_{1} \cdot x + a_{2} \cdot y - z_{i}) \cdot y

\to

0 = a_{0} \cdot \sum y + a_{1} \cdot \sum x \cdot y + a_{2} \cdot \sum y^{2} - \sum y \cdot z_{i}

Jetzt mache dir klar, sämtliche folgende Größen sind bekannt:

n

\sum x

\sum y

\sum z_{i}

\sum x^{2}

\sum x \cdot y

\sum y^{2}

\sum x \cdot z_{i}

\sum y \cdot z_{i}

D . h

. du hast ein lineares Gleichungssystem mit nur noch den unbekannten Größen

a_{0}, a_{1}, a_{2}

und eben die drei Bestimmungsgleichungen:

0 = n \cdot a_{0} + a_{1} \cdot \sum x + a_{2} \cdot \sum y - \sum z_{i}

0 = a_{0} \cdot \sum x + a_{1} \cdot \sum x^{2} + a_{2} \cdot \sum x \cdot y - \sum x \cdot z_{i}

0 = a_{0} \cdot \sum y + a_{1} \cdot \sum x \cdot y + a_{2} \cdot \sum y^{2} - \sum y \cdot z_{i}

Wunderbar - einfach auflösen. Ich habe das mal für dich gemacht. (Aber bitte erwarte nicht, dass ich den ganzen Rechenweg hier in Tastendrücke ausführe)
Nenner=n*sum

x^{2} \cdot \sum y^{2} - \sum x^{2} \cdot \sum y \cdot \sum y - \sum x \cdot \sum x \cdot \sum y^{2} - n \cdot \sum x \cdot y \cdot \sum x \cdot y + 2 \cdot \sum x \cdot \sum y \cdot \sum x \cdot y

a_{0} = (\sum x^{2} \cdot \sum z_{i} \cdot \sum y^{2} - \sum z_{i} \cdot \sum x \cdot y \cdot \sum x \cdot y + \sum x \cdot \sum x \cdot y \cdot \sum y \cdot z_{i} - \sum x \cdot \sum x \cdot z_{i} \cdot \sum y^{2} - \sum x^{2} \cdot \sum y \cdot \sum y \cdot z_{i} + \sum y \cdot \sum x \cdot z_{i} \cdot \sum

x*y)/Nenner

a_{1} = (\sum x \cdot \sum y \cdot \sum y \cdot z_{i} - n \cdot \sum x \cdot y \cdot \sum y \cdot z_{i} - \sum x \cdot \sum z_{i} \cdot \sum y^{2} + \sum z_{i} \cdot \sum y \cdot \sum x \cdot y - \sum x \cdot z_{i} \cdot \sum y \cdot \sum y + n \cdot \sum x \cdot z_{i} \cdot \sum

y^2)/Nenner

a_{2} = (n \cdot \sum x^{2} \cdot \sum y \cdot z_{i} - \sum x^{2} \cdot \sum z_{i} \cdot \sum y - \sum x \cdot \sum x \cdot \sum y \cdot z_{i} - n \cdot \sum x \cdot z_{i} \cdot \sum x \cdot y + \sum x \cdot \sum z_{i} \cdot \sum x \cdot y + \sum x \cdot \sum y \cdot \sum

x*z_i)/Nenner

swali aktiv_icon

09:02 Uhr, 06.01.2017

Vielen Dank für die Hilfe! Das Abtasten der Oberfläche und Generieren der Ausgleichsebene klappt inzwischen über das Script.

swali aktiv_icon

09:02 Uhr, 06.01.2017

EDIT

Diese Frage wurde automatisch geschlossen, da der Fragesteller kein Interesse mehr an der Frage gezeigt hat.

1347253

1345527

	Status: nicht eingeloggt	Noch nicht registriert?