24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1672  |  回复: 7

kz341

新虫 (小有名气)

[求助] 一个关于 quadgk 函数求积分的问题 已有1人参与

小弟想求解一个半无限震荡积分问题:int(x*m.*exp(-x/t2).*sin(x/t1), 0, inf),其中m, t1, t2 为已知常数。
具体程序如下
clear all
t1=12.2*(10^-15);
t2=32*(10^-15);
m=(t1.^2+t2.^2)/(t1*t2.^2);
F=@(x)x*m.*exp(-x/t2).*sin(x/t1);
fplot(F,[0,1.4*(10^-12)]);
[q,errbnd] = quadgk(F,0,inf)
画出的函数形状无误,但积分后得到的结果与实际结果不符,q和errbnd均为0,实际结果q应接近于1.67*(10^-14)
还请高人帮忙看一下是什么原因

非常感谢!!
回复此楼

» 猜你喜欢

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

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

ju5200

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
可能是显示精度的问题 用format Eng 声明一下试试
2楼2014-04-29 12:40:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kz341

新虫 (小有名气)

引用回帖:
2楼: Originally posted by ju5200 at 2014-04-29 12:40:37
可能是显示精度的问题 用format Eng 声明一下试试

用了之后还是0....
q =

    0.00000000000000e+000


errbnd =

    0.00000000000000e+000
还是非常感谢你的回复。
3楼2014-04-29 16:09:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by kz341 at 2014-04-29 16:09:45
用了之后还是0....
q =

    0.00000000000000e+000


errbnd =

    0.00000000000000e+000
还是非常感谢你的回复。...

一个关于 quadgk 函数求积分的问题

根据帮助文件  这个积分函数是近似计算 并不是精确计算 所以才有你的问题
4楼2014-04-29 19:12:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

32100004

铜虫 (初入文坛)

用用其它函数试试
我用的是Mathematica

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : tmp.png
  • 2014-04-30 13:39:02, 12.55 K
5楼2014-04-30 13:39:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kz341

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 32100004 at 2014-04-30 13:39:03
用用其它函数试试
我用的是Mathematica

非常感谢!我发现之前的函数是因为上限不能设成inf,根据震荡函数设到收敛的位置得出的结果差了1/2,量级是对了
6楼2014-05-06 15:48:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kz341

新虫 (小有名气)

引用回帖:
4楼: Originally posted by ju5200 at 2014-04-29 19:12:00
根据帮助文件  这个积分函数是近似计算 并不是精确计算 所以才有你的问题...

感谢回复!
7楼2014-05-06 15:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shanlihong

银虫 (小有名气)

你用的是matlab什么版本?
8楼2014-11-24 19:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kz341 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +9 我不是只因 2026-04-02 9/450 2026-04-03 00:55 by 小徐0109
[考研] 一志愿北京科技大学材料工程085601,求调剂 +11 cdyw 2026-04-02 11/550 2026-04-02 23:06 by 马儿快快地跑
[考研] 求调剂 +5 朔朔话 2026-04-02 6/300 2026-04-02 22:02 by barlinike
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 求调剂22408 288分 +4 new382 2026-04-02 4/200 2026-04-02 16:07 by 氮气气气
[考研] 0832食品科学与工程学硕282调剂 +4 鱼在水中游a 2026-04-02 7/350 2026-04-02 14:12 by baoball
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666!
[考研] 270求调剂 +8 小杰pp 2026-03-31 10/500 2026-04-02 12:57 by yulian1987
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 070300求调剂306分 +5 26要上岸 2026-03-27 5/250 2026-04-01 11:09 by oooqiao
[考研] 299求调剂 +8 嗯嗯嗯嗯2 2026-03-27 8/400 2026-03-31 18:23 by lizhi8172
[考研] 本2一志愿C9-333分,材料科学与工程,求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 277跪求调剂 +8 1915668 2026-03-27 13/650 2026-03-31 14:58 by 王亮_大连医科大
[考研] 286求调剂 +5 丢掉懒惰 2026-03-27 8/400 2026-03-31 11:27 by Delta2012
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 283求调剂 +3 A child 2026-03-28 3/150 2026-03-28 15:41 by ms629
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
信息提示
请填处理意见