Matlab中的多变量积分问题

Question

我只想按如下方式评估积分：

enter image description here

函数rho_delta定义如下：

p = 3;
delta = 0.1;
A_p = 1/integral(@(s) exp(-s.^2),-p,p);
rho_delta = @(t) A_p*exp(-t.^2/delta^2)/delta;

我尝试使用int 函数来评估被积函数的符号表达式，然后将其转换为匿名函数。最后用integral函数计算，还是不行。

Answer 1

下面的代码可以计算这个积分：

I = quadgk(@(t1) arrayfun(@(t1) quadgk(@(s) rho_delta(t1-s),1,1.1).*(E(t1)-quadgk(@(s) rho_delta(t1-s).*E(s),0,1)),t1),0.9,1);