與奇異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;