- 在线时间
- 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( T; S5 C2 C2 N + C G* B: p; Y* W
* t ?/ H5 R4 L, W8 V# C. R5 S
4 I8 I0 m* G/ W2 V& D
| 再求数学高人 |
6 x% I! ]/ U! @
8 V/ p* K2 D. x3 a+ r0 d |
+ O/ M+ b8 t' Y) d, u
- g3 x; s* }4 y+ |/ `: m7 l. R# a7 x! F
! B+ l4 a2 J! C7 Y8 C4 b用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), 3 |3 V4 I3 w- q% N6 B1 W
Min x12-4x1-2x2
6 { D9 k( @' k9 y. C3 p9 Xs.t.不等式组条件如下 - g8 Z1 P, z7 v' k5 B; B
X1+X2<=4
! h. V; _- |% `) e9 w2X1+X2<=5 8 z) R& g) ~. J" X. F0 o- k7 F* e1 a& b
-X1+4X2>=2
! _- Z( J. ~3 e- iX1>=0 X2>=0
; b& q g) Z$ {# U) [$ v p【方法一】: 7 D, ^3 l( O7 U8 c
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 & N9 Z2 [) H$ I1 R
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: \6 \( l; N& T2 Y" ]- {: X
点O x1=0 y= x12=0
& j) J4 i) H, T) ^% r' f. L, y4 H点A x1=1 y= x12=1
* j6 s% Y5 d. H5 p* l) B点B x1=2 y= x12=4 3 X5 B- q; V: G
点C x1=2.5 y= x12=6.25 9 }& h$ ~9 G9 D" V( ]
如图1:
X* ?- ^& j+ f. B9 |$ i* B ^1 f! G! X. W% ^, a
) n6 h4 }: y ~7 Y- u9 } y- O" g用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: & u# a( P; n/ l# j2 D: B
Min Y-4X1-2X2
, P% W% _! n+ Z0 b6 D' [+ `6 p4 m! Xs.t.
6 U) x* b3 `/ Q2 U& [ F/ w' Z$ SX1+X2<=4
g7 a8 ~, S3 j" f1 R2X1+X2<=5
9 Z6 G) U& ?' [7 V- \! {% l-X1+4X2>=2 * V# N/ ~3 n, S [" Y
X1-Y<=0
& v' @+ j3 W, T# r3 y$ w v3X1-Y<=2
, h0 G0 K% _* c5 T; ^. x, V4.5X1-Y<=5
! j. P# n' f3 S3 w7 }X1>=0 X2>=0 ! U5 ~& |) o. U2 r; P% l- l
至此没有问题, 但解下来该怎样解就不知道了.
$ E; r+ _3 }# U9 A$ k5 s6 T% U
, K4 R5 M! _% s8 R& P! S: F: m. t% n【方法二】: . t. l1 y) }! @( }2 F
取近似值的方法不一样,
5 a U& {7 c% \! G/ x) A- Z% d6 rX1=0p1+1P2+2P3+2.5p4
+ O) y: \- {# M2 {: t JY =0p1+1*p2+4*p3+6.25*P4
: D8 N* T5 b7 h/ Q- q8 B1 = p1+p2+p3+p4
8 C& b0 U# k6 P+ y+ w( @5 o& l( y+ O
原规划可表示如下:
" k, z: l" S4 N! {Min Y-4X1-2X2 4 U& V" |8 E' O. S
s.t.
$ ?* f; T3 J/ ~; r. aX1+X2<=4
! q4 |$ h9 y1 k$ ~2X1+X2<=5 a, S/ t) s5 P& r, ^# q
-X1+4X2>=2
! M; S7 f8 y* t9 C0 w; @( @-X1+P2+2P3+2.5P4=0
% O* R9 l1 K4 c8 D$ H( ?-Y+P2+4P3+6.25P4=0 0 Q" B( F, D$ B7 ?( [2 P
P1+P2+P3+P4=1 5 m3 r- X; d, ?: \
Y,X1,X2,P1,P2,P3,P4>=0 - D/ H; n4 u! L4 E
同样的问题,到这里就不知道下面怎么解了.
4 E3 \" y9 J4 r* y2 U4 m, D图2
4 E. Z" H7 X, T" A/ T* o& I 7 ?5 ?7 ?7 D7 p% s6 W( q
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
& ^* o5 I, v( I/ x& U: B: v: F, [0 K! @
" F2 y6 y6 `. Q) c0 D7 `我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 5 n) }3 w" _# y1 P& q
|    |
zan
|