- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
. l) G3 z* g" G) A % k% a( g8 K* N8 z% Y2 A
* l9 m& Y& M# g& {: V, o- E; e( B! X8 D
再求数学高人 | 0 o6 R$ H+ u; z; P/ ~5 _
; R6 x3 r9 M4 `$ L" C# h | ' ~( g1 z# d9 m- q m
7 f, e+ l. `3 ? O1 k# K) b# R+ {8 @7 O1 `5 t7 |: a: L" W
w6 O- p1 N" G/ c, g+ C用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), & h7 H! N% m; i4 f
Min x12-4x1-2x2 ) S0 \5 o. x9 z% U, D7 E p% {* T
s.t.不等式组条件如下
" ^; i. b# \5 qX1+X2<=4 7 X4 }- l$ _; L$ B+ v
2X1+X2<=5 0 K: H4 @1 D+ V8 |
-X1+4X2>=2 ' Z4 m( |2 e" h7 e: H
X1>=0 X2>=0
9 C( O7 |# N# H【方法一】: : J5 \4 y: s: F4 C3 m3 N$ D2 P! M
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 - j' w O) w1 P* `
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: ! Y9 T$ K6 z) b" G, [, d
点O x1=0 y= x12=0
! d/ f8 R" V! I1 g9 T2 T& X点A x1=1 y= x12=1 / Z3 N0 F, x! I" c7 f0 `# k
点B x1=2 y= x12=4 / S) y+ V0 T2 I/ N! n. y0 J
点C x1=2.5 y= x12=6.25
: p2 F0 g7 [* c! X+ S7 _7 X7 _- W如图1:
* w8 Y0 [$ Z% |5 a \ b$ ~* C1 ? 4 ~# |, l5 d5 B0 f" m
* I/ d- l8 X$ R( [' N1 C用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
3 u& `' c( d8 `% a0 s7 FMin Y-4X1-2X2 . o" T! Z7 ]( s- d5 f4 C
s.t.
& g. [# H6 S$ x6 j6 DX1+X2<=4 5 e* k7 ]9 Q6 e& M ~0 X; s
2X1+X2<=5 x: }) i5 J: _9 ^
-X1+4X2>=2 % S1 q7 H/ ~( t/ u/ |
X1-Y<=0 ! ?6 `% z& c9 _, i t2 D
3X1-Y<=2
$ V8 b8 E$ S9 Z: e4.5X1-Y<=5 9 z# c* x- A5 V- {2 u
X1>=0 X2>=0
1 L: M: k& }: L" W至此没有问题, 但解下来该怎样解就不知道了. ' f3 q- E# p& H# t8 j1 ~% Z
$ s: t/ O6 W, u$ m; q# J7 x
【方法二】:
# y# D1 f: ~2 I2 H R$ S取近似值的方法不一样,
# H3 s2 z) @3 n) }, H7 \X1=0p1+1P2+2P3+2.5p4
. l, T4 _, ~4 IY =0p1+1*p2+4*p3+6.25*P4 # K# [3 g6 H/ V6 N$ g/ u; g/ {/ Q$ j
1 = p1+p2+p3+p4
3 R0 b8 _, n; Z/ e' C! [) X. v, H: r" H2 C0 W# y* x
原规划可表示如下: - r) Q# T7 J; x( }
Min Y-4X1-2X2
' H8 S' O- j* h/ r7 ]s.t.
( Y$ U4 D* q# v) L! [& M' Z8 JX1+X2<=4
4 f) m5 W$ w' A% y$ M ~, `2X1+X2<=5
) c X F- z4 g1 j. l-X1+4X2>=2 ; q, Z7 ~4 z6 W
-X1+P2+2P3+2.5P4=0
O- q! c5 l- Z-Y+P2+4P3+6.25P4=0 ' z* A0 e. A D% l4 l; l
P1+P2+P3+P4=1
- F/ x2 t: Y9 u. ?8 aY,X1,X2,P1,P2,P3,P4>=0 3 L0 N; L% H# m; t% }4 @) x
同样的问题,到这里就不知道下面怎么解了.
4 }! J5 ]; M7 ]! h% ~ g8 a2 S# q图2 / J- v* b5 ?5 O; R$ F
1 G; @9 D: y: g: _( {1 n
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! n+ g3 J8 K4 {5 R: s
4 N8 l# ?' ]7 H$ @) Y2 `! V4 t, f( }; C# P
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
5 }+ C6 k5 s. M$ a ~) i |    |
zan
|