24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2233  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料考研求调剂总分280 +29 mkjlz1 2026-04-06 33/1650 2026-04-08 20:14 by 逆水乘风
[考研] 生物与医药273求调剂 +17 荔题南墙 2026-04-05 18/900 2026-04-08 19:12 by 我减肥1
[考研] 085600材料与化工301分求调剂院校 +27 刺痛jk 2026-04-06 28/1400 2026-04-08 16:16 by luoyongfeng
[考研] 求调剂 +14 张zic 2026-04-05 15/750 2026-04-08 16:14 by 一只好果子?
[考研] 求调剂,262机械专硕 +4 嗯yyl 2026-04-08 4/200 2026-04-08 16:12 by 一只好果子?
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 材料调剂 +19 一样YWY 2026-04-02 20/1000 2026-04-08 11:02 by 不吃魚的貓
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 一志愿国科大信工所,英二数二408总分293分求调剂 +3 ilcyuan 2026-04-02 4/200 2026-04-06 16:35 by likeihood
[考研] 第一志愿东南大学物理313,有科研竞赛获奖经历,希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 一志愿C9的化学工程(085602) 340分,感觉校内调剂无望,求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 313求调剂 +3 ~微微凉~ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
信息提示
请填处理意见