极限测试之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+...
... ...（省略）
; V- W5 C: Z: b8 Z+ \. dx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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;

- }+ T1 ~  Q! f: ]- P完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

测试如下代码速度：

1. tic;
   7 r! l# S& ~\" b- x6 R' Q! d
2. s=0;
   : T8 Y! t4 j# S  U) J3 g! i2 n% W
3. for i=0:1000
   \" @% w\" x$ E9 ?1 T5 ~
4. s=s+f(1);% y* l: Q; [5 L( b% [! l
   
5. end# D! J( U0 s% ?! ]\" z1 e; Q6 u
   
6. s
   8 n0 o, H. Y# f\" l
7. toc. r% [% h, x1 d: O+ t
   
8. 2 {) k; }2 B5 R
   
9. s =- b& u) Y! f% X! G! E\" ]) @
   
10. : Z7 l$ P$ ?( [& y* ]2 i' \) `) F) ~
    
11.     67100033
    5 g; c# z0 m1 Z+ N
12. 
    # m+ L4 y\" a; Q, e
13. Elapsed time is 11.902260 seconds.

复制代码

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

8 y* z( W" P4 G. |3 hForcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：
3 M* _' S5 o  n7 o0 z
f(x)=
8 F  f0 `# ~5 s% h9 u: e# Ox+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
6 z" {  o' C. c$ B/ ]x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
! A3 u. A! `, u. `7 h. a8 Rx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 a/ z6 R4 H; f1 P5 ~) K... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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;
- U, _+ u; b/ H; A
8 w4 c2 _) Q1 {+ o% y完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

# C) `" e+ n' Y  {
测试如下代码速度：

1. mvar:
   1 a1 r/ a( f\" b1 A- g' t
2. t0=sys::clock(),# o0 ?# x! d\" ?# L3 F
   
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;5 N# G/ R+ _+ O  N' h
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
. g- Z9 \: u2 e! o/ _! p67100033.
0.64  秒
( |. x' a$ I! S2 L7 r* P: a
: N5 h( x  b; ~( ?7 m) ~＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
* E) E1 J- i8 \$ Q
4 _! D- V4 Z4 y. `7 ]可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。

1 ]+ f# X3 B1 p! n- M5 f＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
) ?: g) s0 q6 n
/ c0 i) t* R- R' |; g- }4 v为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
' o, `5 n' M- j/ c
' W# U' v/ K3 f8 L+ ?7 b目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。

对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

+ E6 N' t( k. f% Q0 H来电来函..欢迎洽购...
& d: G4 c2 \. K1 G9 Y5 G" H) a

wangchen881202

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

支持一下