関数step3では共有メモリー2をnumpyの行列だと思って行列ノルムを算出している importmultiprocessingimportmultiprocessing.sharedctypesimporttimeimportnumpyasnpbyte_buf_size=4096float_buf_size=4096defstep1(out_buf,buf1_ready):f=open("/dev/urandom","rb",buffering=True)# /dev/randomにすると遅くな...
データセット構築の際にNumpyの多次元配列(主に2次元)に識別子や別の配列の情報を付加するしたいときにどのような書き方をすればいいか忘れないための自分用メモ。 importnumpyasnpA=np.array([np.arange(0,6),np.arange(10,16),np.arange(20,26)])B=np.array([100,200,300])C=np.array([np...
decomposition import PCA if __name__ == '__main__': # 行列シートMBook1を探す。このシートには、200バンドのハイパースペクトルデータ用に200の行列オブジェクトが画像スタックとして含まれる msX = op.find_sheet('m', 'MBook1') # 行列シートから200バンドのデータを3Dnumpy配列...
2D numpy arrayのインデックスは数学の行列と同じ定義なので、data[0,0]は中身を表示させた時の左上にあり、[i,j]のiの正方向は下でjの正方向は右です。ところで、dataの可視化の際は慣習的なx軸とy軸の正方向とiとjの正方向を合わせるために、data[0,0]が左下にくるように設定するのがよい...
下記画像で、初見でつまづきそうなポイントに注釈を入れてみました。特に式展開においては、「確率変数ベクトルX、Yに対し、それらが独立なら、AX+BYの分散共分散行列はAΣXAT+BΣYBTになる」という内容がよく出てくるので、おさえておくと良いと思います。(Σは分散共分散行列の意味です。V[...
例えば値を除いて分析した後にこの最大値を取るものを個別に見てみたいとか思っても、消している要素があるのでそれが元のデータ上では何番目になるのかわかりません。 Masked Array そんな時に便利なのがmasked arrayです。masked arrayはNumPyに元から入っている機能です。masked arrayは普通の...
となる。import numpy as np import matplotlib.pyplot as plt def make_ma1(mu, theta1, sigma, t=100): y = [0] e = [np.random.normal(scale=sigma)] for _ in range(t): e.append(np.random.normal(scale=sigma)) y.append(mu+e[-1]+theta1*e[-2]) return y Mu = [0, 2, -2,...
ここで、定義からβは負なので、よりエネルギーが低い軌道はψ1であり、規格化条件c12+c22=1を満たすために12をつけています。 解き方2) numpyを使って import numpy as np a, b = 0, 1 H = np.array([ [a, b], [b, a] ]) ee, cc = np.linalg.eigh(H) print(ee) print(cc.rou...
Pythonで並列処理・分散処理などを簡単に扱ってくれる。 Pythonでよく使われるライブラリとかなり近いインターフェイスを提供している(NumPy、Pandas、Scikit-Learnを中心に、他にもTensorFlow・XGBoostなども)。 必要な場合には、PySparkなんかでやるような、複数のマシンを用意して分散処理でスケ...
import pandas as pd import numpy as np PyCaretとは データの前処理および機械学習モデルのトレーニングを自動化し、ローコード環境でデプロイ可能なライブラリです。 https://pycaret.org/ なお、インストールはpipコマンド一発です。とても簡単。。 pip install pycaret 概要やパイプライン...