24小时热门版块排行榜    

查看: 2416  |  回复: 6

rolita

银虫 (小有名气)

[求助] 数值积分方法求助

各位虫友们,谁做过无穷积分区间的数值模拟吗?我用的变步长辛普森积分法,但是收敛的太慢了,用高斯拉盖尔算出来结果不对,请问有没有收敛快点的计算无穷积分区间或者积分区间很大的数值方法?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Nothingiseasyintheworld,justkeepingrightattitudealltimes.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
各种数值积分方法都试试吧,比如龙贝格法。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-07-18 16:32:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rolita

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 月只蓝 at 2013-07-18 16:32:50
各种数值积分方法都试试吧,比如龙贝格法。

现在主要不是积分方法的问题,因为积分上限太大导致用各种积分方法判断精度时for循环里积分区间以2的N次方细分,这样计算区间和计算量也是如此增长,这样的话就算不出来。现在就想问问有没有不用划分区间的近似计算大积分区间或无穷积分区间的方法 。。求各位虫友们帮帮忙啊
Nothingiseasyintheworld,justkeepingrightattitudealltimes.
3楼2013-07-19 09:00:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★
csgt0: 金币+1, 谢谢 2013-07-19 14:53:17
rolita: 金币+2, 有帮助, 用其它方法做出来了,很谢谢你的热心帮忙,希望以后有机会能多多交流! 2013-07-22 13:52:29
引用回帖:
3楼: Originally posted by rolita at 2013-07-19 09:00:52
现在主要不是积分方法的问题,因为积分上限太大导致用各种积分方法判断精度时for循环里积分区间以2的N次方细分,这样计算区间和计算量也是如此增长,这样的话就算不出来。现在就想问问有没有不用划分区间的近似计算 ...

可不可以先求得原函数,或者对待积分函数作泰勒展开、级数趋近等类似的近似解析表达,展开的项数多一些,截取至阶很高的无穷小,然后不定积分?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2013-07-19 09:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rolita

银虫 (小有名气)

引用回帖:
4楼: Originally posted by 月只蓝 at 2013-07-19 09:08:14
可不可以先求得原函数,或者对待积分函数作泰勒展开、级数趋近等类似的近似解析表达,展开的项数多一些,截取至阶很高的无穷小,然后不定积分?...

被积函数非常复杂,所以不能求得原函数。解析表达式做出来了,现在就要的是数值解,所以比较头疼。被积函数里头有很多以指数表示的复数,如果积分变量变大的话就会变成高震荡函数,之前用了反比的方法把积分区间化成0到1/a的形式还是不行 。。
Nothingiseasyintheworld,justkeepingrightattitudealltimes.
5楼2013-07-22 09:35:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★
csgt0: 金币+2, 应助指数+1, 多谢指导 2013-07-23 11:40:31
1 被积函数 f(x) 在无穷远处是否为震荡性(変符号),还是同号衰减型?
2 有没有估计值, 比如 | f(x) | <1/ x^n 或  | f(x) | < exp(-a*x)之类?
3 先把函数画一下,实部虚部分别画一下,了解函数趋势
4 如果是震荡型,需要小心,有大数相减的误差累计,可能需要增加有效位数,你需要什么精度(5位,10位有效位数)?
6楼2013-07-22 10:01:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rolita

银虫 (小有名气)

引用回帖:
6楼: Originally posted by pippi6 at 2013-07-22 10:01:07
1 被积函数 f(x) 在无穷远处是否为震荡性(変符号),还是同号衰减型?
2 有没有估计值, 比如 | f(x) | <1/ x^n 或  | f(x) | < exp(-a*x)之类?
3 先把函数画一下,实部虚部分别画一下,了解函数趋势
4  ...

在无穷远处应该是震荡型的,因为其欧拉表示形式里指数的i后面所带的数是与积分限一致的。然后想请问下如何估值?增加有效位数是设定成长整型吗?
Nothingiseasyintheworld,justkeepingrightattitudealltimes.
7楼2013-07-22 13:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rolita 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见