我是新手，MATLAB接触不久。有一题不会，求教求教

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
6 m% \) i  f7 k0 v4 Rmain(-xn)
st
. h+ D3 ?* T, l" p4x1-4r1=1;
) ^- c4 _" D- X" ]0 T6 Zx1+s1=1;
7 g: v$ |5 {, q& K" S4xj-x(j-1)-4rj=0;(j=2,3,4......)
4 @4 M3 v0 }  ?) U4xj+x(j-1)+4sj=0;(j=2,3,4......)
. O+ Y& P- c, ~xj,rj,sj>=0;(j=2,3,4......)
3 C' Q2 A& ?) v- f" a3 o, ^; {: e5 S
' U! x" n! ~0 ?5 X6 m请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
其中j，n是下标 。

玖玖

有大神在不？帮忙啊！！！

玖玖

好吧，经过不断地努力我自己做出来了。

玖玖

    syms n
2 b) H0 w  p2 d& un=input('input n please.(press enter)n=')
if(n==1)
; c! Y: H$ Y+ ~    a=[4,-4,0;1,0,1];
    b=[1;1];
' ~6 F. @5 N6 K7 y& J. C! w    v=zeros(1,3);
    c=[-1;0;0];
    [x,f]=linprog(c,[],[],a,b,v)
' [2 l  V- D  D, f$ ^, Belse
8 X4 k  c1 H# F/ U' s( u, u    a=zeros(2*n,3*n);
2 b; m* |( z1 s: G6 ?! q    for j=2:n
                   x=zeros(1,n);
9 n1 C- D$ r% G) c! ?                   r=zeros(1,n);
                   s=zeros(1,n);
                   x(j)=4;
+ F/ u* F& K) w, ~                   x(j-1)=-1;
                   r(j)=-4;
                   a((j+1),=[x,r,s];
; o  S6 T2 M4 H; x% N                   x(j)=4;
                   x(j-1)=1;
                   r(j)=0;
7 w: q8 q5 c7 h1 z/ ^                   s(j)=4;
/ [* b' V7 C: I3 |1 P                   a((j+n),=[x,r,s];
3 x: b! `3 R5 Z$ W    end
% O5 p% u, K" v& W3 F( D9 X               v=zeros(1,3*n);
# Z; x2 Z7 s& ^, z2 y              c1=-ones(n,1);
: p2 y% ?& o6 B6 }              c2=zeros(2*n,1);
              c=[c1;c2];
3 j1 T# c: ]( O, g  B+ O               x=zeros(1,n);
               r=zeros(1,n);
9 i3 i6 P* B3 `) d% U2 m               s=zeros(1,n);
$ d( S) `6 E+ r* X. |- @               x(1)=4;
; y! Q/ ?$ }* P) s               r(1)=-4;
               a(1,=[x,r,s];
               x(1)=1;
( G7 e  N% ~$ r4 E               r(1)=0;
) O( W; v9 q  a4 u0 _               s(1)=1;
% ?5 B/ G! y; d0 q$ @6 `8 }( @               a(2,=[x,r,s];
8 b4 C6 U0 Q0 f  v4 V" g4 q               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
               [x,f]=linprog(c,[],[],a,b,v)
end

江边一粟

楼主，笑脸是什么？求解。。。昨晚就在做这题，到现在还没做出来

qianchuan11

/wx ssssssssssssssssssss