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:
Y la ecuación a resolver es:
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:
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;
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)
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;
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
0