| 查看: 1155 | 回复: 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) |
回复此楼» 猜你喜欢
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有235人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有22人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.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 (硕士)
- 金币: 18111.5
- 红花: 53
- 帖子: 2879
- 在线: 603.7小时
- 虫号: 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