极限测试之Matlab与Forcal编译计算含自变量的简单长表达式

forcal

Matlab 2009a编译了如下简单长函数，多次编译，耗时均大约为3分钟：

f=@(x)x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
/ |9 P8 m5 h- Jx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
3 b7 N& J" Y1 c6 `  X1;

完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

点击文件名下载附件
下载积分: 体力 -2 点

测试如下代码速度：

1. tic;2 m/ M# {- J' ~! r! t8 N/ ]
   
2. s=0;
   * |/ M; J% \. j. J5 e1 [7 u+ G. ^
3. for i=0:1000
   / w* }\" X! U% N
4. s=s+f(1);; r9 [2 Q: J4 S9 x5 ~
   
5. end
   & I\" X. n; `- r& c. F: v/ c: d, Y
6. s2 ^6 K9 T/ x2 t\" a! {\" b
   
7. toc$ M. w! r$ x4 I
   
8. 2 J7 i' a1 O5 C/ c2 F5 R
   
9. s =
   ; U- z8 f* R0 v7 i2 L
10. 
    2 H# ]0 `! G/ ?' Y$ N# V
11.     671000330 N* r# ?! s3 v3 g. L8 e( W6 _  N
    
12. , g3 N+ b  |6 B! {) T$ W) m
    
13. Elapsed time is 11.902260 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

3 i" j. H) W8 f& Y. F& d' L5 jForcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：

f(x)=
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
; U8 B4 \- z0 l1 \, o4 P0 Dx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
( T: M* g& c: y% A6 V... ...（省略）
. M6 J" Y6 d+ E8 H8 M; P# @x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
. n+ H% e( ]/ u2 Q1;
( Z" ~6 t6 m- D( {  M7 @
完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

点击文件名下载附件
下载积分: 体力 -2 点

. V+ ]& p8 ~2 m% N5 w
/ Q$ f- c; }2 j# v) Q2 t测试如下代码速度：

1. mvar:
   ( e2 E- v: I! u& u. Y- q5 s  j
2. t0=sys::clock(),
   * a. Y0 t3 W. K4 D/ `* k% h
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;
   ! L: Z* |0 `2 t
4. [sys::clock()-t0]/1000;

复制代码

结果：
$ w' E- D; W9 r0 x0 b! F67100033.
0.64  秒
$ f$ L2 J" ]. n0 ^
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
9 a& _2 \, T1 g
8 o7 b/ d5 U$ j3 w; @- h可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
2 {+ O) e1 R" R' R
6 M7 t& ^% J1 ?. q4 g7 |# V＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

+ v2 }# Y3 j" j5 {5 Q为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

4 }1 e4 o4 l2 S% A6 }  ~' T7 KForcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
, y0 t  J0 L1 H' v, ?) l+ w
1 u( ^- V/ @% \9 s3 h1 ~6 Y$ V目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。

对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

没来得急看，应该不错，先帮你顶

yunbuhuiku

如果跟导师讲不清楚，那么就把他搞胡涂吧!

* \0 F5 K% j8 @* v. i! T- `% z0 S0 }来电来函..欢迎洽购...
+ Y% b3 k8 L- A# m3 b

wangchen881202

“拆东墙补西墙”之升级版“拆别人墙补自己墙”

支持一下
: b; V6 r! i' h9 O& [' ^