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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
& e% Y# L- B( @- omain(-xn)
& d8 |8 [5 s; [st
4x1-4r1=1;
3 O* z3 _+ L3 l8 U, u2 S7 Q1 Qx1+s1=1;
/ s& N6 n! {8 V! S6 m* z  _. w4xj-x(j-1)-4rj=0;(j=2,3,4......)
9 {# ?) v) D( ~, X! I4xj+x(j-1)+4sj=0;(j=2,3,4......)
xj,rj,sj>=0;(j=2,3,4......)

请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
n=input('input n please.(press enter)n=')
if(n==1)
& X: {  T8 J& @2 r" P8 u    a=[4,-4,0;1,0,1];
" q6 [" D9 d* J6 T, W8 \& `& S, J    b=[1;1];
$ w4 \# X1 b# X* _    v=zeros(1,3);
    c=[-1;0;0];
    [x,f]=linprog(c,[],[],a,b,v)
else
    a=zeros(2*n,3*n);
    for j=2:n
                   x=zeros(1,n);
                   r=zeros(1,n);
                   s=zeros(1,n);
9 j& h9 P: A+ p; Z; M% Z. L                   x(j)=4;
5 ^0 u/ T/ v5 {: D' Z* G0 y                   x(j-1)=-1;
                   r(j)=-4;
8 _( g/ {* w& s  t' m                   a((j+1),=[x,r,s];
* I; J0 P/ {8 O3 Z! W5 `* H                   x(j)=4;
                   x(j-1)=1;
! _3 _6 }4 S  K5 H0 v4 D                   r(j)=0;
1 _" ]1 x) p/ V/ _2 L0 d( L5 |                   s(j)=4;
) W& t; C: U1 \; [6 n                   a((j+n),=[x,r,s];
5 V9 |) s' r# _    end
5 \7 I3 C. ?. q0 z% q" j0 i* n/ E               v=zeros(1,3*n);
              c1=-ones(n,1);
              c2=zeros(2*n,1);
/ q' {5 u/ L% P* ]4 u# S              c=[c1;c2];
8 ?" j2 \5 R( C2 ]& E9 e# k               x=zeros(1,n);
, X5 \) h# R/ K               r=zeros(1,n);
! c, L, P2 \" O  w( n9 W. q               s=zeros(1,n);
' h( F& g' h2 k: P# W, U6 Z! ~               x(1)=4;
: G. X8 O$ V  U( E               r(1)=-4;
" j7 ~) e! Z4 N# C, p               a(1,=[x,r,s];
1 m$ |& Y. P+ c  n+ x" F               x(1)=1;
               r(1)=0;
4 X1 u3 l9 r4 p3 y% u               s(1)=1;
+ h- H2 D! w4 g# C( K* }8 H7 t               a(2,=[x,r,s];
) k- i4 I1 C+ k- G0 c               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
9 R$ ^1 ]2 k: a; P: }               [x,f]=linprog(c,[],[],a,b,v)
, \7 N; [/ S: Y1 v  s) S! l$ D  Zend

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss