在线时间 514 小时 最后登录 2023-12-1 注册时间 2018-7-17 听众数 15 收听数 0 能力 0 分 体力 40255 点 威望 0 点 阅读权限 255 积分 12788 相册 0 日志 0 记录 0 帖子 1419 主题 1178 精华 0 分享 0 好友 15
TA的每日心情 开心 2023-7-31 10:17
签到天数: 198 天
[LV.7]常住居民III
自我介绍 数学中国浅夏
/ I, U) _" @, \) h5 z& E 数学建模算法之优化模型【线性规划问题、非线性规划问题、整数规划问题、二次规划问题】 1. 线性规划问题(LP) - O0 J+ ^+ `) i" b
线性规划问题是要最小化或最大化一个受限于一组有限的线性约束的线性函数。
7 z+ Q7 Q. p T% X3 U( k7 V 9 y# d! p( l; }. }
* s% h- `- O, N" p& Y Matlab 中规定线性规划的标准形式为 & T3 l# q* E+ E
' U& a, Q! e9 Q- x+ _+ x
5 T' e/ r y( X) a2 _* {
第一个式子为目标函数,s.t. 式是约束条件。其中 c 和 x 为 n 维列向量,A、Aeq 为适当维数矩阵,b、beq 为适当维数列向量 5 R/ `+ M. h4 X% f r. w( h0 W g1 \: `
0 l" Y. {; g- `: R+ f- ~
$ t# z# }7 f1 f y" T: _
在 matlab 中,线性规划的函数为 linprog() ,有两种常用形式: * V" R+ U1 A, b# L
X = linprog(f,A,b,Aeq,beq,LB,UB,X0)
0 L( r3 o% m7 V5 n6 c [X,FVAL]=linprog(f,A,b,Aeq,beq,LB,UB,X0) 1 F6 _5 w2 z$ J- C
返回的值 X 是向量 x 的值,FVAL 是目标函数的值,LB 和 UB 分别是变量 x 的下界和上界, 是 x 的初始值。
; t- k( \3 V" h; v, j
p5 f* G1 K4 W# p# L% x: ~
) J4 v* q2 `6 p1 g. } 1.2 应用例子
. H0 S4 c' d7 Q5 ^ 求下列线性规划问题: ! T+ r/ }2 o/ [; b2 ], b" E: U
" h+ p2 ^: i; f1 ?$ m I. `
' Z: [5 J, K4 r1 K5 z4 U, D
依据 Matlab 的标准,默认求解是求最小值,而本例是求的最大值,把 z 的系数变为相反数,即 -1 就好了,同理下面的大于等于号也做同样处理,然后没有上界 UB,下界 LB 为三个变量都为 0,也就是一个全零的矩阵 zeros(3, 1)
' Y" {2 r5 b: d, m& f8 N 2 ~; ^3 d5 {' h/ l' O
8 [' B' z* j* j3 L8 p( e) W 编写一个 .m 文件:
# E; E8 `2 @8 @! | % e2 Y' K, ?9 f
3 Y9 H5 j, W, F; s5 Y* F
1.3 相关问题
- b$ V9 T" K4 S- M 运输问题(产销平衡)
S( L+ X( z+ j- y 指派问题(匈牙利算法) & _- ^* M+ Z- T1 [8 f: \
对偶理论与灵敏度分析
+ a1 f8 I+ U) I# q( S- t* D 投资的收益和风险 3 N2 S& s& i/ P% }' D u$ v
2. 非线性规划问题(NLP) : o& J/ w$ P: s1 V& `
如果目标函数或者约束条件中至少有一个是非线性函数时的最优化问题就叫做非线性规划问题。 1 J2 X: V$ s* n5 C7 F7 t- a
& @: U) b R; o
* z8 k+ @* H8 y3 R# G
# P. @& l! ~9 L- P5 J' e 2.2 非线性规划的基本解法 . S3 s0 h+ j0 S: U* i; ?1 v B
罚函数法 ; T4 g5 L& K( |# [+ M9 n7 r
近似规划法
7 ?9 q4 r" F6 l" C 近似规划法的基本思想;是将问题(3)中的目标函数f(x)和约束条件g(x);h(x)近似为线性函数,并对变量的取值范围加以限制。从而得到一个近似线性规划问题。再用单纯形法求解。把其符合原始条件的最优解作为(3)解的近似。 / y1 ]+ @9 k4 U! X$ H9 M
每得到一个近似解后,都从 [color=rgba(0, 0, 0, 0.75)]都从这点出发,重复以上操作。 2 b0 D' [& \( U6 m! p/ B! F
r" E0 C% E9 M9 b+ V9 ^
. q6 p2 k. t: T4 x$ @/ {
$ g& T8 \& _: [5 G: |) y3 H
. j; ?, F D) J& S+ ` 2.3 相应问题
5 N% A! F+ T( F+ b 无约束问题(一维搜索方法、二次插值法、无约束极值问题的解法) % }) m: \; Y. p
约束极值问题(二次规划、罚函数法)
+ V. y, x' v" e p0 K/ |# C) u. ` 飞行管理问题 ~' L. c. S' Z' m `+ `
3. 整数规划问题(IP) # M$ a" ~ x; Q& I
数学规划中的变量(全部或部分)限制为整数时,称为整数规划,例如,所求的解是机器的台数,人数,车辆船只数等等。 0 x: q" W4 W7 G s
$ M% l: r5 w c4 H8 p3 a+ n
; e; E: \% r0 S8 G/ h 整数规划的分类:
1 D: Y% y: E( } 纯整数规划:全部决策变量只能取整数的线性规划 5 w" f$ Z. R' `; {/ a$ n' s) U
混合整数规划:决策变量中有一部分必须取整数,另一部分可以不取整数的线性规划
. x) i: z( W5 A& M& i 0-1整数规划:决策变量只能取0,1的线性变化
* v0 c/ z8 K5 n4 o 3.1 混合整数规划问题(MIP) * c: {" U, z: u' ~* r' @
混合整数线性规划是整数线性规划模型的一种。 / ^+ e+ h& Z# e, o) {
, B* L8 c7 E# k6 {8 k# T3 c1 F( U 6 H& k# T7 X% \# \' W$ H+ c
整数线性规划模型分类: 6 A. K+ ]# |- U) Y/ F
% G. T) [& u2 I5 }1 S- E% ]: H9 G7 v/ \ # ~; K9 W1 K, Z
若I={0,1},J={1,…,n},即全部的决策变量仅取0或1,称之为0-1规划;
: ?4 O% a q0 \$ ?) \" q! ~5 N) D6 a 若J是{1,2…n}的非空真子集,即仅有部分决策变量要求取整数,称为混合整数线性规划;
! z0 T w; T" h. k3 d 若J={1,2,…n},即全部的决策变量都取整数,称为纯整数线性规划; * S, P% Q5 Z" e) `, D4 Q5 _
常用的整数规划问题解法有:
9 W, l2 J4 C& j$ K 8 w' l) ^' L1 O" ]: m! P. L: b! e
+ ~8 [6 B* y/ B9 V& E% l$ q 分枝定界法:可求纯或混合整数线性规划
! v2 g$ p& O$ x6 }/ c+ A" `; M, b 割平面法:可求纯或混合整数线性规划
7 U1 t" v( i4 o 隐枚举法:用于求解0-1整数规划,有过滤法和分枝法。 9 L$ j9 X5 l" l1 I+ h9 f9 x
匈牙利法:解决指派问题(0-1规划特殊情形) 7 D8 B$ u' a) @) P5 y9 r
蒙特卡罗法:求解各种类型规划 9 M9 z# q! q8 @& k1 F
3.2 常用方法讲解
4 C/ w! p7 k9 e- f% x
+ p! I0 d. g/ I+ z" D" P: M8 b
) L+ p- Y5 J/ k5 q$ t+ y
# Q. _/ t& b" ^4 q5 H: \ % l- V( I# }2 i+ B6 w* W" ?- Y2 y6 j
0 S8 B. X" W8 G+ g5 e- G7 V' P
4. 二次规划问题(Quadratic Programming)
/ V* S( [( b# h4 m
' I) ^' u8 ?4 `4 i! _ 3 K6 C6 x% t9 J& k( X: p% ^
6 s! C6 Q9 W) z* y( t
5. 混合整数二次规划问题(MIQP) 5 j" L$ N5 d) j& u! k
通过上面的定义我们不难看出混合整数二次规划问题本质就是一种混合整数规划问题的一种特例,其中他的目标函数为二次型,其约束条件满足混合整数规划问题。
/ B- O# {6 `2 A$ l, Y
5 q$ D% P T* `
, a4 I) E5 d' T ( O" e( N1 @/ x/ N
zan