matlab中integral2内的嵌套积分答案

【问题标题】：Nested integral within integral2 in matlabmatlab中integral2内的嵌套积分
【发布时间】：2013-08-26 18:34:31
【问题描述】：

我正在尝试对由积分定义的函数进行双积分（使用integral2）。

http://i.imgur.com/gIUsLSw.jpg

这是我目前正在尝试的：

t=linspace(0,1,50);
fun_1= @(v) exp(.071*v)
fun  = @(x,y) exp(0.14*0.00607*integral(@(u)fun_1(u),0,x)).*exp(-(x-y).^2).*exp(0.14*0.00607*integral(@(u)fun_1(u),0,x));
for i=2:length(t)
    for j=i:length(t)
    A(i,j)=integral2(fun,t(i-1),t(i),t(j-1),t(j));
    end
end

我收到错误消息

使用积分时出错（第 86 行）A 和 B 必须是浮点标量。

谁能提供有关如何解决此问题的任何信息。

【问题讨论】：

标签： matlab nested integration integral

【解决方案1】：

给你：

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

可以，但是速度很慢。

【讨论】：