jiguozhao
2009-06-15, 16:13
各位高手：
大家好！
我遇到这样一个问题：
R=0.5;
g=@(theta,ro,l)(ro.*l.*(R-ro.*cos(theta)))./((l.^2+R.^2+ro.^2-2*R.*ro.*(cos(theta))).^2);
G=triplequad(g,0,2.*pi,0,0.5,0,1)
中间出现了很多warning（Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In triplequad>innerintegral at 67
In quad at 71
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51
Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In triplequad>innerintegral at 67
In quad at 80
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51
Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51）
最后结果G =
0/0
左侧G的value项显示NaN
希望有高人给我指点一下
大家好！
我遇到这样一个问题：
R=0.5;
g=@(theta,ro,l)(ro.*l.*(R-ro.*cos(theta)))./((l.^2+R.^2+ro.^2-2*R.*ro.*(cos(theta))).^2);
G=triplequad(g,0,2.*pi,0,0.5,0,1)
中间出现了很多warning（Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In triplequad>innerintegral at 67
In quad at 71
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51
Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In triplequad>innerintegral at 67
In quad at 80
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51
Warning: Infinite or Not-a-Number function value encountered.
> In quad at 103
In dblquad>innerintegral at 82
In quad at 71
In dblquad at 58
In triplequad at 51）
最后结果G =
0/0
左侧G的value项显示NaN
希望有高人给我指点一下