锥规划问题存储格式——CBF格式简介

2023-11-09 17:20

本文主要是介绍锥规划问题存储格式——CBF格式简介,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

锥规划问题存储格式——CBF格式简介

文章目录

  • 锥规划问题存储格式——CBF格式简介
  • 前言
  • 1 问题形式
  • 2 CBF文件结构
  • 3 实例
  • 参考文献


前言

CBF(Conic Benchmark Format),这是一种文件存储格式,整合了混合整数变量在混合圆锥(包括线性锥、二阶锥、半定锥、指数锥和幂锥)上的优化问题实例。该格式在设计时考虑了锥规划问题的基准库,因此侧重于紧凑且易于解析的表示。


1 问题形式

要使用CBF储存格式,首先要将凸优化问题化为CBLIB中使用的原始问题的标准参考形式:
minimize ⁡ x , X c T x + ⟨ C , X ⟩ subject to  A x + F ( X ) − b ∈ K g n g , H ∗ ( x ) − B ∈ S + n G , x ∈ K x n x , X ∈ S + n X , and  x j ∈ Z for  j ∈ I , \begin{aligned} \underset{x, X}{\operatorname{minimize}} & c^T x+\langle C, X\rangle \\ \text { subject to } & A x \quad+\mathcal{F}(X)-b \in \mathcal{K}_g^{n_g}, \\ & \mathcal{H}^*(x) \quad-B \in \mathcal{S}_{+}^{n_G}, \\ & x \in \mathcal{K}_x^{n_x}, X \in \mathcal{S}_{+}^{n_X}, \text { and } x_j \in \mathbb{Z} \text { for } j \in \mathcal{I}, \end{aligned} x,Xminimize subject to cTx+C,XAx+F(X)bKgng,H(x)BS+nG,xKxnx,XS+nX, and xjZ for jI,
其中 F ( X ) \mathcal{F}(X) F(X) H ∗ ( x ) \mathcal{H}^*(x) H(x)定义如下:
F ( X ) = [ ⟨ F 1 , X ⟩ ⋮ ⟨ F n g , X ⟩ ] , H ∗ ( x ) = ∑ j = 1 n x x j H j \mathcal{F}(X)=\left[\begin{array}{c} \left\langle F_1, X\right\rangle \\ \vdots \\ \left\langle F_{n_g}, X\right\rangle \end{array}\right], \quad \mathcal{H}^*(x)=\sum_{j=1}^{n_x} x_j H_j F(X)= F1,XFng,X ,H(x)=j=1nxxjHj
C ∈ S n X , B ∈ S n G , F i ∈ S n X C \in \mathcal{S}^{n_X}, B \in \mathcal{S}^{n_G}, F_i \in \mathcal{S}^{n_X} CSnX,BSnG,FiSnX for i = i= i= 1 , … , n g 1, \ldots, n_g 1,,ng, and H j ∈ S n G H_j \in \mathcal{S}^{n_G} HjSnG for j = 1 , … , n x j=1, \ldots, n_x j=1,,nx.

Note:

  1. 上面的记号与大多数锥规划文献中所用记号类似,故不在此赘述

2 CBF文件结构

CBF文件由一系列信息块列表组成,而每个信息块的结构如下:

key
  1. 第一行信息块是关键字,显示所提供信息的类型;
  2. 第二行(某些关键字包含)是HEADER,说明变量的大小等信息,从而使求解器能够预先分配数据结构;
  3. 其余的行是保存实际信息的BODY。

CBF文件主要由以下四类信息块构成:

在这里插入图片描述

3 实例

考虑如下凸优化问题,有半定锥约束,二次锥在无序子指标上的约束以及两个等式约束。

在这里插入图片描述

注意到,这并不是直接CBLIB中的标准参考形式,所以需要进行相应的转换。

  1. 等式约束可以很容易重写为锥形式的约束 ( g 0 , g 1 ) ∈ { 0 } 2 (g_0, g_1) \in \{0\}^2 g0g1{0}2,方法是移动常量使右侧变为零;

  2. x 1 ≥ x 0 2 + x 2 2 x_1 \geq \sqrt{x_0^2+x_2^2} x1x02+x22 表示 ( x 1 , x 0 , x 2 ) ∈ Q 3 (x_1,x_0,x_2)\in Q^3 (x1,x0,x2)Q3,而我们需要的是 ( x 0 , x 1 , x 2 ) (x_0,x_1,x_2) (x0,x1,x2),所以这里用一个仿射变换:
    [ 0 1 0 1 0 0 0 0 1 ] ( x 0 x 1 x 2 ) = ( x 1 x 0 x 2 ) ∈ Q 3 \left[\begin{array}{lll} 0 & 1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{array}\right] \left(\begin{array}{l} x_0 \\ x_1 \\ x_2 \end{array}\right)= \left(\begin{array}{l} x_1 \\ x_0 \\ x_2 \end{array}\right) \in \mathbf{Q^3} 010100001 x0x1x2 = x1x0x2 Q3

通过上面两步转换,可以得到该问题CBF格式:

在这里插入图片描述

Note:

  1. PSDVAR
    • 1 : 表示只有一个属于半定锥的变量
    • 3:表示这个矩阵变量是 3 × 3 3\times 3 3×3
  2. VAR
    • 3 1: 表示有三个标量变量属于1个锥
    • F 3:表示这个锥是F ,3维的
