matlab——数值积分integral和integral2

1、integral

    

注意：fun 必须使用数组运算符而不是矩阵运算符。例如，使用 .* (times) 而不是 * (mtimes)。

eg：

  fun = @(x) exp(-x.^2).*log(x).^2;
  q = integral(fun,0,Inf)
  q = 1.9475

含有参数：

  fun = @(x,c) 1./(x.^3-2*x-c);
  q = integral(@(x) fun(x,5),0,2)
  q = -0.4605
参考：https://ww2.mathworks.cn/help/matlab/ref/integral.html#btdd9x5

2、integral2

     

eg：

  fun = @(x,y) 1./( sqrt(x + y) .* (1 + x + y).^2 );

    

注意：仅y的上下限含变量（可以在对非矩形区域积分时将 ymin 和ymax 指定为函数句柄（x 的函数）），x的上下限必须是有限或无限的实标量值。

在对非矩形区域积分时，如果 ymin、ymax（或两者）为函数句柄，则可实现最佳性能和精确度。

如果 ymin、ymax（或两者）为无边界函数，请使用 'iterated' 方法。

     

参考：https://ww2.mathworks.cn/help/matlab/ref/integral2.html