极限测试之Matlab与Forcal编译计算含自变量的简单长表达式

forcal

Matlab 2009a编译了如下简单长函数，多次编译，耗时均大约为3分钟：
/ W- V  S7 r. w& L% u
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+...
9 j/ L; s$ V# F" e5 [2 Z; n... ...（省略）
, V8 x& o( f  e6 B; 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+...
8 [% _4 k+ |  B; u1 K1;
- ^& O2 S6 f* j- i. g
完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

& ]8 Q4 [' z0 q! x7 U- D
测试如下代码速度：

1. tic;
   9 U: w) j3 v! a1 p
2. s=0;
   \" f8 g% j( g& w\" ], t, n! H
3. for i=0:1000. {+ q/ u: V) `8 K# Z, ]
   
4. s=s+f(1);& W6 t7 g- @$ P- I! ^5 A$ K
   
5. end3 z. C/ u8 d' r2 x( ?
   
6. s' T! D* J4 \2 L5 p9 `% t
   
7. toc
   ; }  R: f3 X4 x7 F8 u
8. 
   ) y0 N. y8 \7 ?
9. s =' H; k/ `( A$ P9 ~. P5 x
   
10. 
    3 Q3 n. j; U* ]( t7 S( o
11.     671000333 y, n\" P* o  [6 w
    
12. 
    ) d; U' x! s\" R  s3 n! s) x
13. Elapsed time is 11.902260 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
, o# P8 J. T$ T- Y9 I
Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：

- C. c3 ]) Q1 _" Q2 wf(x)=
7 w& X) z, h2 A3 l/ [4 Lx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 Q' W/ u/ R/ f4 i7 }x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 点

: ?6 O# v. y0 c% J; g+ C! M测试如下代码速度：

1. mvar:
   ! f% t: D) ?9 b$ s% [  n/ y6 b' L
2. t0=sys::clock(),
   + F8 E# n+ q  L& M5 J; M
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;. m! h- n8 c6 k) L9 i) P
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
3 j, L% A2 N6 ]" P' O$ z7 i67100033.
# n" k# `6 C$ d0 q5 f1 q* Z2 |/ \0.64  秒

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
: B0 {7 d* g7 a+ Z3 b
可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
% n1 ~# ~/ M% ~* h
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

6 {1 i/ _' o4 [! F2 ^为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？
8 \! X: @  r) X# s+ I- i8 b% j; Y
+ F  C: `( e+ cForcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。

目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。

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

wangchen881202

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

支持一下
/ q& l& L# T+ G; S8 V4 f2 u. _

yunbuhuiku

如果跟导师讲不清楚，那么就把他搞胡涂吧!
7 u& n# W. V. ^" b& T" ]
来电来函..欢迎洽购...

安平池伟

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

jmedwardliu1205

geilivable