23得票4回答
我如何在Julia中高效地计算二次型?

我想计算一个二次形式:x' Q y,使用Julia编程语言。 以下是在不同情况下最有效的计算方式: 没有假设。 Q 是对称的。 x 和 y 相同 (x = y)。 同时满足 Q 对称和 x = y。 我知道Julia有dot()函数。但我想知道它是否比BLAS调用更快。

17得票1回答
如何将二次规划转化为线性规划?

我有一个优化问题,目标函数中有两个变量相乘,使得模型成为二次型。 我目前正在使用zimpl来解析该模型,并使用glpk来求解它。由于它们不支持二次规划,因此我需要将其转换为混合整数线性规划问题(MILP)。 第一个变量是实数,在[0,1]范围内,第二个变量是实数,在0到正无穷的范围内。这个...

14得票4回答
CVXOPT QP求解器:TypeError:'A'必须是具有1000列的'd'矩阵

我试图使用CVXOPT qp求解器来计算支持向量机的拉格朗日乘数def svm(X, Y, c): m = len(X) P = matrix(np.dot(Y, Y.T) * np.dot(X, X.T)) q = matrix(np.ones(m) * ...

11得票4回答
MATLAB:查找缩写版本的矩阵,使矩阵元素之和最小化

我有一个 151-by-151 的矩阵 A。这是一个相关系数矩阵,因此主对角线上有 1,在主对角线上方和下方有重复的值。每行/列代表一个人。 给定整数 n,我将寻求通过剔除一些人来减小矩阵的大小,从而留下一个最小化元素总和的 n-by-n 相关系数矩阵。除了获得缩写矩阵外,我还需要知道应该从...

7得票1回答
CVXPY抛出“SolverError”异常的具体原因是什么?

我正在使用CVXPY(版本1.0)解决一个二次规划(QP),但通常会出现以下异常:SolverError: 求解器“xxx”失败。请尝试另一个求解器。这使得我的程序非常脆弱。我已经尝试了不同的求解器,包括CVXOPT、OSQP、ECOS、ECOS_BB、SCS等。它们都或多或少存在相同的问题。...

7得票4回答
使用Numpy在Python中进行二次规划?

我正在将一些 MATLAB 代码翻译成 Python。有一行让我有些困惑: [q,f_dummy,exitflag, output] = quadprog(H,f,-A,zeros(p*N,1),E,qm,[],[],q0,options); 我查阅了MATLAB的文档,发现quadpr...