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

forcal

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

7 J6 D8 _4 s6 \) w( L" Of=@(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+...
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+...
8 r8 q% {5 F5 y9 ^: @$ U2 r+ `  F... ...（省略）
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+...
: V) e7 X7 ^, j1;
* @# r! l1 X* c
完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

* P- \) v4 ?0 V5 ?% q8 r+ l
测试如下代码速度：

1. tic;
   8 ^0 @0 Z3 S% p1 p, m4 B9 c
2. s=0;
   / R\" F\" t4 l% w6 y1 j
3. for i=0:1000
   9 _0 |\" F6 K/ X5 x, ]& A# ]
4. s=s+f(1);
   / E+ r* \. ?2 @  R! ^! X
5. end
   ) f+ D) Y  H( Z4 ~* J3 o8 R) L
6. s
   $ W1 Z. K& R( [) c$ r* y
7. toc: v- q2 q6 N: B$ U7 z
   
8. 8 g, R% A+ T* h/ J. F% F
   
9. s =
   2 L\" w4 Y' N( t; ~% s$ [& N# Q
10. : j% d0 I* J$ M
    
11.     671000334 p7 M( f1 o# `4 X\" R. H% y
    
12. 9 g% V0 [) d; D/ r
    
13. Elapsed time is 11.902260 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
0 ~7 j* ~$ I! w0 t, h/ a
Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：
/ L( Q. x5 ^0 b' C! {! t
  f1 ]1 N. w$ @. O5 @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+
1 @& V( V# e3 l9 }& [) fx+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+
... ...（省略）
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 W* }; c. F$ S1;
# \5 w/ C7 Q, |) w7 }+ Z
完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

7 q0 m7 O7 _  E2 |" f" V1 J测试如下代码速度：

1. mvar:
   ) I; [! q) e: _7 o3 w2 J
2. t0=sys::clock(),
   , Z8 e: j. Z, V; z4 v  F- _
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;% n5 i4 U/ n3 }1 F2 o5 W) M! K2 x1 ?
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
67100033.
0.64  秒
" k' c' m4 v, ]% k5 t* ?7 H, h
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
( W: u! }* O  j* [3 ^. J3 ?% f
可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
3 V6 _6 |7 x, p! ^! w2 O
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

( \3 d' K$ [, w3 H/ I6 y为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？
  E" O( E+ u: |" A0 ]& k& W1 l
Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。

目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
) x3 [2 H' ~& p. N5 O; @
对于脚本，编译效率和运行效率同样重要。

wangchen881202

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

支持一下

yunbuhuiku

如果跟导师讲不清楚，那么就把他搞胡涂吧!
9 N) B5 P+ g1 w3 O# F+ d
来电来函..欢迎洽购...
6 S; `# m& k0 c/ d4 O5 f' Q2 z9 X7 K

安平池伟

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

jmedwardliu1205

geilivable