我很难找到一个比原始循环更快的解决方案……
生成测试数据
n_values = 10000;
数控= 100;
pre_matrix =兰德(n_values、数控);
深圳=大小(pre_matrix);
J =细胞(1、n_values);
为2 = 1: n_values
n0 =兰迪(数控,1);
J{2} =兰迪(数控1 n0);
结束
%原始循环
抽搐
为我= 1:n_values
pre_matrix (i, J{我})= 0;
结束
toc
%直接索引，但使用循环获取行索引(Rik的建议)
抽搐
我=细胞(大小(J));
为2 = 1: n_values
{二}= ii。*的(大小({2}));
结束
idx = sub2ind(sz, cat(2, I{:})， cat(2, J{:}));
pre_matrix (idx) = 0;
toc
%直接索引，没有显式循环
抽搐
我= cellfun(@(a,b) ones(size(a)))。* b J num2cell (1: n_values),“大学”, 0);
idx = sub2ind(sz, cat(2, I{:})， cat(2, J{:}));
pre_matrix (idx) = 0;
toc
结果:
运行时间为0.017597秒。< -循环
运行时间为0.047797秒。<-循环索引
运行时间为0.119106秒。<- cellfun索引
我猜罪魁祸首可能是循环中缺少的分号?如果你每次迭代都要打印出巨大的矩阵，那肯定需要一些时间!但是循环本身非常快。
