双循环优化

Question

我已经编写了一个我必须控制的代码，如果位置 (x,y)（保存在矩阵垫中）位于以 (posx,posy) 为中心的圆形对象内。如果是这样，该点将获得值 val 否则为零。
我的代码看起来像这样，但事实上它被宣传为不在 matlab 中使用循环。由于我使用的不是 1 个而是 2 个循环，我想知道是否有更有效的方法来解决我的问题。

Mat = zeros(300); %creates my coordinate system with zeros
...
for i =lowlimitx:highlimitx %variable boundary of my object
   for j=lowlimity:highlimity
      helpsqrdstnc = abs(posx-i)^2 + abs(posy-j)^2; %square distance from center
      if helpsqrdstnc < radius^2
         Mat(i,j)= val(helpsqrdstnc); 
       end
   end
end

Answer 1

优化 matlab 代码的常用方法是向量化操作。这是因为内置函数和运算符通常要快得多。对于您的情况，这将为您留下以下代码：

Mat = zeros(300); %creates my coordinate system with zeros
...
xSq = abs(posx-(lowlimitx:highlimitx)).^2;
ySq = abs(posy-(lowlimity:highlimity)).^2;

helpsqrdstnc = bsxfun(@plus,xSq,ySq.'); %bsxfun to do [xSq(1)+ySq(1),xSq(2)+ySq(1),...; xSq(1)+ySq(2),xSq(2)+ySq(2)...; ...]
Mat(helpsqrdstnc < radius^2)= val(helpsqrdstnc(helpsqrdstnc < radius^2)); 

其中helpsqrdstnc 的大小必须与 Mat 相同。此处可能还需要进行 reshape，但如果您获得列向量，您会自己注意到这一点。

这当然假设radius、posx 和posy 是不变的，但阅读问题似乎就是这种情况。但是，我不确切知道 val 的外观，所以我没有设法测试代码。我也认为val(helpsqrdstnc) 很乏味，因为它指的是距离，它不一定是整数。