24小时热门版块排行榜    

查看: 469  |  回复: 4
当前主题已经存档。

abandonkkk

铁虫 (小有名气)

[交流] 【求助】一个方程组的求解

谢谢

matlab程序如下: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的回帖

haixing2008

荣誉版主 (文坛精英)

★ ★
小木虫(金币+0.5):恭喜抢沙发,给个红包
formleaf(金币+1,VIP+0):谢谢 1-22 09:30
会不会这个出问题呀? x0=[1000,0.002,58 52];
平平淡淡才是真!
2楼2010-01-21 21:54:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abandonkkk

铁虫 (小有名气)

bluesine(金币+0):58和52中间差个逗号吧:D 1-23 20:33
引用回帖:
Originally posted by haixing2008 at 2010-1-21 21:54:
会不会这个出问题呀? x0=[1000,0.002,58 52];

请问怎么改 谢谢
3楼2010-01-22 16:16:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

引用回帖:
Originally posted by abandonkkk at 2010-1-22 16:16:

请问怎么改 谢谢

[...]中第三个数据后面差个逗号
板凳要做十年冷文章不发一个字
4楼2010-01-23 20:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

有结果啊
x =

  1.0e+002 *

  Columns 1 through 2

   3.47281491474341   0.00002986498385

  Columns 3 through 4

   0.59789000050529   0.50429589722392
板凳要做十年冷文章不发一个字
5楼2010-01-23 20:39:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abandonkkk 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见