数学建模社区-数学中国
标题:
数学建模：matlab解决线性规划问题
[打印本页]
作者:
zhangtt123
时间:
2020-5-28 15:56
标题:
数学建模：matlab解决线性规划问题
在人们的生产实践中,经常会遇到如何利用现有资源来安排生产,以取得最大经济效
* q, J" t2 B7 @( s8 }
益的问题。此类问题构成了运筹学的一个重要分支一数学规划,而线性规划(Linear
6 Q4 b/ ?: I p; R- A# q
Programming, LP)则是数学规划的一个重要分支。自从1947年G. B. Dantzig 提出求解线
$ I: K, @8 Z: v0 G+ _
性规划的单纯形方法以来,线性规划在理论上趋向成熟,在实用中日益广泛与深入。特别
5 q% F, _) h% f" Y" b( B
是在计算机能处理成千上万个约束条件和决策变量的线性规划问题之后，线性规划的适
: g' Q2 q& o/ r* v& l- o4 P
用领域更为广泛了,已成为现代管理中经常采用的基本方法之一。
0 u2 f# h3 Y$ o& @/ A- I, i0 _; \5 {; c
% q4 q4 x! t- N0 E
例如，给定m个数据点(xi,yi),i=1,2,…,m，拟合一条直线y=ax+b(即确定参数a、b)，使得所有数据点(xi,yi)和拟合直线上对应的点(xi,axi+b)之间距离的最大值rmax最小。也就是说，对整个这组数据点而言，最大绝对偏差r=max{ | yi-y(xi) | }最小。这种准则实际上定义了如下优化问题：
: M0 h' D/ J* ^7 ?
2020-5-28 15:58 上传
下载附件
(6.42 KB)
& X3 T2 j0 S8 x- f' G2 ^" Y* ]
这是一个线性规划问题，具有很广泛的应用。
% [7 @" C( W8 Y/ B8 L
6 g7 @* M) l5 {. t
木匠问题
; h- k2 a0 g& @- K3 s5 u* o
木匠销售桌子和书架的单位净利润分别为25美元和30美元，他希望确定每周生产的桌子数量x和书架数量y。他每周最多有690 张木板和120小时的劳动时间可以利用，如果木板和工时不用于生产桌子和书架，他能够将它们有效地使用在其他方面。他估计，生产一张桌子需要20张木板和5小时劳动时间，生产一个书架需要30张木板和4小时劳动时间。模型为
/ ?5 s% W. r* `2 j# k5 g7 U
2020-5-28 15:58 上传
下载附件
(11.62 KB)
6 h: R& h! H0 B; A! ]0 [
- n$ F- ?9 h, F7 U
木匠问题中的约束所代表的凸集在图中用多边形ABCD表示。请注意，约束所代表的直线有6个交点，但只有四个交点（即A—D）满足所有约束从而属于该凸集。点A—D是该多边形的极点。
) I3 A7 W5 `/ L8 _" I5 i
如果一个线性规划存在最优解，它必然也会出现在约束所形成的凸集的某个极点上，极点上目标函数的值（木匠问题的利润）是
- B' M& S/ `- u/ ^' W4 m
% w+ M0 u$ o# ~! |7 k/ m+ h( g+ G" K
2020-5-28 15:53 上传
下载附件
(6.63 KB)
0 V% d2 |0 S$ \0 w z P
7 o- o1 K4 `! V8 h
- b) e: Y1 l! i5 t2 v
因此，木匠每周应该制作12张桌子和15个支架，每周最大利润为750美元。
7 H6 b- `8 g8 j( c5 G5 r
6 M6 c6 G; t% I4 ~& |
木匠问题(线性规划)的matlab解法
0 n) q+ p5 ?$ L6 f+ ]- T* S, Q
1 V/ C7 I4 H: @' x/ R
线性规划的目标函数可以是求最大值,也可以是求最小值,约束条件的不等号可以是
: b8 G2 k0 ]$ L2 I% W" Z
小于等号也可以是大于等号。为了避免这种形式多样性带来的不便, Matlab中规定线性
) w; f+ C0 @6 u- C, z+ ]1 H
规划的标准形式为
. B; Z/ Y3 \, u8 T! i
2020-5-28 15:55 上传
下载附件
(6.74 KB)
+ H$ V, F; f2 r) d1 ` o; q! g
7 e- M, h( v3 V$ x4 P
式中:f,x ,b ,beq,lb, ub为列向量,其中f称为价值向量,b称为资源向量;A , Aeq为矩阵。
! {; d3 z, @5 @8 N& p
Matlab中求解线性规划的命令为
' J+ ^; W/ e8 l- b% l$ Y$ V
[x,fval] = linprog(f,A,b)
' O+ F b1 j# h* x X
[x,fval] = linprog( f,A,b,Aeq, beq)
/ D) K0 {9 p9 d
[x,fval] = linprog( f ,A, b,Aeq , beq,1b,ub)
7 B9 g K3 R3 H, \
式中:x返回决策向量的取值;fval返回目标函数的最优值;f为价值向量;A和b对应线性
; e, y1 g- h: a4 U/ P8 P2 r6 W
不等式约束;Aeq和beq对应线性等式约束;lb和ub分别对应决策向量的下界向量和上
* _; U" R; }; ^8 F% ~
界向量。
7 Q: C- a- o* n3 C! x
例如,线性规划
+ P! G# d) ^& `8 z
2020-5-28 15:59 上传
下载附件
(2.61 KB)
4 I+ B5 k& N& A9 Y; f# U2 J7 }
的matlab标准型为
' }1 o, c' V. b0 @% z7 f
2020-5-28 15:59 上传
下载附件
(2.96 KB)
& g+ q+ c9 {& f
5 p& `& _( W- v
4 n* _; k3 j. V% i
附上木桶问题的matlab代码
3 ^) e0 e5 q2 ^
2020-5-28 15:56 上传
下载附件
(10.68 KB)
0 S4 c. I% t/ x' N
. w1 n' r( w8 V1 a
————————————————
1 M/ `! J2 D' @
版权声明：本文为CSDN博主「板砖大的砖头」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
+ C6 z7 J9 j4 q
原文链接：https://blog.csdn.net/qq_44953660/article/details/104227831
`4 h7 }- @" l, |' ?
4 J' J, d4 ?* ^4 [( m
) m# C1 J8 u: ~2 x/ H+ d: J \3 d
6 u" f0 C' U) V& A% g3 k: i9 M
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
7 o- o1 K4 `! V8 h