用python求解多元多项式方程组

Question

编辑：我从中得到的参考方程中包含了几个错误。我已经修好了。现在解决方案可能真的有意义了！

当两层流体在地形上流动时，存在一个数 取决于流动的相对大小 速度和流体中的波速。在

它们被称为“超临界”、“次临界”和“临界”（即 前两个我在这里称之为“额外关键”）。在

以下方程式定义了临界点之间的边界线 （h，U0）参数空间中的超临界行为：

我想消除d_1c（即我不在乎它是什么）并找到 这些方程的解。在

简化因素：

• 我只需要给定的答案
• 我不需要精确的解，只需要一个解的概要 曲线，所以这可以通过解析或数值求解。在
• 我只想在区域（h，U0）=（0，0）到（0.5，1）上绘制。在

我想用热情中可用的模块来解决这个问题 分发（numpy，scipy，sympy），但真的不知道该去哪里 开始。它的消除变量d1c真正混淆 我。在

以下是python中的方程式：

def eq1(h, U0, d1c, d0=0.1):
    f = (U0) ** 2 * ((d0 ** 2 / d1c ** 3) + (1 - d0) ** 2 / (1 - d1c - d0) ** 3) - 1
    return f

def eq2(h, U0, d1c, d0=0.1):
    f = 0.5 * (U0) ** 2 * ((d0 ** 2 / d1c ** 2) - (1 - d0) ** 2 / (1 - d1c - d0) ** 2) + d1c + (h - d_0)
    return f

我期望有一个解决方案，它有许多解决方案分支（不是 总是身体上的，但不用担心），看起来很粗糙 像这样：

我如何着手实施这一点？在

Answer 1

让我先来处理d1c的消去。假设您设法将第一个等式转换成d1c = f(U, h, d0)的形式。然后把它代入第二个方程，在U，h和{}之间有一定的关系。在d0固定的情况下，它为两个变量U和{}定义了一个方程，从中可以找到任何给定的h的{}。这似乎就是你所说的解决方案，基于你最后的草图。 坏消息是，从任何一个方程中得到d1c并不容易。好消息是你不需要。在

fsolve可以取一个方程组，比如两个依赖于两个变量的方程，然后给出解。在本例中：fix h（d0已经修复），并将其作为变量U0和d1c作为变量提供给d0。记录U0的值，对下一个值{}重复上述操作，依此类推。在

请注意，与@duffymo的建议相反，我建议使用fsolve，或者至少从它开始，只有当它耗尽蒸汽时才寻找其他解算器。在

一个可能的警告是，对于给定的h，U0您需要一个开始的猜测，fsolve需要一个开始的猜测，并且没有简单的方法告诉它收敛到解决方案分支之一。如果这是一个问题，请查看brentq解算器。在

另一种方法是观察您可以轻松地从系统中消除U0。这样，您将得到h和d1c的一个方程，为h的每个值求解{}，然后使用原始方程中的任一个来计算给定的d1c和{}。在

使用fsolve的示例：

>>> from scipy.optimize import fsolve
>>> def f(x, p):
...   return x**2 -p
... 
>>> fsolve(f, 0.5, args=(2,))
array([ 1.41421356])
>>> 

这里的args=(2,)是用来告诉fsolve的语法，如果f(x,2)=0，你真正想要解决的是0.5的值。在