qcp。R

# 2021年版权,Gurobi优化狗万app足彩,LLC # #这个例子中制定和解决以下简单的QCP模型:#最大化# #话题# x + y + z = 1 # x ^ 2 + y ^ 2 z ^ < = 2(二阶锥)# x ^ 2 < = yz(旋转二阶锥)# x, y, z非负库(Gurobi)库(矩阵)模型< -列表()模型＄一个<-矩阵(c(1,1,1)， nrow=1, byrow=T)模型＄Modelsense <- 'max'模型＄Obj <- c(1,0,0)模型＄RHS <- c(1)模型＄第一个二次约束:x^2 + y^2 - z^2 <= 0 qc1 <- list() qc1＄Qc < - spMatrix (3 3 c(1、2、3),c(1、2、3),c (1.0, 1.0, -1.0)) qc1＄rhs <- 0.0 #第二个二次约束:x^2 - yz <= 0 qc2 <- list() qc2＄Qc < - spMatrix (3 3 c(1、2),c (1,3), c (1.0, -1.0)) qc2＄RHS <- 0.0模型＄Quadcon <- list(qc1, qc2) result <- gurobi(model) print(result＄objval)打印(结果＄x) #清除空间rm(模型，结果)