Einzelnen Beitrag anzeigen

SMO

Registriert seit: 20. Jul 2005
178 Beiträge
 
Delphi XE6 Professional
 
#3

AW: Determinante 3 Punkte Kreis

  Alt 19. Nov 2015, 17:13
Ich verstehe aber nicht warum das so ist, kann mir das einer erklären?
Das ist auch auf den ersten Blick nicht leicht zu verstehen, besonders wenn man noch keine Vektorrechnung hatte.
Die von Sherlock verlinkte Erklärung ist gut, aber nicht besonders anschaulich, finde ich.

Deine "isLeft"-Funktion sollte wohl eher "isRight" heißen, denn sie liefert False für den kleinen Test A(1|1), B(3|2), C(2|3), obwohl dieser Punkt C links von AB liegt. Es kommt darauf an, welche "Richtung" man der Geraden gibt. Interpretiert man sie als "von A nach B", dann liefert die Formel positive Werte wenn C links liegt, interpretiert man die Gerade aber als "von B nach A", dann liefert sie negative wenn C links liegt.

Die Formel ist im Prinzip nichts anderes als ein zweidimensionals Kreuzprodukt (Kreuzprodukt in R²).
Falls du das Kreuzprodukt nicht kennst und auch von Vektoren keine Ahnung hast, müsste ich viel zu weit ausholen... Aber wenn du Grundkenntnisse hast, lass es mich wissen und ich werde es erklären.

Ansonsten hier eine Zerlegung des Codes von mir, die helfen könnte:
Delphi-Quellcode:
function Kreuzprodukt2D(x0, y0, x1, y1: Single): Single; inline;
begin
  Result := x0 * y1 - x1 * y0;
end;

// PunktAufSeite: Bestimmt, auf welcher Seite einer Gerade ein Punkt ist.
// Der zu testende Punkt P ist (px,py), die Gerade geht durch die Punkte A (ax,ay) und B (bx,by)
// Rückgabewert: < 0 Punkt ist rechts, > 0 Punkt ist links, = 0 Punkt ist auf der Geraden
function PunktAufSeite(ax, ay, bx, by, px, py: Single): Single;
begin
  // 2D Kreuzprodukt der Vektoren von A nach B und von A nach P
  // = Kreuzprodukt(B-A, P-A)
  // = (bx-ax) * (py-ay) - (px-ax) * (by-ay)
  Result := Kreuzprodukt2D(bx - ax, by - ay, px - ax, py - ay);
end;
Und hier kannst du interaktiv mit dem 3D-Kreuzprodukt herumspielen.
Da das aber etwas unübersichtlich ist und dein Problem nur 2D ist, habe ich schnell ein kleines GeoGebra-Projekt zusammengeklickt, siehe Anhang. Einfach die Zip-Datei entpacken, Web-App starten, und im Menü rechts oben die .ggb-Datei öffnen (Drag & Drop der Datei ins Browserfenster geht auch). "uxv" ist dabei das Kreuzprodukt, d.h. der Wert der Formel, bevor auf Vorzeichen getestet wird (<0, >0).
Miniaturansicht angehängter Grafiken
punktaufwelcherseitegerade.png  
Angehängte Dateien
Dateityp: zip PunktAufWelcherSeiteGerade.zip (7,6 KB, 10x aufgerufen)

Geändert von SMO (19. Nov 2015 um 17:24 Uhr)
  Mit Zitat antworten Zitat