在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)