Question

我有一些代码需要很长时间才能运行（几个小时），我认为这是因为它在if语句中进行了大量的比较。我希望它运行得更快，有没有人有任何有用的建议来改善运行时？如果有人对减慢代码的速度有不同的看法，那么我可以尝试修复它，我们将不胜感激。

xPI = zeros(1,1783);
argList2 = zeros(1,1783);
aspList2 = zeros(1,1783);
cysList2 = zeros(1,1783);
gluList2 = zeros(1,1783);
hisList2 = zeros(1,1783);
lysList2 = zeros(1,1783);
tyrList2 = zeros(1,1783);

minList= xlsread('20110627.xls','CM19:CM25');
maxList= xlsread('20110627.xls','CN19:CN25');
N = length(pIList);
for i = 1:N
    if (argList(i)>= minList(1) && argList(i) <= maxList(1)) ...
        && (aspList(i)>= minList(2) && aspList(i) <= maxList(2)) ...
        && (cysList(i)>= minList(3) && cysList(i) <= maxList(3)) ...
        && (gluList(i)>= minList(4) && gluList(i) <= maxList(4)) ...
        && (hisList(i)>= minList(5) && hisList(i) <= maxList(5)) ...
        && (lysList(i)>= minList(6) && lysList(i) <= maxList(6)) ...
        && (tyrList(i)>= minList(7) && tyrList(i) <= maxList(7))

        xPI(i) = pIList(i);
        argList2(i) = argList(i);
        aspList2(i) = aspList(i);
        cysList2(i) = cysList(i);
        gluList2(i) = gluList(i);
        hisList2(i) = hisList(i);
        lysList2(i) = lysList(i);
        tyrList2(i) = tyrList(i);
        disp('passed test');
    end
end

Answer 1

您可以尝试对代码进行矢量化;我编写了一些示例数据集，并复制了您在下面执行的一些操作。

matA1 = floor(rand(10)*1000); 
matB1 = floor(rand(10)*1000);

matA2 = zeros(10); 
matB2 = zeros(10);

minList = [10, 20]; 
maxList = [100, 200];

indicesToCopy = ( matA1 >= minList(1) ) & ( matA1 <= maxList(1) ) & ( matB1 >= minList(2) ) & ( matB1 <= maxList(2) );

matA2(indicesToCopy) = matA1(indicesToCopy); 
matB2(indicesToCopy) = matB1(indicesToCopy);

不知道这是否更快，你必须尝试一下。

编辑：
这并不重要，因为你只打了两个电话，但是xlsread非常慢。您可以使用函数的这种变体语法来加速这些调用。

num = xlsread(filename, sheet, 'range', 'basic')

问题是忽略了range参数并且读取了整个工作表，因此您必须正确地对结果进行索引编制。

Answer 2

使用profiler查看哪些行或函数的执行时间最长。

您可以通过vectorizing代码大幅提高执行速度。这意味着使用一次操作整个向量的操作，而不是使用for - 循环来迭代它。类似的东西：

% make a logical vector indicating what you want to include
ii = (argList >= minList(1) & argList  <= maxList(1)) & ...

% use it
argList2(ii) = arglist(ii); % copies over every element where the corresponding ii is 1
...

在Matlab中改进运行时？

2 个答案: