- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
8 [+ u5 V, ^2 J0 T8 }4 O* n9 S" _2 j % K3 [2 y. I9 o0 r3 ~' r! ^
. L1 b6 H, A& p( V5 P( {, p' t
n9 K+ | g. l" f| 再求数学高人 | 7 Q; r3 P- N9 M0 {* S" g+ p- S8 ]
& l$ T$ b0 z, y! A+ ?, I& Z+ F
|
S5 e; s: S3 d; Z3 C& f/ D
3 B: [, }1 C: v) ^; m: _& `1 F" R) Q& Q/ s) b0 e
# T6 }5 r1 b' O" v
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), $ _* k; T1 _) j- o7 O
Min x12-4x1-2x2 ' z- G0 S3 F/ s; R$ @0 v
s.t.不等式组条件如下
f2 V* i2 J5 Z8 k0 L( PX1+X2<=4
w8 r! c" E! N$ j2X1+X2<=5
6 C8 [& n; m0 A0 ?-X1+4X2>=2
/ i* C, d% d" XX1>=0 X2>=0
& l1 \1 c8 V/ M3 u【方法一】:
6 |5 P. r1 o2 O/ R7 Q2 ^( ~) b# Q7 h+ F因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。
' m, X; w& x9 a" ~4 c7 V首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: ) [5 ^4 m* K" z
点O x1=0 y= x12=0
' a' H5 p1 V' X, t2 `/ y ?/ L点A x1=1 y= x12=1
9 G' D$ @" m: p8 G' ^/ ]$ W0 N4 {& E点B x1=2 y= x12=4 4 r) x# d+ N9 R. N+ i/ `: W
点C x1=2.5 y= x12=6.25
% _* }) M2 y% f2 _$ S) x7 L4 W* i如图1:
$ K; _. b" Q% B3 Z' p) a. |* a
2 U! {4 l2 R* h; { }: v$ J: F( z( u& x
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: 2 E' k* Y8 l' T. b0 C
Min Y-4X1-2X2
R* h$ o5 R) d# R! @s.t.
2 G# W$ }: W: T8 Z! ^/ k8 xX1+X2<=4
. T/ M% D( Q7 o' ?% p2X1+X2<=5
0 M4 j8 ^0 L* n5 M* |3 }& C8 f! V-X1+4X2>=2 5 T1 v5 Y+ r" u: G
X1-Y<=0
2 D' W8 }. Y$ p* P2 E: D3X1-Y<=2 % _3 R$ j2 ^# g; d( T
4.5X1-Y<=5
% Z8 H2 ^8 v3 w4 HX1>=0 X2>=0 + m( a: E, d \+ R( k
至此没有问题, 但解下来该怎样解就不知道了.
r1 `- a6 o( d: r" Z N
* d$ ?. j( N# p. |' e【方法二】:
n! }) @$ x9 K2 K8 v4 D4 S取近似值的方法不一样, 0 e, k; t; v& c* w' S7 i" M+ F
X1=0p1+1P2+2P3+2.5p4
& Q% O0 ]* p- M, X4 Z, q y3 b" mY =0p1+1*p2+4*p3+6.25*P4
$ ]) S6 s/ D- ~' M2 X# ?1 = p1+p2+p3+p4
; }$ Y0 l; n" ^9 m0 t: V4 q/ _# U3 w; B" W3 l
原规划可表示如下:
$ Q! T8 D2 m0 s( ], D# ~Min Y-4X1-2X2
- U( N7 v9 u+ Z. p G& is.t. 7 m( j: ]! v4 q' C& {( d
X1+X2<=4 , _4 }$ b4 }% |) h! T& `7 b* Y
2X1+X2<=5
( P7 B, q; X$ f: U Z' t; Q-X1+4X2>=2
+ [% |8 O* l, t: R-X1+P2+2P3+2.5P4=0 6 |6 \7 _6 M1 M8 \
-Y+P2+4P3+6.25P4=0 * M" ]# ^2 X" Z1 W$ O$ k
P1+P2+P3+P4=1
& Y: g- \, T. H( s: \4 ZY,X1,X2,P1,P2,P3,P4>=0
6 L$ V& M3 U3 h, V同样的问题,到这里就不知道下面怎么解了. , z- e& [) M0 c
图2
9 A' h( Q/ N9 f4 t , h4 f& A8 F' _* Z# L9 U0 e9 X
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!
' X5 r; N/ }! t2 W6 `/ @ b+ x+ L$ \# E% Q% K
$ B- n* f" x& R9 D2 r我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
/ v6 b8 t X8 z: } |    |
zan
|