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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
main(-xn)
" V# `. l% S1 i" D& ]4 O2 D" Nst
4x1-4r1=1;
7 c6 z; D+ \: a. A) |x1+s1=1;
' A. M) |, ~1 c  J- B4 T4xj-x(j-1)-4rj=0;(j=2,3,4......)
4xj+x(j-1)+4sj=0;(j=2,3,4......)
xj,rj,sj>=0;(j=2,3,4......)
; L6 ?3 s7 G, w, a; L
) g0 \/ p& b6 N* d0 q请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
* `0 y3 u: M: r+ V2 j" x2 k其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
0 A- @% @4 g) M0 Y5 w$ ?n=input('input n please.(press enter)n=')
, H: i' X& D" J; Tif(n==1)
    a=[4,-4,0;1,0,1];
# d/ [" D" K- h; U+ [5 U    b=[1;1];
    v=zeros(1,3);
    c=[-1;0;0];
( }7 u* ]- a9 L  N! F; p    [x,f]=linprog(c,[],[],a,b,v)
7 t  E7 ]* Y9 A! b% \% }else
# `# \6 _+ n+ i  k- L+ E4 X    a=zeros(2*n,3*n);
7 Q4 S1 t4 D6 Z5 A5 G* e    for j=2:n
" V+ h% Y0 L4 d5 P; G4 G                   x=zeros(1,n);
$ x; W) ]0 b8 S                   r=zeros(1,n);
                   s=zeros(1,n);
9 b' T' A, N7 M                   x(j)=4;
( A/ p  Z& [5 _8 t. P  D                   x(j-1)=-1;
8 g4 P) K$ j4 t/ E! X                   r(j)=-4;
$ N' C2 n# ?$ K* r, R5 u                   a((j+1),=[x,r,s];
& |1 d) u0 H- a5 K! H" s) c; F3 ^                   x(j)=4;
                   x(j-1)=1;
                   r(j)=0;
                   s(j)=4;
' w$ {5 n  W: ]. O" k( P                   a((j+n),=[x,r,s];
    end
3 j3 m- k" G/ G4 T+ o               v=zeros(1,3*n);
              c1=-ones(n,1);
' G, i) \5 A: ?; q              c2=zeros(2*n,1);
$ b6 |1 y! O: k5 v: L& D$ V7 @              c=[c1;c2];
' p& N0 W- D0 d& n' r. K6 F               x=zeros(1,n);
! w% f; Y$ i/ h" {2 }, e/ b               r=zeros(1,n);
" t  O( g% M. V               s=zeros(1,n);
3 G: J5 S0 P/ W               x(1)=4;
% Y7 _5 f9 ^  c               r(1)=-4;
( U1 y3 a7 u9 m' \: K% Y& a+ Q               a(1,=[x,r,s];
! \! W! v" q" N( s               x(1)=1;
( U( i( x0 d2 p/ X9 e8 t6 e               r(1)=0;
               s(1)=1;
               a(2,=[x,r,s];
, P3 V! z0 k. N$ q  F* x8 v+ g               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
  J+ r2 V! b' n& p. D" j* Z               [x,f]=linprog(c,[],[],a,b,v)
1 b( f: I& w0 ]- send

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss