Решение уравнения в гекко - PullRequest
Новая
       17

Решение уравнения в гекко

2 голосов
/ 18 октября 2019

Это часть моего кода: 

import gekko as GEKKO
for i in range(0,4,1):
    m = GEKKO()
    x = m.Var()
    e_i1 = R_n[0,0] * A_slice[i,0] + R_n[1, 0] * A_slice[i, 1]
    e_i2 = R_n[0, 1] * A_slice[i, 1] + R_n[1, 1] * A_slice[i, 1]
    r = get_determinante(R_n)
    d = P[0,i]
    m.Equation(e_i1*R_n[0,1]*x + e_i1*R_n[1,1]*A_slice[i,4] - e_i2*R_n[0,0]*x + e_i2*R_n[1,0]*A_slice[i,4] == r * d)
    m.solve()
    print(x)

однако он не печатает x.

Какие-либо причины?

1 Ответ

2 голосов
/ 18 октября 2019

Попробуйте напечатать значение x как: 

print(x.value)

Переменная x является типом переменной Gekko, а результаты загружаются в x.value в виде массива. Чтобы получить доступ к значению, вы можете напечатать x.value или x.value[0], чтобы получить только первый элемент массива.

Вы создаете новую модель Гекко, переменные и уравнения в каждом цикле. Вы также можете объявить массив значений x с помощью x=m.Array(m.Var,4) и решить все уравнения и переменные одновременно без цикла. Это может помочь ускорить ваш код. Если вы сделаете это, вам нужно будет напечатать каждое значение x как: 

print([x[i].value[0] for i in range(4)])

Вот отдельный пример с функциями массива: 

from gekko import GEKKO
m = GEKKO()
# variable array dimension
n = 3 # rows
p = 2 # columns
# create array
x = m.Array(m.Var,(n,p))
for i in range(n):
   for j in range(p):
      x[i,j].value = 2.0
      x[i,j].lower = -10.0
      x[i,j].upper = 10.0
# create parameter
y = m.Param(value = 1.0)
# sum columns
z = [None]*p
for j in range(p):
   z[j] = m.Intermediate(sum([x[i,j] for i in range(n)]))
# objective
m.Obj(sum([z[j]**2 + y for j in range(p)]))
# minimize objective
m.solve()
print(x)
...