| 查看: 218 | 回复: 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个未知数 请教怎么修改 我会送上金币的 谢谢！！ |
回复此楼» 猜你喜欢
这年头没有找到涵评专家，还有中面上的可能吗
已经有15人回复
-
26/27博士推荐
已经有3人回复
-
找博士生导师
已经有4人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
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
