-
最近の投稿
最近のコメント
Archive
カテゴリー
メタ情報
「未分類」カテゴリーアーカイブ
よく使う VS Code のキーボードショートカット
ターミナルにフォーカス Ctrl + @ ターミナルパネルを閉じる Ctrl + @ + @ エディタにフォーカス Ctrl + 1 エディタのタブ移動 Ctrl + Tab サイドバーを閉じる Ctrl + B &n … 続きを読む
ガウス関数の FWHM を Mathematica の数式処理で求める
g[x_, mu_, sigma_] := 1/(Sqrt[2 *Pi]*sigma)*Exp[-(x – mu)^2/(2*sigma^2)] x = mu で最大値を取るので,その値の 1/2 になる x … 続きを読む
フィッティング Scipy optimize curve_fit
|
1 2 3 4 5 6 7 8 9 10 |
import numpy as np from scipy.optimize import curve_fit param_initial = np.array((10., 100., 0.2, 0.0)) # initial guess param_bounds = ((0.0, -np.inf, 0.1, -np.pi), (np.inf, np.inf, 0.25, np.pi)) # bounds for parameter popt, pcov = curve_fit(fit_func, x, y, p0=param_initial, sigma=err, bounds=param_bounds) print(popt) # optimized parameters print(np.sqrt(np.diag(pcov))) # standard error |
パラメータの値をある範囲内に指定するときは，bounds option を使う． 値の与え方は， bounds = ((min1, max1), (min … 続きを読む
Mathematica で高次の項を無視する
https://stackoverflow.com/questions/21517020/in-mathematica-how-can-i-cut-off-the-high-order-terms-in-a-polyno … 続きを読む
