- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
' W9 y' A& W& E. n c
+ `9 Y1 H( m4 V! ]5 j+ T' o, l7 T" n8 m9 X4 x
' B- t8 D- Z+ E: ^- x, M0 j" \* b| 再求数学高人 |
( H& b L; ?$ [! \) J3 E& G( f& o/ j, h" H6 P3 z* {
|
$ B+ p" o5 S- q8 i# ^* R" L& V9 W) d& O5 }1 |5 T! Q
, ^- Z$ ?4 L( o% [
5 \7 k; f2 D1 K a7 @- ]" u
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
( P( g& Z, x+ b" q( F) ~4 { zMin x12-4x1-2x2
% G+ Z/ p9 o3 h* k5 s' t" s) n! m$ j% cs.t.不等式组条件如下
+ [. l% @5 S, B2 U: [3 ZX1+X2<=4
, T- s, o m/ c. |2X1+X2<=5 4 }9 G/ y4 `# c. U
-X1+4X2>=2 3 E2 V y- x8 x8 V4 _$ P8 e
X1>=0 X2>=0 . O* X- F* z6 M) p3 |4 O1 m$ w
【方法一】: 9 y8 B2 `2 r% N9 `- J
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 $ \( \3 b# b4 L& f+ f$ U
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 2 P$ T, i6 C' a0 k' n7 w; j
点O x1=0 y= x12=0
7 h0 `, X4 {9 }- `) Q8 i点A x1=1 y= x12=1 " L1 @: z! |8 @
点B x1=2 y= x12=4
7 F v7 t0 o4 o点C x1=2.5 y= x12=6.25
: f, o8 Q" a( a- y如图1: 1 j7 A- t2 b0 A- e* P1 d
: U* F3 n/ U* X2 Q# m& s
/ V3 m8 T9 m+ O1 u3 w5 m& A( S* v
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
5 i" B) D9 O& SMin Y-4X1-2X2 " O3 C {5 C; _2 z. x
s.t.
% ]9 a$ B( h8 S8 AX1+X2<=4 ; k0 i; I8 {5 |3 c }: Q
2X1+X2<=5
. }# w+ t: S8 Y5 F( p& G( q& {-X1+4X2>=2 8 O) e$ ^; X7 }% g6 U
X1-Y<=0
. X* S) t4 ~( S* Z6 S1 A/ |! L3X1-Y<=2
) W$ w. Z4 i. i* g9 G$ \ F, m) E1 H4.5X1-Y<=5
: y; {( {% K) h# mX1>=0 X2>=0
# z6 J3 I2 x6 j1 B4 C+ ]. Z至此没有问题, 但解下来该怎样解就不知道了.
* ^6 _3 p1 E0 l: h
5 I+ k! m6 v: k9 P, a( ]7 ~# p9 [【方法二】: 1 W& H5 }7 A* m9 N1 D1 h
取近似值的方法不一样,
% `. U2 U4 S$ s& l) V4 QX1=0p1+1P2+2P3+2.5p4 2 D' a, z8 J" D& J
Y =0p1+1*p2+4*p3+6.25*P4 4 o, e/ R! b2 Q, f
1 = p1+p2+p3+p4
, q2 P( v- N% H. t9 H1 R
% G1 L' L1 u4 n E. x) N原规划可表示如下: . P1 C1 s+ _4 ?% S% D( m! k' t8 R
Min Y-4X1-2X2 9 U' A- Y X6 F) n: T$ @( {
s.t. ( _: Z: D; ^% ]1 J
X1+X2<=4 2 R; _6 a& y- _& N! [- l8 h/ V
2X1+X2<=5
# E4 b. N4 I7 q5 ?; V6 y" P% M-X1+4X2>=2
4 f, u# N! E2 c-X1+P2+2P3+2.5P4=0 % y4 ]3 R7 v; j, l" ^* k
-Y+P2+4P3+6.25P4=0
8 [3 x0 v0 G* t4 A) uP1+P2+P3+P4=1 ; H2 i) \- S. c) \
Y,X1,X2,P1,P2,P3,P4>=0
. D) b8 ^& N' ^6 |& y$ u0 Y9 P同样的问题,到这里就不知道下面怎么解了. 4 p x# v8 A' E7 G' Q
图2 Z: A( A- w; E; ]
, _5 g. o4 x2 ~ e. I) K. r4 b0 c
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
: Z: d( i6 {% y% B( D
/ G+ M! i. ?# r6 y2 q
! U6 R. {* H9 k我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 ) C# Z( ^- {! v5 m1 q# {9 z' @+ u
|    |
zan
|