Punkt in Pyramide: unterbestimmtes LGS

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

crysel

Forum-Newbie


	Beiträge: 6

	Anmeldedatum: 15.08.10

	Wohnort: ---

	Version: ---

Verfasst am: 16.09.2010, 12:47 Titel: Punkt in Pyramide: unterbestimmtes LGS

Hi,

ich muss nachweisen können, ob ein Punkt Teil einer Pyramide (rechteckige Grundfläche) ist. Die Grundfläche wird von A,B,C,D aufgespannt, die Spitze ist S. Im Folgenden sind das die Vektoren von 0 bis zu dem jeweiligen Punkt.
Ich muss also folgendes LGS lösen, wobei der Punkt (x,y,z) gegeben ist und nur die Parameter a,b,c,d,s fehlen:

(x,y,z) = a*A + b*B + c*C + d*D + s*S

Dies liefert mir 3 Gleichungen für 5 Unbekannte. Eine zusätzliche Gleichung ist:
a + b +c + d + s = 1 (damit der Punkt in der Pyramide liegt).

Somit also 4 Gleichungen für 5 Unbekannte. Eine weitere Bedingung ist, dass alle Parameter a...s > 0 sind. Dies bekomme ich aber nicht als 5. Gleichung entsprechend eingebaut, sodass es bei einem unterbestimmten LGS bleibt.

Code:

% Pyramide aufstellen

A = [0 0 0];
B = [10 0 0];
C = [0 10 0];
D = [10 10 0];
S = [5 5 10];

% Test-Punkt

X = [5 5 5 1];

% Matrix bilden
M = [A' B' C' D' S'];
M(4,:) = [1 1 1 1 1];

% Lösen mit Backslash-Operator

L = M\X'

Funktion ohne Link?

Das Problem ist, dass mit der Backslash-Operator irgendein Ergebnis liefert, das die 4 Gleichungen erfüllt, aber häufig nicht die letzte Bedingung (a...s > 0), also er sagt mir quasi, dass der Punkt nicht in der Pyramide liegt, obwohl er drin liegt!

Ich bräuchte entweder eine 5. Gleichung, um ein vollwertiges LGS zu bekommen oder eine andere Methode zur Berechnung der Parameter, die die genannte Bedingung in irgendeiner Weise erfüllt.

Thomas84

Forum-Meister


	Beiträge: 546

	Anmeldedatum: 10.02.10

	Wohnort: ---

	Version: ---

Verfasst am: 17.09.2010, 07:09 Titel:

Entweder du unterteilst die Pyramide in 2 Tetraeder und schaust ob P in einem von den beiden ist, oder du rechnest symbolisch, oder du verwendest lsqnonneg oder fmincon.
Ich glaub die erste Variante ist die einfachste.

crysel

Themenstarter

Forum-Newbie


	Beiträge: 6

	Anmeldedatum: 15.08.10

	Wohnort: ---

	Version: ---

Verfasst am: 17.09.2010, 10:06 Titel:

Vielen Dank! Hätte ich eigentlich auch selbst drauf kommen können. Very Happy

Die erste Methode mit 2 Tetraedern scheint zu funktionieren (nach ausgiebig getesteten Punkten). Very Happy


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.