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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
main(-xn)
st
4 N) T4 i7 l4 j4x1-4r1=1;
x1+s1=1;
9 r# h& }/ k4 i* e' ^- G4xj-x(j-1)-4rj=0;(j=2,3,4......)
4xj+x(j-1)+4sj=0;(j=2,3,4......)
4 n0 P8 |( A' O  c; Qxj,rj,sj>=0;(j=2,3,4......)
. c0 i9 ]9 t! f9 f& O7 V" w
& G3 h/ v' k) p/ H) z4 v2 t请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
/ j9 W/ o1 ]  l0 V# H! Bn=input('input n please.(press enter)n=')
- n) g- r' h1 f& E! g( Rif(n==1)
% z( W. V8 F1 i! c) o& j    a=[4,-4,0;1,0,1];
4 b' w# ?& q  ?7 Q# E+ ^  t# j. a    b=[1;1];
9 c1 \2 w, e9 Z3 b- M    v=zeros(1,3);
    c=[-1;0;0];
5 G' g8 T3 M3 d& h) ^  B    [x,f]=linprog(c,[],[],a,b,v)
. {9 _9 z7 {5 jelse
7 Y  A" f- l8 T, G- h    a=zeros(2*n,3*n);
    for j=2:n
                   x=zeros(1,n);
: T: J$ r( |3 E; A. O                   r=zeros(1,n);
0 A5 m/ |: `: Z' F3 g                   s=zeros(1,n);
                   x(j)=4;
8 w/ g. x) f0 B* u% K2 K                   x(j-1)=-1;
                   r(j)=-4;
& ?( a5 w- f/ g' ?$ F+ t                   a((j+1),=[x,r,s];
) ^( I* ^# K3 H& U, p2 y. y7 }                   x(j)=4;
                   x(j-1)=1;
* H5 X, j0 O3 r6 V* j' a2 l8 N+ |4 v                   r(j)=0;
                   s(j)=4;
                   a((j+n),=[x,r,s];
) V9 u3 S9 K+ F0 Z, E5 j  b    end
3 I, n- c1 T8 g6 f# [               v=zeros(1,3*n);
) v, H5 I8 `0 b              c1=-ones(n,1);
% I+ K: l, B2 f0 z* _0 [              c2=zeros(2*n,1);
              c=[c1;c2];
' y: o$ {% u, h. z! X! ?               x=zeros(1,n);
               r=zeros(1,n);
               s=zeros(1,n);
               x(1)=4;
               r(1)=-4;
, {! X7 r5 |3 O1 ?: G9 q4 \               a(1,=[x,r,s];
$ S$ D8 _% J, ^# \! E- E               x(1)=1;
3 z+ n: `- r/ D8 V9 C               r(1)=0;
  L# O9 s0 c# q               s(1)=1;
               a(2,=[x,r,s];
               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
               [x,f]=linprog(c,[],[],a,b,v)
2 Z- a* c. L) Q' M) M1 Kend

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss