- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
( z4 ~5 ]5 q& N0 |8 g( X9 L# V4 q; w
; ^- W A: B7 w' W Y' r- v( }3 S- M5 m' Q
( ?& x& r3 u \& l0 j$ q| 再求数学高人 | 7 [ _) ^# T7 c) f" {! R3 S8 w
) G3 O. C; v/ O8 x# H |
' H) g( m' _0 X" o9 R+ `, k- h# I8 h
5 a3 B7 Z; E! d( p; _/ y* q8 t) U
3 [7 z. _* D7 e/ L3 a用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
$ t6 Q. p0 P5 O; jMin x12-4x1-2x2
- Y3 G# E8 T" q6 R) Z% l: ps.t.不等式组条件如下
) |0 w, x3 z/ Z0 |& d sX1+X2<=4
4 Y2 \' I8 x$ ~- g w: M8 K2X1+X2<=5 1 P4 d$ B! c% T) W- L
-X1+4X2>=2
1 ^2 N7 W" [6 x& P' _3 Y( kX1>=0 X2>=0
X% z7 s! A4 q y/ z7 f2 t/ w【方法一】:
& r( {/ _2 c6 E/ ?因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
- l1 v% l$ H% i- t$ |" V& I首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 4 s4 x5 W; [; \, w5 N* y
点O x1=0 y= x12=0 6 Y1 ~* C$ v# p$ M# s: S; l1 `
点A x1=1 y= x12=1
9 U7 Z3 L; [, F. `; y点B x1=2 y= x12=4
; {' E+ R0 B1 _点C x1=2.5 y= x12=6.25 & V4 k( v: w3 {% U4 X" d8 G0 e. M
如图1: ) ~9 @3 I6 l8 P+ l- ?$ F
$ ]' s5 A) L1 A: A- W3 d Y
R! ]7 X$ Y4 N$ G0 f; b' `: u
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
% j: u9 N7 _$ k1 r5 RMin Y-4X1-2X2
$ l% g% }9 [- d9 |+ Ns.t.
. I& ~% @3 R [2 N& O0 ^X1+X2<=4 1 y7 _! `( L, T1 E) E
2X1+X2<=5 5 l* F$ l8 W- X6 \, h9 a. I9 } S" P
-X1+4X2>=2 * G# y( H' K- Q5 e p- g/ U
X1-Y<=0 ( y# k: U3 T. I( G
3X1-Y<=2
9 l( l* r8 ^: V: W- \# _3 u4.5X1-Y<=5
$ Q& B- \4 g2 [2 z% mX1>=0 X2>=0 $ I; J. S4 p" R( N+ A7 V
至此没有问题, 但解下来该怎样解就不知道了. ; h! P: j$ X Z: d& W& W, u
: H7 `' Y! n. Y
【方法二】:
) L& v3 K( ]5 [, s, p! P* I8 c+ T取近似值的方法不一样, ! N$ k; _' K7 d, t
X1=0p1+1P2+2P3+2.5p4 7 j* d4 p0 c8 r z0 n! ]( R
Y =0p1+1*p2+4*p3+6.25*P4 4 @2 w2 U6 [* o) {4 `# o9 g
1 = p1+p2+p3+p4 0 i H' w4 z+ G: h0 W8 j, d8 m7 e6 v
* n% ~: _0 G0 S
原规划可表示如下:
# V' T; O( j$ j, w4 \" M# c% u2 w/ QMin Y-4X1-2X2 1 Y6 m! y6 c! q5 `; \ f, s2 c
s.t.
4 |2 s- K2 [( T4 [ N( RX1+X2<=4
) h# B$ u/ c1 x5 H, M! j, s7 M2X1+X2<=5 : F2 D2 ?8 G- p! V$ _2 V9 F' M
-X1+4X2>=2
' `$ A( S$ ^, D* k# ^8 g-X1+P2+2P3+2.5P4=0
1 b( j" D a- J$ @$ r: G3 {-Y+P2+4P3+6.25P4=0 ' N( e/ Q6 A7 L% K: Z, z. w
P1+P2+P3+P4=1
6 _% d; Y+ Q! n. d" q* hY,X1,X2,P1,P2,P3,P4>=0 ) v. w2 C( F) P. I. D
同样的问题,到这里就不知道下面怎么解了. % o( t8 f1 i6 H7 A2 H: o
图2
& ~) _' s1 Y T( i 0 I) B$ P" L1 {5 D6 R
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! ) g! z- Y* Y/ I) Z" c: }6 u
+ u4 C7 ?9 b% F9 _- y, v4 k5 @/ C: s5 J% e0 ], x+ U
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 & W' m( p" n! i" U9 a
|    |
zan
|