- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
4 H- J x$ t: D# r4 W) p- S/ ~ + K! x& }. K% ]2 q4 w, I
+ h$ L8 O3 ~: q$ K8 b7 G4 L8 I" |: Y# A7 ^# t
| 再求数学高人 | 7 E, v G, E" R0 a
4 d* f( c+ S6 Z2 h | ) y: C$ i# U: [- N" b) I& U
: T3 F& e& b% ^& k' ?) Y5 t& G* A6 s" @! O2 F
$ Z5 p; D [% R* Z! Z! p9 Y; h用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
+ n8 \) b+ N! |5 o+ t6 _Min x12-4x1-2x2
7 J; b6 |8 M. Cs.t.不等式组条件如下 . w% `0 u! n* E [, N
X1+X2<=4
1 G% Y5 P# S$ u+ d2 z1 h- K U2 H V2X1+X2<=5 : P# J. H' g! b# N" G
-X1+4X2>=2
) H) U- r* S; s8 s8 IX1>=0 X2>=0 , ~0 x4 N3 v% M
【方法一】: 2 o/ g9 r- j0 k# R8 m, t; I
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 / z; c- C& t1 A V; I; `4 u
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 1 O/ E5 x# f, Z
点O x1=0 y= x12=0 ! {9 x6 ]1 [' M1 s6 P" F7 p# F
点A x1=1 y= x12=1
1 D# Y) C' U/ c( x" m) V. e5 ~点B x1=2 y= x12=4
, u" o0 E: `0 Q点C x1=2.5 y= x12=6.25
* K+ J H: m! i1 N8 J6 K如图1: ! y0 G( \1 ~# I/ {( |$ h0 p
2 j5 q4 m+ D- ]9 G+ h" m, ?% I v2 u- z1 C7 O9 X
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: * L2 z+ j$ E5 Z$ K9 B8 ~7 A: L
Min Y-4X1-2X2
+ ]" ^( ^8 P# K4 Fs.t.
6 c0 s, n2 X5 F* A3 [X1+X2<=4
4 J2 d* u* j8 t' |2X1+X2<=5
, h% F0 z+ E0 J- P1 J/ l-X1+4X2>=2
+ S( H! S, T( |9 tX1-Y<=0 3 u2 H6 N; k6 M2 V
3X1-Y<=2 ' M2 u, Y. I7 d) G# A; u3 Y
4.5X1-Y<=5
& Y0 q- U7 ]0 aX1>=0 X2>=0 6 v- o% v5 ^# h1 o8 \
至此没有问题, 但解下来该怎样解就不知道了.
' k2 b! e( \; b& y6 R# z & u9 D: o4 {* p4 Q9 ^
【方法二】: 0 j5 Z" e* a( e' r" j
取近似值的方法不一样,
+ S% H7 w: ~, GX1=0p1+1P2+2P3+2.5p4 $ E' X9 Z9 ?* V+ S) E6 H# x9 U
Y =0p1+1*p2+4*p3+6.25*P4 5 }6 ^0 g0 u8 _" f6 I* A& ^
1 = p1+p2+p3+p4 1 H) k( O/ x5 h
& x1 b6 u+ K% R$ x& Y7 D( h, |
原规划可表示如下: # L: N! Y5 m. d- G' Y7 I% L
Min Y-4X1-2X2
2 H% o7 ?/ x; S: G9 js.t.
+ ^0 J& `5 i( e# MX1+X2<=4
- G4 i% z9 Z/ k0 B$ v2X1+X2<=5 ( z/ W8 w/ p' e; Q- G( |
-X1+4X2>=2 & [6 _8 ~' k1 M p. T$ V
-X1+P2+2P3+2.5P4=0 4 U" \' k- m) `7 P$ U
-Y+P2+4P3+6.25P4=0 # [8 k3 v0 S6 K/ R2 S7 t
P1+P2+P3+P4=1
) @+ Z4 T# [% M2 q" p4 B# OY,X1,X2,P1,P2,P3,P4>=0 : T* r4 ^; n4 C3 ]9 |4 Y. ]8 |0 r) x6 ]
同样的问题,到这里就不知道下面怎么解了.
0 e {! O2 |- }% {8 t图2
6 O. j) J \$ g4 s& E+ \
0 z/ g5 H7 B0 k, m3 P8 ?; ]原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! * u6 x1 n. m. p+ g. v: @
, I R% A( F6 q' I, g( n
" g/ m# ?6 \" u& l: W6 y我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法 1 M1 p! E2 g$ f, ^: n9 I
|    |
zan
|