Sum[Cos(rn(1)-Sin(rn(1.2)*x^(y/rn(2))+Cos(rn(1)-Sin(rn(1.2)*y^(x/rn (2)))))),0] T7 ?: ]& s9 t1 F& K
}; ! Q0 y# M9 P0 i: D
[clock()-t]/1000;
复制代码
结果:# d8 g" i1 T0 }( G
1008606.64947441" X% E* F I- p! q2 M! N
0.6417 N6 q$ x7 l5 z" W7 }" ^: E
. \8 o) c+ d* k$ S( P7 rForcal比Matlab稍慢些。) E5 J- b' P+ G M) t. U
" d# x m1 F; h
---------- 3 V; F% A5 W# B' }$ c5 a/ V $ }- d$ x9 Y7 j% y6 u M+ L% U再看循环效率。! t# g% {- k; s# h# m5 k6 |$ S
while{x<=1, //while循环算法; & \ b: h' i( @% w. T! Z! s! A
y=1, 8 K3 r K. @* ` c" @
while{y<=2, # N/ A' T+ N% i. Y
s=s+cos(1-sin(1.2*x^(y/2)+cos(1-sin(1.2*y^(x/2))))), 1 g5 d" X8 T1 s6 G# e
y=y+0.0009 3 N' m! ?9 O7 }- }) \! f
}, " L1 n2 R/ B9 ?
x=x+0.0011 / j; T6 o! F) \' y
}, 8 O" {$ K1 G1 }. t
s;6 W" W) X2 z7 E
[sys::clock()-t]/1000;
复制代码
结果:) W. i& [3 Z. _. V6 P& w' c
1008606.649474410 b- Y5 k0 q B
0.734 //时间,秒 % v m8 r; l; ~) J, H I) _ $ H: r ]5 e& v# Y1 W7 }6 u b# b+ B我很奇怪,在这个例子中,matlab的JIT加速器为什么没有起作用?, z O' Z+ k x0 U
: O1 M8 Y8 r4 j/ b+ ?2 t- X
-------! i5 _" _6 H0 k* i. k, @4 T