QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2245|回复: 5
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大
玖玖        

1

主题

8

听众

66

积分

升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    跳转到指定楼层
    1#
    发表于 2013-7-15 16:44 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    对于如下线性规划问题(有3n个决策变量(x,r,s)和2n个约束)
    , C6 }9 X6 m& O) x$ t0 ^: Lmain(-xn)
    & Z* H' W3 u. b) ~st 5 y8 D, i+ S; t$ [/ l
    4x1-4r1=1;& Q. p* j1 r5 p$ Q! w0 C7 \
    x1+s1=1;5 X: G" a. z& f# G7 S. {
    4xj-x(j-1)-4rj=0;(j=2,3,4......)2 s/ b+ r4 \* m) Q
    4xj+x(j-1)+4sj=0;(j=2,3,4......)
    0 U1 y3 w+ ^2 T7 x3 }xj,rj,sj>=0;(j=2,3,4......)# x& t: E6 r, B, ~
    & J1 a5 W' f9 c$ u! R4 a0 F2 T
    请分别对n的不同取值(如n=2,10,50等),利用MATLAB软件求解上述规划。9 u) m" @# O8 a6 n' B
    其中j,n是下标 。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信
    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    回复

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    回复

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

        syms n. \/ V, O* R8 u) B& n+ }! _2 n7 g: F
    n=input('input n please.(press enter)n=')
    ' w# c6 x9 N1 dif(n==1)" A. W# e* i' B1 k
        a=[4,-4,0;1,0,1];
    5 @2 o  C' v2 \6 [    b=[1;1];
    3 Y( W/ D! a8 f2 r- I    v=zeros(1,3);! T0 J0 d! i8 i/ J9 u7 v  Z
        c=[-1;0;0];' f5 O  f+ x  [) u6 A8 B: @
        [x,f]=linprog(c,[],[],a,b,v). J% J/ I: q1 Y$ o
    else! @5 \2 @% A* ~' H/ U0 l
        a=zeros(2*n,3*n);6 ?" Z. o; H( [0 t8 E
        for j=2:n$ w# h9 S7 H& G
                       x=zeros(1,n);! w1 a/ q; W( `) D* h
                       r=zeros(1,n);( |% [+ H  m: o- G( z4 }# }/ ~
                       s=zeros(1,n);
    ; R1 }; `# i. ~6 x. M                   x(j)=4;
    & d8 B9 J1 Y& n2 s                   x(j-1)=-1;7 L) W+ f% x6 Y$ _& \$ v. d- e
                       r(j)=-4;0 `2 H. C8 G7 g* |/ o
                       a((j+1),=[x,r,s];$ x* `# O' ~& T. N& X& @5 B, r( T; i/ q+ [
                       x(j)=4;# z8 z, K3 O6 Y
                       x(j-1)=1;' ~! f5 U9 e% t, o% k
                       r(j)=0;) x% z9 ?4 @' @# g" O3 N
                       s(j)=4;$ Y/ K% m0 `, L0 o% C* j
                       a((j+n),=[x,r,s];0 _( F& n6 x2 L8 `$ Q$ m
        end
    6 y+ j. p% _* a* U* S               v=zeros(1,3*n);4 B; \+ _, H! @: x- j2 ~4 F4 B
                  c1=-ones(n,1);, e7 _) J* D' {* z. f. ~
                  c2=zeros(2*n,1);
    + m8 W  R6 L. K3 w0 c              c=[c1;c2];
    ; b9 k% @+ r" V               x=zeros(1,n);' m. M' S7 n7 X
                   r=zeros(1,n);
    . F  f; c* p6 N               s=zeros(1,n);
    " Q2 B; D. I* h) b' p# r               x(1)=4;
    " D( S: G1 K! I, X( Q, ^& I               r(1)=-4;
    # r8 Q' p& ?+ ^. ^6 T" Q4 q               a(1,=[x,r,s];
    9 V! |( ?) h: k3 b4 I3 _5 p" e: i$ Q               x(1)=1;7 N) h  e5 F' [
                   r(1)=0;5 j2 D! l( ]8 @# R7 e6 _( x
                   s(1)=1;
    " X5 l' j1 Z+ i& f8 u" t               a(2,=[x,r,s];
    * _( q/ q2 t# |& K               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
    1 x3 t2 V. V% m3 r! o, {: n  q1 E  [* g               [x,f]=linprog(c,[],[],a,b,v)9 e4 X( |0 t/ N$ m0 N
    end
    回复

    使用道具 举报

    0

    主题

    7

    听众

    37

    积分

    升级  33.68%

  • TA的每日心情
    开心
    2014-1-29 21:12
  • 签到天数: 11 天

    [LV.3]偶尔看看II

    自我介绍
    一名医学生,对数学建模很感兴趣,目前正在参加培训,还望各位大神多多指教
    回复

    使用道具 举报

    0

    主题

    11

    听众

    23

    积分

    升级  18.95%

  • TA的每日心情
    擦汗
    2014-9-14 19:04
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱数学建模
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-11 23:49 , Processed in 0.457398 second(s), 77 queries .

    回顶部