| 查看: 1223 | 回复: 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) |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有169人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有3人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有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 (硕士)
- 金币: 17808.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 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
