24小时热门版块排行榜    

查看: 478  |  回复: 2

chaoweihui

新虫 (初入文坛)

[求助] 新手求助非线性方程组解

我今天试着变了一个程序解非线性方程组,可是计算结果不对,初值我给的是正确解,请高手帮我指导一下
function xfsolve
clear all;clc

x0=[1 0.8 1.5 0.3 5 6666 1000 2345]';x1=fsolve(@abc,x0)

function f=abc(x)
R=8.314;T=1000;
f(1)=19.3+R*T*log(x(1)/(x(1)+x(2)+x(3)+x(4)+x(5)))+x(6)+4*x(7);  %x(1)_x(5)五个组分,x(6)=λc,x(7)=λh,x(8)=λo
f(2)=-192.6+R*T*log(x(2)/(x(1)+x(2)+x(3)+x(4)+x(5)))+x(8)+2*x(7);
f(3)=-200.6+R*T*log(x(3)/(x(1)+x(2)+x(3)+x(4)+x(5)))+x(6)+x(8);
f(4)=-395.9+R*T*log(x(4)/(x(1)+x(2)+x(3)+x(4)+x(5)))+x(6)+2*x(8);
f(5)=R*T*log(x(5)/(x(1)+x(2)+x(3)+x(4)+x(5)))+2*x(7);
f(6)=x(1)+x(4)+x(3)-2;
f(7)=4*x(1)+2*x(2)+2*x(5)-14;
f(8)=x(3)+2*x(4)+x(5)-3;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

zhg19

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
chaoweihui: 金币+1, 有帮助 2013-10-25 10:41:41
计算结果不对主要和初值设定以及迭代次数限定有关,LZ可尝试:
1、更换初值,使之更接近正确值。
2、增大函数赋值的最大次数options.MaxFunEvals和迭代次数上限options.MaxIter
2楼2013-10-19 00:15:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


xzhdty: 金币+1, 谢谢参与 2013-10-27 08:25:53
下面结果参考下:
x1: 2.00001110466504
x2: 1.32537151762595E-13
x3: 2.46246539795621E-13
x4: 1.08824143288337E-26
x5: 3.00001586674963
x6: -895.270393284408
x7: 2123.50255603622
x8: 255852.51453193
3楼2013-10-26 23:21:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chaoweihui 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 enl5e5vjt5 2026-06-12 3/150 2026-06-13 19:53 by 1js5n6upxm
[公派出国] +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[公派出国] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:33 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[公派出国] 售SCI文章,我:8.O.551O.5.4,科目齐全,可十急 +3 i7b870kxb0 2026-06-12 4/200 2026-06-13 16:07 by 999yhm2znn
[教师之家] 售SCI文章,我:8.O.551O.5.4,科目齐全,可十急 +3 i7b870kxb0 2026-06-12 3/150 2026-06-13 15:50 by 999yhm2znn
[考博] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 uaz7enpl3v 2026-06-12 4/200 2026-06-13 15:25 by 999yhm2znn
[硕博家园] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 iatmx9xeph 2026-06-12 4/200 2026-06-13 13:22 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[教师之家] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 uaz7enpl3v 2026-06-12 4/200 2026-06-13 09:59 by ZYNi吖
[公派出国] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 uaz7enpl3v 2026-06-12 3/150 2026-06-13 09:22 by t4hz8376wy
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[硕博家园] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 iatmx9xeph 2026-06-12 3/150 2026-06-13 05:27 by fjz4noejhx
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +6 yexuqing 2026-06-07 7/350 2026-06-12 16:46 by wwwkkk83
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
信息提示
请填处理意见