与奇异MATLAB

Question

我想如下

我知道如何写在MATLAB二重积分来执行功能的双重积分二重积分。正如你所看到的，这个函数不容易计算，因为变量是分母。假设我称之为integrad fun1。

我写了下面的代码来解决它

gamma2min=@(gamma1)gamma1; 

I= integral2(fun1,0,+inf,gamma2min,+inf); 

我收到以下错误

最小步长接近X = 0.160784到达。可能存在奇点，或者容差对于这个问题可能太严格。

我也明白问题是什么，但我不知道如何解决它。

感谢所有帮助

Answer 1

Matlab的文档中说，“integral2执行最好的时候，奇点积分边界上” [src。因此，由于您的积分区域是在伽马空间简单的矩形，我打破它分为两个区域对奇点（gamma1 == gamma2）和金额：

I1 = integral2(fun1, 0  , gamma2min , gamma2min , +Inf); 
I2 = integral2(fun1,gamma2min , +Inf , gamma2min , +Inf); 
I = I1 + I2;