2009-01-15から1日間の記事一覧

行列方程式を解く(2)

昨日の問題設定で行列方程式Ax=bを解いてみる。 numpyモジュールのarrayで行列を作り、solveメソッドがあるのでとりあえずそれを使う。計算速度を計測する。読み込むデータはあらかじめ作成しておく。 ソースは以下の通り。 # (2009/01/15) # solving matrix…