24小时热门版块排行榜    

查看: 1812  |  回复: 6
【悬赏金币】回答本帖问题,作者rolita将赠送您 3 个金币

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 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 大学老师 +10 考研一路顺风 2024-11-13 12/600 2024-11-15 14:01 by 阿荣喝酒
[精细化工] 同一个反应相同的反应条件,是不是反应结果相差不大? +5 青霉素 2024-11-11 5/250 2024-11-15 13:43 by zyqchem
[硕博家园] 实验进度 +9 麻辣烫。 2024-11-11 14/700 2024-11-15 13:14 by hanweiguo
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[论文投稿] 职称论文 +9 xiaoqi_66 2024-11-09 11/550 2024-11-15 10:39 by 15035121292
[硕博家园] 和导师撕破脸了 +16 圆满123456 2024-11-13 20/1000 2024-11-15 10:09 by 小杰1572
[基金申请] 求助 +4 Enenenene 2024-11-15 4/200 2024-11-15 09:36 by 榨菜拌青椒
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[论文投稿] 二审审一年的佛系编辑 +10 thefoxrain 2024-11-09 15/750 2024-11-12 19:27 by lide966
[论文投稿] 投稿求助 40+3 liukun11 2024-11-12 4/200 2024-11-12 19:13 by liukun11
[硕博家园] 同步辐射 +4 小张要加油努力 2024-11-10 7/350 2024-11-11 13:47 by mpdfwxgui
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见