24小时热门版块排行榜    

查看: 2901  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lgycjpcqu

金虫 (正式写手)

[求助] 求助 MATLAB解方程组-fslove

求解x,y,z的数值,我用MATLAB 的fslove函数求解,解出来结果为
>> x=fsolve(@myfun1,[0.8 0.9 0.05],optimset('Display','off'))
x =
   0.800000000000194   0.900000000000932   0.050000000000305
这个结果似乎就是在初始值附件变动了一下,我改变初始值,结果也发生变化,是不是解法有问题,下面是我写的程序
function q=myfun1(p)
K1=10^(-12.43);
K2=10^(19.4);
K3=5.6*10^(-20);
Ccu=0.24;
Coh=0.75;
Csa=0.26;
x(1)=p(1);
x(2)=p(2);
x(3)=p(3);
q(1)=((1/Coh)*(x(1)/((1-x(1)-x(2))*(1-x(1)-3*x(2)-2*x(3)))))-K1;
q(2)=((1/(Csa*Coh^3))*(x(2)/((1-x(2)-x(3))*(1-x(1)-x(2))*(1-x(1)-3*x(2)-2*x(3))^3)))-K2;
q(3)=(Ccu*Coh^2*(1-x(2)-x(3))*(1-x(1)-3*x(2)-2*x(3))^2)-K3;
%x=fsolve(@myfun1,[0.8 0.9 0.05],optimset('Display','off'));
要是该方程解出来就是这个这个结果,就说明我方程有问题,我就去重新考虑数学模型,现在求高手帮忙看看这些方程的求解是否有问题,多谢!或者哪位高手能够有其它方法准确的求解这个方程,还请指点一下。

方程组及参数

[ Last edited by lgycjpcqu on 2012-7-16 at 16:13 ]
回复此楼

» 本帖@通知

» 猜你喜欢

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

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

emspire529

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
dbb627: 金币+1, 欢迎交流 2012-07-17 19:17:50
lgycjpcqu: 金币+5 2012-07-18 09:19:34
这个应该是化学反应动力学里的吧?据我的初步经验不要用这种,应该考虑一下反应引起的刚性问题。可以减小就算步长或者换成刚性求解器。
7楼2012-07-17 09:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

zt_chem

木虫 (正式写手)

非线性方程组求值时,求解结果与初值有关
2楼2012-07-16 17:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgycjpcqu

金虫 (正式写手)

引用回帖:
2楼: Originally posted by zt_chem at 2012-07-16 17:00:58
非线性方程组求值时,求解结果与初值有关

但是每次求解的结果只在初始值附近波动一个极小的数值,这样正常吗?而且一直随初始值改变,如何确定哪一个是方程的解?我水平比较低,还望您给详细解释下,谢谢!
3楼2012-07-16 17:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lgycjpcqu: 金币+5 2012-07-18 09:20:05
第二个式子似乎有问题,没法满足。
4楼2012-07-16 18:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
信息提示
请填处理意见