了解矢量化答案

【问题标题】：Understanding Vectorization了解矢量化
【发布时间】：2013-06-13 15:24:26
【问题描述】：

我正在尝试使一些代码更高效，并且想知道如何对其进行矢量化：

%==========================================================================
% MinutesInDayTable.m
% 
% Creates a table identifying every minute as an integer from 0 to 2400.

minuteTableInDay=zeros(24*60,1);
k=1;
for i=1:24
    for j=1:60
        minuteTableInDay(k) = ((i-1)*100+(j-1))*100;
        k=k+1;
    end
end
%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ end of code ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

【问题讨论】：

请注意，您的代码和描述不完全匹配。

标签： matlab optimization vectorization

【解决方案1】：

这是另一个：

minuteTableInDay = (datestr(0:1/24/60:1, 'HHMMSS')-48) * 10.^(5:-1:0)'

玩得开心:)

【讨论】：

对于那些试图阅读它的人来说，48 与小时或天无关，而是'0' 的字符代码，所以“正常”（无聊）的方式是@ 987654323@
@DennisJaheruddin：当然，num2str 应该是str2num。请注意，“无聊的解决方案”要慢很多（str2num 不是内置的）。

【解决方案2】：

这很简单，向量化这个计算的一般方法是：

[I, J] = meshgrid(1:24, 1:60);                      %// 2-D arrays for I and J
minuteTableInDay = ((I - 1) * 100 + (J - 1)) * 100; %// Compute all values at once
minuteTableInDay = minuteTableInDay(:);             %// Convert to a column vector

向量化的关键是为 MATLAB 提供所有值，以便可以一次对所有元素执行操作。

在您的情况下，您的代码包含基本的算术函数，因此您可以进一步减少代码（类似于 Oleg 的建议）：

minuteTableInDay = bsxfun(@plus, (0:59)', (0:23) * 100) * 100;

【讨论】：

猜猜它实际上是相同的解决方案，但我会这样做：[H, M] = meshgrid(0:100:2300,0:59);minuteTableInDay=H(:)+M(:);

【解决方案3】：

您可以创建从0,...,115900,120000,...,235900 开始的时间：

out = bsxfun(@plus, (0:59)' , 0:23)*100;
out = out(:);

如果你希望它是100,...,115900,120000,...,240000，添加这一行：

out = [out(2:end); 240000]

【讨论】：

+1。不过，与原始代码相比，您会少 100 倍。
@EitanT 谢谢，一开始我没明白他也想要秒。