QQ登录

只需要一步,快速开始

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

极限测试之Matlab与Forcal动态生成函数的效率

[复制链接]
字体大小: 正常 放大
forcal 实名认证       

45

主题

3

听众

282

积分

升级  91%

  • TA的每日心情
    难过
    2012-8-27 18:22
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2011-8-3 07:40 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    程序动态生成一些函数,然后调用这些函数,体现了脚本的优势。
      j( ^0 M7 Q5 ]4 P2 W
    + Q1 C7 |2 r  E6 m7 y' z) L) W* T2 b以下Forcal程序动态编译了函数 f(x,y)=x+y ,然后执行该函数,共进行了10000次。
    1. !using("sys","string");  {$ ^* @) Q; |9 T5 H3 r
    2. mvar:
    3. , f9 ]2 m\\" j; b$ C0 J. f9 o1 m
    4. t0=clock();
    5. . C0 [+ w3 ]- z. a6 I3 }/ Q, b
    6. main(:pf,a,i,s)=/ Y+ t3 y. V8 M) t7 ]. X
    7.   oo{pf=fcfor()},         //申请保存Forcal表达式句柄的FcData数据指针) @, e. o\\" a, r, Q. e, s
    8.   s=0, i=0,
    9. * U& \9 c% i% b
    10.   while{++i<=10000,6 C7 y1 i5 l\\" `, N0 H9 O/ r
    11.     oo{comfors(pf,ws("f(x,y)=")+ws("x+y"),2,100,0,0)},    //将Forcal动态字符串编译为全局的实数表达式& Q# [+ e\\" q% J9 x5 v) |* I% t
    12.     a=HFor("f"),          //获得函数f的句柄1 A7 m4 ?+ B- P6 K
    13.     s=s+a[1,0]            //动态调用函数f
    14. 4 B3 S+ j5 D  |& s: X9 a, s: g
    15.   },0 c& V# V, n, D2 B- A5 v- ?: O& c
    16.   s;/ L7 [% E5 C, Q5 l* c4 Q5 B8 m  Z
    17. [clock()-t0]/1000;
    结果:
    + M* Q' w4 v* A% z10000.
      i- Y/ J, C4 d8 [  ?4 E/ C3 s0.359  秒
    2 w7 ~$ `* A# o/ A  U: J3 M) S$ M9 m$ [; P) D" ]( I6 {  R6 a% Y" j; L: j4 p
    matlab完成类似的功能参考:http://forum.simwe.com/archiver/tid-983780.html. P- \- i% k2 [
    * i  p, x( b; Z# h# p, @5 z
    由于matlab的函数必须存为磁盘文件才能使用(这应该是一个很大的缺点),故动态编译运行的效率必然很低。
    # z& x6 `# X# j. S( Z2 f1 I
    / w6 I1 C) d6 `$ V2 k程序动态生成函数的一个应用是:数据拟合中,常常要为数据寻找一个最佳的拟合公式,此时可按一定的规则生成一些函数,然后检查这些函数哪个最佳的。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-4-11 20:17 , Processed in 0.393745 second(s), 50 queries .

    回顶部