我在尝试使用solve.QP来解决二次规划的最优问题,希望最小化目标函数Q(x)=X^TDX,D是10X10的矩阵,X是10维待求权重
有两个约束:
1.权重之和等于1
2.权重均大于0
Dmat是系数矩阵
Dmat <-C3 #C3是10X10的系数矩阵
Amat是约束矩阵,约束A^TX=1
A <- matrix(c(1,1,1,1,1,1,1,1,1,1),ncol=1)
bvec <- c(1)
Amat <- t(A)
然后我运行Solve.QP
sol <- solve.QP(Dmat, Amat, bvec)
给出了错误
Error in rep(0, q) : invalid 'times' argument
求指教
我不知道我哪里做错了
相似问题