Question

from scipy.optimize import fsolve

x1=input("P1")
y1=input("rho1")
v1=input("velocity1")


def eqn(x):                              #three jump condition equations

    f1=(x[1]*x[2])-(y1*v1)
    f2=x[0]+(0.5*(y1**2)*(v1**2)/x[1])-x1-(0.5*y1*v1*v1)
    f3=(0.5*(y1**2)*(v1**2)/(x[1]**2))+(2.5*(x[0]/x[1]))-(0.5*v1*v1)-(2.5*(x1/y1))

    return(f1,f2,f3)

y=[10,11,12]         #initial guess
print fsolve(eqn,y)

`我必须解决三个非线性方程组，也称为冲击波中的跳跃条件。我被要求使用scipy.optimize这样做，我们采用三个已知参数并使用它们来查找其余三个参数。但是我没有得到正确的结果，并且根据初始猜测值而变化很大。

FAQ

Answer 1

除了Python特有的细节之外，你的方程似乎是正确的。解决方案非常简单：将所有输入转换为浮动。

x1=float(input("P1"))
y1=float(input("rho1"))
v1=float(input("velocity1"))

问题是x1/y1表达式中的f3是用整数除法（截断结果）而非浮点除法

这意味着您可能使用整数值进行测试。给原始代码浮点格式输入（例如5.0, 14.0, 3.0）也按预期工作。

冲击波跳跃条件

1 个答案: