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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
( S0 p7 w: ~" D; l. V1 ^# Zmain(-xn)
9 l1 A+ u' V1 P0 vst
- x  \& p* g( }- l2 ^3 {# {4x1-4r1=1;
- ?% j$ z2 Z) H$ I- }' X- Sx1+s1=1;
! K$ Y9 _0 Y: {: ]4 S) E. c- M8 {4xj-x(j-1)-4rj=0;(j=2,3,4......)
! r* P" A9 v! P4xj+x(j-1)+4sj=0;(j=2,3,4......)
xj,rj,sj>=0;(j=2,3,4......)
4 v8 k, n2 i' {
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
: V' I% `  B8 j3 }其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
n=input('input n please.(press enter)n=')
if(n==1)
% y- ~  l- R( n* W    a=[4,-4,0;1,0,1];
  b' b# P$ v4 X$ A+ K7 X& g8 O    b=[1;1];
    v=zeros(1,3);
    c=[-1;0;0];
& Q9 Z; i! O8 e    [x,f]=linprog(c,[],[],a,b,v)
  n( d% K- r' i: Z# [$ h8 o3 h  nelse
% Q  c4 l$ k) p# S) O+ v    a=zeros(2*n,3*n);
    for j=2:n
/ L5 ^  Y/ t5 f1 F5 k3 W3 Y4 O$ W                   x=zeros(1,n);
                   r=zeros(1,n);
                   s=zeros(1,n);
: |' R$ k7 d+ P8 s! X                   x(j)=4;
) A. K8 R5 W( D  w5 n/ g# o                   x(j-1)=-1;
4 \3 m# x: O" r8 s# N5 n6 w/ l                   r(j)=-4;
- q9 @6 \2 [) p% J" P6 Y2 ~# o                   a((j+1),=[x,r,s];
0 P; V' d& e6 \) q                   x(j)=4;
) B" Z% I, `& J" p5 j' J1 [                   x(j-1)=1;
+ O% ^/ Y6 D! {2 U# V* Z4 _                   r(j)=0;
8 i, @) B3 p3 z2 r                   s(j)=4;
                   a((j+n),=[x,r,s];
7 g! O0 b. @* E, I    end
               v=zeros(1,3*n);
              c1=-ones(n,1);
! p3 F9 x. ^2 |3 K              c2=zeros(2*n,1);
' h+ _6 ]$ D/ b, b+ r* ?              c=[c1;c2];
               x=zeros(1,n);
               r=zeros(1,n);
               s=zeros(1,n);
               x(1)=4;
               r(1)=-4;
4 q; S# ?" |  s/ ~" {* n) x               a(1,=[x,r,s];
               x(1)=1;
               r(1)=0;
: M7 i8 q  M' p' r6 ^               s(1)=1;
               a(2,=[x,r,s];
2 Y5 E* H8 O4 o2 A/ q  A0 B               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
7 R& G  m; f6 W  x! P0 O; t" j               [x,f]=linprog(c,[],[],a,b,v)
end

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss