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

玖玖

对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
- h, Z) @3 R& W/ ?. Nmain(-xn)
+ l' r6 p& ~5 Q% ?4 B- E0 ]st
4x1-4r1=1;
x1+s1=1;
9 |9 B3 U+ T0 _7 e; [. M7 }4xj-x(j-1)-4rj=0;(j=2,3,4......)
- ^$ @0 s+ u' n9 S. H4xj+x(j-1)+4sj=0;(j=2,3,4......)
xj,rj,sj>=0;(j=2,3,4......)

! B9 U+ X0 K7 \/ i: J/ m请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
  g9 N6 t) Y+ o4 W0 h其中j，n是下标 。

玖玖

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

玖玖

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

玖玖

    syms n
" J$ [% p( C7 Y1 t: g4 Hn=input('input n please.(press enter)n=')
5 T0 O; i$ h( K1 y2 r  Q5 t, ~if(n==1)
* O$ K% d8 R+ k% w9 g    a=[4,-4,0;1,0,1];
0 M: e- y8 q- T7 U, l    b=[1;1];
    v=zeros(1,3);
    c=[-1;0;0];
! _" \5 ?0 ^1 q) u    [x,f]=linprog(c,[],[],a,b,v)
: A7 V* D- I  k! {! G. ]2 Qelse
    a=zeros(2*n,3*n);
    for j=2:n
) t# T4 G2 R& h; w                   x=zeros(1,n);
1 X! i! F6 a1 Y+ p& E/ N5 `                   r=zeros(1,n);
                   s=zeros(1,n);
8 P: t( x: P- v, X, I6 x                   x(j)=4;
0 W4 ~$ H/ \) Q2 Q' |% |: f4 G                   x(j-1)=-1;
9 O) ?- `& C: J5 p5 J# e                   r(j)=-4;
+ r5 E! I9 v2 _/ V, [                   a((j+1),=[x,r,s];
                   x(j)=4;
) ^& O# J7 B) y0 h  ^  x                   x(j-1)=1;
. ^4 Q8 ^5 q$ T( g0 [                   r(j)=0;
0 \: o3 y6 I+ N+ |8 D( G/ @/ G  Z* {                   s(j)=4;
! Q2 B, @' U0 O4 ^8 N                   a((j+n),=[x,r,s];
    end
4 _/ r3 T* {! _% d7 L4 }* h9 W               v=zeros(1,3*n);
0 f* F& g4 t8 ~+ @: s- I" U              c1=-ones(n,1);
- y- t' S/ u5 Q: Q3 g  X              c2=zeros(2*n,1);
6 [$ I, y5 P1 B' x/ [9 S2 v9 T              c=[c1;c2];
               x=zeros(1,n);
% z! A# j  H5 p0 u               r=zeros(1,n);
7 M) G2 F7 ^, K+ I" {  u* M7 B               s=zeros(1,n);
5 {% f0 ]( Z3 I) f/ L               x(1)=4;
               r(1)=-4;
               a(1,=[x,r,s];
               x(1)=1;
               r(1)=0;
: ~1 Y4 Q# K* I0 u, t               s(1)=1;
               a(2,=[x,r,s];
, @0 `! B3 U* Y: k               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
               [x,f]=linprog(c,[],[],a,b,v)
" ?7 X% W8 H# A) K; S9 B0 E5 }end

江边一粟

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

qianchuan11

/wx ssssssssssssssssssss