方程组求解，计算不出来

综合讨论区
1

请教各位大神，这个方程组怎么求解，Matlab求解出来x,y,z有10组解，julia怎么算不出来？

using SymPy
@syms x y z e₁ e₂  c₁ c₂ c₃  c₅ s₁ s₂ s₃ g₁ g₂
fx = x*(x - 1)*(c₁ + c₂*y*z - c₂*y - c₂*z - e₁ - e₂*y*z + e₂*y + 
e₂*z - g₁*y - g₂*z - s₁*y)
fy = y*(y - 1)*(c₃ + g₁*x - g₁ + s₁*x)
fz = -z*(z - 1)*(-c₅ + s₂ + s₃)
solve([fx, fy, fz], [x, y, z])
2

这个求解使用Maltab 的solve ,很容易算出来结果啊，julia怎么就算不出来呢？

3

你试过 Python 的 sympy 吗？SymPy.jl 实际是用 Python 的 sympy 求解的。我试了一下 Python，也是半天没解出来。

import sympy

x, y, z, e1, e2, c1, c2, c3, c5, s1, s2, s3, g1, g2 = sympy.symbols("x y z e1 e2  c1 c2 c3  c5 s1 s2 s3 g1 g2")

fx = x*(x - 1)*(c1 + c2*y*z - c2*y - c2*z - e1 - e2*y*z + e2*y + 
e2*z - g1*y - g2*z - s1*y)
fy = y*(y - 1)*(c3 + g1*x - g1 + s1*x)
fz = -z*(z - 1)*(-c5 + s2 + s3)
sympy.solve([fx, fy, fz], [x, y, z])
4

Python也算不出来，好奇怪啊，手动算都能很快算出来，为啥python算不出来呢

5

遇到过类似问题， sympy 处理稍复杂的计算就会卡住，但同样方程用 MATLAB 就可以很快搞定。

对于这个例子，由于已经因式分解，分类讨论能降低问题计算难度。

using SymPy
@syms x y z e₁ e₂  c₁ c₂ c₃  c₅ s₁ s₂ s₃ g₁ g₂
fxs = [x, x-1, c₁ + c₂*y*z - c₂*y - c₂*z - e₁ - e₂*y*z + e₂*y + e₂*z - g₁*y - g₂*z - s₁*y]
fys = [y, y-1, c₃ + g₁*x - g₁ + s₁*x]
fzs = [z, z-1, -c₅ + s₂ + s₃]
sols = []
for fx in fxs, fy in fys, fz in fzs
    sol = solve([fx, fy, fz], [x, y, z]; set=true)
    isempty(sol[1]) || push!(sols, sol)
end
sols

求解结果

10-element Vector{Any}:
 (Sym[x, y, z], PyObject {(0, 0, 0)})
 (Sym[x, y, z], PyObject {(0, 0, 1)})
 (Sym[x, y, z], PyObject {(0, 1, 0)})
 (Sym[x, y, z], PyObject {(0, 1, 1)})
 (Sym[x, y, z], PyObject {(1, 0, 0)})
 (Sym[x, y, z], PyObject {(1, 0, 1)})
 (Sym[x, y, z], PyObject {(1, 1, 0)})
 (Sym[x, y, z], PyObject {(1, 1, 1)})
 (Sym[x, y, z], PyObject {((-c₃ + g₁)/(g₁ + s₁), (c₁ - e₁)/(c₂ - e₂ + g₁ + s₁), 0)})
 (Sym[x, y, z], PyObject {((-c₃ + g₁)/(g₁ + s₁), (c₁ - c₂ - e₁ + e₂ - g₂)/(g₁ + s₁), 1)})

由于要求解一般的公式，x, y, z 以外的符号可以取任意值，第三条式子 fz = -z*(z - 1)*(-c₅ + s₂ + s₃) 可以直接解出 z=0,1。这 3*3*3 条等式里绝大多数是平凡的，只有两条稍微要点计算，但带入 z 后也只是简单的线性方程。

2 个赞