极限测试之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+...
3 W  `4 J. [4 `  T0 P+ `9 Ax+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
7 k7 q7 Q8 R% _6 y' ~8 x& C! M- A2 I... ...（省略）
5 O: _% Z2 l8 @9 ?# ~+ j0 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+...
" F1 y7 Y$ z; U% ^/ i3 A5 y1;

3 t  }& g9 w! ]" ]: x. g完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

4 e. K# f3 o( C+ Z+ J4 M
测试如下代码速度：

1. tic;
   4 o* W9 E$ {8 `- P; J1 m
2. s=0;
   7 N  y% s\" ~' v. O
3. for i=0:1000) [8 n/ R+ P$ Z4 l+ J\" W
   
4. s=s+f(1);
   4 g8 h+ R/ y, P# h3 V1 T
5. end4 s! I6 M6 K; l+ W% o! `
   
6. s: `7 ]/ E: E* S0 O0 V1 `  o
   
7. toc) q7 k5 X! L2 {
   
8. ; v, W3 k+ m# V1 ]
   
9. s =
   ' x9 P$ V) N  @5 F2 c
10. & y% ~. V* I) G
    
11.     67100033
    : m: N\" s, e2 i8 U
12. 
    1 e& _/ f( l5 W, ^# @( n* @1 W4 _\" g
13. Elapsed time is 11.902260 seconds.

复制代码

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

3 H8 c+ H  {4 i0 ~' m4 l& T2 fForcal 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+
8 P  ?" s, G" y2 G- 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+
- H8 x6 _4 w% A; R' C/ yx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

$ L% R  z6 O0 f1 Y9 J- U* M9 k7 g4 q
测试如下代码速度：

1. mvar:  W  w$ z, y! {$ r& n
   
2. t0=sys::clock(),
   ( m3 T: ?& b. L4 ^6 S0 D( U
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;3 H3 }6 N7 i3 _+ k- P8 @/ O/ }2 ?
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
% M  k4 ?% A8 @1 J1 [: ^) t' Y# R67100033.
0.64  秒

: G( @( \0 p: W4 V, N# a＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

  x3 x, k! g% V" w* c可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
! u  F& j, m8 I, p) j
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
$ R* E$ ?0 T5 u' C
为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
) @! F: ~5 N9 W6 b
. j1 A0 N/ Y: N) ?' Y目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
# p' |2 T2 Y, O: _# {% S
对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

如果跟导师讲不清楚，那么就把他搞胡涂吧!
5 @; U& \% q, g
# p$ ]" \' |: i0 R# v5 a/ y来电来函..欢迎洽购...

wangchen881202

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

  w/ n: H0 W: c! I- ^支持一下