- 在线时间
- 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& _) Y" Z3 r, [3 m; P
9 E& Z* X; U! {, L0 R& ]' I
" g7 ?6 ?6 \# o P
; w% L# i& d& ]6 V6 D再求数学高人 | ! l9 H/ B' M, x, I0 c* G+ A
! B |' N- Y4 U |
( }( P$ Z; f4 i
$ F( }3 |- }- c0 o5 h b
1 h' o/ b2 P5 x- K4 Z! ]' L$ d3 h% r8 L' Y9 [
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
- b. @7 d) f, o6 {% H& k7 B QMin x12-4x1-2x2 2 b- p! R& |1 Z" i; h3 @
s.t.不等式组条件如下
4 w9 s) N/ z5 L, ~4 gX1+X2<=4
: L. s4 @ v9 f# g2X1+X2<=5
1 g- R: v# B# H4 u0 t( C) K2 Q' c-X1+4X2>=2 * U# A) U; R5 j, J
X1>=0 X2>=0
- i/ E% z, M8 J6 U1 c6 A- b9 A【方法一】: 1 l& r9 L$ q2 S( Q9 m* M" C/ {; S
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
! A4 \& ?+ y2 y: o% ^& a+ ]首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: % p4 d4 V0 ?' B u- R% v
点O x1=0 y= x12=0
- n# e0 j, O: s- X点A x1=1 y= x12=1 ; x2 E1 l- y% ~$ a. V* \
点B x1=2 y= x12=4
+ n/ G* i9 A6 B( |3 d$ {点C x1=2.5 y= x12=6.25
' P9 O# ^) L0 y如图1:
0 z | X5 q+ ^0 ] " i5 _& D- R! j- _, n. E( m! _. P
& Y' W$ D9 C" S% u用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
4 I6 e0 z% `* L- X: KMin Y-4X1-2X2
* M* f/ O! ]& P/ g4 [. c$ D, ms.t. 8 y: L. m& S- m; K' l+ N8 u
X1+X2<=4 . V ?2 E5 s. k3 Q
2X1+X2<=5 " k6 i! k, ^0 d' {
-X1+4X2>=2
. B8 j3 R% N: v7 i T- bX1-Y<=0
0 t& n0 T4 g) J& Z3X1-Y<=2 2 O& n+ h' d8 P* M& q T, W y
4.5X1-Y<=5
8 J/ t, }) ~2 y: Z% cX1>=0 X2>=0
- K- r3 l, ?0 U9 @7 g, y4 M7 V至此没有问题, 但解下来该怎样解就不知道了.
& { o T# n& Z y7 ?/ M9 x 1 [. M( E% |# S' T6 ^
【方法二】:
7 a. ]( \- y s( w5 ?3 g* C1 K6 b& I取近似值的方法不一样, ( q7 z6 G# {$ t2 z; \1 ~" ~
X1=0p1+1P2+2P3+2.5p4
( y2 L9 e) k+ x# v5 X: aY =0p1+1*p2+4*p3+6.25*P4
+ i! R8 D1 k, E1 = p1+p2+p3+p4
9 i% s5 }4 y1 k
`6 N F% U: `/ J原规划可表示如下: 2 O- i$ i8 A$ {% K$ H6 f
Min Y-4X1-2X2
6 _- A) J, C: n5 v& es.t.
/ C- P: P( w L1 ?/ Z' OX1+X2<=4
X6 f( ?: L0 }0 F: x2X1+X2<=5
% U8 ?9 T$ C0 L. ^0 X) w-X1+4X2>=2
) c! r: j3 P I5 [& ~-X1+P2+2P3+2.5P4=0
8 a# I& Z9 Y3 T) I( s+ q# R-Y+P2+4P3+6.25P4=0 8 h a J$ O5 `2 ~. l
P1+P2+P3+P4=1
) a" X6 e9 \" I q0 j& |Y,X1,X2,P1,P2,P3,P4>=0 ! L9 M4 D1 n! C/ {- A# a
同样的问题,到这里就不知道下面怎么解了. ; }7 c: e K" }9 h D$ }4 M3 \
图2 4 ^3 k/ K' ]5 l+ f7 z
" q5 x8 S5 w4 f. j7 R9 F! q
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! & Y# J# T3 C) Z" n* O
) ~# v3 x1 a! ?
7 ]. |5 ~3 z$ S4 ?我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 3 x- p, j$ p1 b s. z
|    |
zan
|