| 查看: 1281 | 回复: 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) |
回复此楼» 猜你喜欢
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有4人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
独孤神宇
版主 (知名作家)
- 应助: 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 (硕士)
- 金币: 17827.5
- 红花: 53
- 帖子: 2880
- 在线: 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