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

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

新規登録して質問してみよう
ただいま回答率
87.20%
NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

解決済

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

tallman0202
tallman0202

総合スコア6

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

2回答

0評価

0クリップ

458閲覧

投稿2020/07/07 15:10

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

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

試したコード

python

#ライブラリをインポート 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ページの「注目」タブのフィードに表示されやすくなります。

気になる質問をクリップする

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

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

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

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

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

まだ回答がついていません

会員登録して回答してみよう

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

ただいまの回答率
87.20%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問

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

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。