数学建模社区-数学中国

标题: [求助]一道优化的题目，除了库恩塔克的解法，有些不明白，欢迎指教 [打印本页]
作者: 徐帆    时间: 2003-1-15 11:03
标题: [求助]一道优化的题目，除了库恩塔克的解法，有些不明白，欢迎指教

7 V7 d! m3 L' Y& d- A: D: y" b) x1 W @6 ?4 O' `! E7 W- p/ O9 { x" F+ } N3 ~; t0 m* R2 j2 B4 p/ |; R$ S$ v. x& t1 j: y! x5 h# d2 |9 x! Y% ]3 U, S. h8 M4 l! `5 J
再求数学高人
% O* o: s5 \2 s/ U/ \: g
; h- o- ?+ Z4 g8 P4 Y8 `# l" ?

用分段线性函数近似求解非线性规划(不知道标准数学表达法是怎样的)，

% _8 q! S5 G& b$ e

Min x12-4x1-2x2

3 R. ?. q! W3 f7 x+ @$ b

s.t.不等式组条件如下

7 E/ Z" j3 c1 e

X1+X2<=4

1 D! Z, m: m, Y5 {

2X1+X2<=5

/ ^, P+ f/ X5 t

-X1+4X2>=2

- l- T3 ^# a( H5 Y. f' a

X1>=0 X2>=0

: j) M" F( ~8 ~4 `

【方法一】:

3 |' @+ H: \" c% q4 i

因目标函数中变量是可分离的，所以用替代变量y代替二次方的变量x12，将原目标表示成一次函数，以便于用线性（单纯形）方法求解。

, |! ^5 f2 G* D9 q! w

首先给定一个取值的区间，设0≤x1≤2.5, 选择一系列x1的值，比如0，1，2 和2.5, 计算对应的函数值y如下:

$ c5 r7 n; V% m

点O x1=0 y= x12=0

% \9 p" a: k1 {0 w+ K

点A x1=1 y= x12=1

$ V! D2 v- L4 k/ w! ^, `

点B x1=2 y= x12=4

, T" c9 P2 V; o! ?7 R

点C x1=2.5 y= x12=6.25

* A$ `9 O6 q4 ~; _9 l. A

如图1:

$ f- ^1 E! M0 [1 w+ z

! _0 ]/ k0 v: f6 y2 C6 t

0 k. t$ J# z2 }

用分段的线段OA, AB, BC近似代替y=x2, 原规划表达为线性规划如下:

4 } O. b* T+ L, Z

Min Y-4X1-2X2

, M7 j, X/ w6 b

s.t.

5 Y" ~: R$ w- b: A0 Y

X1+X2<=4

7 o( G. N$ V0 t( v5 S0 N

2X1+X2<=5

( }( J0 q/ `( }

-X1+4X2>=2

/ k3 O* i9 k8 Q9 U6 n

X1-Y<=0

% m8 ^, y' }' p3 C- @

3X1-Y<=2

! D& H1 x/ Y1 v9 T4 X( C0 F2 U

4.5X1-Y<=5

3 P( R. r, m: z5 i# d

X1>=0 X2>=0

) e9 t0 M* e3 d. m7 f

至此没有问题, 但解下来该怎样解就不知道了.

( k+ I+ _) U9 ?1 ^

! {# f+ {! n0 h2 ]

【方法二】:

7 k. G2 U1 S5 a6 X* w- a# J4 F

取近似值的方法不一样，

' z$ p; `6 r% E: x, n. u

X1=0p1+1P2+2P3+2.5p4

% V5 I& O7 [ K' R

Y =0p1+1*p2+4*p3+6.25*P4

# @. }5 X/ f( E( {

1 = p1+p2+p3+p4

1 I# C$ Y4 U: J( c

( y. {9 O. |! L6 V! w8 r

原规划可表示如下:

9 e. ?+ {. P" D" V) _" ~1 ~+ {

Min Y-4X1-2X2

* t# H% R/ u6 e- l6 @1 o% P

s.t.

. i: ~& k' \: r

X1+X2<=4

) @: H6 Z2 @% _) ]. b

2X1+X2<=5

! S, k" S+ d4 _# _2 @

-X1+4X2>=2

3 ^9 `4 w( D) _+ n" l* d

-X1+P2+2P3+2.5P4=0

5 h: X# d/ l1 X4 X: |

-Y+P2+4P3+6.25P4=0

4 G0 G: F+ {- c6 X- l

P1+P2+P3+P4=1

# U9 ]" p2 I, _4 P3 [% \

Y,X1,X2,P1,P2,P3,P4>=0

8 U, I: K1 x6 [# a1 J

同样的问题,到这里就不知道下面怎么解了.

" i, J' q! {1 e* U) b: h6 _

图2

6 J- \8 g8 r* T- Y4 i

" F( W. s6 k6 O5 l2 G2 v! t+ L

原规划用Kuhn-Tucker方法可解, 但老师非要求用上面方法解, 所以请不吝赐教, 谢谢!

- v3 Z _" Y, W, {& | j0 T3 Q$ d3 E1 b. i% \1 [" P! f 6 [0 O2 N# b- e( E

我用kuhn tucker解得答案为 min=-9， x1=1， x2=3,验算觉得答案是对的，可还是不是老师要求的方法

( W& V8 U! k/ v. j \

作者: madio    时间: 2003-1-16 14:55
你不断地加密分点不就可以使得近似成的线性规划的可行域与原先非线性规划的可行域不断接近吗？细分到一定程度就可以认为差不多了呀！这样相当于增加了约束条件的个数来换取目标函数的简化！
作者: 徐帆    时间: 2005-12-29 18:04
这些我都明白，就是到了最后，不知道该怎么解这个不等式组了
作者: madio    时间: 2005-12-29 22:37

晕，大型线性不等式的求解问题至今好像没有什么好的办法呀！

作者: 徐帆    时间: 2005-12-29 22:45

不用大型的，就一般的，比如这一题，算三个点就好了，

" _8 t3 t. K4 K! N( c2 ^! u9 w. t

如下了，该怎么解这个不等式，我在国内是英文专业，现在上物流master，碰到这种玩意实在不知道从哪里下手了，谢谢拉

3 X$ _/ F, c. v9 \. `6 R3 [

Min Y-4X1-2X2

* x8 f- F, |+ _, P$ O8 j/ V3 ^

s.t.

# z7 H; {6 G" t: t k

X1+X2<=4

( j; ` k& m: c$ q# b" ?( h- F7 q

2X1+X2<=5

. ^5 K" l, K* H7 Z0 g, h. _

-X1+4X2>=2

1 G0 I. }6 W2 S8 Q& n9 K- S

X1-Y<=0

5 O# }1 }3 n& s6 {4 @

3X1-Y<=2

) _2 r- ?7 H# g0 s/ d

4.5X1-Y<=5

; N$ J9 q0 G! I8 D

X1>=0 X2>=0

[em06][em06]
作者: madio    时间: 2005-12-29 22:55
好在是三维的向量不等式，你把每个不等式看成三维空间中的半平面，将这些半空间画出图来然后做交集就可以了！
作者: chenyugf    时间: 2006-3-13 15:58
唉,都是高手,就我不是
作者: yjshappy    时间: 2006-4-23 02:23

没看明白哦。。。。。。

作者: youshanglanyu    时间: 2007-1-14 23:34
还没看懂



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5