QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2209|回复: 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个约束)
    " a6 ~3 c8 p, X, Vmain(-xn)
    : v6 Y. G$ J, o8 X# F+ qst + m+ L( T* }2 I8 R5 B5 g: r6 w$ P
    4x1-4r1=1;/ v5 N# w2 J4 I  F$ n& e/ [. N
    x1+s1=1;/ T0 {3 K" C* ?  L4 x5 f( c
    4xj-x(j-1)-4rj=0;(j=2,3,4......)
    * Q# ?! ^$ B1 u1 T8 M- E3 f# H4xj+x(j-1)+4sj=0;(j=2,3,4......)
    8 x" q! A+ s! t8 Q3 Y( Gxj,rj,sj>=0;(j=2,3,4......)0 C, c& ?9 O/ `6 @
    9 t) h* W  r) O- w
    请分别对n的不同取值(如n=2,10,50等),利用MATLAB软件求解上述规划。
    * Y$ b+ C. K! {) u' K/ ?其中j,n是下标 。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信

    0

    主题

    11

    听众

    23

    积分

    升级  18.95%

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

    [LV.3]偶尔看看II

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

    使用道具 举报

    0

    主题

    7

    听众

    37

    积分

    升级  33.68%

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

    [LV.3]偶尔看看II

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

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

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

    [LV.4]偶尔看看III

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

    群组Matlab讨论组

        syms n! h0 n  b/ V$ Q* C/ ^) O. v8 f" P$ w
    n=input('input n please.(press enter)n=')
    5 F: L% ?- R* Z' Tif(n==1). N9 D5 p* y( W4 s
        a=[4,-4,0;1,0,1];( I1 r$ Z/ B* m& {5 O. Q
        b=[1;1];* R+ ^* F' P+ ]6 Q4 X; P
        v=zeros(1,3);
    : U$ k0 Y; @! c+ ^, c    c=[-1;0;0];0 P/ M4 ]# P& l" r' Y
        [x,f]=linprog(c,[],[],a,b,v)2 f( E: o1 f2 G% S# ^4 O
    else3 ?2 U4 w( _9 |" n
        a=zeros(2*n,3*n);2 {# O" v6 L9 D0 _3 u3 ^
        for j=2:n. B6 Z! I6 ?  v! T0 g) [
                       x=zeros(1,n);
    . C, F& w: _! U% c* D7 S( ]                   r=zeros(1,n);
    & M9 d7 x  z9 i$ N4 A+ U  N                   s=zeros(1,n);
    ' y, Q. Q: p  C5 V( l" U( I. z                   x(j)=4;2 v8 d/ E( e  T# G* i' I8 |
                       x(j-1)=-1;# @$ t& h& f4 k* F
                       r(j)=-4;
    7 O/ `/ A( t. m5 X3 X! o                   a((j+1),=[x,r,s];8 d) }! i. B& w! B
                       x(j)=4;
    " y' c8 p; i- n$ c                   x(j-1)=1;: o1 D3 [0 W( k- s! x
                       r(j)=0;) B/ B1 |% e' Y- o1 h
                       s(j)=4;
    & K0 {, a& n5 T; V                   a((j+n),=[x,r,s];
    5 Y5 z: c0 v6 l1 Q    end
    $ i5 }5 ^7 e1 k. |! a# d; d' k               v=zeros(1,3*n);
    5 O& _, e; m% x" C7 ?- _              c1=-ones(n,1);
    9 K& X. H, U+ x3 S8 G- I" j              c2=zeros(2*n,1);
    0 O! R  B# V6 \              c=[c1;c2];
    & C7 P6 i0 p* A- e- x               x=zeros(1,n);7 a# E2 `& |4 x& Q
                   r=zeros(1,n);5 \% M2 }4 j& e5 N5 V$ I8 `
                   s=zeros(1,n);/ R3 y" z) s+ [# \
                   x(1)=4;' w! s; S6 i4 D* Q
                   r(1)=-4;
    ( Q2 ?5 ^( C% b0 T" s               a(1,=[x,r,s];" A" K5 \: `8 i6 R2 @* h% |
                   x(1)=1;* \' p  u/ o% a' N' q  J5 J
                   r(1)=0;6 E# g/ R6 i# K$ U
                   s(1)=1;
    9 G0 K, m- ~7 G# Y$ e) e0 h               a(2,=[x,r,s];
    + c5 P; V( l/ c9 K4 E# H( i               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];* r( X4 Q  Q# w+ Y) {
                   [x,f]=linprog(c,[],[],a,b,v)
    0 ]% w3 y9 Z6 K, K: r% Aend
    回复

    使用道具 举报

    玖玖        

    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讨论组

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2025-10-28 01:06 , Processed in 0.516915 second(s), 77 queries .

    回顶部