| 查看: 1320 | 回复: 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客户端 |
回复此楼» 猜你喜欢
时间戳今天，20号变了
已经有5人回复
-
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2019-03-28 15:14:18
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1569
- 散金: 50
- 帖子: 432
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
3楼2019-03-28 15:32:37
訾zhou
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1569
- 散金: 50
- 帖子: 432
- 在线: 38.2小时
- 虫号: 12877988
- 注册: 2018-11-16
- 性别: GG
- 专业: 化工系统工程
4楼2019-03-28 15:33:11
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
5楼2019-03-29 13:10:28
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
6楼2019-03-29 21:14:36
40