数学建模社区-数学中国

标题: 我是新手，MATLAB接触不久。有一题不会，求教求教 [打印本页]
作者: 玖玖    时间: 2013-7-15 16:44
标题: 我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）$ l5 I5 a# m: p3 B- V+ C
main(-xn)
2 U/ Q0 x7 u: V3 d* xst " b% m* s2 z9 V' F% @
4x1-4r1=1;
! [4 a9 J( g" F9 S" @* D7 {  ]x1+s1=1;
; o# w/ _; x. W& T2 h" }4xj-x(j-1)-4rj=0;(j=2,3,4......)4 W, {' O/ w4 y- g# X' x
4xj+x(j-1)+4sj=0;(j=2,3,4......)
. v# ~% N# G& p6 uxj,rj,sj>=0;(j=2,3,4......)
$ j( H# z- ], S2 }% \+ e% `1 e: z3 d6 n( p, P- a- G
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。8 B; Q( M& A6 \: \8 z& u
其中j，n是下标 。
作者: 玖玖    时间: 2013-7-15 20:29
有大神在不？帮忙啊！！！
作者: 玖玖    时间: 2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者: 玖玖    时间: 2013-7-16 14:57
    syms n
  |8 T" f' a9 [, q5 L& U3 F& w4 wn=input('input n please.(press enter)n=')& u& ~) H# c0 ?% ^* l- H+ Z
if(n==1)/ ]$ y5 \5 e# k3 Q7 g% j
    a=[4,-4,0;1,0,1];& r/ X, b! V3 P4 e
    b=[1;1];
$ }' [) c/ D" _4 b% C6 I3 G& }    v=zeros(1,3);: J# ?# k/ `" I0 H6 w
    c=[-1;0;0];( o. |4 r% l. U/ b; a
    [x,f]=linprog(c,[],[],a,b,v)
0 p! p' h3 R- m  `else
, c1 q* h4 l5 P# w+ M    a=zeros(2*n,3*n);& w4 n9 c* y' t* p+ C# \
    for j=2:n
3 d6 J. ^7 r; M- {8 ~& `                   x=zeros(1,n);
+ u& \+ T2 w+ P" ^" G                   r=zeros(1,n);
5 Y# @3 ^3 Z/ y! r! T0 ]- u) h                   s=zeros(1,n);
) O( Y& b- l, M                   x(j)=4;
! {" q" [) z" ?- H/ N                   x(j-1)=-1;0 V+ t; K- D7 R% B
                   r(j)=-4;
% c$ j+ g1 J: a4 `# b: Y. P                   a((j+1),=[x,r,s];
, N% o5 |- @# x1 Q) @                   x(j)=4;7 e5 n5 V3 @( n6 C
                   x(j-1)=1;+ e! J+ |& O6 {
                   r(j)=0;
3 L, k+ t) e4 K7 C* f0 u                   s(j)=4;
+ {) K( z/ Z7 G  P; L; @                   a((j+n),=[x,r,s];/ h& Q" E$ ~1 G" |& `, Q
    end* A8 H0 t: O# \( x+ n. {. V
               v=zeros(1,3*n);
4 u$ i7 u- Q' p              c1=-ones(n,1);
" r2 ]! f1 Z2 a+ `% g' x2 e              c2=zeros(2*n,1);
/ a! U0 f( F7 X9 {2 R+ P              c=[c1;c2];
9 N" C5 N, Y' ?/ Z! E" @/ h               x=zeros(1,n);0 p2 U3 T# R$ x$ ~
               r=zeros(1,n);" u# J3 r6 n" b3 c
               s=zeros(1,n);
6 l; G7 s; I, Y               x(1)=4;# c2 p3 h, U9 s9 t4 t7 L  ?
               r(1)=-4;
3 z6 o  E) N2 k* ?               a(1,=[x,r,s];( X  \  n! @0 G! y% r; b  p
               x(1)=1;! R" \. b, d! {. c7 p7 J
               r(1)=0;1 ~5 W+ j1 _) d5 B
               s(1)=1;
" [! k) @0 P. ~4 w; l' e               a(2,=[x,r,s];/ g5 m5 c, l) \) C2 q7 H% X8 `
               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];" a% g' }0 G* z8 ~7 S
               [x,f]=linprog(c,[],[],a,b,v)
5 F% Q. ]8 v8 N& Oend
作者: 江边一粟    时间: 2013-8-15 11:37
楼主，笑脸是什么？求解。。。昨晚就在做这题，到现在还没做出来
作者: qianchuan11    时间: 2014-8-11 20:30
/wx ssssssssssssssssssss



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