Verfasst am: 09.11.2022, 16:12
Titel: Hyperbeln als Implicit Function: Schnittpunkte finden
Hallo zusammen,
ich benötige mal wieder eure Hilfe.
Ich habe mehrere Hyperbeln, welche als ImplicitFunctionLine vorliegen. (Erstellt mit fimplicit)
Die XY-Werte kann ich erfolgreich auslesen, allerdings sind diese immer unterschiedlich fein skaliert.
Nun schneiden sich alle Hyperbeln miteinander.
Ich möchte in einem bestimmten Bereich immer zwei Hyperbeln miteinander vergleichen und gucken an welchem Punkt sie sich schneiden. Für diesen Punkt interessiert ich dann der X und Y Wert.
Gerne würde ich den Abstand mit least mean square bestimmen, ist aber nur eine idee.
Ich verzweifle aktuell daran mit den FunctionLine Objekts umzugehen und diese zu verwenden. Mein Ansatz war, die XY-Werte der Lines zu bekommen, was auch funktioniert hat. Allerdings sind diese immer nur in dem sichtbaren Bereich mit unterschiedlicher Auflösung definiert und so ist es mir z.b. nicht möglich diese mit find auf einen Schnittpunkt zu untersuchen.
bitte poste deinen Code. Dann kann man viel besser nachvollziehen, was du machst.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
diese zeichne ich dann in meinen Plot. das Speichern der XY Daten war ein Ansatz um den Abstand zu bestimmen, funktioniert aber (wie oben beschrieben) nicht.
Jetzt bin ich wieder an dem Punkt, wo ich den Abstand der Hyperbeln zueinander bestimmen möchte, um entweder einen Schnittpunkt(e) oder ebenhalt den kleinsten Abstand zu finden. Am liebsten wie oben beschrieben.
ich würde es mit
lsqnonlin
versuchen.
Solltest du damit nicht zurecht kommen, gib bitte für alle beteiligten Variablen Beispielwerte vor. Das erleichtert das Ausführen und Experimentieren mit dem Code ungemein.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Wenn es nur um zwei Hyperbel den Fall ist, würde es vielleicht interessant, die Gleichung für jede Kurve in die Form umzuwandeln:
(x-x0)^2/a^2 - (y-y0)^2/b^2=1
und
(u-u1)^2/a^2 - (v-v1)^2/b^2=1
u,v sind Cartesian neuen Koordinaten, aus x,y mit dem Gl. berechnet (,wo der Winkel theta1 zwischen der Achse der ersten Hyperbel und der zweiten Hyperbel ist)
u=x*cos(theta1)+y*sin(theta1)
v=-x*sin(theta1)+y*cos(theta1)
in dem Fall, theta1=arctan( (yfocus2-yfocus1)/(xfocus2-xfocus1) )
Ich habe es nicht zum Test durchgeführen. Wenn es möglich für mich zu probieren , bitte die Daten senden
FP
ich würde es mit
lsqnonlin
versuchen.
Solltest du damit nicht zurecht kommen, gib bitte für alle beteiligten Variablen Beispielwerte vor. Das erleichtert das Ausführen und Experimentieren mit dem Code ungemein.
Grüße,
Harald
Hallo Harald,
danke für deine Antwort. Ich habe leider nicht die Optimization Toolbox zur Verfügung und kann diese auch nicht anschaffen.
Hier mal der gesamte Code mit entsprechenden Werten:
(Die Werte sind erstmal ohne Einheit, da es vorerst nur um den TDOA Algorithmus geht)
plot(e1(1),e1(2),'Marker','*','LineStyle','none') xlim([-1,11]);
ylim([-1,11]);
hold on
plot(e2(1),e2(2),'Marker','*','LineStyle','none') plot(e3(1),e3(2),'Marker','*','LineStyle','none') plot(e4(1),e4(2),'Marker','*','LineStyle','none')
Meine Idee wäre es, die Funktionen nach Y umzustellen und somit in meinem interessanten Bereich über einen entsprechenden X Vektor zu berechnen.
Allerdings scheitere ich da gerade sehr dran!
du kannst auch fminsearch nutzen. Allerdings geht da der Least Squares-Charakter verloren und es besteht eine höhere Gefahr, in ein lokales Minimum zu laufen.
Ich habe leider nicht die Optimization Toolbox zur Verfügung und kann diese auch nicht anschaffen.
Aus meiner Sicht ist es eine Kosten-Nutzen-Abwägung. In diesem Fall ist es mit überschaubarem Aufwand möglich, die Nutzung der Toolbox zu umgehen. In anderen Fällen wird das nicht so leicht möglich sein, und dann muss man (bzw. der Vorgesetzte) sich entscheiden, ob man (sehr) viel Zeit investiert oder Geld.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Einstellungen und Berechtigungen
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
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.
RSS