24小时热门版块排行榜    

查看: 1811  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 二甲胺的使用 20+4 太阳谷 2024-11-14 12/600 2024-11-15 11:45 by luolinfeng
[基金申请] 广东省自然科学基金-面上项目消息有了么? +4 pachang 2024-11-14 6/300 2024-11-15 09:39 by zhoulei7170
[基金申请] 求助 +4 Enenenene 2024-11-15 4/200 2024-11-15 09:36 by 榨菜拌青椒
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+3 猪小耍 2024-11-13 12/600 2024-11-15 08:41 by 北京莱茵润色
[基金申请] 博后资助名单出来了 +11 Shxjjxjkx 2024-11-14 13/650 2024-11-15 08:41 by homojojo
[论文投稿] 找不到审稿人 +5 吉祥如意141913 2024-11-13 5/250 2024-11-15 07:09 by 半生梦君
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[考博] 随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利 +5 Ricoch4t 2024-11-13 14/700 2024-11-15 01:23 by Ricoch4t
[论文投稿] chemical science和advanced science哪个好一点 +5 yly150 2024-11-12 5/250 2024-11-15 01:09 by ca0yan9
[教师之家] 大学老师 +9 考研一路顺风 2024-11-13 11/550 2024-11-14 21:47 by EndNoted
[教师之家] 北大教授何怀宏曾如此描述他的同行 +11 zju2000 2024-11-09 11/550 2024-11-14 18:14 by frks
[精细化工] 同一个反应相同的反应条件,是不是反应结果相差不大? +4 青霉素 2024-11-11 4/200 2024-11-14 17:33 by zinedine
[基金申请] 变性人申请基金咋算性别? +20 leicahfs 2024-11-11 20/1000 2024-11-14 15:00 by wanghuawei
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[找工作] 咨询一下江西的高校待遇,人文氛围怎么样? +5 akslis2024 2024-11-09 5/250 2024-11-14 13:53 by 啄木鸟、
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 核心期刊那个阶段最难 +3 公共茅斯 2024-11-11 3/150 2024-11-13 15:27 by 13333105938
[论文投稿] 论文接受后,编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见