版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

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

回复此楼

» 猜你喜欢

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