24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 3981  |  回复: 7

1612524436

新虫 (正式写手)

[求助] 使用MATLAB计算定积分方程组,求系数已有1人参与

已知函数f1=exp[-q/(8.314*8800*t)],f2=exp[-q/(8.314*7800*t)],y1=int(f1,t,0,0.1),y2=int(f1,t,0,0.08),y1=2*y2,求q。请问大神用MATLAB应该如何编写,多谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
sol=fzero(@(q)quadl(exp(-q/(8.314*8800*t)),t,0,0.1)-2*quadl(exp(-q/(8.314*7800*t)),t,0,0.08),1)
数值计算
2楼2018-10-17 10:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

sol=fzero(@(q)quadl(exp(-q./(8.314*8800*t)),t,0,0.1)-2*quadl(exp(-q./(8.314*7800*t)),t,0,0.08),1)
数值计算
3楼2018-10-17 10:27:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1612524436

新虫 (正式写手)

引用回帖:
3楼: Originally posted by 独孤神宇 at 2018-10-17 10:27:29
sol=fzero(@(q)quadl(exp(-q./(8.314*8800*t)),t,0,0.1)-2*quadl(exp(-q./(8.314*7800*t)),t,0,0.08),1)

报错了呀老哥,完整的应该是啥样的呢?
4楼2018-10-17 10:36:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

引用回帖:
4楼: Originally posted by 1612524436 at 2018-10-17 10:36:47
报错了呀老哥,完整的应该是啥样的呢?...

错误信息是什么?我检查一下
数值计算
5楼2018-10-17 11:08:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1612524436

新虫 (正式写手)

引用回帖:
5楼: Originally posted by 独孤神宇 at 2018-10-17 11:08:51
错误信息是什么?我检查一下...

错误使用 fzero (line 287)
FZERO 无法继续,因为用户提供的 function_handle ==>
@(q)quadl(exp(-q./(8.314*8800*t)),t,0,0.1)-2*quadl(exp(-q./(8.314*7800*t)),t,0,0.08)
失败,出现下面的错误。

如果 FUN 为 MATLAB 对象,则它必须具有 feval 方法。
6楼2018-10-17 11:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 1612524436 at 2018-10-17 11:13:17
错误使用 fzero (line 287)
FZERO 无法继续,因为用户提供的 function_handle ==>
@(q)quadl(exp(-q./(8.314*8800*t)),t,0,0.1)-2*quadl(exp(-q./(8.314*7800*t)),t,0,0.08)
失败,出现下面的错误。

...

sol=fzero(@(q)quadl(@(t)(exp(-q./(8.314*8800*t))),0,0.1)-2*quadl(@(t)(exp(-q./(8.314*7800*t))),0,0.08),5000000)
数值计算
7楼2018-10-17 11:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1612524436: 金币+10, ★★★★★最佳答案 2018-10-17 14:48:17
参考这个结果

Parameter q;
Function int(exp(-q/(8.314*8800*t)),t=0,0.1)-2*int(exp(-q/(8.314*7800*t)),t=0,0.08)=0;


====== 结果 ======

迭代数: 21
计算用时(时:分:秒:微秒): 00:00:00:184
计算结束原因: 达到收敛判断标准
优化算法: 通用全局优化算法(UGO1)
目标函数值(最小): 0
q: 3853554.16362925

====== 计算结束 ======
数值计算
8楼2018-10-17 11:26:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1612524436 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 北邮女博导并没有强制学生必须参与项目,说了学生可以不再参与项目 +38 苏东坡二世 2024-04-10 59/2950 2024-04-16 17:05 by dujianchao
[基金申请] 基金和生小孩 +26 Ausy 2024-04-15 26/1300 2024-04-16 14:23 by dudujing
[论文投稿] with efitor 越久是不是越容易拒稿。我的已经一个多月了 +5 lizhengke06 2024-04-14 5/250 2024-04-15 18:33 by jonewore
[考研] 273求调剂 +5 Late婉安 2024-04-15 7/350 2024-04-15 13:01 by Late婉安
[考研] 304求调剂 +14 gsugsysy 2024-04-10 14/700 2024-04-14 22:40 by 永字号
[硕博家园] 双9 调剂求捞一把 +12 随鱼而安. 2024-04-10 15/750 2024-04-14 21:15 by lizh5064
[考研] 314求调剂 +18 阿洛法 2024-04-12 20/1000 2024-04-13 13:40 by kxxmaqq
[考研] 290 生物与医药求调剂 +7 GPX4 2024-04-11 9/450 2024-04-13 09:40 by lincunhui
[考研] 290 生物与医药专硕(086000)求调剂 +5 GPX4 2024-04-12 5/250 2024-04-13 07:48 by coco1981
[考研] 北信科第二批调剂,工科相关 +4 biser 科研人 2024-04-11 6/300 2024-04-12 21:26 by biser 科研人
[考研] 266求调剂 +3 ~酥糖 2024-04-12 15/750 2024-04-12 21:22 by 赵燕高仙兰
[考研] 280求调剂 +3 黑皮冰棒 2024-04-12 6/300 2024-04-12 16:18 by 黑皮冰棒
[考研] 322求调剂 +3 努力进步man 2024-04-12 6/300 2024-04-12 15:24 by 1145075130
[硕博家园] 求购镍包碳化钨,镍包氧化铝,镍包碳化硅 +3 于于于与于 2024-04-11 3/150 2024-04-12 14:52 by 不发光就发疯
[考研] 双9 +5 随鱼而安. 2024-04-11 5/250 2024-04-11 16:40 by 公瑾逍遥
[考研] 282求调剂 +4 王子豪桑 2024-04-10 5/250 2024-04-11 10:41 by 阿杰鲁呀
[论文投稿] 关于转投问题 5+3 xyn@990905 2024-04-09 7/350 2024-04-11 09:10 by bobvan
[考研] 329求调剂 +12 瞧瞧乔儿 2024-04-10 17/850 2024-04-11 00:24 by xuxiaodai
[考研] 350求调剂 +4 Midori0103 2024-04-10 5/250 2024-04-10 22:39 by 15081537629
[考研] 320求调剂 +3 安尧0108 2024-04-09 3/150 2024-04-10 06:38 by MS.
信息提示
请填处理意见