介绍
每次经过 for 和 while 循环时,这些循环都会递增数据结构体的大小,这会对性能和内存的使用产生不利影响。反复重新调整数组大小往往需要 MATLAB花费额外的时间来寻找更大的连续内存块,然后将数组移入这些块中。通常,您可以通过预分配数组所需的最大空间量来缩短代码的执行时间。
解决方法
下面的代码显示了创建标量变量 x,然后在 for 循环中逐步增加 x 大小所需的时间量。
最低性能
tic
x = 0;
for k = 2:1000000
x = [x x(k-1) + 5];
end
toc
时间需要数分钟
中等性能
tic
x = 0;
for k = 2:1000000
x(k) = x(k-1) + 5;
end
toc
Elapsed time is 0.301528 seconds.
如果您为 x 预分配一个 1×1,000,000 的内存块并将其初始化为零,则代码的运行速度更快,这是因为无需反复为不断增长的数据结构体重新分配内存。
最高性能
tic
x = zeros(1, 1000000);
for k = 2:1000000
x(k) = x(k-1) + 5;
end
toc
Elapsed time is 0.011938 seconds.
对要初始化的数组类型使用适当的预分配函数:
对于数值数组,使用 zeros
对于字符数组,使用 cell
预分配非双精度类型的矩阵
当您预分配一个内存块来存储除 double 外的某类型的矩阵时,避免使用以下方法
A = int8(zeros(100));
该语句预分配了一个 100×100 的 int8 矩阵,方法是先创建一个由 double 值组成的满矩阵,然后将每个元素转换为 int8。以 int8 值的形式创建数组可节省时间和内存。例如:
A = zeros(100, 'int8');
数组 行和列数组性能比较
数组采用行存储的性能略微高于列,但是当有运算时候,列性能高于行性能,推荐将数组以列形式存储。
tic
t1 = zeros(100000000,1);
toc
tic
t2 = zeros(1,100000000);
toc
时间已过 0.000473 秒。
时间已过 0.000180 秒。
总结
建议使用预先分配存储空间,如果不知道初始数据大小,可以采用第二种解决方案。或者预先分配一个超大的数组,在存储时候记录存储个数,最后将多余的数据删除即可。