2017-09-22 6 views
1

私はJupyter QtConsole上で動作するPython v.3.6を使用しています。私は、予測を調査スコアと結びつける個人プロジェクトのために、Sympyを使ってデータセット上の線形代数を試みています。線形システムを解くsympy問題

本質的には、N = 14の線形方程式とM = 5の未知数を持つ拡張行列を設定し、システムを解決しようとしています。私の問題は、私は私の拡大行列にsolve_linear_systemコマンドを使用するとき、私は私のコードのいずれかの出力を得ることはありませんということです。

import sympy 
from sympy import * 
from sympy import Matrix, solve_linear_system 
from sympy.abc import x, y, z, u, v 
system = Matrix(((1,1,-1,0,0,1),(1,1,-1,0,0,2),(0,0,-1,0,-1,3), 
(0,0,-1,0,-1,2),(0,0,0,1,0,1),(1,0,1,1,-1,2),(0,0,-1,0,-1,2),(1,0,1,0,0,1), 
(1,1,1,0,1,3),(1,1,1,0,0,2),(-1,1,0,0,-1,3),(1,-1,-1,-1,0,2),(-1,1,1,1,-1,3), 
(0,-1,0,0,0,2))) 
solve_linear_system(system, x, y, z, u, v) 
>> 

誰かが問題とどのような状況を改善するかもしれないかを説明することはできますか?私は他の行列を試しましたが、それはそれらと一緒に使えるように見えるので、私はSympy todoを求めているか、それとも方法ですか? ありがとうございます。

答えて

1

理由は、拡張システムの解決策がないためです。
(おそらくあまりにも多くの制約、あなたは余分方程式の一部を排除することによって、それを緩和しようとすることができ)

あなたはしばらくの間、あなたの行列を凝視する場合は、例えば、互換性のない方程式があることがわかります、行2 0 & 3:(0,0、-1,0、-1,3)、(0,0、-1,0、-1,2)、または行0と1:(1,1、 0,0,1)、(1,1、-1,0,0,2)である。冗長なものもあります。

+0

私は7つの方程式(最後の値が2と等しくないもの)を試してみましたが、それ以上の出力はありません。それ以上の情報損失が心配です。しかし、それが唯一の解決策であれば、私はそれをさらに調べることを検討することができますか? – Jonesn11

+0

行列を少し見ていると、行2と行3のような互換性のない方程式があることがわかります。 '(0,0、-1,0、-1,3)、(0,0 、-1,0、-1,2) 'または行0と1:'(1,1、-1,0,0,1)、(1,1、-1,0,0,2) ' 。冗長なものもあります。 –

+0

正しいことです。しかし、私は2と等しいものを除いたすべての方程式でも試してみましたが、それでも解を見つけることはできませんか? – Jonesn11

関連する問題