推 kaige:解決了~~非常感謝:) 11/15 12:57
※ 引述《kaige (kaige)》之銘言：
: 請問...
: 想要求解下列非線性方程
: 試過fsolve 但常出現方程式解不出來或者超過最大求解次數的錯誤訊息
: 方程式:
: y1-x1-x4*x3=0
: y2-x2+2*x4*x2=0
: y3-x3-x4*x1=0
: x2^2-x1*x3=0
: 有什麼改善fsolve的辦法嗎?
: 或者其他解法??
: 懇請解答
可以直接用solve求解析解
y1~y3我隨便帶值
y1=rand;
y2=rand;
y3=rand;
syms x1 x2 x3 x4
eq1=y1-x1-x4*x3;
eq2=y2-x2+2*x4*x2;
eq3=y3-x3-x4*x1;
eq4=x2^2-x1*x3;
[x1,x2,x3,x4]=solve(eq1,eq2,eq3,eq4);
x1=eval(x1)
x2=eval(x2)
x3=eval(x3)
x4=eval(x4)
有4組解
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.167.168.239