Matlab - Problemas al resolver ecuación de elipse

   
Vista:

Problemas al resolver ecuación de elipse

Publicado por aol (2 intervenciones) el 01/04/2014 12:26:01
Hola a todos,

Estoy intentando dibujar una ecuación no lineal.

Mi parámetros de entrada son:
x1 = 109.0912;
y1 = 154.8336;
x2 = 198.0985;
y2 = 162.3950;
x3 = 205.7735;
y3 = 250.0836;
x4 = 116.7662;
y4 = 242.5223;
a = 135.7204;
b = 114.1723;
x = x1:0.1:x3;

Y la ecuación a resolver es:
syms y real
fun = (((y3-y1)*x-(x3-x1)*y-x1*y3+x3*y1)/sqrt((y3-y1)^2+(x3-x1)^2))^2/a^2+(((y4-y2)*x-(x4-x2)*y-x2*y4+x4*y2)/sqrt((y4-y2)^2+(x4-x2)^2))^2/b^2-1;
y = solve(fun,y);
plot(x,y)

Me da un error cuando lee la función de la elipse. ¿Por qué podría ser? ¿Se os ocurre otra manera de resolverla?

El error que me muestra es:
Error using mupadmex
Error in MuPAD command: Not a square matrix.
[(Dom::Matrix(Dom::ExpressionField()))::_power]

Error in sym/privBinaryOp (line 1694)
Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});

Error in sym/mpower (line 194)
B = privBinaryOp(A, p, 'symobj::mpower');

Error in centro_morf (line 158)
fun =
(((y3-y1)*x-(x3-x1)*y-x1*y3+x3*y1)/sqrt((y3-y1)^2+(x3-x1)^2))^2/a^2+(((y4-y2)*x-(x4-x2)*y-x2*y4+x4*y2)/sqrt((y4-y2)^2+(x4-x2)^2))^2/b^2-1;
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder