質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

88.93%

【python】numpyとscipyで変数を含んだ行列式を計算したい

解決済

回答 2

投稿

  • 評価
  • クリップ 0
  • VIEW 220

tallman0202

score 6

numpyとscipyで変数を含んだ行列式を作り、その変数を求めるような計算をしたいのですがうまくいきません。
例えば、
[x, 3]
[4, -1]
のような行列Mの行列式det[M]が0に等しいとすると
det[M]=0
x*(-1)-3*4=0
x=-12
のようにxを求めることができます。
これを実施したいのですがやり方がわかりません。

ご教示いただけますと幸いです。
宜しくお願い致します。

試したコード

#ライブラリをインポート
import numpy as np
import numpy.linalg as LA
import scipy as sc
from scipy import optimize

#行列式det[M]を定義
def M(x):
    M = np.array([[x, 3],[4, -1]])
    return LA.det(M)

#def[M]=0として解く
print(optimize.fsolve(M, 0))

上記コードのエラーメッセージ

TypeError                                 Traceback (most recent call last)
<ipython-input-37-4d75ddba5d6c> in <module>
      8     return LA.det(M)
      9 
---> 10 print(optimize.fsolve(M, 0))

~\Anaconda3\lib\site-packages\scipy\optimize\minpack.py in fsolve(func, x0, args, fprime, full_output, col_deriv, xtol, maxfev, band, epsfcn, factor, diag)
    145                'diag': diag}
    146 
--> 147     res = _root_hybr(func, x0, args, jac=fprime, **options)
    148     if full_output:
    149         x = res['x']

~\Anaconda3\lib\site-packages\scipy\optimize\minpack.py in _root_hybr(func, x0, args, jac, col_deriv, xtol, maxfev, band, eps, factor, diag, **unknown_options)
    211     if not isinstance(args, tuple):
    212         args = (args,)
--> 213     shape, dtype = _check_func('fsolve', 'func', func, x0, args, n, (n,))
    214     if epsfcn is None:
    215         epsfcn = finfo(dtype).eps

~\Anaconda3\lib\site-packages\scipy\optimize\minpack.py in _check_func(checker, argname, thefunc, x0, args, numinputs, output_shape)
     24 def _check_func(checker, argname, thefunc, x0, args, numinputs,
     25                 output_shape=None):
---> 26     res = atleast_1d(thefunc(*((x0[:numinputs],) + args)))
     27     if (output_shape is not None) and (shape(res) != output_shape):
     28         if (output_shape[0] != 1):

<ipython-input-37-4d75ddba5d6c> in M(x)
      6 def M(x):
      7     M = np.array([[x, 3],[4, -1]])
----> 8     return LA.det(M)
      9 
     10 print(optimize.fsolve(M, 0))

~\Anaconda3\lib\site-packages\numpy\linalg\linalg.py in det(a)
   2091     t, result_t = _commonType(a)
   2092     signature = 'D->D' if isComplexType(t) else 'd->d'
-> 2093     r = _umath_linalg.det(a, signature=signature)
   2094     r = r.astype(result_t, copy=False)
   2095     return r

TypeError: No loop matching the specified signature and casting
was found for ufunc det
  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 2

checkベストアンサー

+1

xをx[0]としてください。

#ライブラリをインポート
import numpy as np
import numpy.linalg as LA
import scipy as sc
from scipy import optimize

#行列式det[M]を定義
def M(x):
    M = np.array([[x[0], 3],[4, -1]]) #ここを変更
    return LA.det(M)

#def[M]=0として解く
print(optimize.fsolve(M, 0))
[-12.]

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/07/08 09:34

    ご回答いただき、有難うございます。
    0という初期値を与えてあげるのが重要だったのでしょうか。
    大変助かりました、有難うございました。

    キャンセル

  • 2020/07/08 14:05 編集

    xは単なる実数値で与えられるのではなく、配列(array)で与えられるためです。
    これは、optimize.fsolveが1変数関数だけでなく、多変数関数に対応しているための仕様です。
    今回は1変数なのでx[0]と一個だけ要素をとってやります。
    x[0]の0は初期値ではなく、添字です。fsolve(M,0)の0が初期値になります。

    キャンセル

  • 2020/07/08 17:22

    添え字だったのですね、理解しました。
    ご丁寧に有難うございました。

    キャンセル

+1

numpy, scipy は数値計算ライブラリなので、計算対象はすべて具体的な数値が与えられている必要があります。
記号計算を行いたい場合は sympy を使いましょう。

3.2. Sympy : Python での代数計算 — Scipy lecture notes

import sympy as sp

x = sp.Symbol("x")
M = sp.Matrix([[x, 3], [4, -1]])
det = M.det()
print(det)  # -x - 12

追記

すいません。質問の意図を勘違いしてました。
行列式を求めるというよりは、方程式 det(M) = 0 を解くということだったのですね。

import sympy as sp

x = sp.Symbol("x")
M = sp.Matrix([[x, 3], [4, -1]])


ret = sp.solve(M.det())


print(ret)  # -12

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/07/08 09:31

    ご回答いただき有難うございます。
    追記までいただき、助かります。
    実はこの後scipyの特殊関数を組み込んだ行列式に変更したかったため、sympyではなくscipyで計算したいという質問でございました。
    質問が曖昧で申し訳ありません。

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 88.93%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る