| 查看: 212 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】新手弱问 一个方程组的求解
|
|||
|
程序如下：function agmd clear all clc x0=[1000,0.002,58 52]; x=fsolve(@NonlinEqs,x0) function p=P(T) p=exp(23.1964-3816.44/(-46.13+(T+273.15))); function f=NonlinEqs(x) Q=x(1); J=x(2); Thm=x(3); Tgf=x(4); Tc=50; Th=60; f(1)=Q-808.4027*(Tgf-Tc); f(2)=Q-1645.88424*(Th-Thm); f(3)=Q/0.38-98.39476*(Thm-Tgf)+2350*J; f(4)=J-23.38876*log((101325-P(Tgf))/(101325-P(Thm)))/(Thm+Tgf+546.3); 运行没有反应 4个方程 4个未知数 请教怎么修改 我会送上金币的 谢谢！！ |
回复此楼» 猜你喜欢
292求调剂
已经有17人回复
-
266分，求材料相关专业调剂
已经有15人回复
-
化工求调剂
已经有4人回复
-
283分材料与化工求调剂
已经有23人回复
-
320分，材料与化工专业，求调剂
已经有23人回复
-
293求调剂
已经有6人回复
-
334求调剂
已经有9人回复
-
esi高被引论文是不是能对中标有所加分和帮助呢
已经有6人回复
-
求调剂机会
已经有5人回复
-
0856，269分求调剂
已经有17人回复
chxmright
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2084.1
- 散金: 10
- 红花: 3
- 帖子: 358
- 在线: 117.7小时
- 虫号: 879714
- 注册: 2009-10-21
- 性别: GG
- 专业: 高等教育学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):谢谢应助！欢迎常来仿真模拟版。 1-27 19:33
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):谢谢应助！欢迎常来仿真模拟版。 1-27 19:33
|
我把你的程序拷在我的Matlab中运行了下，可以运行得到了结果： Optimization terminated: first-order optimality is less than options.TolFun. x = 347.2815 0.0030 59.7890 50.4296 前面那行说明的是循环成功，你要验算就重新把这个计算结果代到非线性方程中，检验一下。 有个小错误，就是clear all后面换行，将clc换到下一行。 |
2楼2010-01-27 17:20:34
