matplotlibを使って昨日書いた連成振動をグラフ化してみた。 #!/usr/bin/python # -*- coding:utf-8 -*- import numpy as np from scipy import stats import matplotlib.pyplot as plt #-- initialization m = 1.0 k = 1.0 kappa = 1.0 dt = 0.01 omega1 = …

行列計算と正規分布が使えるようになったので、バネマスの連成振動を状態空間モデルで解いてみることにした。・2質点系 ・減衰、外力は働かないが、観測時に観測誤差が入るとするmatplotlibとかで出力できるとかっこいいのだが、とりあえずはファイル出力す…

SciPyを使って正規乱数N(mu,sigma)を発生させるには、SciPyのstatsを利用する。 #!/usr/bin/python # -*- coding:utf-8 -*- >>>import numpy as np >>>from scipy import stats #-- 正規分布N(0,1) >>>mu = 0 >>>sigma = 1 >>>print stats.norm.rvs(0,1) [[ …

pythonのnumpyの行列で、等差級数な行列(ベクトル)の作り方メモ。1.等間隔(間隔指定)なベクトル #!/usr/bin/python # -*- coding: utf-8 -*- >>>import numpy as np >>>print np.arange(2,3,0.1) [ 2. 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9] 2.0から3.0まで…