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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
main(-xn)
7 c( o; R7 t9 W; vst
; Q, f/ D) x/ @9 v5 D% f& E% E5 q4x1-4r1=1;
x1+s1=1;
4xj-x(j-1)-4rj=0;(j=2,3,4......)
4xj+x(j-1)+4sj=0;(j=2,3,4......)
# x0 X5 R% @2 ^6 J( D- B" Y3 e6 wxj,rj,sj>=0;(j=2,3,4......)
& l0 _0 M5 m! j3 R$ n7 X2 b
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
, c+ c1 R$ r: L. ?+ ~% p2 {$ d其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
n=input('input n please.(press enter)n=')
8 `5 V# x/ o  d- b4 A) uif(n==1)
    a=[4,-4,0;1,0,1];
% j) B' k) p  [, u( B    b=[1;1];
    v=zeros(1,3);
    c=[-1;0;0];
5 T+ ^. T( z' D% G& ~: M    [x,f]=linprog(c,[],[],a,b,v)
else
    a=zeros(2*n,3*n);
    for j=2:n
                   x=zeros(1,n);
8 _7 r. d0 \2 Q* v+ s                   r=zeros(1,n);
4 D- J; @; k* c* n3 S( c  G! x                   s=zeros(1,n);
                   x(j)=4;
                   x(j-1)=-1;
1 x9 W! b. h9 F% }9 c, U( ?                   r(j)=-4;
( N1 a- Y- S- J$ q$ f+ E. t7 B) W                   a((j+1),=[x,r,s];
) t. E) C* h# [                   x(j)=4;
. g7 \! D) H3 q" _) t                   x(j-1)=1;
7 L4 W! S& d0 v% k                   r(j)=0;
) \9 w5 J' q5 U0 [( {; C8 J                   s(j)=4;
: M# @8 G0 t5 {# P+ b: j9 `                   a((j+n),=[x,r,s];
    end
               v=zeros(1,3*n);
              c1=-ones(n,1);
              c2=zeros(2*n,1);
              c=[c1;c2];
8 ~! {/ N! J' c/ m: K               x=zeros(1,n);
               r=zeros(1,n);
               s=zeros(1,n);
3 ^6 `: P" n) Y4 ]. {# \; R               x(1)=4;
, \( L, K6 {6 }7 ?               r(1)=-4;
0 q0 B2 p/ z. {5 `+ \; S) {0 n% k               a(1,=[x,r,s];
6 T3 Z/ D+ N3 f# E4 ?3 _               x(1)=1;
4 f* L* Y+ i, O6 U3 [" G$ Y/ ~1 t3 d               r(1)=0;
               s(1)=1;
* W! ]% y1 p" e6 f, \: [+ b0 S. ?               a(2,=[x,r,s];
& M. S. G5 A9 L; J* W* s  w               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
               [x,f]=linprog(c,[],[],a,b,v)
3 a- ?0 O  D8 @9 `( ^; m. X; Gend

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss