nonlinsolve в SymPy дает мнимые решения, когда символы реальны - PullRequest
Новая
       50

nonlinsolve в SymPy дает мнимые решения, когда символы реальны

0 голосов
/ 18 апреля 2020

Все мои символы были созданы как настоящие (используя real=True). Если я использую nonlinsolve для решения x_1^2+1, я правильно не получаю никаких решений. Однако более сложная система уравнений приводит к мнимым решениям. Я неправильно понимаю, как работает nonlinsolve 

# Create real symbols x_1,x_2,...,x_18
syms = [Symbol('x_{}'.format(i), real=True) for i in range(1,19)]
for i,s in enumerate(syms):
    exec('x_{} = syms[{}]'.format(i+1, i))

# This code respects the real attribute
test = [Eq(eq) for eq in [x_1*x_1+1]]
print(nonlinsolve(test, x_1))

# This code violates the real attribute and gives imaginary solutions
test = [Eq(eq) for eq in [36*x_1 + 15*x_18**3 - 50*x_18**2 + 6*x_18 - 4, 15*x_18**3 - 50*x_18**2 + 114*x_18 + 36*x_2 - 40, -15*x_18**3 + 50*x_18**2 - 54*x_18 + 12*x_3 + 28, x_4 + 1, -x_18 + x_5, x_6 - 1, -15*x_18**3 + 50*x_18**2 - 42*x_18 + 36*x_7 + 4, 15*x_18**3 - 50*x_18**2 + 6*x_18 + 36*x_8 - 4, 15*x_18**3 - 50*x_18**2 + 54*x_18 + 12*x_9 - 16, x_10 - 3*x_18 + 2, x_11 - 1, x_12 + x_18 - 1, 36*x_13 + 15*x_18**3 - 50*x_18**2 + 114*x_18 - 76, x_14 + 1, 18*x_15 - 15*x_18**3 + 50*x_18**2 - 24*x_18 + 4, 6*x_16 + 15*x_18**3 - 50*x_18**2 + 54*x_18 - 22, x_17 + 2*x_18 - 1, 15*x_18**4 - 20*x_18**3 + 14*x_18**2 + 8*x_18 - 8]]
print(nonlinsolve(test, syms))

1 Ответ

0 голосов
/ 18 апреля 2020

Документально подтверждено, что nonlinsolve возвращает как реальные, так и комплексные решения: https://docs.sympy.org/latest/modules/solvers/solveset.html#nonlinsolve

Поскольку nonlinsolve является частью новых решателей, он работает с использованием наборов и предназначен для игнорировать любые предположения относительно входных символов. Вы можете получить только реальные решения, пересекаясь с набором действительных чисел: 

In [10]: nonlinsolve(test, syms) & Reals ** 18                                                                                    
Out[10]: 
⎧⎛2   √10  5   2⋅√10                √10       2   4⋅√10  2   √10                3⋅√10         √10  8   2⋅√10        4   √10     
⎨⎜─ - ───, ─ - ─────, -4 + √10, -1, ───, 1, - ─ + ─────, ─ - ───, 3 - √10, -2 + ─────, 1, 1 - ───, ─ - ─────, -1, - ─ + ───, 7 -
⎩⎝3    15  3     3                   5        3     15   3    15                  5            5   3     3          3    3      

            2⋅√10  √10⎞  ⎛√10   2  5   2⋅√10                -√10        4⋅√10   2  √10   2                3⋅√10     √10      2⋅√
 2⋅√10, 1 - ─────, ───⎟, ⎜─── + ─, ─ + ─────, -4 - √10, -1, ─────, 1, - ───── - ─, ─── + ─, 3 + √10, -2 - ─────, 1, ─── + 1, ───
              5     5 ⎠  ⎝ 15   3  3     3                    5           15    3   15   3                  5        5         3

10   8        4   √10                 2⋅√10  -√10 ⎞⎫
── + ─, -1, - ─ - ───, 2⋅√10 + 7, 1 + ─────, ─────⎟⎬
     3        3    3                    5      5  ⎠⎭
...