- 数学建模与数学规划:方法、案例及编程实战(Python+COPT/Gurobi实现)
- 刘兴禄 赖克凡 杉数求解器COPT团队主编
- 595字
- 2024-11-28 16:21:10
1.2.6 二阶锥规划
二阶锥规划(Second-Order Cone Programming,SOCP)的目标函数为线性表达式,约束包含二阶锥约束,是一种非常特殊的非线性优化模型。在给出其一般形式之前,我们需要了解一下什么是(凸)锥,什么是二阶锥。
·锥(Cone):对于一个向量空间Rn与它的一个子集C,如果子集C中的任意一点x与任意正数α的积αx仍然属于子集C,则称C为一个锥。若C中任意两点x与y,以及任意两个正数α与β,都有αx+βy∈C,则C为凸锥。
·二阶锥(Second-order Cone):以二范数定义的锥被称为二阶锥。在k维空间中,标准的二阶锥数学定义为式(1.5)。图1.2为三维空间中的二阶锥示意图。
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_22.jpg?sign=1739384811-YltcvZslYawsAPo4zYK6czcAQ1MyRVrf-0-331f43d939fdb4c0f669e1437b3c400b)
需要说明的是,符号‖·‖2表示向量的二范数(L2 norm)。列向量x=[x1,x2,…,xn]T的二范数定义为:
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_23.jpg?sign=1739384811-i6tT0Bddk8MJiYac43Hg7aKCjau8HPMN-0-d4ed19cfe3cdac423a3934c642c0b267)
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_24.jpg?sign=1739384811-l4rdlag9rMhl7lqWgTmd1YurEjmbo0rl-0-5d53c1a7564cf2c20847423dceb33a11)
图1.2 三维空间中的二阶锥示意图
形如式(1.6)的约束即为二阶锥约束(Second-order Cone Constraint)。
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_25.jpg?sign=1739384811-ZAvbAIdjlBGmXAE1K0nStXrl8GllhSYx-0-3ef9f69471c0a34896a251eee099fcae)
其中,A∈Rk×n,表示系数矩阵;x∈Rn×1,为列向量,是决策变量;b∈Rk×1,为列向量;c∈Rn×1,为列向量;d为常数。
二阶锥规划的一般形式[16]如下:
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_26.jpg?sign=1739384811-6rp7IxfNJtgqLh5jSwdJLSPHiNdzz8q3-0-3b3b1b45a8ad6a104bf787127c0f5070)
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_27.jpg?sign=1739384811-HskRkTrylum47n3zvDwA9FhpRL26lRzC-0-65ade12aae3526c583fc4db948a01142)
其中,f∈Rn×1,为列向量;x∈Rn×1,为列向量,表示连续型决策变量;;
;ci∈Rn×1,为列向量;di∈R;F∈Rh×n,表示约束系数矩阵;g∈Rh×1,为列向量,表示右端常数。
当ci=0(∀i=1,…,m)时,SOCP可以等价转换为QCQP。当Ai=O(即零矩阵)(∀i=1,…,m)时,SOCP退化为LP。
下面给出两个二阶锥规划的简单例子。
【例1.1】
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_30.jpg?sign=1739384811-9ZfQoInzZuOjcF1itcyA7uds1ITIet4j-0-64d65cb17f124ae4a50d86e3ae909af3)
若表示成紧凑的矩阵形式,则
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_31.jpg?sign=1739384811-aLJfJf32dyZN4tWxFXBNxQY0QYPu3pEy-0-5155e479a40ae08b8deb20e7d5cbf67a)
【例1.2】
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_32.jpg?sign=1739384811-Th0AAJj4JuYWcCYZnHSZep6BikF5WA2C-0-90e189d41e0b8b57fff9c3e19690d546)
若表示成紧凑的矩阵形式,则
![](https://epubservercos.yuewen.com/0DD641/31155568907421606/epubprivate/OEBPS/Images/txt002_33.jpg?sign=1739384811-nUG7eN8w7f10RkxI8wLqbfi7NUPGfGGm-0-24692f9cc5dfcb50c65e02d0dbe8bda7)