Julia言語
久しぶりの投稿です。今回は非因果モデリングのためのJuliaパッケージModelingToolkit.jlを使用してタンクモデルを作成してみます。 ModelingToolkit.jlについて ModelingToolkitはJuliaネイティブのモデリングツールで、方程式ベースでモデルを記述できます…
今回はJuliaの強力な微分方程式数値計算パッケージDifferentialEquations.jlを使用してRLC回路を解いていきます。 パッケージの詳細は以下です。 diffeq.sciml.ai また関連する記事は以下になります。 otepipi.hatenablog.com otepipi.hatenablog.com RLC回…
今回もJuliaの数式処理パッケージSymbolics.jlの事例紹介になります。 過去の記事は以下になります。 otepipi.hatenablog.com otepipi.hatenablog.com また、記事内容は以下書籍のプログラムを参考にしています。 MATLABとOctaveによる科学技術計算丸善出版A…
今回も前回に引き続き数式処理システムSymbolics.jlの解説です。今回はこのパッケージを用いて簡単な多変数最適化の文章題を解いてみたいと思います。 問題は以下書籍からそのまま抜き出しました。 数理モデリング入門 ―ファイブ・ステップ法― 原著第4版作者…
Juliaの数式処理パッケージSymbolics.jlを試してみます。 github.com ドキュメントは↓です symbolics.juliasymbolics.org パッケージの宣言 using Symbolics 変数の宣言 @variables x \begin{equation} \left[ \begin{array}{c} x \ \end{array} \right] \en…
今回もJuliaとmodelicaの連携記事になります。 関連記事は以下です。 otepipi.hatenablog.com otepipi.hatenablog.com 今回は制御工学の題材としてメジャーな倒立振子をmodelicaでモデリングし、Juliaでアニメーションを行います。 倒立振子のモデル 倒立振…
今回もこれまでに引き続いてJuliaとmodelicaの連携を実践してみます。 今回は「魚の体重」をJuliaとmodelicaを使ってPI制御してみます。 魚の体重の数理モデル 魚の体重の増加量を表す数理モデルとしてフォン・ベルタンフィーモデルが知られています。 ここ…
※今回の内容は以下論文を再現しただけになります https://www.researchgate.net/publication/330809787_OMJulia_An_OpenModelica_API_for_Julia-Modelica_Interactionwww.researchgate.net 前回に引き続きmodelicaとJuliaを繋ぐOMJulia関連の記事になります…
今回は久々にmodeica関連の記事になります。 JuliaでOpen Modelicaを操作するパッケージ「OMJulia」を用いてopen modelicaのモデルのパラメーターを最適化します。 関連する記事は以下になります。 otepipi.hatenablog.com 今回のモデル 今回モデリングする…
2021年初投稿です。今回は婚活アプリのシステムモデリングをしてみましたが、Jupyterで作成したものをそのままブログに貼り付けてみました。こちらのほうが手数少なく記事が作成できそうですが、見栄えがどうでしょうか……？ Jupyterからの投稿手順については…
前回の続きでセルオートマトンについて書いていきます。 otepipi.hatenablog.com 今回はセルオートマトンの分類についてになります。 以下書籍を参考にしています。セルオートマトンに関する和書は数少ないですが、その中でも特に詳細に記述されている用に思…
今回はシステムの離散モデルの1種であるセルオートマトンをJuliaで描画してみたいと思います。 最近、『東京大学工学教程 システム工学』シリーズを読む機会がありました。 東京大学工学教程 システム工学 システム理論I作者:大橋 弘忠,古田 一雄発売日: 201…
今回も書籍『A Biologist's Guide to Mathematical Modeling in Ecology and Evolution』 で学んだ内容を復習します。内容としては、計算モデルの収束判定です。
これまでいくつかの計算モデルを静的にプロットしてきました。 今回はJuliaのパッケージInteract.jlを使って プロットをgifにしてみます。
今回も書籍『A Biologist's Guide to Mathematical Modeling in Ecology and Evolution』 から学習した内容の備忘録になります。 計算モデルを解析する際に使用されるCobweb Plotsを作成してみます。
今回も書籍『A Biologist's Guide to Mathematical Modeling in Ecology and Evolution』 から学習した内容の備忘録になります。 ネズミが生殖により増えたり、銀行に預けたお金が複利で増えたり、人気商品だということが広まるとさらに人気になったり…… 「…
今回はJuliaで分岐図を作成してみたいと思います。 分岐図についてですが例として以下のようなものになります。 これはロジスティック関数の分岐図でパラメーターrを増やしていくと値が発散してカオス減少が確認できるものになります。 ja.wikipedia.org 今…
唐突ですがJuliaを使って信頼性工学の勉強をしていきます。 今回参考にした書籍は以下です。 入門 信頼性工学(第2版)：確率・統計の信頼性への適用作者:福井 泰好発売日: 2016/07/30メディア: 単行本（ソフトカバー） ワイブル分布とは ワイブル分布は1939年…
久々にOpen Modelicaに手を出してみました。 今回はJuliaを使ってOpen Modelicaを操作してシミュレーションする試みです。 今年になってOpen ModlicaのJulia APIが開発されために可能になりました。 www.openmodelica.org 手順 パッケージOMJuliaを使用しま…
久々にDesign Structure Matrixの記事になります。 今までDesign Structure Matrixの描画にはGadfly.jlパッケージのヒートマップ機能を使用していました。 例としては以下のようなものになります。 ただこの描画方法ですとDSMの要素の有無を色でしか表現でき…
Design Structure Matrix ( DSM )関係の関数を集めたJuliaパッケージDesignStructureMatrix.jlを公開しました。 github.com 現時点の内容としては DSMの描画 DSMのシーケンシング DSMのクラスタリング が利用できます。詳細はREADMEを御覧ください。 今回、J…
今まではDesign Structure Matrixの描画にSeaborn.jlを使用してきました。 otepipi.hatenablog.com ただSeaborn.jlはPython環境がないと使用できない点、日本語に対応していない点が不満でしたので、描画をGadfly.jlに乗り換えたいと思います。 Gadfly.jlで…
Design Structure Matrixのシーケンシング関連の記事です。今回の内容は以下の記事をJuliaで自動化したものになります。 otepipi.hatenablog.com 可到達行列を並び替えて下三角行列に近い形にします。 並べ替える前の可到達行列は以下のような形になります。…
DSM(DesignStructureMatrix)のシーケンシングに必要なDSMの可到達行列をJuliaのプログラムで求めてみます。 可到達行列とは 以下のスライドの説明がわかりやすいと思います。 http://ise.me.tut.ac.jp/lectures/decision-making/DSS12-2.pdf DSMの隣接行列A…
前回の記事の続きです。 otepipi.hatenablog.com 前回はJuliaでDSMのクラスタリングのみしましたが、今回はそのクラスタリング結果に基づいてDSMを並べ替えるところまでやります。 実装 実装したコードは以下のようになりました。 using LinearAlgebra using…
今回はDSM(Design Structure Matrix)のクラスタリングコードをJuliaで実装してみます。 クラスタリングアルゴリズムについては下記記事を参照してください。 otepipi.hatenablog.com 実際には以下の記事の擬似コードをJuliaで実装し直したものになります。 o…
前回の記事で行列を並べ替える置換行列を勉強しましたが、これをDSM(Design Structure Matrix)に応用してみます。 otepipi.hatenablog.com 例題 以下のDSMを対角線上に非ゼロ要素が集まるように行列演算で並べ替えます。 また、今回のDSMの隣接行列は以下の…
今回はJuliaのパッケージMetis.jlでグラフを分割してみます。 github.com 以下の隣接行列からなるグラフを２つに分割してみます。 A=[0 1 0 0 1 1 0; 0 0 0 1 0 0 1; 0 1 0 1 0 0 1; 0 1 1 0 1 0 1; 0 0 0 1 0 1 0; 1 0 0 0 1 0 0; 0 1 1 1 0 0 0] Metis.par…
前々回Juliaでランダムネットワークを生成しましたが、今回はスケールフリー・ネットワークを生成してみます。 otepipi.hatenablog.com スケールフリーネットワークとは ノードの次数分布がべき則で近似できるネットワークのことです。 一方、以前紹介したラ…
この記事は以下の大著を大いに参考にしています。 いつか書評記事を書きたいですね。 ネットワーク科学: ひと・もの・ことの関係性をデータから解き明かす新しいアプローチ作者: アルバート・ラズロバラバシ,Albert‐L´aszl´o Barab´asi,池田裕一,井上寛康,谷…
