极限测试之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+...
2 Y  z$ }$ M+ K5 bx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
# B9 f& V4 |& N6 `* m) 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+...
... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 K7 Q' U6 {3 C& j1;
" P; E% @1 H# m, [1 x5 ~. V
完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

" v- s( ~  q- P/ }) I: f* s
测试如下代码速度：

1. tic;
   / r# ?3 w8 J4 G% R: G/ {
2. s=0;
   # m2 c# x. A+ v
3. for i=0:1000
   ( ]; p0 P1 B1 o9 `- k
4. s=s+f(1);
   * W\" f  E1 M: C% U& ^
5. end
   ( O1 N6 u' W0 G+ f; \, h
6. s: z2 {( e2 g8 J; ]) Z
   
7. toc7 a/ _' v1 f- L0 U
   
8. 
   7 D\" e* g* L1 ?6 l# k3 J
9. s =
   - Z6 k% |  c! M, u6 n; K$ Y' O
10. & g' }6 z2 U9 n  j# U2 q7 Q! F: i6 w
    
11.     67100033
    8 O% q$ Z1 m, g
12. 
    * Y$ w4 r\" ~' m8 s3 g  _
13. Elapsed time is 11.902260 seconds.

复制代码

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

' b" V7 P- m/ t3 G" U! jForcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：

$ q9 ]3 B3 f- l0 X/ tf(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+
- x4 x2 i. K9 r+ @0 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+
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
) [4 O6 T7 B% @1 B9 a) f% g... ...（省略）
/ w  z) F5 ]  x" R! Hx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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;
3 m* K; D: W4 a
3 N+ O0 b8 b; C* O完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

测试如下代码速度：

1. mvar:
   % y+ p5 U# N/ ]
2. t0=sys::clock(),
   # ~6 a  l, y\" o4 V: o: B\" z
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;
   7 K. N4 z+ [/ U1 x0 Y
4. [sys::clock()-t0]/1000;

复制代码

结果：
" h  m( h8 R. v- L& [67100033.
0.64  秒

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

0 t" v: B- i6 C5 y; V3 j可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
- X$ y5 O# N: {( L. ^; h
3 [# v2 n, e4 J3 W) u/ F＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

4 ~' t: p# ^! Q7 ~为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？
" P9 L) y/ P0 v2 Y$ `# V7 O) d4 u
# `, B5 g- w/ x1 VForcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
' x" J/ E5 i7 Y  _; J' @
# m. z7 B* [. r; I% A目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
7 p8 C* ^) R$ _% D  _
对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

来电来函..欢迎洽购...

wangchen881202

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

支持一下