motivation:
1.要实现“符号几何规划”的功能,使用MATLAB自带的fmincon函数来实现
2.工程的代码是通过Python来完成的,需要Python调用MATLAB实现
思路:
1.使用Python直接生成约束条件,目标函数,初始值等函数
2.使用Python来调取MATLAB接口,执行MATLAB的.m文件
3.MATLAB的.m文件中,将输出结果写入txt文档
4.Python读取txt文档结果
1.MATLAB部分
首先,简要介绍下MATLAB中自带的fmincon函数。fmincon函数是用来解决“非线性约束规划”,本文需要用到的是它解决“符合几何规划”的能力。【几何规划】的定义是,目标函数和约束条件均为系数为正的多项式,而【符号几何规划】在此基础上拓展,允许符号为负。
fmincon函数的约束形式如下:
其中,c(x)表示非线性多项式的不等式约束部分,注意,不等式需要严格按照“小于等于0”的格式(可以进行一些简单的不等式变换,写成右边为0的形式),ceq(x)表示非线性多项式的等式约束部分,形式同上。(b和beq是线性约束部分的参数,写成向量形式,A 和 Aeq 是线性约束部分的参数,写成矩阵形式,本文的问题场景不涉及“线性约束”条件)。lb和ub分别为x变量的上下限。在MATLAB中使用fmincon,传入参数如下:x=fmincon( objfunction, x0, A, b, Aeq, beq, lb, ub, st)
在本文的问题背景下,目标函数的.m文件代码如下