24小时热门版块排行榜    

查看: 1532  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

ju5200

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
可能是显示精度的问题 用format Eng 声明一下试试
2楼2014-04-29 12:40:37
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见