Problem mit dem Definieren einer Funktion
octavo
Forum-Newbie
Beiträge: 7
Anmeldedatum: 11.06.12
Wohnort: ---
Version: ---
Verfasst am : 22.06.2012, 10:18
Titel : Problem mit dem Definieren einer Funktion
Und noch ein Problem, dann sollte allerdings mein Programm endlich fertig sein
Ich möchte eine Funktion definieren aus Werten die mir Matlab vorher berechnet hat:
k= z1*x.^2 + z2*y.^2 + b*x + c*y + e;
mein k soll allerdings zusätzlich gleich null sein (also ein Kegelschnitt!)
wenn ich den dann plotten will mit plot(k) weiß ich nicht, wie ich die Null da einbinde.
Wisst ihr, wie ich das hinbekomme?
Gruß
Christian
octavo
Themenstarter
Forum-Newbie
Beiträge: 7
Anmeldedatum: 11.06.12
Wohnort: ---
Version: ---
Verfasst am : 23.06.2012, 08:33
Titel :
Hey,
ich weiß nicht ob das jetzt tatsächlich so funktioniert wie es soll (da ich im Moment kein Matlab zur Hand hab!), aber falls noch jemand das Problem hat:
Mit
syms x y;
lhs = z1*x^2 + z2*y^2 + b*x + c*y + e;
ezplot(lhs);
Sollte das eigentlich lösbar sein
octavo
Themenstarter
Forum-Newbie
Beiträge: 7
Anmeldedatum: 11.06.12
Wohnort: ---
Version: ---
Verfasst am : 25.06.2012, 09:57
Titel :
Hey,
mein Programm funktioniert trotzdem nicht. Kann mir jemand helfen?
Bis zur "Kegelschnittgleichung funktioniert (glaube ich) alles. Dann sagt er mir:
??? Error using ==> char
Cell elements must be character arrays.
Error in ==> ezplot at 158
fmsg = char(f);
Error in ==> Untitled2 at 52
ezplot(lhs);
Code:
%
% Darstellung von Kegelschnitten mit 5 Punkten
%
% Punkte eingeben:
%
x1=input ( 'Bitte ersten Kegelpunkt eingeben: ') ;
x2=input ( 'Bitte zweiten Kegelpunkt eingeben: ') ;
x3=input ( 'Bitte dritten Kegelpunkt eingeben: ') ;
x4=input ( 'Bitte vierten Kegelpunkt eingeben: ') ;
x5=input ( 'Bitte fünften Kegelpunkt eingeben: ') ;
%
% Plückermethode
%
C1 = cross ( x1',x2') ;
C2 = cross ( x3',x4') ;
D1 = cross ( x1',x3') ;
D2 = cross ( x2',x4') ;
%
A1 = C1*C2';
A2 = C2*C1';
B1 = D1*D2';
B2 = D2*D1';
%
A = A1 + A2;
B = B1 + B2;
%
l = x5*B*x5';
m = -x5*B*x5';
%
% C ist Kegelschnitt-Matrix
%
C = l*A + m*B;
%
% Hauptachsen-Transformation
%
D=[ C( 1 ,1 ) C( 1 ,2 ) ; C( 1 ,2 ) C( 2 ,2 ) ] ;
values = eig ( D) ;
z1 = values( 1 ) ;
z2 = values( 2 ) ;
%
[ V,D] = eig ( D) ;
%
a=[ 2 *C( 3 ,1 ) 2 *C( 3 ,2 ) ] ;
b=a*V( :,1 ) ;
c=a*V( :,2 ) ;
%
% Kegelschnitt-Gleichung
%
e=C( 3 ,3 ) ;
syms x y;
lhs = 'z1*x^2 + z2*y^2 + b*x + c*y + e ';
ezplot ( lhs) ;
und mit fprintf hat er auch so seine Probleme. Dabei wurde mir der Tipp aus einem anderen Thread in diesem Forum gegeben. Ich bin verzweifelt.
Code:
%
% Unterscheidung in:
% Parabel
% Kreis
% Ellipse
% Hyperbel
%
if ( ( z1 == 0 ) && ( z2 ~= 0 ) )
plot ( lhs) ;
fprintf ( 'Parabel ') ;
fprintf ( 'Scheitelpunkt = ( %g, %g) ) ', ( -b/( 2 *z1) , -e+( b*b/4 *z1) ) ;
elseif ( ( z2 == 0 ) && ( z1 ~= 0 ) )
plot ( lhs) ;
disp ( 'Parabel ') ;
fprintf ( 'Scheitelpunkt = ( %g, %g) ) ', ( -c/( 2 *z2) , -e+( c*c/4 *z2) ) ;
elseif ( z1 == z2)
plot ( lhs) ;
fprintf ( 'Kreis ')
fprintf ( 'Mittelpunkt = ( %g, %g) ) ', ( '-b/( 2 *z2) ,-( c/2 *z2) ) ;
fprintf ( 'Radius = ( %g) ) ' ( sqrt ( -e+( b*b/4 *z1) +( c*c/4 *z2) ) ) ;
elseif ( ( z1*z2 > 1 ) && ( ( -e + ( b*b/2 *z1) ) >0 ) )
plot ( lhs) ;
fprintf ( 'Ellipse ')
fprintf ( 'Mittelpunkt = ( %g, %g) ) ', ( '-b/( 2 *z2) ,-( c/2 *z2) ] ) ;
fprintf ( 'Halbachse a = ( %g) ) ', ( z1/( -e+( b*b/4 *z1) +( c*c/4 *z2) ) ) ;
fprintf ( 'Halbachse b = ( %g) ) ', ( z2/( -e+( b*b/4 *z1) +( c*c/4 *z2) ) ) ;
elseif ( ( z1*z2 < 1 ) && ( ( b*b/2 *z1 + c*c/2 *z2) ~=0 ) )
plot ( lhs) ;
fprintf ( 'Hyperbel ')
fprintf ( 'Asymptoten = ( %g) ) ', ( y + b/2 *z1 '= ' ( c/d) *( x + d/2 *z2) ) ;
else
fprintf ( 'Der Kegelschnitt besitzt keine Achse ') ;
end
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
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.