- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
/ e+ k# w, \+ c/ _
4 B) d% x$ j. G
9 p! f J. Q' E# I9 ~
/ i# y9 h3 L6 M# ^4 p+ R& {| 再求数学高人 |
9 V9 Q1 {$ R1 C' U7 }. k) d# m
' K: C% A; V8 g! N |
. v$ q# U5 Q& h9 A+ q; ]; p$ w5 c8 E, q; r
$ f t) _4 t3 r" O; c6 U2 \* M5 `
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
& D6 B" v; n9 ^- k2 Q1 K$ nMin x12-4x1-2x2 . b2 Q& F; C0 Y' ]# c
s.t.不等式组条件如下 $ v% ~, P4 ], [, E5 g
X1+X2<=4 ; k' Q* ^( I" f4 s. `
2X1+X2<=5
. I, o8 u2 n% i1 I, n$ g-X1+4X2>=2 1 _: H! Q' W- L) S/ k; H" V
X1>=0 X2>=0
& Y* ]3 h& J ?8 t& y" {【方法一】:
1 ~3 x5 F( b! p因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 8 V J& C( L- I. _+ V l/ y# A
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: ) _) l& o( b( T& O
点O x1=0 y= x12=0 ' x" _5 U E0 w
点A x1=1 y= x12=1 2 o5 a& x4 K4 ^; ~4 Z
点B x1=2 y= x12=4
9 o& b. j- Y- b* r* \点C x1=2.5 y= x12=6.25 " ^% J y9 y. ^0 d5 m/ W3 h
如图1: , @, T& U3 ]/ u5 E
3 q& \; |9 O) O* r( n; p1 w
4 z/ p' `% y% W1 I% p
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: + v6 o q8 u8 ]! z: N
Min Y-4X1-2X2
! B' Q! ^# g( ^: S( v* D8 fs.t. ' L, Q$ z; u0 k* K! n! h
X1+X2<=4 3 l6 ~( \" q+ U; P, W& y
2X1+X2<=5 " ?( d; C$ b0 x. _( `* c- ]
-X1+4X2>=2 ( H$ t# {9 i. v' O$ g% W
X1-Y<=0 l* o7 a8 s' o* k# `
3X1-Y<=2 7 N t& M y) a; S
4.5X1-Y<=5
5 v. X; t# r% j3 y: G! BX1>=0 X2>=0 9 \, r( D" q# e a
至此没有问题, 但解下来该怎样解就不知道了.
0 D1 a7 x4 \! E4 Y1 v+ B
# c; p. C6 p3 |1 |【方法二】:
1 M. G; A- G* v! ^ N. s0 d p取近似值的方法不一样, 2 x+ O' k5 c3 U) V# N, n( K8 |
X1=0p1+1P2+2P3+2.5p4
' i5 U1 w# M& R. _ R; Y3 M3 yY =0p1+1*p2+4*p3+6.25*P4 ' j; v! Y/ D8 Y% r6 D) k; ^7 o
1 = p1+p2+p3+p4 " O: X; d A. A" w0 N, V& E
- S- O- _: d: \, L; \原规划可表示如下:
9 @) q: ?0 c, N2 VMin Y-4X1-2X2 ( d# H5 A% p) j
s.t.
7 S/ e; C. _* @9 S: ~2 p5 qX1+X2<=4 4 a2 f/ q* B& u# H) p! `
2X1+X2<=5 + A. n$ u7 b4 F6 }
-X1+4X2>=2
; @/ _0 _' T; C-X1+P2+2P3+2.5P4=0 8 j/ M$ E1 F# a8 I9 Q: y- B9 r
-Y+P2+4P3+6.25P4=0 - A" z' y! C( L/ X. b
P1+P2+P3+P4=1
/ c$ e' l3 s% VY,X1,X2,P1,P2,P3,P4>=0 " ?9 d. X- d5 z) p% S' X% Y
同样的问题,到这里就不知道下面怎么解了. , L" x- R6 f5 F9 J$ }% [
图2 Q9 o; E6 e6 Q5 p
0 ~0 R2 F" z9 D4 s原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
- a( |0 b6 `; C" ?4 T8 N/ w% Z8 { ~: {% @" B* y* d% B
# A |! O8 D- d: e+ }- u8 h
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 4 c' o' m/ `0 Q- N5 A
|    |
zan
|