matlab中的稀疏索引

Question

我有一个很长的代码，其中充满了以下“if”，matlab 编辑器给了我如下建议：

这种稀疏索引表达式可能很慢

mt = rand(200,200);
[c r] = size(mt);
T = sparse(r*c,2);

for i = 1:c 
    for j = 1:r 
        if(ind(j,i)==1)
            templat = template + 1;
            T((i-1)*r+j,2)=100000;
        end
    end;
end;

有什么方法可以让我更快地编写代码并按照 matlab 的建议进行操作吗？ （代码可能无法运行，因为我只是选择了几行并试图显示问题）

Answer 1

嵌套的for's 和if 等价于以下矢量化代码：

[jj, ii] = find(ind==1); %// jj is rows, ii is columns
templat = template + numel(ii);
T((ii-1)*r+jj,2) = 10000;