- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
4 I7 X0 ?+ R% H+ e% u, S
0 w+ ]/ G8 W" a4 Y. p) B# b
E Q# e; K# Z- p7 @1 g2 A
* w ~1 h8 F- R0 p7 ?| 再求数学高人 | * u; j( {9 @" Z. N; s* \4 R. g
8 L H9 i2 l3 f7 S# t
|
N3 u& t) O/ n$ v
0 t# f' O4 \6 s, ~/ m4 E
: P/ A+ v a8 R+ N! @% k( T8 q0 R
0 O: ]4 i. @1 y) n8 m! U用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), % Q# @9 X: ^' {7 }% Q% u0 b7 H
Min x12-4x1-2x2
* a- V+ E% W/ b! E. @* ps.t.不等式组条件如下 : i" q1 t7 J: x- ?& P
X1+X2<=4
+ l, K: {: a1 O& w4 Q; D0 c2X1+X2<=5
( D" j$ K% Y ^. w-X1+4X2>=2 ( b8 e' v/ L) y' k6 E b
X1>=0 X2>=0 5 }* ]2 T) r D V" E, f9 Q
【方法一】:
0 H- N; U0 K6 T- f因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 % ^7 F% T! t0 `2 h
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 3 X, \; z6 A3 A# u8 k9 k" z8 N
点O x1=0 y= x12=0 ) j5 z# o0 U# G+ S* u9 F
点A x1=1 y= x12=1 # C% W8 F7 E, i- C* F
点B x1=2 y= x12=4 3 _& l: H+ z. t- {
点C x1=2.5 y= x12=6.25 9 Q0 Y( \# |- N) W; J0 Y, n
如图1: " O5 l& |; k# W3 L3 ^
% W' {1 M5 p7 f/ z) K
2 v& M3 b: B+ K: v% _1 N& ^. `) x用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
$ O/ t( R; v8 r( H* `% {" d! o9 t" kMin Y-4X1-2X2
: Q, P2 S3 Q$ \2 {8 fs.t. 0 c: Y- y( ~ E1 E5 |( u# q E
X1+X2<=4
@; _) e( Y7 M2X1+X2<=5 1 o b0 m2 L4 o! G8 X6 x1 c
-X1+4X2>=2
+ c0 z* m3 q; gX1-Y<=0
( \" f: c& m7 H& K; @/ ~ n8 J3X1-Y<=2
6 Q2 \/ U" P( S, S5 r4.5X1-Y<=5
5 g1 a& k, Y+ G% L- {2 H" [X1>=0 X2>=0 & ?$ C) V8 {. k+ G- a- W
至此没有问题, 但解下来该怎样解就不知道了.
4 e4 V& c/ [, S( {; S% [6 y# E+ P ' t h/ ^5 N4 x5 M8 y3 _' b
【方法二】:
1 _+ k$ r# k2 b2 E6 R8 F# c取近似值的方法不一样, 6 S1 n6 a( T: S" D) G O
X1=0p1+1P2+2P3+2.5p4 * R% n# N, G9 k1 J9 L1 m
Y =0p1+1*p2+4*p3+6.25*P4
0 l& }" s' m. E, V" L/ P1 = p1+p2+p3+p4
' C ~& R* L4 [9 q, X: A) r
: J9 G" N/ z( Q- V' v原规划可表示如下:
) D7 Z4 {* B5 P, n2 sMin Y-4X1-2X2 ! e8 A# \! u! @7 g: M
s.t. ; I9 T+ e0 B6 \5 z T9 w
X1+X2<=4 ! D) E: m2 R7 I/ |: Q
2X1+X2<=5 ) f3 j4 h8 P" T
-X1+4X2>=2 " k. p. d. q$ O5 V) H P
-X1+P2+2P3+2.5P4=0
# M: U# E. }4 J-Y+P2+4P3+6.25P4=0 5 e4 Z5 t0 M- I; N# `
P1+P2+P3+P4=1 2 Y/ a/ t8 x1 u* F: @- t6 `
Y,X1,X2,P1,P2,P3,P4>=0 # u' x5 A/ |' n* g1 y
同样的问题,到这里就不知道下面怎么解了.
7 ^) v/ x1 A2 M4 y+ _# c k( e图2
- F, z/ }- k$ y% F# s4 b% F
7 r9 @% H, K) O' o原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! : B, k `$ s( y% w8 l) ~1 j l
! u( H( W( e) o3 @9 }3 u5 g: F1 g( O: U
9 d0 X$ b8 T. v7 D. J我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 9 d7 z# H5 r) C( j) H( z) T
|    |
zan
|