本人小白,因為論文的原因要計算一個比較復雜的方程組如下:
0.073796 = x1*(0.012686*y2 + 0.013109*y3 + 0.012871*y4 + 0.022445*y5)
0.601394 = x2*(0.012686*y1 + 0.02059*y3 + 0.034696*y4 + 0.011424*y5)
0.239644 = x3*(0.013109*y1 + 0.02059*y2 + 0.029915*y4 + 0.009379*y5)
0.075492 = x4*(0.012871*y1 + 0.034696*y2 + 0.029915*y3 + 0.011717*y5)
0.009675 = x5*(0.022445*y1 + 0.011424*y2 + 0.009379*y3 + 0.011717*y4)
0.073796 = y1*(0.012686*x2 + 0.013109*x3 + 0.012871*x4 + 0.022445*x5)
0.601394 = y2*(0.012686*x1 + 0.02059*x3 + 0.034696*x4 + 0.011424*x5)
0.239644 = y3*(0.013109*x1 + 0.02059*x2 + 0.029915*x4 + 0.009379*x5)
0.075492 = y4*(0.012871*x1 + 0.034696*x2 + 0.029915*x3 + 0.011717*x5)
0.009675 = y5*(0.022445*x1 + 0.011424*x2 + 0.009379*x3 + 0.011717*x4)
想求出所有的x和y,我用的下面這個方法

但是程式一直沒有反應,也不報錯,是我哪里寫錯了么

如果有知道應該怎么算的,麻煩指點一二,謝謝謝謝
uj5u.com熱心網友回復:

貼出代碼來吧
uj5u.com熱心網友回復:
import sympy as sy
x1,x2,x3,x4,x5,y1,y2,y3,y4,y5 = sy.symbols('x1,x2,x3,x4,x5,y1,y2,y3,y4,y5')
result = sy.nonlinsolve([x1*(0.012686*y2 + 0.013109*y3 + 0.012871*y4 + 0.022445*y5)-0.073796,
x2*(0.012686*y1 + 0.02059*y3 + 0.034696*y4 + 0.011424*y5)-0.601394,
x3*(0.013109*y1 + 0.02059*y2 + 0.029915*y4 + 0.009379*y5)-0.239644,
x3*(0.013109*y1 + 0.02059*y2 + 0.029915*y4 + 0.009379*y5)-0.239644,
x4*(0.012871*y1 + 0.034696*y2 + 0.029915*y3 + 0.011717*y5)-0.075492,
x5*(0.022445*y1 + 0.011424*y2 + 0.009379*y3 + 0.011717*y4)-0.009675,
y1*(0.012686*x2 + 0.013109*x3 + 0.012871*x4 + 0.022445*x5)-0.073796,
y2*(0.012686*x1 + 0.02059*x3 + 0.034696*x4 + 0.011424*x5)-0.601394,
y3*(0.013109*x1 + 0.02059*x2 + 0.029915*x4 + 0.009379*x5)-0.239644,
y4*(0.012871*x1 + 0.034696*x2 + 0.029915*x3 + 0.011717*x5)-0.075492,
y5*(0.022445*x1 + 0.011424*x2 + 0.009379*x3 + 0.011717*x4)-0.009675],[x1,x2,x3,x4,x5,y1,y2,y3,y4,y5])
print(result)
uj5u.com熱心網友回復:
有解嗎,不是無解的題吧轉載請註明出處,本文鏈接:https://www.uj5u.com/qita/54181.html
