数学建模社区-数学中国

标题: 我是新手，MATLAB接触不久。有一题不会，求教求教 [打印本页]
作者: 玖玖    时间: 2013-7-15 16:44
标题: 我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）2 q8 L/ a0 @+ D1 l7 M6 G( @  T
main(-xn)
$ M4 f7 G0 L9 m5 tst 8 r% {8 c5 K+ q
4x1-4r1=1;
% {3 h! |/ O  k2 K: W: O) Bx1+s1=1;; T  K4 J; W% t& \2 C
4xj-x(j-1)-4rj=0;(j=2,3,4......)
) ^4 }, C" U4 `/ [! P7 R# A# ~4xj+x(j-1)+4sj=0;(j=2,3,4......)* J0 U7 ^& P' i: a# V. P! y4 f0 k
xj,rj,sj>=0;(j=2,3,4......)* T5 b" M' Z/ u. I2 S
, z; {$ o1 w( |  x- z
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
7 ?+ N& ], ~6 v& V4 C9 i; t7 |其中j，n是下标 。
作者: 玖玖    时间: 2013-7-15 20:29
有大神在不？帮忙啊！！！
作者: 玖玖    时间: 2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者: 玖玖    时间: 2013-7-16 14:57
    syms n3 u/ x4 w$ }  i
n=input('input n please.(press enter)n='): c1 Z: ^2 T* J1 O# a( }1 b
if(n==1)
3 h9 c  A- e/ ]$ W5 {% L- I    a=[4,-4,0;1,0,1];: r: G# T) {. x
    b=[1;1];
( O- Q7 V4 @  o. u0 O+ Z    v=zeros(1,3);; W$ S  ~6 w: E5 m& J
    c=[-1;0;0];
" A5 e4 [+ h, B6 _+ Z    [x,f]=linprog(c,[],[],a,b,v)
8 e/ A* p* X- n) o1 U/ B& eelse/ c, ^, Q* y9 L) Y: s" Q  g/ {
    a=zeros(2*n,3*n);
% Y- A( R( R: U: n% D    for j=2:n
  {3 y+ v9 |5 z                   x=zeros(1,n);
# O/ B1 y$ `7 U  {3 F! p                   r=zeros(1,n);
) {3 y% |: h7 u2 Y                   s=zeros(1,n);  E3 z" Z- H  B% H4 G, e% g+ x
                   x(j)=4;) w# y+ _) S! O8 n
                   x(j-1)=-1;
# ]1 ]3 O4 I8 c% V& K: v( R8 R                   r(j)=-4;) ^+ D2 p% W1 j. G% Y& Z/ s
                   a((j+1),=[x,r,s];
' D, s9 t7 m6 N, _1 i                   x(j)=4;
; d1 X" [0 A  ^. Q" o                   x(j-1)=1;3 j. ]5 o6 N5 @$ s" B+ X' F$ s
                   r(j)=0;% f) k7 W; Z; O3 }5 ]- r5 l9 p1 d/ N
                   s(j)=4;
1 p7 O) b" \  q                   a((j+n),=[x,r,s];
  f7 u9 r* _5 }9 ~    end
" Q4 T5 U. o$ ?2 x3 J: q/ c               v=zeros(1,3*n);
$ I  i' z' y, `) [* g              c1=-ones(n,1);
, x; B( i& D: _) D# ]) K              c2=zeros(2*n,1);
. ^0 v) M4 x7 }              c=[c1;c2];
7 W. r! Q/ z3 q& L4 B               x=zeros(1,n);
6 \& W+ C: n  v& x7 \               r=zeros(1,n);
, m5 E7 Z9 j' |& }% H' v$ c               s=zeros(1,n);& y2 Q; o8 H9 Z: b, ^
               x(1)=4;6 @: a* X7 A) _+ L: L7 l
               r(1)=-4;
) {' o9 i' ~* g  z# e% j4 _               a(1,=[x,r,s];
: i2 j* P# c7 A0 y- T$ g# m0 o               x(1)=1;
, t. K2 i: z" f5 M, E; c5 M               r(1)=0;
, N% _5 N# @  ?0 K+ u               s(1)=1;
- |: i2 U, b. g; m( F               a(2,=[x,r,s];
; c6 d& S9 W" f3 E, [8 P$ p5 q5 V               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
/ ^: B5 }4 k" o' e7 ~               [x,f]=linprog(c,[],[],a,b,v)3 {( L- T' E& ~' c
end
作者: 江边一粟    时间: 2013-8-15 11:37
楼主，笑脸是什么？求解。。。昨晚就在做这题，到现在还没做出来
作者: qianchuan11    时间: 2014-8-11 20:30
/wx ssssssssssssssssssss



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5