- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
) F# U7 H8 A' A0 K, ~8 l
' V$ H6 N# Z; r+ _) h' l4 `2 P& d" C1 F+ T1 h& H( _+ w7 n( F
8 s! J9 c) ]0 c6 u# p! r. V4 k2 Y
| 再求数学高人 | 6 ?3 S( ]% T! O# _$ p, _- V& A% X
/ W' n2 |# V; r a2 O2 W+ V
|
, z: W4 O3 W9 }
1 d7 _$ F5 J* n: O: F% B7 d( r* u1 B* x5 m/ C- f
* x% _; E @" ^/ f: p H1 X
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), 4 z( R# k9 s: e) q) K* e" w$ L
Min x12-4x1-2x2 " ^% e# K+ @. K. V6 T& W: a' j* S
s.t.不等式组条件如下 * T5 q$ [6 m2 v
X1+X2<=4
! W* \; p# b8 M" ?' S1 [/ S2X1+X2<=5
5 ?2 a* ]* j& |& l, [; h8 Y-X1+4X2>=2 p1 B O% R ]& r, z4 [" V
X1>=0 X2>=0 4 U2 c+ F' L; }9 y7 }# G: p
【方法一】: " E' u9 n# w9 s# i9 f( I+ K, ~" I
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
P5 T% p9 t) Q! @( d6 t' s首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下:
2 u; J0 v ^' p' E* B' N ~0 \点O x1=0 y= x12=0 0 Y" W. Q" l3 Y
点A x1=1 y= x12=1
8 ?: |8 I. B+ k9 |9 ]点B x1=2 y= x12=4 - m) U. i' A. Q7 {; ~; `! @
点C x1=2.5 y= x12=6.25
$ q% R! w8 [4 U3 U如图1: , r" R. C+ F8 V# r' @" j
6 {6 J0 S$ P' j: U9 t8 z
* J. Q9 N, R% ^
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
6 n8 z+ `) j* m( M/ t$ j# A8 f) R4 Z( aMin Y-4X1-2X2 4 ^7 X; Q ]7 V! l# \# s5 g
s.t. : ~, Q+ J0 X2 L6 B5 m7 E* Y
X1+X2<=4
$ |0 o3 F# A" O0 q# X2 U2X1+X2<=5 : R8 \, u! r, }% X1 g7 H L
-X1+4X2>=2 8 W) e* ~2 s* ^8 Y U1 s
X1-Y<=0
. T' y- q$ q E5 E3X1-Y<=2
+ ]! J x% g, i9 E' [. y4.5X1-Y<=5
$ y4 c- U* Z) @' K lX1>=0 X2>=0
' H' J8 a9 ?4 | q- E B至此没有问题, 但解下来该怎样解就不知道了. 7 j+ I, R. ]& \3 f/ s- D
; U: C; G0 o$ H5 c: ^8 ?; o1 M/ ?【方法二】:
# `; d9 I- `( p' F& v& u# J取近似值的方法不一样,
& Y9 G3 c/ ?* i& P! uX1=0p1+1P2+2P3+2.5p4
( K/ H' L# Z* L# @/ N5 T5 KY =0p1+1*p2+4*p3+6.25*P4 . I9 {) O d2 a* v, A5 M
1 = p1+p2+p3+p4
' w) M, x( |6 Z0 q4 O# ~; A
& r2 v% y, R9 W1 t" ^原规划可表示如下: % u# L# \: B) H" w: s. Z
Min Y-4X1-2X2 * {. K" U% W( I" |, f/ m
s.t. 9 [ y8 L* d& }0 l4 }% Y
X1+X2<=4 1 O8 R- V( ?7 u `3 x2 Y
2X1+X2<=5 * z" {# k. n! s4 U" N6 v
-X1+4X2>=2
0 [5 s( k4 G" k$ A3 D& U-X1+P2+2P3+2.5P4=0 ' o* } a* k, M( p" c
-Y+P2+4P3+6.25P4=0 2 Q+ B6 w1 k$ I( ]2 D
P1+P2+P3+P4=1
2 U8 m! t" W: v# f% wY,X1,X2,P1,P2,P3,P4>=0
5 y" R9 y \8 Y同样的问题,到这里就不知道下面怎么解了. 8 [/ u3 Y7 r; B* g
图2 I8 G( \6 Q% L. M% C X/ ]
/ s2 H* F: T ~- I3 M/ H
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! " E$ C4 A- }* J/ X7 p& {8 k* P1 Y& V
+ H' F2 ^" R+ I9 [ Q3 @3 a6 ~- L* K1 r# E! F" D
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
" e6 _8 d0 P/ Y/ l* u6 z |    |
zan
|