Блог

Я пытаюсь максимизировать служебную функцию с двумя ограничениями. Но я не уверен, где мне нужно исправить код. Я попытался сделать размеры матриц одинаковыми и указанными в коде

Я определил глобальные параметры, я указал ограничения

 theta1 = 1
theta2 = 0.1
big_gamma1 = 1
big_gamma2= 0.5
gamma1 = 0.5
gamma2 = 0.5
z2 = 0.7*np.ones(5)
z1 = np.linspace(0.01, 4, 5)
xi_1 = 1/z1
xi1 = xi_1.transpose
print(z1)

z2 = 0.07*np.ones((1,5))
print(z2)
z1 = np.linspace(0.01, 4, 5)
xi_1 = 1/z1
xi1 = xi_1.transpose

r1 = cp.Variable((5,1))
r2 = cp.Variable((5,1))
f1 = cp.Variable((1,5))
f2 = cp.Variable((1,5))

constraints = [f1   f2 ==1,
              f1@r1   f2@r2 == 1]

obj = cp.Problem(cp.Minimize((f1@r1   f2@xi1)@z1   f2@(r2 - xi1)@(z2)), [f1   f2 ==1,
              f1*r1   f2*r2 == 1])
print("obj is DCP:", obj.is_dcp())
  

Ошибка в коде.

 [[0.07 0.07 0.07 0.07 0.07]]
---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-21-b52595a9b4b9> in <module>
     13               f1@r1   f2@r2 == 1]
     14 
---> 15 obj = cp.Problem(cp.Minimize((f1@r1   f2@xi1)@z1   f2@(r2 - xi1)@(z2)), [f1   f2 ==1,
     16               f1*r1   f2*r2 == 1])
     17 print("obj is DCP:", obj.is_dcp())

~/opt/anaconda3/lib/python3.7/site-packages/cvxpy/expressions/expression.py in cast_op(self, other)
     44         """A wrapped binary operator that can handle non-Expression arguments.
     45         """
---> 46         other = self.cast_to_const(other)
     47         return binary_op(self, other)
     48     return cast_op

~/opt/anaconda3/lib/python3.7/site-packages/cvxpy/expressions/expression.py in cast_to_const(expr)
    465                         "Combine Expressions using atoms such as bmat, hstack, and vstack."
    466                     )
--> 467         return expr if isinstance(expr, Expression) else cvxtypes.constant()(expr)
    468 
    469     @staticmethod

~/opt/anaconda3/lib/python3.7/site-packages/cvxpy/expressions/constants/constant.py in __init__(self, value)
     43             self._sparse = True
     44         else:
---> 45             self._value = intf.DEFAULT_INTF.const_to_matrix(value)
     46             self._sparse = False
     47         self._imag = None

~/opt/anaconda3/lib/python3.7/site-packages/cvxpy/interface/numpy_interface/ndarray_interface.py in const_to_matrix(self, value, convert_scalars)
     48             return result
     49         else:
---> 50             return result.astype(numpy.float64)
     51 
     52     # Return an identity matrix.

TypeError: float() argument must be a string or a number, not 'builtin_function_or_method'