82563341 发表于 2012-3-21 09:17:55

含奇异点的数值积分

请问一下大家含有奇异点的积分该如何处理呢？

function y = fx1(m,m1,R,r0,t)

FF = @(x)cos(lam(R,m)*x).*cos(lam(R,m1)*x)./(x+r0).^t;
y = quadgk(FF,0,R);

其中t是分母的次数，可取1，2，3

当r0=0，x=0时该积分有奇异。

这是对圆环扩展到整圆，遇到这样一个积分。现在奇异点始终会影响积分结果，请大家帮忙看看，怎么能够较好地解决这个问题，谢谢。
页: [1]
查看完整版本: 含奇异点的数值积分