24小时热门版块排行榜    

CyRhmU.jpeg
查看: 196  |  回复: 1
当前主题已经存档。

abandonkkk

铁虫 (小有名气)

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

程序如下: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个未知数 请教怎么修改
我会送上金币的 谢谢!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chxmright

木虫 (正式写手)

★ ★ ★
小木虫(金币+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换到下一行。
一个人静静地看看天~。~
2楼2010-01-27 17:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abandonkkk 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见