- 在线时间
- 0 小时
- 最后登录
- 2005-12-31
- 注册时间
- 2003-1-15
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 60 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
 |
|
' \' y3 c2 H" P+ I6 p# C
3 {# [; B1 c" _3 L: o/ N( k. C( U7 e; j3 ?% |* T/ i- f+ S
& p2 Y3 k. Y8 D! Z2 P& V) }
| 再求数学高人 |
/ } H( i( x4 \3 b% `
" w/ z5 J! G+ b1 [5 y/ U2 g9 y | 9 X- v8 d! U5 F
6 m9 [; }9 Z0 {: P$ Z; W* _% b) e
: ~9 B1 [' e! A' P* w: H" @; V1 M% h# l2 T; Z
用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的), " G4 N1 g* v! T9 x+ p( P0 Y
Min x12-4x1-2x2 # v2 F7 b6 @2 |1 m% g4 h
s.t.不等式组条件如下
; H9 R! Q0 _6 e! `' GX1+X2<=4 + H4 Q% y9 e U( i( v1 _
2X1+X2<=5 # U! B) t6 ^7 c
-X1+4X2>=2 0 F: f. P9 f8 e' _) Q
X1>=0 X2>=0
3 }% X- u* C6 m4 i【方法一】: ( Z, H* ^, i7 d; E9 c
因目标函数中变量是可分离的,所以用替代变量y代替二次方的变量x12,将原目标表示成一次函数,以便于用线性(单纯形)方法求解。 8 ^# B _5 S1 Q
首先给定一个取值的区间,设0≤x1≤2.5, 选择一系列x1的值,比如0,1,2 和2.5, 计算对应的函数值y如下:
( C% _1 G' r( B. u/ T0 H点O x1=0 y= x12=0 * d. @( W# J# ^" }8 o
点A x1=1 y= x12=1
1 H9 s3 A1 @) e' z! @- j点B x1=2 y= x12=4 7 e0 ~! ~( Z$ ^1 U* g
点C x1=2.5 y= x12=6.25 5 m& C9 q, E) V3 m" C8 v
如图1:
& h4 c* ], Q! t! c2 R! T J4 c9 B. e5 G, T, p/ e0 k% N0 B
( r! l* A$ ]* c( S/ X: s用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下: % |) @; m. G4 h; B( k* t
Min Y-4X1-2X2 1 V/ h. o! F. h0 r9 i+ x
s.t. # q" n- V! u0 k1 B! F6 [# {
X1+X2<=4
- W* Y( Z! d. P4 Q2 a2X1+X2<=5 4 A2 A1 r: e. v4 u/ @. v
-X1+4X2>=2
. l- j. [9 A+ t; @$ e9 u' j! IX1-Y<=0
+ |2 \* e9 J. W6 r3X1-Y<=2
8 \7 _1 ~$ Y1 n: a4.5X1-Y<=5 6 H7 ?; F$ ?4 T; P, x
X1>=0 X2>=0
7 X0 p. n& M) D% r1 L/ D" L至此没有问题, 但解下来该怎样解就不知道了. 4 L2 W8 r5 y( Z
8 A& E, q" R! J【方法二】: 2 b4 q" Z$ T3 g! b8 n3 X
取近似值的方法不一样, * u1 n! J2 J/ q! l% C9 E0 s
X1=0p1+1P2+2P3+2.5p4 ( u) R# q2 M, q
Y =0p1+1*p2+4*p3+6.25*P4
0 B& [& B1 w8 i$ X; M& e1 = p1+p2+p3+p4 : N. ^' E" r* o
" T; \6 p @" p6 P0 m, A, u3 t原规划可表示如下: ; G0 H* y, N# [
Min Y-4X1-2X2
4 N# E" q, t. \: A2 O* M7 c0 Cs.t.
4 O; Z5 J8 A' l: C/ xX1+X2<=4 : ?. K+ i2 T+ M- I5 e# R7 O
2X1+X2<=5 2 c/ K0 A5 D/ `' }: G8 g
-X1+4X2>=2
( o+ C* X8 n' d7 I-X1+P2+2P3+2.5P4=0 " z o. g# d) B2 v9 L2 @2 {0 M
-Y+P2+4P3+6.25P4=0
2 i7 l9 P0 o9 W# p$ M# yP1+P2+P3+P4=1 / S& @& X! S+ m7 m
Y,X1,X2,P1,P2,P3,P4>=0
6 `- V; t; }6 m( x; h( M: ^同样的问题,到这里就不知道下面怎么解了. ( l. I4 q, R" z
图2 ! B/ _9 l! m2 c6 D4 r, h
/ g: E8 w) e, d2 h' \, G' b
原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢! / U8 h' w" \% I3 ~
0 I4 o B; A0 w- t* G8 a7 j0 R$ ^0 i; v1 I# p# b
我用kuhn tucker解得答案为 min=-9, x1=1, x2=3,验算觉得答案是对的,可还是不是老师要求的方法
% y) n+ b K! u( G9 i |    |
zan
|