24小时热门版块排行榜

返回列表

当前主题已经存档。

abandonkkk

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 26.3
帖子: 119
在线: 25.8小时
虫号: 428985
注册: 2007-08-04
专业: 生物化工与食品化工

[交流] 【求助】新手弱问一个方程组的求解

程序如下：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个未知数请教怎么修改
我会送上金币的谢谢！！

回复此楼

» 猜你喜欢

实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
假如你的研究生提出不合理要求已经有6人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复
26申博（荧光探针方向，有机合成）已经有4人回复

1楼 2010-01-20 18:34:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chxmright

木虫 (正式写手)

应助: 3 (幼儿园)
金币: 2084.1
散金: 10
红花: 3
帖子: 358
在线: 117.7小时
虫号: 879714
注册: 2009-10-21
性别: GG
专业: 高等教育学

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):谢谢应助！欢迎常来仿真模拟版。 1-27 19:33

引用回帖:

Originally posted by abandonkkk at 2010-01-20 18:34:30:
程序如下：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) ...

我把你的程序拷在我的Matlab中运行了下，可以运行得到了结果：

Optimization terminated: first-order optimality is less than options.TolFun.

x =

347.2815 0.0030 59.7890 50.4296

前面那行说明的是循环成功，你要验算就重新把这个计算结果代到非线性方程中，检验一下。

有个小错误，就是clear all后面换行，将clc换到下一行。

赞一下(21人)

回复此楼

一个人静静地看看天~。~

2楼2010-01-27 17:20:34