24小时热门版块排行榜    

查看: 499  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 8/400 2026-08-15 11:53 by KxMI1BYBxWX1
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 11:41 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +6 Ldrop2023 2026-08-13 6/300 2026-08-15 10:09 by tuanggou
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见