极限测试之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+...
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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;

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

2011-7-30 12:06 上传

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

6 I6 e7 s) i9 O7 H) C+ m( y" G
测试如下代码速度：

1. tic;
   \" v& H6 E( s7 o$ _, d5 u
2. s=0;1 y1 x9 l1 h4 _. K2 q3 F
   
3. for i=0:10004 G- r7 K  p! E, J! L
   
4. s=s+f(1);
   2 I) X! C- v9 V% r; {
5. end
   . p) n! a6 a2 ]: T# q' t
6. s( G' I$ |& h! C
   
7. toc
   2 b  h: q  i4 E
8. . F/ T- p# K! C# J
   
9. s =) R  G9 Q/ P5 K
   
10. 6 K+ t9 X3 R* t  G
    
11.     67100033/ @  y4 p0 O/ M% t% J
    
12. $ r6 o! Y4 `% [( r\" u$ Q4 z; ~
    
13. Elapsed time is 11.902260 seconds.

复制代码

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

Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：
% S( |* v1 ]0 ^# x0 b/ P: s% R7 b1 _
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+
( c' T  O3 z3 R' b3 b1 Vx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
5 o  F% C+ }. a/ |... ...（省略）
+ q# M, N% u( h* i0 y6 v: gx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
/ @5 J- z8 S4 o# S8 O1;

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

2011-7-30 12:06 上传

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

测试如下代码速度：

1. mvar:9 u* y4 H4 I6 ^. N; T$ T  V
   
2. t0=sys::clock(),
   . w& r# q  e$ N' q+ q* K
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;# |9 o9 P7 j8 V$ K
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
  Z5 h3 \5 K2 P3 f0 v67100033.
/ x7 E  C  _8 u2 h1 R$ d: n0.64  秒
# n; A5 ?) V1 `! D
% Q: v. o; t/ F0 y, b/ u0 _$ X＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。

% D8 o- `% d  w8 ?0 F/ [( U' d$ ~7 T＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？
* j' I1 h  h7 [) A  C
Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。

) w2 N( N$ {* Q" e* I( q5 Z目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。

* {0 x7 i- A0 j* I6 ^7 b$ _- @对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

8 R+ T" s7 u3 U% n8 h5 f来电来函..欢迎洽购...
8 ]( |, r/ F- [! Z6 k; B" }3 X

wangchen881202

“拆东墙补西墙”之升级版“拆别人墙补自己墙”
/ X$ \' a8 N5 d( {
4 r$ }! l$ N# E; N& C  o; m- T7 m支持一下