| 查看: 1377 | 回复: 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客户端 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
5楼2019-03-29 13:10:28
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2019-03-28 15:14:18
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1579
- 散金: 50
- 帖子: 434
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
3楼2019-03-28 15:32:37
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1579
- 散金: 50
- 帖子: 434
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
4楼2019-03-28 15:33:11
