作者boffin (新的开始)
看板MATLAB
标题[运算] 三元二次方程式 求解
时间Thu Aug 16 21:26:00 2012
参考了上面的例子後仍无法了解该如何进行
希望可以得到一些解答
我是用fsolve去解 因为2009找不到有关solve的说明
式子在F中 应该只需要前三个就可以求得解了
但同时需要符合下面两个式子
此外所求得的解要为正值且 x(1) < n1
x(2) < n2
x(3) < n3
试了很多初始值代入都得不到合理的答案
请问该如何修改或是用什麽方式来求解
感谢!!
==============================================
[x] = fsolve(@myfun,x0,options)
function F = myfun(x)
n1=7.1638e-04;
n2=4.1505e-09;
n3=5.3156e-12;
P1=3.1316e-02;
P2=9.7300e-03;
P2=7.4428e-03;
R=0.082;
T=298;
K2=1e5;
K3=5.7e3;
Va=1;
F = [ n1-P1*(x(1)/((x(1)+x(2)+x(3))))-x(1);
n2-P2*(x(2)/(x(1)+x(2)+x(3)))-x(2);
n3-P3*(x(3)/(x(1)+x(2)+x(3)))-x(3);
x(2)-(n2-x(2))*R*T/Va*K2*x(1)*18/1000;
x(3)-(n3-x(3))*R*T/Va*K3*x(1)*18/1000];
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 111.243.43.27
※ 编辑: boffin 来自: 111.243.43.27 (08/16 21:27)