声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5534|回复: 10

[综合讨论] 用quadl算积分时的积分限的取值

[复制链接]
发表于 2008-6-12 20:51 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
被积函数为:besselh(0,k*l*x)
积分区间是:(0,inf),我近似取inf=20;
可是由于被积函数随自变量的增大而迅速衰减,因此本来inf不用取太大,且取大和取小的差异应该很小,但是计算结果取显示当inf=40是inf=20时的10000倍,不知道是怎么回事?请教各位高手。
计算过程中出现以下warning:
Warning: Maximum function count exceeded; singularity likely.。

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2008-6-12 21:10 | 显示全部楼层
积分奇异,对besselh积分的时候刨去零点
 楼主| 发表于 2008-6-14 10:12 | 显示全部楼层

回复 2楼 的帖子

是刚用matlab编程的,有很多不懂,谢谢指导。
再请教一下:用quad和quadl算同一个积分时,为什么结果有不同?
另外:交换quadl的积分上下限为什么结果不是互为相反数?而是两种结果相差很多。(结果为复数)
谢谢
发表于 2008-6-14 12:42 | 显示全部楼层

回复 3楼 的帖子

quadl精度应该更高点

另外,被积函数奇异了,结果当然不准,后面的那些问题也就不是问题了
你可以用简单的被积函数,验证你那2个问题
 楼主| 发表于 2008-6-14 14:12 | 显示全部楼层

回复 4楼 的帖子

我试了一下:
我的被积函数是:
function y=myfun(x)
rr =0.40000000000000;
k =1.84799567858223;
hr=besselh(0,-k*rr*cosh(x));
y=hr.*(-cosh(x));
分别有quad和quadl计算得到的结果为:
>> a=quadl('myfun',0,10)
Warning: Maximum function count exceeded; singularity likely.
> In quadl at 98

a =

-89.3828 +33.0058i

>> a=quad('myfun',0,10)
Warning: Maximum function count exceeded; singularity likely.
> In quad at 92

a =

  91.7643 + 9.4596i

既然差不多,只是一个比一个精度差点,为什么结果会相差这么多呢?
如果没有办法让warning不出现,那么哪个计算结果更可信一点呢?谢谢各位高手了。
发表于 2008-6-14 14:49 | 显示全部楼层
发表于 2008-6-14 15:24 | 显示全部楼层
既然差不多,只是一个比一个精度差点,为什么结果会相差这么多呢?

那是因为被积函数奇异了

参考2楼的做法
 楼主| 发表于 2008-6-14 19:15 | 显示全部楼层

回复 7楼 的帖子

>> a=quadl('myfun',0.001,10)
Warning: Maximum function count exceeded; singularity likely.
> In quadl at 98

因为besselh在综量等于零时是无穷大,综量为零是奇点故挖去‘0’改为‘0.001’可是结果还是有问题,再次请教各位。
附计算结果:
a =

-81.3031 +45.3588i

>> a=quad('myfun',0.001,10)
Warning: Maximum function count exceeded; singularity likely.
> In quad at 92

a =

  91.0831 -16.0628i
发表于 2008-6-27 10:14 | 显示全部楼层
那试用quadgk函数什么样?
发表于 2010-1-27 12:59 | 显示全部楼层
原帖由 happy 于 2008-6-12 21:10 发表
积分奇异,对besselh积分的时候刨去零点


这个具体怎么操作啊?
发表于 2010-1-28 00:50 | 显示全部楼层

回复 8楼 Caroli 的帖子

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 18:49 , Processed in 0.062525 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表