QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2242|回复: 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个约束)8 {# d" G; b" n* M9 r) l/ Z
    main(-xn)  e  x* I# _; k; @$ o
    st
    ) X0 g  Q* B! m4x1-4r1=1;
    " m2 _# U0 M: r# }2 ix1+s1=1;: A/ i7 d' r) W* r3 m3 L
    4xj-x(j-1)-4rj=0;(j=2,3,4......)
    # e7 Y- N1 ~* j* r4xj+x(j-1)+4sj=0;(j=2,3,4......)
      S* A2 \; m8 W  a* N  u& @xj,rj,sj>=0;(j=2,3,4......)1 J( O8 M' }! L. E
    ( x; b$ G( K6 s/ ~$ @  ]; |
    请分别对n的不同取值(如n=2,10,50等),利用MATLAB软件求解上述规划。
    7 C9 o( g6 i& F4 ]0 R0 v其中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
    - y# _$ g% t. b$ v) h  z' [( m8 tn=input('input n please.(press enter)n=')
    & Y7 w( Q, K. B: U7 m% \* p: c/ Oif(n==1)* Q7 ]  Q) h* u" K5 s' N) h5 r
        a=[4,-4,0;1,0,1];/ N" s4 ]: K& p% n5 C1 B$ i3 n
        b=[1;1];
      A, r; P  O& Q  d' {7 @# t    v=zeros(1,3);4 _0 N6 {* y: A/ h% l% W
        c=[-1;0;0];
    / \; Y5 N2 i5 o; k    [x,f]=linprog(c,[],[],a,b,v): ?) H& A% R' R1 F/ e
    else7 J) K8 B% W, I4 ]% s
        a=zeros(2*n,3*n);- A5 d  Q* i5 t& i" Q/ x) W. V
        for j=2:n" n: h( M9 z3 U9 P' P9 b8 O' @
                       x=zeros(1,n);
    2 t# G! O, l3 t  r                   r=zeros(1,n);7 t$ X1 d1 O) s
                       s=zeros(1,n);
    ' R* |7 M" m4 C, K% _' }                   x(j)=4;6 }2 ^9 n$ }- G% K% w% ]7 p. @5 O
                       x(j-1)=-1;7 C2 {2 \$ U) S1 g8 m8 q, f) g
                       r(j)=-4;5 Y2 K) D, w& g* W# }" s* |% k5 \
                       a((j+1),=[x,r,s];
    9 `, L  t& x/ j! o! g                   x(j)=4;, v5 Y+ {8 Y$ Y+ x4 n
                       x(j-1)=1;
    " U7 a9 N9 a% e  O                   r(j)=0;0 M9 k5 B# Y3 ]" R- [
                       s(j)=4;
    4 f  C" Q3 r8 i- y                   a((j+n),=[x,r,s];
    ' Y# o/ P- Q. H    end# e7 j$ r! D, `
                   v=zeros(1,3*n);
    8 w: q1 B* C% S- I5 i1 i              c1=-ones(n,1);
    % R6 `6 T* g7 V; `% X              c2=zeros(2*n,1);) F9 l, ]9 W# f
                  c=[c1;c2];
    ! v* R$ N1 Z3 @: j$ G8 z               x=zeros(1,n);
    8 f2 ^7 d3 T9 p& u               r=zeros(1,n);
    3 ]/ ]6 Q% O8 e- v! T- L; j               s=zeros(1,n);. x% U' J4 j$ @8 Y& u
                   x(1)=4;
    ) j: g' \  d+ _4 _               r(1)=-4;
    # R7 Y6 h0 I8 C7 j# J! I               a(1,=[x,r,s];
    - ^3 g% y+ C5 _3 X( X9 l               x(1)=1;9 ?+ V4 x  ?4 }
                   r(1)=0;" F* t: ]" v+ D' |4 [, k6 M
                   s(1)=1;8 ]6 g2 ^# O3 C5 n7 S/ s) o) ~
                   a(2,=[x,r,s];
    3 C9 w8 @7 N. O1 x2 p               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
    & W  K( O0 P2 Y5 C3 A4 {1 X               [x,f]=linprog(c,[],[],a,b,v)' C) l9 _4 |; L0 z0 R
    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 06:28 , Processed in 0.437268 second(s), 77 queries .

    回顶部