%Hallo,
%ich habe Folgendes Problem:
%Ich muss eigentlich nur noch beweisen, das der Ausdruck
>>(2*XM(1))^(2/3)+XM(2)^(2/3)-3^(2/3) %mit syms t und
>>XM =
-1/2*cos(t)*(-4+(3*sin(t)^2+1)^(1/6))
-sin(t)*(-1+(3*sin(t)^2+1)^(1/6));
%gleich null ergibt (simple usw. liefern leider keine Vereinfachung), aber %unabhängig von t muss Null rauskommen(klappt auch bei jedem eingestzten %t nur allgemein vereinfacht Matlab nicht)
%Danke im vorraus für eure %Vorschläge[color=darkred][/color][size=18][/size]
% Ortsvektor der Ellipse
vec_x=[a*cos(t);b*sin(t)]
% Tangentenvektor der Ellipse
vec_t=simplify(diff(vec_x,t)./sqrt(dot(diff(vec_x,t),diff(vec_x,t))))
% Normalenvektor der Ellipse
vec_n=simplify(diff(vec_t,t)./sqrt(dot(diff(vec_t,t),diff(vec_t,t))))
% Kruemmung der Ellipse
k=simplify((diff(vec_x(1),t).*diff(vec_x(2),t,2)-diff(vec_x(2),t).*diff(vec_x(1),t,2))/(sqrt(dot(diff(vec_x,t),diff(vec_x,t)))).^3)
% Evolute der Ellipse
vec_ev=simplify(vec_x+(1/k)*vec_n)
% Test, hier sollte 0 rauskommen
maple('simplify',(vec_ev(1)/(b*c))^(2/3)+(vec_ev(2)/(a*c))^(2/3)-1)
Vielleicht bekommt noch jemand diesen Test für allgemeine a und b hin?
_________________
Viele Grüße
nschlange
"Chuck Norris ejakuliert fluessigen Stahl!"
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.