| 查看: 1338 | 回复: 4 | ||
100121114金虫 (小有名气)
|
[求助]
matlab求解物理函数 已有2人参与
|
|
需要用matlab计算，求T下的值，谢谢大佬。 clear; q=1; k=8.6174e-5; BN=0.45; A=110; VR=45; JF=IF/S; IF=20; s=112; S=(112*25.4/1000)^2; f=JF*(BN+k*T/q*ln(JF/(A*T^2)))-A*T^2*e^(-q*BN/(k*T)*VR); ans=solve(f==0,T) |
回复此楼» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有673人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2020-07-09 16:17:00
100121114
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1496.6
- 散金: 30
- 帖子: 66
- 在线: 21小时
- 虫号: 3202917
- 注册: 2014-05-13
- 专业: 半导体科学与信息器件
3楼2020-07-09 16:51:08
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2020-07-09 20:48:10
100121114(独孤神宇代发): 金币+50, 鼓励交流 2020-07-10 08:24:25
感谢参与，应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2020-07-09 20:48:10
100121114(独孤神宇代发): 金币+50, 鼓励交流 2020-07-10 08:24:25
|
题主的问题是没有按软件函数的要求编程。完整的代码为 clear; syms T q=1; k=8.6174e-5; BN=0.45; A=110; VR=45; IF=20; S=(112*25.4/1000)^2; JF=IF/S; f=JF*(BN+k*T/q*log(JF/(A*T^2)))-A*T^2*exp(-q*BN/(k*T)*VR); ans=vpasolve(f==0,[T]) 这里，用vpasolve计算比solve更精确。 计算结果如下。  123.jpg |
4楼2020-07-09 19:07:25
100121114
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1496.6
- 散金: 30
- 帖子: 66
- 在线: 21小时
- 虫号: 3202917
- 注册: 2014-05-13
- 专业: 半导体科学与信息器件
5楼2020-07-09 21:44:58
