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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
' Y, q( c. z0 tmain(-xn)
3 K( Q  f8 T0 X8 Q* N, a/ j/ Zst
4x1-4r1=1;
x1+s1=1;
4xj-x(j-1)-4rj=0;(j=2,3,4......)
, k  B5 u4 g8 {$ c) n% l4xj+x(j-1)+4sj=0;(j=2,3,4......)
0 t/ f; X% r$ p5 vxj,rj,sj>=0;(j=2,3,4......)
& V# D' U% v3 [( a8 ^/ ?  k7 {
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
* s" O. D& L' Wn=input('input n please.(press enter)n=')
1 x4 e' O6 }0 C, h; ?+ Aif(n==1)
7 R7 [* g. r0 L    a=[4,-4,0;1,0,1];
# S- q4 h8 ~& n2 _    b=[1;1];
8 |" {8 @: Q6 ?1 |    v=zeros(1,3);
    c=[-1;0;0];
    [x,f]=linprog(c,[],[],a,b,v)
: j, |3 M, F& V' D8 z: `4 yelse
, ?5 v6 u9 V3 ?) I    a=zeros(2*n,3*n);
( [" h; k8 K& a5 n2 K- K2 D    for j=2:n
4 A) O/ \) L2 g; a& D                   x=zeros(1,n);
                   r=zeros(1,n);
& _3 e" Q* X0 r3 l0 i- C' v                   s=zeros(1,n);
                   x(j)=4;
                   x(j-1)=-1;
8 j- g, f) k2 S/ Y# j& a                   r(j)=-4;
                   a((j+1),=[x,r,s];
4 Z" K' C, t, x* R3 @% w) e' X- M                   x(j)=4;
4 G) `  u% }3 \: S+ [$ K, n                   x(j-1)=1;
                   r(j)=0;
                   s(j)=4;
                   a((j+n),=[x,r,s];
% |) e% f. V" w3 Q% _    end
( i% f4 a1 ]* U. p8 G3 h               v=zeros(1,3*n);
, }" ~  v3 s+ ~- P              c1=-ones(n,1);
              c2=zeros(2*n,1);
              c=[c1;c2];
/ o9 N5 B5 b+ Y$ P* }: S0 W               x=zeros(1,n);
               r=zeros(1,n);
               s=zeros(1,n);
  D. S  C# y6 P6 u7 G6 f% o% A               x(1)=4;
( T% ?9 q5 {, ^! t( {               r(1)=-4;
4 v/ ?5 n1 N) h! V, l# F, g& G               a(1,=[x,r,s];
2 M9 |1 {+ N. R0 F% N               x(1)=1;
7 j4 x' _# a* f1 F- q% }; ]( `               r(1)=0;
9 E1 B4 O! U' @               s(1)=1;
5 M3 a" j6 c6 l1 q5 m; y6 L+ C               a(2,=[x,r,s];
) Y5 f0 P! T* Z+ z4 H               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
6 m, W* x4 }: ^  l               [x,f]=linprog(c,[],[],a,b,v)
end

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss