Как реализовать Python код для решения программы quadrati c?

Question

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

Я хочу минимизировать (1 / 4b) [(Π1-s) ' K (Π1-s) + (Π'1-t) 'K (Π'1-t)] - tr (KΠ) относительно Π.

Π - матрица nxn, а 1 обозначает все свой вектор. Также s и t - векторы размерности n, а b - фиксированный скаляр. Таким образом, единственное количество, которое изменяется, это Π, и для этой матрицы у нас есть ограничение, что все записи суммируют до 1.

Как бы я это сделал? Или, если это не правильное место, где я должен спросить это?

Answer 1

Сначала вам нужно express ваше уравнение в виде python кода. Необработанный python не так хорош в чистом вычислении чисел, поэтому вам стоит подумать о такой библиотеке, как Numpy, чтобы сделать тяжелую работу за вас.

Как только вы это сделаете. Вы можете попробовать использовать один из оптимизаторов , которые поставляются с scikit-learn package

Если домен ird странный (например, не непрерывный), попробуйте использовать HyperOpt пакет