- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
; G- g" E4 J8 C/ i1 B' | ' g8 u$ w% c* I$ v+ ?
* P( ^: t% P/ G: h
( T/ h) g5 P& z0 b| 再求数学高人 | % i# t2 E4 X9 u! F- r/ e
6 Z; r' a" u: f6 d* v | " l8 v1 W1 I4 I" L
" b% I4 A8 i, `% c1 A5 |
9 _' K) d9 q& b2 W, T) @3 y# A4 G7 U
" O& g) w" X5 w2 m# B
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), \. ^( l$ Q' ^
Min x12-4x1-2x2 1 H9 p3 E" d) v7 a, U
s.t.不等式组条件如下
5 c9 U! Y5 k$ Q3 X# y) RX1+X2<=4 x0 Y5 o @7 Z' L/ t
2X1+X2<=5
$ T% o. {5 r' Z$ C: |-X1+4X2>=2
9 ^- k1 g3 O1 AX1>=0 X2>=0 6 @/ D# A H; B% c5 x j5 y
【方法一】:
5 }4 q5 k+ c" F8 t' _; o& K1 i2 ?因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
( U# [5 c2 z, l5 l# M& C. _" F8 _首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 9 M; Z. ? O# r; t% Q3 E
点O x1=0 y= x12=0
+ R7 x6 M! ]& k5 h点A x1=1 y= x12=1
2 l* e# o2 V6 D点B x1=2 y= x12=4 3 M0 Z- }7 R# x
点C x1=2.5 y= x12=6.25
5 J4 b7 _: }* z8 t$ ]( Y. V如图1: ( _! ]1 v0 s' L1 d. q: h
" c* C& o+ W$ N6 |& V" t$ B9 L" L
' K+ Y2 v: \) y/ X8 N, B2 b
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
! v; O* B- x9 Q; NMin Y-4X1-2X2 ! o" V% M+ V8 L; `( C! ?
s.t.
* Z! O* m' a+ W+ p) H( L+ WX1+X2<=4
% w7 w3 g2 L! J! y, u2X1+X2<=5 5 z7 p' N- {; N; K \# ?
-X1+4X2>=2
c' u2 x, I: U3 d# `4 SX1-Y<=0
9 t9 g7 r: t' m" e5 E( o. @* p5 \) T1 O; p7 k3X1-Y<=2 6 Z6 [' U2 @4 U. b* x6 W$ j1 H
4.5X1-Y<=5
F2 x* @0 w; C' J9 ?4 W& JX1>=0 X2>=0
1 K- b1 f2 [3 y/ J至此没有问题, 但解下来该怎样解就不知道了. ) f/ R& c' b8 ]6 H$ r- x5 Y# V
/ ~$ f6 n! {& b }# V【方法二】:
) K4 H) R/ M& i+ K4 k7 w取近似值的方法不一样, 8 I# }; [! c. x9 x& B
X1=0p1+1P2+2P3+2.5p4
: Y X& N- b) iY =0p1+1*p2+4*p3+6.25*P4 9 u6 U8 _0 \/ z* U
1 = p1+p2+p3+p4
. b; d/ y% [; L& b. W. m+ K
2 b, |' A' h ]- e6 h原规划可表示如下:
- T* {$ S; g# B J: lMin Y-4X1-2X2 5 d! Z) J" ?! N5 M6 D
s.t. 7 M4 p: m' H+ V6 U% ]0 N+ ]* K: a
X1+X2<=4 p) b) v O1 L9 k; U- C6 D& C
2X1+X2<=5 0 K! b, z( G3 P! m$ Q
-X1+4X2>=2
! E! M5 E% L# g" d/ L4 X+ y-X1+P2+2P3+2.5P4=0 & b9 V/ e, }. A- h7 l! T A
-Y+P2+4P3+6.25P4=0 / r- P. T; V% ?( n7 X
P1+P2+P3+P4=1 1 c- e1 @% t0 ?1 A7 ~
Y,X1,X2,P1,P2,P3,P4>=0 & B6 K0 c; p+ P. _5 d( R
同样的问题,到这里就不知道下面怎么解了.
4 x2 p# q) j( X. ?8 N图2
- W# p* G n% `! M Z+ B ' P5 z# z/ g! B9 ?! l
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
4 |8 F9 A: u$ h/ m9 l9 b) ?" f- \$ \% h
$ J( o9 p* ]6 I* [3 i3 B! J
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
+ O$ w7 _# Z+ r3 P% \. ^% k |    |
zan
|