| 查看: 1074 | 回复: 5 | |||
訾zhou金虫 (正式写手)
|
[交流]
MATLAB做图出现了问题 已有2人参与
|
|
MATLAB小白来请教问题，为啥我下面的这一段程序不能运行，solve函数应该怎么传递变量啊？ Tc=304.2;Pc=7.376;w=0.225;R=8.314; a=0.42748*R^2*Tc^2.5/Pc；b=0.08664*R*Tc/Pc; syms x for T=linspace(0,2500, 2500) for P=linspace(0, 100, 2500) v = solve('R*T/(x-b)-a /(T^0.5)/x/(x+b)-P'); vmax = max(v); vmin = min(v); if (vmax > b && abs(T*(-((R*v^2)/p + ((R*b + a/(2*T^(3/2)))*v)/p - (a*b)/(2*T^(3/2)*p))/(- 3*v^2 + (2*R*T*v)/p + (b^2*p - a/T^(1/2) + R*T*b)/p))/vmax -1)<0.0001)||(v_min>b &&abs(T*(-((R*v^2)/p + ((R*b + a/(2*T^(3/2)))*v)/p - (a*b)/(2*T^(3/2)*p))/(- 3*v^2 + (2*R*T*v)/p + (b^2*p - a/T^(1/2) + R*T*b)/p))/vmin-1)<0.0001) end plot(P,T) end end 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
不自信的我
已经有10人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2019-03-28 15:14:18
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1549
- 散金: 50
- 帖子: 428
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
3楼2019-03-28 15:32:37
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1549
- 散金: 50
- 帖子: 428
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
4楼2019-03-28 15:33:11
5楼2019-03-29 13:10:28
6楼2019-03-29 21:14:36