QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3299|回复: 0
打印 上一主题 下一主题

[其他资源] 数学建模:matlab解决线性规划问题

[复制链接]
字体大小: 正常 放大

326

主题

32

听众

1万

积分

  • TA的每日心情
    慵懒
    2020-7-12 09:52
  • 签到天数: 116 天

    [LV.6]常住居民II

    管理员

    群组2018教师培训(呼和浩

    群组2017-05-04 量化投资实

    群组2017“草原杯”夏令营

    群组2018美赛冲刺培训

    群组2017 田老师国赛冲刺课

    跳转到指定楼层
    1#
    发表于 2020-5-28 15:56 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    在人们的生产实践中,经常会遇到如何利用现有资源来安排生产,以取得最大经济效# D* S4 `+ {' }. q# p! o8 S* s
    益的问题。此类问题构成了运筹学的一个重要分支一数学规划,而线性规划(Linear
    $ N+ _" O; X; C, H2 ZProgramming, LP)则是数学规划的一个重要分支。自从1947年G. B. Dantzig 提出求解线0 ?* u, X% F9 a- K1 Z( y: S
    性规划的单纯形方法以来,线性规划在理论上趋向成熟,在实用中日益广泛与深入。特别
    ) O. N& C6 r1 k是在计算机能处理成千上万个约束条件和决策变量的线性规划问题之后,线性规划的适
    9 S- q& V! S4 @$ W- b4 g. u用领域更为广泛了,已成为现代管理中经常采用的基本方法之一。
    ) t4 l, \' E  c1 w5 }( i8 y
    ; x* i" {2 G$ W例如,给定m个数据点(xi,yi),i=1,2,…,m,拟合一条直线y=ax+b(即确定参数a、b),使得所有数据点(xi,yi)和拟合直线上对应的点(xi,axi+b)之间距离的最大值rmax最小。也就是说,对整个这组数据点而言,最大绝对偏差r=max{ | yi-y(xi) | }最小。这种准则实际上定义了如下优化问题:: @* {# Y* N& c. k( k) t
    4.png # ?# n# W  B+ O- r4 b" Z1 N+ j
    这是一个线性规划问题,具有很广泛的应用。7 @7 X) T: G( p% _) ^. a
    7 C# D& |3 b( z2 h3 _
    木匠问题7 N4 Z* w  O( `6 Y
    木匠销售桌子和书架的单位净利润分别为25美元和30美元,他希望确定每周生产的桌子数量x和书架数量y。他每周最多有690 张木板和120小时的劳动时间可以利用,如果木板和工时不用于生产桌子和书架,他能够将它们有效地使用在其他方面。他估计,生产一张桌子需要20张木板和5小时劳动时间,生产一个书架需要30张木板和4小时劳动时间。模型为
    8 _1 R: m7 a8 a# L8 ?/ E 5.png ; g3 [1 f7 q$ h3 k' ^

    ' }6 E, @# K) A4 }, {木匠问题中的约束所代表的凸集在图中用多边形ABCD表示。请注意,约束所代表的直线有6个交点,但只有四个交点(即A—D)满足所有约束从而属于该凸集。点A—D是该多边形的极点。
    6 Y; |) S9 p" m6 g. V; Y# X如果一个线性规划存在最优解,它必然也会出现在约束所形成的凸集的某个极点上,极点上目标函数的值(木匠问题的利润)是
      Z& Q4 `$ b' ~! K$ X( o* d' ]! M. j2 z! e$ _+ C
    1.png
    - I. v8 L/ m+ n) Z
      V* y! O7 \" I7 |# F7 y( A8 C1 O- K* O9 r9 F. D' v4 B
    因此,木匠每周应该制作12张桌子和15个支架,每周最大利润为750美元。- x! h( t2 w7 D
      I2 F' P, ?) e( Q5 s3 g- y) A
    木匠问题(线性规划)的matlab解法
    7 Z% @( x( r; q/ v0 o9 T7 L+ m1 ]( L! c5 l
    线性规划的目标函数可以是求最大值,也可以是求最小值,约束条件的不等号可以是6 K8 N* a; _; \) {
    小于等号也可以是大于等号。为了避免这种形式多样性带来的不便, Matlab中规定线性
    . S' |2 Q: ~5 P+ n4 _( p规划的标准形式为& A, }% `5 a( w) b2 h( ~$ C
    2.png

    . V: e9 S0 S& J( S5 @* @1 L4 O* Q: C! W. |9 A
    式中:f,x ,b ,beq,lb, ub为列向量,其中f称为价值向量,b称为资源向量;A , Aeq为矩阵。
    ; D& C( ]* g) C* D8 K; A- XMatlab中求解线性规划的命令为: t6 l, }9 N  ]& H( X+ L2 {  z
    [x,fval] = linprog(f,A,b), C# O+ I4 M) ]; a; E$ T6 `0 B
    [x,fval] = linprog( f,A,b,Aeq, beq)
    * G! |  C4 l+ V4 z% \[x,fval] = linprog( f ,A, b,Aeq , beq,1b,ub)1 u$ t$ j3 v9 P& i
    式中:x返回决策向量的取值;fval返回目标函数的最优值;f为价值向量;A和b对应线性
    ( C3 t+ T) r# h. l5 h( S/ M0 V不等式约束;Aeq和beq对应线性等式约束;lb和ub分别对应决策向量的下界向量和上
    ! V8 u9 n0 D  F. \界向量。
    0 i- S4 D. g' ~. o- P例如,线性规划
    $ c  w' |* e! `  \ 6.png
    / h3 c, P7 f" ~5 g的matlab标准型为1 ?0 }$ f, U7 N* g
    7.png * o3 D4 z, w- F0 M) m: Z: o

    / P  e, X* i) K5 J3 b: L, B1 Z2 J7 \6 P& a6 z" L
    附上木桶问题的matlab代码5 F0 |6 R* u: m
    3.png 1 v1 E( x$ ~6 h0 W7 K/ p4 S

    & B+ N  K' J( A+ r+ x% Y! i# V# j————————————————
    5 v" ?4 ]2 K: d% I2 y5 ?* e版权声明:本文为CSDN博主「板砖大的砖头」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    + U+ t0 U: |) ?1 T. |8 ^+ L- Z8 \& x原文链接:https://blog.csdn.net/qq_44953660/article/details/104227831( t; V6 l8 z) P. _7 Q
    7 o( W7 s" D6 O! _5 B

    7 y! y3 g* O2 N$ s: s2 X9 t1 f5 ~
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-13 06:35 , Processed in 0.426097 second(s), 54 queries .

    回顶部