Matlab：for向量上的循環。奇怪的速度行爲？

Question

我在matlab上用for循環播放了一下，我知道經常可以避開它們，在這種情況下，速度要快得多。但是，如果我真的想通過矢量V的所有元素，我做了小測試：

n=50000000; 
V =1:n; 

s1 = 0; 
tic 
for x=V 
    s1 = s1+x; 
end 
toc 

s2 = 0; 
tic 
for ind=1:numel(V) 
    s2 = s2+V(ind); 
end 
toc 

S1等S2是相等的（正常的），但它需要24.63秒的第一個環和僅0.48秒爲第二個。

我對這些數字有點驚訝。它是否已知？有沒有人有任何解釋？

Answer 1

這可能是由內存分配造成的。情況1

​​

有聲明，創建副本的V。我們爲什麼知道？如果在循環內修改V，則x不會受到影響：它仍將貫穿原始的V值。

另外，在殼體2中的語句，

for ind=1:numel(V)

實際上並不創造矢量1:numel(V)。從help for,

由於從未創建索引向量，所以在for語句中出現冒號表達式時，長循環的內存效率更高。

事實上，不需要分配內存可能會增加速度，至少部分是因爲速度的增加。

要測試這個，讓我們將for ind=1:numel(V)更改爲for ind=[1:numel(V)]。這將強制創建矢量1:numel(V)。然後，運行時間應該與案例1相似，或者確實稍大一點，因爲我們仍然需要使用V(ind)索引到V。

這是我的電腦上運行時間：

% Case 1 
n=50000000; 
V =1:n; 
s1 = 0; 
tic 
for x=V 
    s1 = s1+x; 
end 
toc 

% Case 2 
s2 = 0; 
tic 
for ind=1:numel(V) 
    s2 = s2+V(ind); 
end 
toc 

% Case 3 
s3 = 0; 
tic 
for ind=[1:numel(V)] 
    s3 = s3+V(ind); 
end 
toc 

結果：

Elapsed time is 0.610825 seconds. 
Elapsed time is 0.182983 seconds. 
Elapsed time is 0.831321 seconds.