- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
- V, e( K7 a, Y- u, k3 d 7 `* N* B, y. e" ]; L6 `
/ j* ^8 J0 p/ h$ i G+ n9 o" I! Q5 z% E# Y* e) x7 p- \% _$ i
| 再求数学高人 |
3 K1 c7 x. }2 D. C
2 [+ {: u; c# ] k | + E$ U( G3 h7 o8 ]* t5 s
! u8 N8 N" o# `, } _% g/ ^% q
3 F, q3 P8 U* g+ @4 |2 B) C$ d. ~3 g) _/ P; M( B
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
: |' p# \ w+ b# B( W1 bMin x12-4x1-2x2 3 D$ k: {! K! ]- i" Z( i
s.t.不等式组条件如下 / w- q1 y- {. z4 D5 p
X1+X2<=4 / L" D* K( P- [# J
2X1+X2<=5
+ c* l' L* J3 H# J$ X& F/ A5 U4 `-X1+4X2>=2
) T$ W3 |$ x( zX1>=0 X2>=0 + U' \& ^$ S4 Y6 i( V8 K: g
【方法一】: 6 x/ `0 ~" q N! x9 J% ]3 N
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 ( B5 i7 p7 B3 K* j' d$ J& W
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: + z: g) j/ T' i! ]+ N6 e
点O x1=0 y= x12=0
5 w. k1 i9 ?; f- y4 z点A x1=1 y= x12=1
1 I+ Q4 e. \" ?5 ]! W2 z j点B x1=2 y= x12=4
) B2 i+ z! L4 [2 k点C x1=2.5 y= x12=6.25 4 L8 ]6 w1 m/ x- k( p0 m
如图1: * `" s# p. ?' F% b6 Y% P9 l4 W
( C3 ^! J8 i; V; s3 u
) l# ^. N$ f) P" }用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: 2 L5 F% `; w. I
Min Y-4X1-2X2
! O& h- G. d/ Y( L) u3 \s.t.
. P( Q$ ^( O& W" I4 e, H2 c3 kX1+X2<=4 1 x0 v$ E: l$ }3 K
2X1+X2<=5 - v7 F+ I& k$ M: h |% {8 S8 L
-X1+4X2>=2
3 V# E6 R" Y( CX1-Y<=0 % T c. L2 ?3 E+ C+ B r2 n
3X1-Y<=2
4 n1 o) ?; L: x7 v8 t2 c& r4.5X1-Y<=5 . M( r1 w/ m& B- t
X1>=0 X2>=0
: {6 D- m# ]6 L( t至此没有问题, 但解下来该怎样解就不知道了. , T" T" w2 @- y, W4 v' t3 X
1 _* p. ^6 S9 Y: }! X9 x8 p' K4 m【方法二】: 6 f/ j; T: C) W/ }
取近似值的方法不一样,
+ U* t& _; M2 k: I; e) [8 @& }4 ~X1=0p1+1P2+2P3+2.5p4 $ P9 { u" B5 L" M
Y =0p1+1*p2+4*p3+6.25*P4 . ?6 b/ Z- c( {' L. X+ ^
1 = p1+p2+p3+p4 9 {. N, a* z0 ]$ Z, B! e7 M
* w$ R- x. n# L) w: T' t原规划可表示如下:
7 R- y5 a5 T; g9 ^! P" nMin Y-4X1-2X2
/ e- E8 r4 z9 G8 j3 Os.t.
3 I0 j; e4 Q0 H$ R9 t) VX1+X2<=4
) n- r Y$ u8 r2 m) J: _6 s A2X1+X2<=5
Z# j, r: A7 h-X1+4X2>=2 6 o: v1 X; W$ U. r5 H7 b
-X1+P2+2P3+2.5P4=0
# W! b* p$ Q/ [; P2 j-Y+P2+4P3+6.25P4=0 % J5 f t5 o4 ^/ Y! G( R; b- u
P1+P2+P3+P4=1
8 f+ a! r; P; G, O# q- ]* U& lY,X1,X2,P1,P2,P3,P4>=0
" d5 N4 t. x: S* c1 Q0 n同样的问题,到这里就不知道下面怎么解了. & k5 I+ [( t5 |. U& ^, ^, k( u
图2
+ X, g5 p% [1 x- ]) ^0 `4 R9 p( b + Y6 s+ ^: x- Z. J6 J- y" _
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
* R) |" K$ k5 {+ R) H8 p3 S: |2 ~( o5 R
: |) Q2 X$ K- f$ G" d9 e2 ]我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
& S: K9 [7 i( ]) ^% V A$ ` V |    |
zan
|