[運算] 三元二次方程式 求解
參考了上面的例子後仍無法了解該如何進行
希望可以得到一些解答
我是用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)
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 2 篇)：