- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
5 b; a5 v; N: E) W; j - W$ I9 }" k3 ? a: \1 k. r
! k+ u$ w0 x- [" u6 j {# z5 I- w2 [# E4 A
| 再求数学高人 | . k5 H% v- X8 U2 b
4 E! Z% V; R& \ P9 P |
5 u- j' S$ K7 O$ K% b. N" I7 p- c5 x. ?
6 {! p* D" T' F) N* z
: T$ D+ ~& R, Y9 H4 R# e用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的),
3 Y p4 ]* F6 n& s1 nMin x12-4x1-2x2 ( ]9 \: H- a7 j4 B
s.t.不等式组条件如下
+ W U' W, ]- P* R! ZX1+X2<=4 ; F( X; S, g6 U
2X1+X2<=5
/ Q; ^' ^6 M/ |0 z$ x/ P7 F-X1+4X2>=2 + Y& |) U- c/ i9 k8 z; Q- j
X1>=0 X2>=0
/ d. {. }1 J6 U" ~- O【方法一】: ) G( @7 u7 t2 r# e7 O
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 ; \& E. A6 Z' b% |% m
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下: 1 b$ C# s, n4 ^
点O x1=0 y= x12=0
' _! Y6 E2 {6 k% B3 \0 ~' S点A x1=1 y= x12=1
0 B; [( w1 `7 }, w% s6 l点B x1=2 y= x12=4
4 z& O2 }& r( C. N5 N" M( U点C x1=2.5 y= x12=6.25 7 F$ i# X! X7 W3 V) E/ u
如图1:
7 a" K9 K7 T# D: @ + I% o% v' g3 T8 t+ U! P
: i2 o. A) f1 Q$ M5 e( B% c
用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: + O7 G( V6 q% b' X) \& b3 [
Min Y-4X1-2X2 * {. u. d1 j9 p! K2 P. Q4 Y
s.t.
* j9 w; s4 O& f( @X1+X2<=4
6 _" H6 k/ U3 R4 {. s+ e: E3 I- P2X1+X2<=5
* B& L/ y% e0 g; ~- b: I-X1+4X2>=2 1 i& m5 l2 O6 }$ X9 h7 l, G
X1-Y<=0
# P& m( ^3 I: F& G. h" y3X1-Y<=2
0 U! t5 r* w7 ^; ~2 `4.5X1-Y<=5 ) D! C$ H H/ I& X, U" Y
X1>=0 X2>=0
& R6 k' S, q9 @" R+ K. u0 p至此没有问题, 但解下来该怎样解就不知道了. - Q; H& s- c; W
. G+ z4 v/ j O4 J }; Z2 i4 d【方法二】: % ^4 p, H8 i4 P
取近似值的方法不一样, 8 k7 b- G9 G4 m0 a" G( [
X1=0p1+1P2+2P3+2.5p4
6 H i5 I9 V6 M6 ~, [% q2 h% dY =0p1+1*p2+4*p3+6.25*P4 1 s: F2 t+ B7 H/ x3 q% S7 N
1 = p1+p2+p3+p4
8 \/ n& W+ ^2 o$ F; o$ g2 E1 f, m& \- f$ F# W9 S9 I3 P' k
原规划可表示如下:
1 Z0 O0 o5 ?3 OMin Y-4X1-2X2
* x: o& j$ b; j. Ys.t.
/ f C( K& e I. G1 w" fX1+X2<=4
/ {+ j3 v) z I( M p; X+ n9 J2X1+X2<=5
7 i, }1 n( k% {) v-X1+4X2>=2
& ?1 B) v: ~# }" `2 x( [-X1+P2+2P3+2.5P4=0 3 F, a: H( u: D+ z, g4 U; F. D
-Y+P2+4P3+6.25P4=0
$ O/ E) A3 k# Z5 y7 _, kP1+P2+P3+P4=1
, [ C X7 B }5 ]Y,X1,X2,P1,P2,P3,P4>=0
! M% G; C: l1 t( C( ?) |同样的问题,到这里就不知道下面怎么解了. ; X; c i m; V5 G# H2 }6 ~/ x' V
图2 ) v( M5 i( v% ^& f3 Y% Y y1 u
/ D4 b3 L* n- X9 Q原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! e X: i1 Y' p6 u
, y) a! ]* e9 p0 z! y
, E! v) Q8 ]) y4 n4 l v6 x我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
, \; {4 e0 A/ {% Y |    |
zan
|