Solve wont return obvious solution
Show older comments
I am trying to solve a symbolic system of nonlinear equations using
solve()
but it will not return obvious solution.
I try to run
sol = solve(K*x0)
with K and x0 being as follows:
x0 =
x1
y1
z1
x2
y2
z2
x3
y3
z3
K =
[ 21216.21*(x1 - 1.0*x2)^2, 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), -21216.21*(x1 - 1.0*x2)^2, -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), 0, 0, 0]
[ 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), 21216.21*(y1 - 1.0*y2)^2, 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), -21216.21*(y1 - 1.0*y2)^2, -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), 0, 0, 0]
[ 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), 21216.21*(z1 - 1.0*z2)^2, -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), -21216.21*(z1 - 1.0*z2)^2, 0, 0, 0]
[ -21216.21*(x1 - 1.0*x2)^2, -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), 21216.21*(x1 - 1.0*x2)^2 + 21216.21*(x2 - 1.0*x3)^2, 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2) + 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2) + 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), -21216.21*(x2 - 1.0*x3)^2, -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3)]
[ -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), -21216.21*(y1 - 1.0*y2)^2, -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2) + 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), 21216.21*(y1 - 1.0*y2)^2 + 21216.21*(y2 - 1.0*y3)^2, 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2) + 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), -21216.21*(y2 - 1.0*y3)^2, -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3)]
[ -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), -21216.21*(z1 - 1.0*z2)^2, 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2) + 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2) + 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), 21216.21*(z1 - 1.0*z2)^2 + 21216.21*(z2 - 1.0*z3)^2, -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), -21216.21*(z2 - 1.0*z3)^2]
[ 0, 0, 0, -21216.21*(x2 - 1.0*x3)^2, -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), 21216.21*(x2 - 1.0*x3)^2, 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3)]
[ 0, 0, 0, -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), -21216.21*(y2 - 1.0*y3)^2, -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), 21216.21*(y2 - 1.0*y3)^2, 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3)]
[ 0, 0, 0, -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), -21216.21*(z2 - 1.0*z3)^2, 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), 21216.21*(z2 - 1.0*z3)^2]
But matlab will not return the 0 vector as a solution, which is obvious here.
Answers (0)
Categories
Find more on Systems of Nonlinear Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!