在matlab中通过改变参数函数值加速积分的双循环[重复]

Question

我问过一个关于在循环中包含积分的问题。我收到了答案，但速度慢得令人难以忍受。我想知道是否有人可以让它更快，我想删除循环并从循环中删除函数定义。

l=3;
t=linspace(0,1,365);
fun3= @(v) integral(@(v)exp(.071*v),0,v,'ArrayValued',true);
for i=2:length(t)
  for j=i:length(t)
      xx=t(i);
      yy=t(j);
      fun  = @(x,y) exp(0.14*0.00607*fun3(yy)).*exp(-(x-y).^2/l).*exp(0.14*0.00607*fun3(xx));
      y(i,j)=integral2(fun,t(i-1),t(i),t(j-1),t(j));
  end
end

Answer 1

是的，for 循环在 Matlab 中是出了名的慢。我通常使用函数arrayfun来避免使用它们。

例子：

[i,j] = meshgrid(1:length(t));
y = arrayfun(@(i,j)integral2(fun,t(i-1),t(i),t(j-1),t(j)),i,j)