- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
. @! u) V+ N( U9 _
! b" u/ G3 ~1 e( s$ q U# ^+ i
@# c5 P, a; K, e" g1 Q
2 d! F. n7 I9 m9 r N( Z: F! H| 再求数学高人 |
# ~" J* b' {2 R$ Y. q
, P: p t% ]) i+ y | ) Q5 I2 w* U% c( e
. I$ K* _1 h/ N7 G% C2 P3 h! P6 b- |$ l4 Q
3 d7 H+ Q" x5 s8 ]: Q: R5 Q6 Q# x0 ]* \2 X
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
% V6 A' ^1 V( b2 ~) h# FMin x12-4x1-2x2 . [; b9 y& i& c) t+ q4 @" D' `
s.t.不等式组条件如下
; B7 g% G) {, n, I& xX1+X2<=4 9 \- t/ i5 e( F5 T% k) b3 p
2X1+X2<=5
, x. j5 b7 U* Z f-X1+4X2>=2 & Q B9 {+ w# x- I1 T
X1>=0 X2>=0
; n& a5 T& H6 m& `) _. r: X【方法一】: 9 T/ v5 ~7 @4 G7 X: ~1 j
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 ( [: W5 N, Y3 k' m% M6 u3 o) @# t
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 9 W0 ]( L+ e( Y6 o
点O x1=0 y= x12=0
- g! k% k1 q$ a; h点A x1=1 y= x12=1
7 k' Y# P* H! K$ Y4 |点B x1=2 y= x12=4
. {5 E( F& W W$ S( F点C x1=2.5 y= x12=6.25 3 h' |% J" A) k6 K
如图1: g; t4 q' Q6 _# J- \ f
2 h0 G+ e. m5 N; }9 @6 b5 {) ]/ t- s! v% }2 \. G' V
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: 4 q7 t6 y8 U4 g7 o$ e% I9 i) M
Min Y-4X1-2X2
# o. S8 x+ K9 i$ Is.t. 9 ` ]; `! w7 X3 F' }8 {' ^
X1+X2<=4 9 }/ ~+ s" U4 W# g7 ?5 S) v! r- g
2X1+X2<=5
/ v/ v& }% H! z V8 {' L7 G- L( p-X1+4X2>=2 1 D0 }9 U8 B0 l; T3 k
X1-Y<=0
! _; t5 l, D1 X1 i' l3X1-Y<=2
' d' k8 u+ k" D1 o" M4.5X1-Y<=5
1 n1 g' K; ^) U2 l, aX1>=0 X2>=0 7 M5 b$ C: D7 j& [& m4 l# r" P) H
至此没有问题, 但解下来该怎样解就不知道了. : a& {3 a) E( p' h
9 d) b+ Y" O2 g* O【方法二】: / [" b9 u# h4 y6 k$ b- Q( a) G' y; q
取近似值的方法不一样, * V. c4 p; A9 F4 M
X1=0p1+1P2+2P3+2.5p4
0 v% c. |" R/ y. |2 O+ m2 LY =0p1+1*p2+4*p3+6.25*P4 ) m4 T e; n9 o' N1 f' E
1 = p1+p2+p3+p4
& Y& y( |) F# w5 I) S8 }1 ]0 u* _+ F! E$ X7 z5 {; N& k
原规划可表示如下:
7 C4 ]4 P! v" Q3 F& q- `' nMin Y-4X1-2X2
$ X; ~0 D( n7 \5 W6 os.t. * `: N" M9 P, `8 P4 o
X1+X2<=4 ' Q# h8 m5 M+ A& X8 h. |
2X1+X2<=5 : u& j) N" S' R" M1 K( U
-X1+4X2>=2
; i4 L: b" D6 R6 e! c# `" q-X1+P2+2P3+2.5P4=0 % g9 h3 P6 U& H" W4 I& v
-Y+P2+4P3+6.25P4=0 {% p' s7 r9 N5 u, B( C+ L
P1+P2+P3+P4=1
+ C3 n4 G4 M0 ^" A# v7 Z1 d, DY,X1,X2,P1,P2,P3,P4>=0
$ f) x4 ^. |2 x) j同样的问题,到这里就不知道下面怎么解了. * B4 X+ U6 A; _
图2 3 c. N1 n* K* z% Z. c
) ]8 Z" h: C( P# i8 m" U$ W( |. _
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
/ Y( ], h* V+ b0 G" h2 N |* \3 F( F# @8 j; l
& E6 W, j+ E# ~ M& ]8 F& d
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
* s2 _, j j! p |    |
zan
|