u2 = exp(i.*zeros(8192,8192));
for ite1 = 1:800
   for ite2 = 1:900
      u = exp(i.*sqrt(xy_2d- z_2d(ite1,ite2)));
      u2 = u2.*u;
   end
end

这个程序越来越慢，请问是什么原因？请问有什么办法加速或者让他匀速运行吗？
※ 修改:·novicer 于 Sep  9 10:01:02 2021 修改本文·[FROM: 143.198.64.*]
※ 来源:·水木社区 http://www.mysmth.net·[FROM: 49.66.217.*]

修改:novicer FROM 143.198.64.*
FROM 49.66.217.*

整个u2提到循环外，怎么保存u呢？u本身是一个800*900的矩阵，如果要保存，循环之后好像是8192*8192*800*900个元素


【 在 zszqzzzf 的大作中提到: 】
: u2整个提取到循环外，可以用prod，对矩阵求累积。
: 双重循环计算也是慢的。如果你能做到消灭两个for，改成矩阵运算，速度可以提高百倍。
: cellfun翻一翻看看能不能就取代ite1和ite2，……
: ...................
--
FROM 143.198.231.*

如果内存足够大，能不能把后一个矩阵变成（8192*8192） * （800*900）的二维矩阵，从而不用for循环？


【 在 shaolimin 的大作中提到: 】
: 2维变1维是为了减少一个循环。
: 循环里面的i确实可以提到循环外。
:
--
FROM 49.66.217.*

基本搞定了，还有个小问题。之前for循环里有一句判断语句，if b(ite)~=0，也就是:
for ite1 = 1:size1*size2  
    if b(ite)~=0
       sum_u2 = sum_u2 + sqrt(1+((x2d-xsource2d(ite1))./zin).^2 +((y2d-ysource2d(ite1))./zin).^2);
    end
end

现在改为：
b_temp = repmat(b,[size1,1 size2]);
b_temp2 = permute(b_temp,[1 3 2]);
ind=find(b_temp2~=0);
[i1,i2,i3]=ind2sub(size(b_temp2),ind);
。。。
sum_u2 = ... x2d_temp(i1,i2,i3) ...
出现错误提示：超出了程序允许的最大变量值。请问这是什么问题？非常感谢！



【 在 shaolimin 的大作中提到: 】
: 那就对了，说明除去循环是成功的。你再检查一下其他语句。
:
--
FROM 192.155.87.*