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

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

  
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Schnittline zweier Objekte
 

StupidFox
Forum-Newbie

Forum-Newbie

Beiträge: 4
Anmeldedatum: 03.11.19
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 05.11.2019, 19:27     Titel: Schnittline zweier Objekte
  Antworten mit Zitat      
Guten Abend!

Ich versuche gerade, einem Matlab Plot zweier Objekte (Cone und Plane) eine Schnittlinie hinzuzufügen. Leider scheine ich irgendwas falsch zu machen. Zunächst mal das Skript, welches mir die dreidimensionalen Objekte erzeugt:

Code:
R = 2;
hold on;

% Plane
xPlane = @(u,v) u;
yPlane = @(u,v) v;
zPlane = @(u,v) ((3/2) * R) - (3/4) * u;
fsurf(xPlane,yPlane,zPlane, [-2*R,2*R]);

% Cone
xCone = @(s,t) (R - 1/3 * t) .* cos(s);
yCone = @(s,t) (R - 1/3 * t) .* sin(s);
zCone = @(s,t) t;
fsurf(xCone, yCone, zCone, [0, 2*pi, 0, 3*R]);

hold off;


Das Funktioniert super. Nun möchte ich allerdings die Schnittkante der beiden Objekte (Ellipse) farblich hervorheben; also einen eigenen Plot3 dafür erstellen. Dafür wird dann ein bisschen Mathe gebraucht und man landet schließlich bei folgender Formel für y in Abhängigkeit der x-Werte:



Jetzt habe ich schon diverse Dinge versucht, das als Schnittlinie in meinen Plot zu integrieren, aber ich komme auf keinen grünen Zweig. Folgendes zum Beispiel ist falsch:

Code:

R = 2;
hold on;

% Plane
xPlane = @(u,v) u;
yPlane = @(u,v) v;
zPlane = @(u,v) ((3/2) * R) - (3/4) * u;
fsurf(xPlane,yPlane,zPlane, [-2*R,2*R]);

% Cone
xCone = @(s,t) (R - 1/3 * t) .* cos(s);
yCone = @(s,t) (R - 1/3 * t) .* sin(s);
zCone = @(s,t) t;
fsurf(xCone, yCone, zCone, [0, 2*pi, 0, 3*R]);

% Intersection line
xIntersection = linspace(-R, R);
yIntersection1 = + sqrt(1/4 * R.^2 + 1/4 * R * xIntersection + 15/16 * xIntersection.^2);
yIntersection2 = - sqrt(1/4 * R.^2 + 1/4 * R * xIntersection + 15/16 * xIntersection.^2);
zIntersection = 3/2 .* R - 3/4 .* xIntersection; % PLANE!
plot3(xIntersection,yIntersection1,zIntersection,'linewidth',5,'col','y');
plot3(xIntersection,yIntersection2,zIntersection,'linewidth',5,'col','y');

hold off;
 


Eventuell kann mir jemand einen Tipp geben. Vielen Dank Smile
Private Nachricht senden Benutzer-Profile anzeigen


Neues Thema eröffnen Neue Antwort erstellen



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 | goMatlab RSS Button 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.