- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
H0 y5 W% ^2 k f
7 @8 ]5 n# Z' T# ?
* E. {4 c# Y% C o! m7 p6 K' m
1 ^9 `2 k! H. j# k) s
| 再求数学高人 |
+ s9 b8 @& v6 F2 p1 B: U P2 x7 F1 \2 ?+ a& i/ A# {
| 8 T* B, J0 J* t$ ?! r
; F0 Z2 A+ S% I( \3 u; V3 {) v6 Y. `6 A. Y) s. [3 T
# ~+ o# E# A4 U5 ^2 @! h% \
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), 7 z8 r2 |3 c5 s2 d8 U
Min x12-4x1-2x2 3 b( s6 j5 ~9 s7 g1 O
s.t.不等式组条件如下
W. ~( C! W/ y5 V4 fX1+X2<=4 " a# G# a/ c, b. A/ E0 L
2X1+X2<=5
}3 l! W' q+ `-X1+4X2>=2 3 {1 _# ~+ P9 l( \
X1>=0 X2>=0
- _" O" w1 B0 C9 f5 q( ], E: l【方法一】: ; W& r5 K; g$ ?( P4 T% S. r: ?
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 2 _8 ^2 T) {' Q2 m$ ~! e
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: ( `, z! ^! @% d# i) y, D
点O x1=0 y= x12=0 * `2 {3 R( U" C8 Z0 @* w1 h% v
点A x1=1 y= x12=1
" `* C9 e) F6 w' ]$ ~! |点B x1=2 y= x12=4 % w0 V8 i& i4 \0 G: H% m: A
点C x1=2.5 y= x12=6.25
( O4 O) R1 Y, L* C7 s. W如图1: . [6 @. P/ ?. _0 }
0 u" ?+ l N+ X
: p: c. V. A5 n; l% l4 }用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:
3 S/ p$ |) X+ x! ~0 k" LMin Y-4X1-2X2
. y: r; ]' a. c0 b# Ws.t.
/ v3 u& @- Q! \X1+X2<=4 % k h; O; ?- o" j! |7 X4 @; U4 z
2X1+X2<=5 # A O0 L% L$ c5 S
-X1+4X2>=2
+ Q3 V8 Y* V$ \) R( g* cX1-Y<=0
4 a4 C5 ], X H0 n" \% h3X1-Y<=2 0 J6 Q% ]! O' ~3 \7 s, x! R' o
4.5X1-Y<=5 9 x% S8 z p# r8 r( F* @) U
X1>=0 X2>=0
' ~* `! k2 }" j5 v' f/ ?至此没有问题, 但解下来该怎样解就不知道了.
2 C% F( H# ~+ C& @* n
- M0 w, l' c% l: F! n【方法二】: 1 N" O- Q0 c$ Y f+ t Z
取近似值的方法不一样,
( n+ b' |& |" yX1=0p1+1P2+2P3+2.5p4 5 t: z9 X8 D1 x* ^/ `8 b1 w
Y =0p1+1*p2+4*p3+6.25*P4
6 k' M! B& s3 ?* K5 c" i1 = p1+p2+p3+p4 0 b9 m; ~/ v4 _
" y. b1 H& ?9 p3 N: z2 I: I/ X原规划可表示如下:
% t0 W P" S W/ UMin Y-4X1-2X2 9 e0 N+ N# T# t
s.t.
# @+ o+ v V% Z2 j. DX1+X2<=4
; S! S. D1 ?2 k5 L2 z; ^. m: o2X1+X2<=5
7 h" f$ s5 V* H( q# ]-X1+4X2>=2 " I/ A2 a7 q( u+ ]. S
-X1+P2+2P3+2.5P4=0
" W! B! g* ^" p5 |-Y+P2+4P3+6.25P4=0 9 O8 v; N3 J8 n% m
P1+P2+P3+P4=1
" D& Q! i# D# ?# j( i) a9 NY,X1,X2,P1,P2,P3,P4>=0 0 e; n" x" ^: Q) H8 _) t9 q9 g
同样的问题,到这里就不知道下面怎么解了.
% m. p2 _8 @. a* c" A图2 1 a1 O% d% e% ~
) K/ U+ n6 @: d, L
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! 8 O. O) H) O: q( W8 \
/ ], z( t$ {2 ?+ O# K. R
: l# R/ _& \" a$ Y" ^我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 9 w2 O0 f; I' a6 Y3 O T4 {3 t
|    |
zan
|