- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
% {$ e2 M8 X. P, C% Q
8 v2 ?6 E! B4 S+ O2 |" r2 A( z; g( _# r: |9 K
* D4 J, c/ v0 w( M| 再求数学高人 |
5 W8 q1 M3 O% ~& m7 u. J
" ~- @( l7 U$ L, I f. q4 @ | # @! y- b% u8 K) I" ~
/ W) P& c- K/ T2 h8 b6 ^& H& `
) k" |4 e( a1 w7 K1 S" I
- U" i6 q* P) \; j; i用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), 4 `; {$ ?$ s/ }
Min x12-4x1-2x2
6 \5 T; e# B5 g. c7 f# g3 o6 as.t.不等式组条件如下
/ `% Z% F" ^) z9 |9 k% M; t/ CX1+X2<=4 4 Y9 x1 g3 ^8 |' D$ p5 j- _/ @
2X1+X2<=5
/ V6 z, }5 j0 ] L; c-X1+4X2>=2
* F+ Z& U2 h4 ^4 s& w& OX1>=0 X2>=0
) V' i; G# o! L0 `5 e+ e【方法一】: 8 w6 h2 q) M6 c* s
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
0 I6 n7 @! B1 f首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下:
0 ~7 d6 d# [2 C7 \1 e点O x1=0 y= x12=0
5 G, s9 m6 J$ ^" v点A x1=1 y= x12=1
7 y* \$ x# l# m6 T点B x1=2 y= x12=4
5 F2 ?& e) t$ k+ p点C x1=2.5 y= x12=6.25
0 U! \5 v1 l. L9 n. ^如图1: . s* n" s' b: V" o
) k6 F f( ^! i! ?8 m8 h( O
( P- o" X5 m& [7 e, h用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: $ C8 J2 q/ b8 X
Min Y-4X1-2X2 Y4 M7 i, b6 ~/ l( ?: z b* L
s.t.
, p; S0 J( A: d. z. HX1+X2<=4
9 |4 a! ~+ X" z% U2 s2X1+X2<=5
" G) y5 f$ M) l6 U; q-X1+4X2>=2 & o% k4 e& L7 }% @
X1-Y<=0 $ `( o# D7 O0 N2 X
3X1-Y<=2
- ^' a/ S. `& @ N! Q7 ~/ W6 U4.5X1-Y<=5 7 h; a/ c# G/ P6 R9 j2 l* p. X
X1>=0 X2>=0 : C, R0 g4 Z( ~. F
至此没有问题, 但解下来该怎样解就不知道了.
% m9 O: k$ \! L; z9 [
& O" U' s8 g; ?. s+ j* F2 D【方法二】:
' S7 B$ `5 Z9 q i0 E取近似值的方法不一样, 5 i8 `0 e7 `1 l+ A. V
X1=0p1+1P2+2P3+2.5p4
( Y6 ^, h7 g JY =0p1+1*p2+4*p3+6.25*P4 ( M X# T) Z! D0 T8 |
1 = p1+p2+p3+p4
$ H8 N' ]" n& D$ P% M8 O$ x4 B6 V( L9 v# W; Z
原规划可表示如下: " u7 m+ s9 G+ |0 |) S
Min Y-4X1-2X2 , N4 r$ z0 x T2 ]
s.t. 3 }. O& g U& E$ n" z
X1+X2<=4 ) {$ z' O v$ V, ]2 f
2X1+X2<=5
+ W. `( K9 ^, M-X1+4X2>=2 % M- c5 U+ p3 \- X* A
-X1+P2+2P3+2.5P4=0
! Y' W' ]# ?: ^ g- s-Y+P2+4P3+6.25P4=0
" C+ _2 Q5 J% i4 Z& a' }P1+P2+P3+P4=1 # R9 N" {, M, u3 t
Y,X1,X2,P1,P2,P3,P4>=0
9 G5 m" Z f: }/ D9 X" ], t同样的问题,到这里就不知道下面怎么解了.
: g( i/ ]& ^( @图2 / [3 I& O) f5 J7 U% m' ?
. F' h/ G# H- m: m) O
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
& p* _6 C" {. h7 ]0 O% q; Q
) U( Y% B7 r& B6 D0 ]* i- X: h8 f4 i/ w1 a- w3 m+ C. M; h, |
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
( ^$ i; f4 G, _5 ]7 R1 g( c) h |    |
zan
|