QQ登录

只需要一步,快速开始

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

【转载】加速matlab运行的三个境界

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

1341

主题

738

听众

2万

积分

数学中国总编辑

  • TA的每日心情

    2016-11-18 10:46
  • 签到天数: 206 天

    [LV.7]常住居民III

    超级版主

    社区QQ达人 邮箱绑定达人 元老勋章 发帖功臣 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组2011年第一期数学建模

    群组第一期sas基础实训课堂

    群组第二届数模基础实训

    群组2012第二期MCM/ICM优秀

    群组MCM优秀论文解析专题

    跳转到指定楼层
    1#
    发表于 2010-1-25 22:39 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    本帖最后由 厚积薄发 于 2010-1-25 22:54 编辑
    ! P2 _6 Z; p8 F  t- ?8 |4 M$ R( ~+ T7 ^1 G
    加速matlab运行的三重境界

    加速matlab运行的三重境界6 {) _0 v+ e8 J! k# w/ t2 P
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    ! [) f& N3 A% S%%%%%%%%%%%%%%%%%%%%%%%
    $ C; V, h6 O, u' a
    一、 遵守Performance Acceleration的规则
    - d" {. n- P" i8 D+ q# j( ?1 l. o二、 遵守三条规则" l" {/ w2 N2 ?; t) `0 ]+ T1 c" L
    三、 绝招
    * u" ^* l7 m) }9 w( z" w%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6 B4 P+ [  e2 d- g7 [
    %%%%%%%%%%%%%%%%%%%%%%%

    一、 遵守Performance Acceleration的规则; ?  N  ]& p) D/ X# I
    关于什么是“Performance Acceleration”请参阅matlab的帮助文件。我只简要的将, r7 {8 [7 L) {& J+ Z
    其规则总结如下7条:& M, ?3 G) n+ E6 _0 q
    1、只有使用以下数据类型,matlab才会对其加速:% l; T2 w# _2 r# u! i
    logical,char,int8,uint8,int16,uint16,int32,uint32,double
    3 C9 ]. C* A) B8 p
    而语句中如果使用了非以上的数据类型则不会加速,如:numeric,cell,structu2 H  @* y8 M: k
    re,single,
    . k# T2 n2 Y5 R) c4 mfunction handle,java classes,user classes,int64,uint64* ~3 u1 c7 w3 |
    2
    matlab不会对超过三维的数组进行加速。
    " {8 d" A3 p4 S
    3、当使用for循环时,只有遵守以下规则才会被加速:afor循环的范围只用标量值3 h. V; H. p# _4 S6 E! p$ I( T
    来表示;/ }9 L! z- _& g/ L
    bfor循环内部的每一条语句都要满足上面的两条规则,即只使用支持加速的数" n' Z" i. ~5 H, R. S/ O
    据类型,只使用8 S! l9 I% P/ ~% @2 C4 u
    三维以下的数组;c、循环内只调用了内建函数(build-in function)。0 Q/ @' O! S. f" ~8 j  _' W
    4、当使用ifelseifwhileswitch时,其条件测试语句中只使用了标量值时,将
    # b/ n  t! P0 ^# p  j9 p加速运行。
    0 e6 L/ y' A% _4 @) M- C: O5、不要在一行中写入多条操作,这样会减慢运行速度。即不要有这样的语句:
    1 K6 L6 r8 @/ F1 n; Tx = a.name; for k=1:10000, sin(A(k)), end;
    * p% H* ?- R% H: R6
    、当某条操作改变了原来变量的数据类型或形状(大小,维数)时将会减慢运行速
      |# W* A5 G5 q" g5 ^+ |- R度。% P: u* i( P) R0 T
    7、应该这样使用复常量x = 7 + 2i,而不应该这样使用:x = 7 + 2*i,后者会降低
    . m$ A/ T2 s* W# `6 u; R$ N' |9 F运行速度。
    5 @- p5 d. i8 L7 h%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    ! {1 |2 Z  R- t0 m( B  g%%%%%%%%%%%%%%%%%%%%%%%
    - m$ w7 W% ^: ^- v
    二、 遵守三条规则
      {5 b7 V% t/ i; Q' c1、尽量避免使用循环,MATLAB的文档中写到“MATLAB is a matrix language, whic) M% L( p+ U7 N  D8 c
    h means it is designed
    : k3 b. i$ N( S+ {. s4 y) d: jfor vector and matrix operations. You can often speed up your M-file c
    2 h1 Y7 n5 ?' ]$ T7 e$ wode by using. t# q' ?+ j) N4 O- {. ]: f( r# S
    vectorizing algorithms that take advantage of this design. Vectorizati% a6 p+ i- e0 I+ F; j
    on means converting: m5 H# j2 W6 U, I2 j
    for and while loops to equivalent vector or matrix operations.”
    。改进
    - {9 U" s1 k8 u# Y% T% L& t
    这样的状况有两种方法:
    2 ~9 A5 C4 @: T* s! Ba、尽量用向量化的运算来代替循环操作。如将下面的程序:
    " ]4 ^4 \: T" R( Vi=0;' e8 W; a, \2 Q6 Z8 ]1 \+ I$ G! K
    for t = 0:.01:10' C0 q% G( E: l' l0 S, _
    i = i+1;, h% T+ C/ x: ]' z* Q& i
    y(i) = sin(t);
    1 P+ T1 P  t9 h9 D9 W) dend
    7 z  q! M' Q' i8 a/ m; E+ c
    替换为:
    + ~8 N5 p0 I6 Y( E; \- Qt = 0:.01:10;
    - ^) t& I' ?" t6 Ey = sin(t);1 }8 M. m; r( D/ N
    速度将会大大加快。最常用的使用vectorizing技术的函数有:Alldiffi0 \4 k+ `( r* X
    permute
    permute
    # S3 E+ z. d6 a" b' w
    reshapesqueezeanyfindlogicalprodshiftdimsub2indcums0 M; _! G9 Z) b
    um
    ind2sub

    ) P1 {9 d& p8 s: j0 F  O  v" @ndgridrepmatsortsum 等。
    $ M: w7 C8 S# N% D+ o1 `& X' ]请注意matlan文档中还有这样一句补充:“Before taking the time to
    3 `* V. c3 v# B# d4 O6 f: gvectorize your code, read the section on Performance Acceleration.8 g; N( c7 R. E  g
    You may be able to
    : o% @' T1 Q# f6 Pspeed up your program by just as much using the MATLAB JIT Accelera8 V; }8 x" v2 V" w
    tor instead of
    6 Y/ Z0 J( X+ Q! ?; W/ H, Avectorizing.”
    。何去何从,自己把握。
    # O+ [) r3 |* o
    b、在必须使用多重循环时下,如果两个循环执行的次数不同,则在循环的外环执
    * s& i8 T" U' W1 D行循环次数少的,. N) K, p2 `8 ~, m) h1 i0 ^: |
    内环执行循环次数多的。这样可以显著提高速度。
    ! q$ A! W; ]( S2a、预分配矩阵空间,即事先确定变量的大小,维数。这一类的函数有zeroson
    % \1 P9 g4 T, a- V/ Jes
    cellstruct

    . N6 m, w4 T4 i+ t0 v- H7 crepmat等。, B6 Y+ L* d0 }  Y0 ?- f
    b、当要预分配一个非double型变量时使用repmat函数以加速,如将以下代码:
    + `$ I9 L3 I4 s" n+ ]& gA = int8(zeros(100));
    2 Q3 b; ~. ~2 z  U1 T7 z5 W
    换成:
    7 B# S6 u0 N" O9 t3 Z' b* F" mA = repmat(int8(0), 100, 100);7 R+ n6 J( V& r5 G
    c
    、当需要扩充一个变量的大小、维数时使用repmat函数。
    " l2 t3 m( ]. g# S" m4 Q3a、优先使用matlab内建函数,将耗时的循环编写进MEX-File中以获得加速。1 R1 d1 `  W- [* e
    b、使用Functions而不是Scripts 4 D5 `' i8 v% L2 T  t
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 e9 e4 ]/ j& }" J) J( o
    %%%%%%%%%%%%%%%%%%%%%%%
    0 Y. I* ]$ e3 }" L' o6 {! b# j
    三、 绝招& B* v" k, `) a* g' J
    你也许觉得下面两条是屁话,但有时候它真的是解决问题的最好方法。$ |* x8 L7 T2 I' b/ W
    1、改用更有效的算法* b: M/ r5 v; Y  d
    2、采用Mex技术,或者利用matlab提供的工具将程序转化为C语言、Fortran语言。

    % e3 H/ ?$ p, l& e- G7 F7 Y

    # I" y4 u: K6 f# Z1 v  来源:iLoveMatlab
    9 c0 ], P$ D' M- j! F  Z$ c" A, }8 F9 ^9 q
    9 b+ y7 h+ T1 r8 y0 K
    在这里本人推荐使用‘绝招’。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    XINGQIBIN 实名认证       

    0

    主题

    5

    听众

    94

    积分

    升级  93.68%

    该用户从未签到

    群组数模应用

    回复

    使用道具 举报

    HSinB 实名认证       

    0

    主题

    5

    听众

    1651

    积分

    升级  65.1%

  • TA的每日心情
    开心
    2012-2-9 17:19
  • 签到天数: 2 天

    [LV.1]初来乍到

    新人进步奖 发帖功臣

    回复

    使用道具 举报

    D-lin 实名认证       

    22

    主题

    5

    听众

    465

    积分

    升级  55%

    该用户从未签到

    自我介绍
    人称玉树临风胜潘安,一支梨花压海棠

    新人进步奖

    回复

    使用道具 举报

    leo12ok 实名认证       

    3

    主题

    4

    听众

    556

    积分

    升级  85.33%

  • TA的每日心情
    擦汗
    2012-6-14 16:31
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    200 字节以内

    不支持自定义 Discuz! 代码

    新人进步奖

    回复

    使用道具 举报

    0

    主题

    3

    听众

    95

    积分

    升级  94.74%

    该用户从未签到

    自我介绍
    嘿嘿~  可爱的我~
    回复

    使用道具 举报

    isnowfy 实名认证       

    6

    主题

    4

    听众

    473

    积分

    升级  57.67%

  • TA的每日心情
    开心
    2013-3-7 02:00
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    200 字节以内
    不支持自定义 Discuz! 代码

    新人进步奖 最具活力勋章 发帖功臣

    群组Matlab讨论组

    mark。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
    回复

    使用道具 举报

    27

    主题

    4

    听众

    196

    积分

    升级  48%

    该用户从未签到

    收藏了~~~对现在建模很有用……我上次写的程序就运行了N久都没出结果……郁闷……
    回复

    使用道具 举报

    schnee        

    0

    主题

    4

    听众

    241

    积分

    升级  70.5%

  • TA的每日心情
    开心
    2013-9-3 12:40
  • 签到天数: 18 天

    [LV.4]偶尔看看III

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2025-8-20 04:28 , Processed in 0.821652 second(s), 99 queries .

    回顶部