2011-11-01 35 views
4

我是新来的matlab,所以我不知道所有的快捷键matlab已使代码更有效和更快。我一直在matlab中一起攻击一些作业,并着重于完成作业而不是效率。现在我发现我花的时间比实际编码时多。下面是一个需要永久完成的嵌套for循环的头痛。有没有一种更快或有效的编码方式,没有太多的前提?更高效的Matlab代码请

for i = 1:ysize 
for j = 1:xsize 
    MArr = zeros(windowSize^2, 2, 2); 
    for i2 = i - floor(windowSize/2): i + floor(windowSize/2) 
     if i2 > 0 && i2 < ysize + 1 
      for j2 = j - floor(windowSize/2): j + floor(windowSize/2) 
       if j2 > 0 && j2 < xsize + 1 
        mat = weight*[mappedGX(i2,j2)^2, mappedGX(i2,j2)*mappedGY(i2,j2); mappedGX(i2,j2)*mappedGY(i2,j2), mappedGY(i2,j2)^2]; 
        for i3 = 1:2 
         for j3 = 1:2 
          MArr(windowSize*(j2-(j - floor(windowSize/2))+1) + (i2-(i - floor(windowSize/2)) + 1),i3,j3) = mat(i3,j3); 
         end 
        end 
       end 
      end 
     end 
    end 
    Msum = zeros(2,2); 
    for k = size(MArr) 
     for i2 = 1:2 
      for j2 = 1:2 
       Msum = Msum + MArr(k,i2,j2); 
      end 
     end 
    end 
    R(i,j) = det(Msum) - alpha*(trace(Msum)^2); 
    R = -1 * R; 
end 
end 
+9

请告诉我们它正在做什么的要点,以便我们了解您的代码正在做什么。看到名为“i,i2,i3”的变量绝不是一个好兆头.... – whatsisname

回答

8

而不是循环,使用冒号。例如:

    for i3 = 1:2 
         for j3 = 1:2 
          MArr(windowSize*(j2-(j - floor(windowSize/2))+1) + (i2-(i - floor(windowSize/2)) + 1),i3,j3) = mat(i3,j3); 
         end 
        end 

可以作为写:

MArr(windowSize*(j2-(j-floor(windowSize/2))+1)+(i2-(i-floor(windowSize/2))+1),:,:)=mat; 

后你会发现在那里可以做到这一点所有的地方,学会使用索引,而不是循环的,例如,

i2 = i - floor(windowSize/2): i + floor(windowSize/2); 
i2=i2(i2>0 && i2<ysize+1); 
j2 = j - floor(windowSize/2): j + floor(windowSize/2); 
j2=j2(j2>0 && j2<xsize+1); 
mat = weight*[mappedGX(i2,j2)^2, mappedGX(i2,j2)*mappedGY(i2,j2); 

(对于高级用户请注意:如果mappedGX是矩阵,并且i2/j2不代表矩形子矩阵,则最后一行可能不起作用。在这种情况下,您将需要sub2ind()