截屏2023-08-11 11.38.58

Note:

  1. CON

    • 5 2:表示有5个约束,这五个约束属于两个锥
    • L= 2:表示第一个锥是L=,2维的(前两个约束属于这个锥)
    • Q 3:表示第二锥是Q,3维的(后三个约束属于Q这个锥)
  2. OBJFCOORD

    • 5 :表示目标函数中F矩阵有五个非零元

    • 第1列说明是第几个F矩阵

    • 第2、3列说明非零元在矩阵中的位置

    • 第4列给出非零元的值

  3. OBJACOORD

    1. 1:表示标量变量只有一个系数非零
    2. 1 1.0:表示 x 1 x_1 x1的系数为1.0
  4. CBF格式只存非零系数

  5. 对称矩阵只存下三角

截屏2023-08-11 11.39.05

Note:

  1. FCOORD
    • 9:表示约束中的F矩阵有9个非零元
    • 第1、2列:说明F矩阵的位置
    • 第3、4列:每个F矩阵中非零元的位置
    • 第5列:非零元的值
  2. ACOORD
    • 6:表示约束中A矩阵有6个非零元
    • 其余表示类似
  3. BCOORD
    1. 2:表示约束中b向量有2个非零元

截屏2023-08-11 14.21.50

参考文献

  1. CBF-format4
  2. CBLIB:2014

这篇关于锥规划问题存储格式——CBF格式简介的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/377509

相关文章

使用python生成固定格式序号的方法详解

《使用python生成固定格式序号的方法详解》这篇文章主要为大家详细介绍了如何使用python生成固定格式序号,文中的示例代码讲解详细,具有一定的借鉴价值,有需要的小伙伴可以参考一下... 目录生成结果验证完整生成代码扩展说明1. 保存到文本文件2. 转换为jsON格式3. 处理特殊序号格式(如带圈数字)4

IDEA和GIT关于文件中LF和CRLF问题及解决

《IDEA和GIT关于文件中LF和CRLF问题及解决》文章总结:因IDEA默认使用CRLF换行符导致Shell脚本在Linux运行报错,需在编辑器和Git中统一为LF,通过调整Git的core.aut... 目录问题描述问题思考解决过程总结问题描述项目软件安装shell脚本上git仓库管理,但拉取后,上l

idea npm install很慢问题及解决(nodejs)

《ideanpminstall很慢问题及解决(nodejs)》npm安装速度慢可通过配置国内镜像源(如淘宝)、清理缓存及切换工具解决,建议设置全局镜像(npmconfigsetregistryht... 目录idea npm install很慢(nodejs)配置国内镜像源清理缓存总结idea npm in

pycharm跑python项目易出错的问题总结

《pycharm跑python项目易出错的问题总结》:本文主要介绍pycharm跑python项目易出错问题的相关资料,当你在PyCharm中运行Python程序时遇到报错,可以按照以下步骤进行排... 1. 一定不要在pycharm终端里面创建环境安装别人的项目子模块等,有可能出现的问题就是你不报错都安装

idea突然报错Malformed \uxxxx encoding问题及解决

《idea突然报错Malformeduxxxxencoding问题及解决》Maven项目在切换Git分支时报错,提示project元素为描述符根元素,解决方法:删除Maven仓库中的resolv... 目www.chinasem.cn录问题解决方式总结问题idea 上的 maven China编程项目突然报错,是

Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题

《Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题》在爬虫工程里,“HTTPS”是绕不开的话题,HTTPS为传输加密提供保护,同时也给爬虫带来证书校验、... 目录一、核心问题与优先级检查(先问三件事)二、基础示例:requests 与证书处理三、高并发选型:

前端导出Excel文件出现乱码或文件损坏问题的解决办法

《前端导出Excel文件出现乱码或文件损坏问题的解决办法》在现代网页应用程序中,前端有时需要与后端进行数据交互,包括下载文件,:本文主要介绍前端导出Excel文件出现乱码或文件损坏问题的解决办法,... 目录1. 检查后端返回的数据格式2. 前端正确处理二进制数据方案 1:直接下载(推荐)方案 2:手动构造

Python绘制TSP、VRP问题求解结果图全过程

《Python绘制TSP、VRP问题求解结果图全过程》本文介绍用Python绘制TSP和VRP问题的静态与动态结果图,静态图展示路径,动态图通过matplotlib.animation模块实现动画效果... 目录一、静态图二、动态图总结【代码】python绘制TSP、VRP问题求解结果图(包含静态图与动态图

MyBatis/MyBatis-Plus同事务循环调用存储过程获取主键重复问题分析及解决

《MyBatis/MyBatis-Plus同事务循环调用存储过程获取主键重复问题分析及解决》MyBatis默认开启一级缓存,同一事务中循环调用查询方法时会重复使用缓存数据,导致获取的序列主键值均为1,... 目录问题原因解决办法如果是存储过程总结问题myBATis有如下代码获取序列作为主键IdMappe

Java Docx4j类库简介及使用示例详解

《JavaDocx4j类库简介及使用示例详解》Docx4j是一个强大而灵活的Java库,非常适合需要自动化生成、处理、转换MicrosoftOffice文档的服务器端或后端应用,本文给大家介绍Jav... 目录1.简介2.安装与依赖3.基础用法示例3.1 创建一个新 DOCX 并添加内容3.2 读取一个已存