時系列データ分析全般 http://elsur.jpn.org/202004MMM/chap5.pdf ランダムウォーク Random Walks Have Never Been Funnier | by Vladimir Ilievski | Towards Data Science 【Python】株価をランダムウォークでシミュレートしてみる - turtlechanのブログ …

作成中ですが、暫定で投稿。 ちなみに、xgboostでは欠損値のnanはそのまま取り扱えるようです。 最適化 試行錯誤的に、動く条件を探索しました。 import xgboost as xgb def objective_xgb(trial): # if trial.number == 0 : learning_rate = trial.suggest_…

XGBoost、LightGBM、アンサンブル学習（Voting regressor）についての情報整理 kaggleはこちら。 www.kaggle.com XGBoost Python: XGBoost を使ってみる - CUBE SUGAR CONTAINER XGBoostハイパーパラメータチューニング Optunaを使ったxgboostの設定方法 - Q…

Optunaを使ってハイパーパラメータを調整します。kaggleにまとめました。 www.kaggle.com Optunaの基本 github.com optuna.readthedocs.io 公式のチュートリアルは下。 import ... # Define an objective function to be minimized. def objective(trial): #…

OpenFOAMでchtMultiRegionFoam+viewFactorの導入について検討したのでそのまとめです。 viewFactorモデルの考え方 形態係数に基づいて輻射を計算します。OpenFOAMではfluid側で新しく輻射の熱流束に関する変数qrを追加して計算します。qrの境界条件は、fluid…

OpenFOAMの結果の可視化・グラフ化についてす。 chtMultiRegionFoamで実施したときにうまくいかなかったので。 $ postProcess -func singleGraph 通常コマンドで実行すると、こんな感じのエラーが出まくります。 Time = 15000 --> FOAM Warning : From funct…

自作の物性の追加します。 前準備（方針の確認） まずはリンクの構成を理解するところから。 $FOAM_SRC/thermophysicalModels/solidSpecie/transport ここに物性。一例として、exponentialSolidTransport.Hのリンク先を確認。 OpenFOAM: API Guide: src/ther…

多項式で近似する場合 numpyの機能で実行できます。 import numpy as np x = np.arange(0, 1000, 1) y_target = xxx ndim = 3 # Coeffs=np.polyfit(x,y_target,3) # plt.plot( x, np.poly1d(Coeffs)(x) ) 一般の関数で近似する場合 scipyのoptimize, curve_f…

chtMultiFegionFoam中で固体（solid）の物性の温度依存性を反映させるのが目的。熱伝導率（transport）、比熱（thermo）の温度依存性を考慮できます。thermophysicalPropertiesが設定ファイル。まずは、thermoTypeで関数の指定。 thermoType { type heSolidT…

chtMultiRegionFoamで輻射を入れる方法についてまとめ。 基本的な考え方 OpenFOAMで輻射を計算する場合、輻射に関する変数（G, IDefault, qr）が追加され、その変数を合わせて解くことになります。 輻射のモデルとしては、P1, viewFactor, FvDOMなどを選択可…

penguinitis.g1.xrea.com 春日氏作成のpythonプログラムblockMesh.pyを使ってみました。ので、簡単なメモ（OpenFOAMも久しぶりなので）。 $ wget http://penguinitis.g1.xrea.com/study/OpenFOAM/blockMeshPython/blockMeshPython-v1912-20200612.tar.gz $ t…

今更ですが。基本的には下サイト通りにやれば、エラーも出ずにうまくいくはずです（私はまずはちゃんとサイトの通りやらなかったので、以下でエラーが出ています。） docs.microsoft.com 管理者権限でpowershellを開いて実行します。 PS C:\WINDOWS\system32…

前回。 matsuo-san.hatenablog.com WSL2になったので、ソフトもxming→VcXsrvに変更しました。 kb.seeck.jp VcXsrvのインストール 以下からソフトウェアをダウンロードします。 sourceforge.net ダウンロードしたら、ポチポチと設定してインストール（あんま…

インターネットを探していたら、BiomassGasificationFoamなるものがあるとのことでインストールしようと思ったのですが、、、 https://iopscience.iop.org/article/10.1088/1742-6596/530/1/012015/pdf 注意しよう！！ 論文中のリンクをたどって開こうとする…

OSSのFEMソフトウェアFrontISTRのDocker利用を通じて、Dockerについての理解を深める回です。 Docker imageをインストールする 以下のFrontISTR公式サイトのダウンロードページに行くと、各バージョンごとのFrontISTRが用意されています。 https://www.front…

Docker Hello-world $ docker run hello-world のコマンドの実行で、hello-worldのimageのダウンロード＋コンテナの展開がおこなわれるようです。 $ docker ps を実行しても特に何も見れないですが、hello-worldは実行後に停止させているからのよう。 $ dock…

WSL2をインストールしたので、Dockerのインストールにトライします。 www.docker.com ページからインストーラをダウンロード・実行。"Enable WSL 2 Windows Features"にチェックを入れるのを忘れない。 qiita.com $ docker ps The command 'docker' could no…

実験装置からの出力データがRS-232Cだったので、勉強ベースで情報のまとめ。 シリアル通信とは？ RS-232はシリアル通信の規格とのこと。 シリアル通信の基礎知識 - シリアル通信とは？RS-232C / RS-422 / RS-485の違いはなんだろう？ - | コンテック RS-232 …

試験結果のデノイズのためにフーリエ解析に基づくローパスフィルタをかけます。 まずは、周波数解析。fftpackを使います。fsには周波数を入れます。 import scipy.fftpack as spfft X = spfft.fft(data) amp = [np.sqrt(c.real ** 2 + c.imag ** 2) for c in…

日本語交じりのファイルをnp.loadtxtで読み込みます。 data = np.loadtxt('data.CSV',delimiter=',',skiprows=1, encoding="utf-8_sig") エラーが発生 UnicodeDecodeError: 'cp932' codec can't decode byte 0xef in position 0: illegal multibyte sequence…

久しぶりに操作すると忘れるのでメモです。 $ conda info -e $ conda activate xxx $ pip install yyy 以上