极限测试之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+...
& T+ e: e& L$ M. B, 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+...
... ...（省略）
. V  y- W0 k1 a- j: e6 x; h' N- Zx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 t: Z+ E, a! I* i, m5 W1 Q
测试如下代码速度：

1. tic;% |$ A* X6 O5 M% m( e( G
   
2. s=0;4 M( R- R) F! p3 ^6 r5 e. ^
   
3. for i=0:1000& O3 U0 L1 |$ G\" ?) [; Z6 b( a
   
4. s=s+f(1);
   + [! j& v. r: ^3 y
5. end
   - u. n( W6 P% [0 P- W2 B\" G
6. s
   1 L! c7 n4 ~7 C$ @  A( y
7. toc( V0 F. |) z' `( f
   
8. 
   9 l# `: a! F' S4 [7 X$ F# G  @
9. s =
     S& [+ ]1 q+ \- q/ W( ]
10. 
    4 L; |1 I9 ]1 A) |3 B! G
11.     67100033! b8 |. S& S$ |( M2 `& [8 S
    
12. 
    5 f+ }; U3 D6 B
13. Elapsed time is 11.902260 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
% F! c! X( }7 _& |  H# }( A/ ]/ A
Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：

0 u( x5 Y. j3 O/ mf(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+
9 [0 @. q$ M" {3 Tx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
# F/ k( }; U, U5 N' nx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
... ...（省略）
# ]2 i  t# m; `; f/ w3 U: h6 A8 Sx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
$ D- {! J' V/ _7 G. w" i2 i6 t1;
+ i. D/ F! s' Z0 X' o0 Z7 V
完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

测试如下代码速度：

1. mvar:$ Y. A) S- Z' {, d
   
2. t0=sys::clock(),2 F6 w\" z) V% n9 M
   
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;2 H9 ~1 t) w. K7 o6 o! @
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
67100033.
* n( Y) w# ~2 y% }- w0.64  秒

+ w6 a, Q- c2 m# m# J6 U＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
8 G4 }0 O/ d8 Q1 X
可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。

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

为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

! r  Z6 o( t( }/ eForcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
- L9 b8 }7 T$ Q# v& D! ]
目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。

, @1 |2 v) I* P对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

来电来函..欢迎洽购...
' L" G' b+ q- z& X* r

wangchen881202

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

2 `- |9 c. x: T3 Y" D$ ^支持一下
  _  U! [: {. t; T0 x