| 查看: 1117 | 回复: 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客户端 |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1554
- 散金: 50
- 帖子: 429
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
3楼2019-03-28 15:32:37
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2019-03-28 15:14:18
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1554
- 散金: 50
- 帖子: 429
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
4楼2019-03-28 15:33:11
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17803.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
5楼2019-03-29 13:10:28
10