24小时热门版块排行榜    

查看: 1647  |  回复: 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

新虫 (小有名气)

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

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +8 天下熯 2026-02-28 8/400 2026-03-02 00:15 by 暮雨星晴
[考研] 材料学硕318求调剂 +11 February_Feb 2026-03-01 13/650 2026-03-01 23:53 by ccp273206157
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 13/650 2026-03-01 23:32 by L135790
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-01 5/250 2026-03-01 23:31 by L135790
[考研] 292求调剂 +6 yhk_819 2026-02-28 6/300 2026-03-01 23:23 by 向上的胖东
[考研] 0856材料与化工,270求调剂 +6 YXCT 2026-03-01 6/300 2026-03-01 23:21 by 向上的胖东
[考研] 275求调剂 +3 明远求学 2026-03-01 3/150 2026-03-01 22:29 by 刘兵
[考研] 材料类求调剂 +10 wana_kiko 2026-02-28 12/600 2026-03-01 22:10 by 海嵙Y
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 6/300 2026-03-01 21:24 by cgyzqwn
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 12/600 2026-03-01 19:48 by 好好好1233
[考研] 材料学调剂 +9 提神豆沙包 2026-02-28 11/550 2026-03-01 18:15 by ms629
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 化工专硕348,一志愿985求调剂 +5 弗格个 2026-02-28 8/400 2026-03-01 17:25 by sunny81
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
信息提示
请填处理意见