Sympy Как установить значение каждого символа в списке - PullRequest
Новая
       41

Sympy Как установить значение каждого символа в списке

0 голосов
/ 03 мая 2019

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

a = symbols('a0:8')

b = symbols('b0:8')

c = symbols('c0:8')

d = symbols('d0:8')

и затем я создал 33 уравнения, которые состоят только из этих символов. Когда я пытаюсь решить для всех этих символов линейно. Я просто продолжаю получать пустой набор как результат.

Я собрал все символы в списке как: 

symbollist=[]
for i in a:
      symbollist.append(i)
for i in b:
      symbollist.append(i)
for i in c:
      symbollist.append(i)
for i in d:
      symbollist.append(i)
print(symbollist)

вывод этого: 

[a0, a1, a2, a3, a4, a5, a6, a7, b0, b1, b2, b3, b4, b5, b6, b7, c0, c1, c2, c3, c4, c5, c6, c7, d0, d1, d2, d3, d4, d5, d6, d7]

и мой полный список уравнений: 

[a0 - 3.0, a1 - 3.7, a2 - 3.9, a3 - 4.2, a4 - 5.7, a5 - 6.6, a6 - 7.1, a7 - 6.7, b0 - b1 + 2*c0 + 3*d0, b1 - b2 + 6*c1 + 27*d1, b2 - b3 + 2*c2 + 3*d2, b3 - b4 + 2*c3 + 3*d3, b4 - b5 + 2*c4 + 3*d4, b5 - b6 + 4*c5 + 12*d5, b6 - b7 + 6*c6 + 27*d6, 2*c0 - 2*c1 + 6*d0, 2*c1 - 2*c2 + 18*d1, 2*c2 - 2*c3 + 6*d2, 2*c3 - 2*c4 + 6*d3, 2*c4 - 2*c5 + 6*d4, 2*c5 - 2*c6 + 12*d5, 2*c6 - 2*c7 + 18*d6, a0 - a1 + b0 + c0 + d0, a1 - a2 + 3*b1 + 9*c1 + 27*d1, a2 - a3 + b2 + c2 + d2, a3 - a4 + b3 + c3 + d3, a4 - a5 + b4 + c4 + d4, a5 - a6 + 2*b5 + 4*c5 + 8*d5, a6 - a7 + 3*b6 + 9*c6 + 27*d6, 2*c0, 2*c7, b0 - 1, b7 + 0.67]

Можете ли вы сказать мне, как решить эти уравнения и получить значения каждого a, b, c, d?

1 Ответ

0 голосов
/ 06 мая 2019

Поможет следующее: 

from sympy import *
var('a:8 b:8 c:8 d:8')
eq=[a0 - 3.0, a1 - 3.7, a2 - 3.9, a3 - 4.2, a4 - 5.7, a5 - 6.6, a6 - 7.1, a7 - 6.7, b0 - b1 + 2*c0 + 3*d0, b1 - b2 + 6*c1 + 27*d1, b2 - b3 + 2*c2 + 3*d2, b3 - b4 + 2*c3 + 3*d3, b4 - b5 + 2*c4 + 3*d4, b5 - b6 + 4*c5 + 12*d5, b6 - b7 + 6*c6 + 27*d6, 2*c0 - 2*c1 + 6*d0, 2*c1 - 2*c2 + 18*d1, 2*c2 - 2*c3 + 6*d2, 2*c3 - 2*c4 + 6*d3, 2*c4 - 2*c5 + 6*d4, 2*c5 - 2*c6 + 12*d5, 2*c6 - 2*c7 + 18*d6, a0 - a1 + b0 + c0 + d0, a1 - a2 + 3*b1 + 9*c1 + 27*d1, a2 - a3 + b2 + c2 + d2, a3 - a4 + b3 + c3 + d3, a4 - a5 + b4 + c4 + d4, a5 - a6 + 2*b5 + 4*c5 + 8*d5, a6 - a7 + 3*b6 + 9*c6 + 27*d6, 2*c0, 2*c7, b0 - 1, b7 + 0.67]
reps={}
do=True
while do:
  do=False
  for i in eq:
   i=i.subs(reps)
   f=i.free_symbols
   if len(f)==1:
    reps[f.pop()]=solve(i)[0]
    do=True
for i in ordered(reps):
  print(i,reps[i].n(2))
for i in eq:
  if i.subs(reps).free_symbols:print(i)

Распечатка покажет решения для всех переменных, кроме d7, что должно означать, что уравнения переопределены.

...