matlab解决有约束的线性规划问题

最新推荐文章于 2024-07-01 14:44:24 发布

此方家的空腹

最新推荐文章于 2024-07-01 14:44:24 发布

阅读量8.7k

点赞数 6

分类专栏：数学建模与matlab

本文链接：https://blog.csdn.net/qq_42138454/article/details/100136396

版权

数学建模与matlab 专栏收录该内容

12 篇文章 15 订阅

订阅专栏

在《matlab解决有约束的二次规划问题》中，提到了线性规划问题（linear programming），而线性规划问题实际与二次规划问题十分相似，现在在matlab中实现求解

这是用于求解问题的linprog（）函数文档

同样使用一个例子说明这个函数的调用方法

【例】求解如下线性规划问题

【分析】对比帮助文档，需要将目标函数和约束条件分别表示成这样的格式：

表示出这个式子中的各个矩阵即可

1.目标函数的表示

与线性代数课中所学的一样，用一个列矩阵f表示表示各个一次项的系数，即 $f=\begin{bmatrix} -5\\ -4\\ -6 \end{bmatrix}$

2.约束条件的表示

由于约束条件中不含等式，所以Aeq与beq全为空 Aeq=[] , beq=[]

A是前三个约束条件的左边部分，表示为 $A=\begin{bmatrix} 1& -1 & 1\\ 3& 2& 4\\ 3& 2& 0 \end{bmatrix}$

b是前三个约束条件的右边部分，表示为 $b=\begin{bmatrix} 20\\ 42\\ 30 \end{bmatrix}$

lb表示变量的下限 $lb=\begin{bmatrix} 0\\ 0\\ 0 \end{bmatrix}$ ,使用这个式子表示第四个约束条件

3.带入程序中进行计算

查看linprog的帮助文档，根据调用格式使用linprog函数

>> f=[-5;-4;-6];
>> A=[1 -1 1;3 2 4;3 2 0];
>> b=[20;42;30];
>> lb=[0;0;0];
>> [x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb)

最后一行不打分号，方便显示结果

产生结果如下：

Optimal solution found.


x =

         0
   15.0000
    3.0000


fval =

   -78


exitflag =

     1


output = 

  包含以下字段的 struct:

         iterations: 3
    constrviolation: 0
            message: 'Optimal solution found.'
          algorithm: 'dual-simplex'
      firstorderopt: 1.7764e-15


lambda = 

  包含以下字段的 struct:

      lower: [3×1 double]
      upper: [3×1 double]
      eqlin: []
    ineqlin: [3×1 double]