- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
" T' f5 c& D* j, W+ H5 I; D1 o+ A4 v5 j
5 p* t( P# U. r
$ f/ j( A) f9 Q" z
b+ }7 o' @2 {, x0 b: U| 再求数学高人 | 1 C5 p8 t. @( \
. B1 O* H0 f$ r, d$ p- |6 q. q$ N |
4 l* k7 y- O( U6 n1 V+ @
6 i: I8 e+ W5 G9 |3 |, i4 M6 i
. [5 T5 i$ ?6 X4 p
3 A; t7 V( s, v7 ~' e用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), 3 m! o7 S' o/ `: Q) y0 e
Min x12-4x1-2x2 3 \, i& w3 e/ m. i/ a2 g* }5 s& i
s.t.不等式组条件如下
2 L' T1 c4 x9 C5 U, D9 I+ |5 j) S% DX1+X2<=4
* B8 [; f) M2 h1 {2X1+X2<=5 $ e" a* l9 G! u- c: r6 [
-X1+4X2>=2
# U( _; V8 q7 R) J+ A- r, CX1>=0 X2>=0 8 j- z: ?0 I/ C' t; N
【方法一】:
7 c: e9 m3 P- p( a5 D" t5 v7 D因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 " _& s1 P$ i# Z) W
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: + A/ |6 i8 g0 R
点O x1=0 y= x12=0
! W# ]0 }$ d" ^! n7 j点A x1=1 y= x12=1
. B/ `: i, I# H8 m/ H' [& E点B x1=2 y= x12=4
$ C2 M: }! U. h点C x1=2.5 y= x12=6.25 . D% h A( A+ A' q" ?. |1 ]
如图1:
! R. K4 Y; j; z9 d6 M- {
/ n1 Q' a7 ^3 L$ t, K9 k# x p4 Y
6 m5 P2 ?# W: H- t$ t用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: 7 D! z z0 z* e% H
Min Y-4X1-2X2 + f) \! p; p, K/ {
s.t.
* U' g/ y' c% E3 I0 `X1+X2<=4
8 n* j# O3 s2 W* a2X1+X2<=5
) d) V2 m( h5 f" j: |4 c-X1+4X2>=2 - ?" h) R5 B# S5 c- G# j" t
X1-Y<=0 * M7 i# S7 s$ J$ I# e0 N7 p
3X1-Y<=2
2 R- f- ~2 z) D. X `+ t4 f( G4.5X1-Y<=5
) K- W8 k; w g. V8 U# o9 `X1>=0 X2>=0 8 r/ S5 W. o' t6 o7 O
至此没有问题, 但解下来该怎样解就不知道了.
5 E# h ], d9 x( S( X" Q
- K) X: f9 U% }1 P1 j9 d- P, n【方法二】: \8 i: d; w' \0 T7 h) w
取近似值的方法不一样,
( u- G n% t) p" p' i4 R/ OX1=0p1+1P2+2P3+2.5p4
- \9 X9 ]( L& B6 P& Q* {+ ~ u+ JY =0p1+1*p2+4*p3+6.25*P4 8 p) U) \" j& h
1 = p1+p2+p3+p4
) ^( b) y) i+ {7 i
! T& V" B$ U0 I3 Y原规划可表示如下:
; _! |, m c* ]4 RMin Y-4X1-2X2 + w& U1 y* F5 l' y. x4 I5 E8 x
s.t.
+ G# @: r6 G9 ^8 r, w& N! `X1+X2<=4
/ |# F; N0 S4 I# Y" I5 R, X2X1+X2<=5 0 _. Y0 S0 H& t8 [. { j. E
-X1+4X2>=2 3 \! R5 }* T0 ~9 }, J" V2 w' v! ]
-X1+P2+2P3+2.5P4=0 + H1 q: Q& a+ J
-Y+P2+4P3+6.25P4=0 ) D6 X5 o7 [+ w9 N' {( d: V" X
P1+P2+P3+P4=1
" }! Z! R5 [4 t; QY,X1,X2,P1,P2,P3,P4>=0
2 A: @& L( [* K" M9 P, f' c同样的问题,到这里就不知道下面怎么解了.
% B- t/ K" w: N1 w图2
' C% L: z% W5 W# z! e 1 b k' \2 n" Y- {3 T
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! 1 I% ~1 _% g9 J& o5 M6 y: t
4 W4 ?7 J; v, ^4 x$ T# Y+ R6 u' O8 M- _2 [1 W
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
1 W8 } Q0 T$ g) K! b# U |    |
zan
|