Solve wont return obvious solution
显示 更早的评论
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.
回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Systems of Nonlinear Equations 的更多信息
2016-6-5
2016-6-5
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
