MATLAB的线性规划工具

Algorithm 同时被 2 个专栏收录
3 篇文章 0 订阅
11 篇文章 0 订阅

MATLAB的线性规划工具

问题描述

某次考试,总分100分.分三种题型:

单选题13题,每题3分,共39分.
多选题10题,每题3.1分,共31分.
判断题12题,每题2.5分,共30分.

交卷之后,只知道分数为79.3,问:单选题、多选题、判断题分别做错了多少?

MATLAB meshz辅助求解

假设单选题、多选题、判断题做错数量分别为x,y,z.则可列方程:
3 x + 3.1 y + 2.5 z = 100 − 79.3 3x + 3.1y + 2.5z = 100-79.3 3x+3.1y+2.5z=10079.3
3 x + 3.1 y + 2.5 z = 20.7 3x + 3.1y + 2.5z = 20.7 3x+3.1y+2.5z=20.7

上述函数是一个空间平面,如果不理解为什么是空间平面可以参考直线与平面方程的几何表达.

最直接的办法就是把平面绘制出来,找到上面的正整数解,脚本如下:

s = linspace(0, 10, 11);
t = linspace(0, 10, 11);
[ss, tt] = meshgrid (s, t);
x = tt;
y = ss;
z = -(3.*x + 3.1.*y  - 20.70)/2.5;
pbaspect ([1 1 1]);

meshz(x, y, z);

patch([0 10 10 0], [0 0 10 10], [0 0 0 0]) 

camproj('perspective')
xlabel('x');
ylabel('y');
set( gca, 'ydir', 'reverse' )

对应图像:

在这里插入图片描述

绿色平面之上的部分就是要搜索的区域.
脚本中x,y变量使用的都是正整数,所以只要确保z也是正整数就可以了,可以查看z的数值:

在这里插入图片描述

显然z只能取1,对应的x,y分别为4与2,所以这个问题的解是: x = 4 , y = 2 , z = 1 x = 4, y = 2, z = 1 x=4,y=2,z=1

MATLAB intlinprog函数求解

显然,上述的方法很难让人满意,因为这个解是靠肉眼观察出来的.
实事上MATLAB有专门函数求解,就是intlinprog().

上述问题可以写成下列条件:

x 1 , x 2 , x 3 x_1, x_2, x_3 x1,x2,x3为整数

x 1 ≥ 0 x_1 \geq 0 x10

x 2 ≥ 0 x_2 \geq 0 x20

x 3 ≥ 0 x_3 \geq 0 x30

3 x 1 + 3.1 x 2 + 2.5 x 3 = 20.7 3x_1 + 3.1x_2 + 2.5x_3 = 20.7 3x1+3.1x2+2.5x3=20.7

3 x 1 + 3.1 x 2 + 2.5 x 3 3x_1 + 3.1x_2 + 2.5x_3 3x1+3.1x2+2.5x3的最小值.

编制脚本:

Aeq = [3,3.1,2.5];
beq = 20.7;
intcon = 1:3;
f = [3;3.1;2.5];

A = [
    -1,0,0;
    0,-1,0;
    0,0,-1;
    ];

b = [0;0;0];

x = intlinprog(f,intcon,A,b,Aeq,beq)

在这里插入图片描述

常见线性规划问题

甲乙两种溶液,要兑成另一个产品.
但有成本限制,不成超过15,甲的单价是3,乙的是5.
甲乙都含有害成分,对它浓度的要求不能超过2%,对于此成分,甲的浓度为5%,乙的浓度为1%.
现在想把这个产品尽可能地做重,问应该分别使用多少的甲与乙?

直观的观察是,要花更少的钱,做去更多的产品,只能尽量使用便宜的甲,但甲的有害成分又比较高,一味的选用甲也不行,这是问题所在.

这是一个典型的线性规划问题.假设甲乙两种溶液的质量分别是x,y,可以列出下面的约束不等式:

3 x + 5 y ≤ 15      ( 1 ) 3x + 5y \leq 15 \ \ \ \ (1) 3x+5y15    (1)

( 0.05 x + 0.01 y ) / ( x + y ) < 0.02 (0.05x + 0.01y)/(x+y) < 0.02 (0.05x+0.01y)/(x+y)<0.02

化简后:

3 x − y ≤ 0      ( 2 ) 3x-y \leq 0 \ \ \ \ (2) 3xy0    (2)

所求为令目标函数x+y最大.

两个不等式边界直线为AB与CG,目标函数直线与EF平行:

在这里插入图片描述

显然G是所求的点,
求解线性方程得AB与CG交点为, ( 5 6 , 5 2 ) (\frac{5}{6}, \frac{5}{2}) (65,25)
最大值为 10 3 \frac{10}{3} 310

linprog求解

编制脚本:

A = [
    3,5;
    3,-1;
    ];
b = [15;0];
f = [-1;-1];
x = linprog(f,A,b)

在这里插入图片描述

结果与手工计算相同.

参考

  • 2
    点赞
  • 5
    评论
  • 17
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

打赏
文章很值,打赏犒劳作者一下
相关推荐
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页

打赏

键盘上的疯兔

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值