他モデリング全般について。唐突ですがJuliaを使って信頼性工学の勉強をしていきます。今回参考にした書籍は以下です。Modelicaに手を出してみました。Modelicaを操作してシミュレーションする試みです。APIが開発されために可能になりました。Matrixの記事になります。Matrixの描画にはGadfly.jlパッケージのヒートマップ機能を使用していました。例としては以下のようなものになります。)関係の関数を集めたJuliaパッケージDesignStructureMatrix.jlを公開しました。詳細はREADMEを御覧ください。Matrixを描画する。Matrixの描画にSeaborn.jlを使用してきました。Seaborn.jlはPython環境がないと使用できない点、日本語に対応していない点が不満でしたので、描画をGadfly.jlに乗り換えたいと思います。Matrixを描画する。Matrixのシーケンシング関連の記事です。今回の内容は以下の記事をJuliaで自動化したものになります。可到達行列を並び替えて下三角行列に近い形にします。並べ替える前の可到達行列は以下のような形になります。Julia言語DSM(DesignStructureMatrix)のシーケンシングに必要なDSMの可到達行列をJuliaのプログラムで求めてみます。以下のスライドの説明がわかりやすいと思います。前回の記事の続きです。前回はJuliaでDSMのクラスタリングのみしましたが、今回はそのクラスタリング結果に基づいてDSMを並べ替えるところまでやります。実装したコードは以下のようになりました。using...DSMをJuliaでクラスタリングし、自動で整列させる2019-05-19DSMをJuliaのクラスタリングプログラムで整理する。Matrix)のクラスタリングコードをJuliaで実装してみます。クラスタリングアルゴリズムについては下記記事を参照してください。実際には以下の記事の擬似コードをJuliaで実装し直したものになります。o...DSMをJuliaのクラスタリングプログラムで整理する。Matrix)に応用してみます。以下のDSMを対角線上に非ゼロ要素が集まるように行列演算で並べ替えます。今回はJuliaのパッケージMetis.jlでグラフを分割してみます。以下の隣接行列からなるグラフを2つに分割してみます。前々回Juliaでランダムネットワークを生成しましたが、今回はスケールフリー・ネットワークを生成してみます。ノードの次数分布がべき則で近似できるネットワークのことです。グラフ理論この記事は以下の大著を大いに参考にしています。いつか書評記事を書きたいですね。Matrix)をMatlabのheatmapで描画し、さらにネットワーク可視化を行いましたが、今回はJulia言語でそれをやってみます。wikipediaのDSM例を描画します。Octave今回はOctaveとJuliaの行列演算速度を比較してみます。×n行列とn×1列の積をとります。この計算速度は一般的にO(n2)になるとされています。以下の記事に使用した計算モデルに間違いがあったため修正と、教訓のために記事を残しておきます。上記記事では以下の例題を扱っておりました。Julia言語これまで何回かJuliaでPID制御をする記事を書きましたが、コードの正しさを確認する検定を行っていなかったので今回はそれを確認します。Enhancedを使い、今回はJuliaを実行してみます。前提として、Juliaがシステム環境変数Pathに登録されている必要があります。Octave最近巷ではDockerが流行していることもあり、私も勉強のためシミュレーション環境をDockerに移行してみたいと考えています。導入に向けての記録を備忘録代わりにとっておきます。制御工学この記事中のプログラムコードには誤りがあります。詳細は下記記事を御覧ください。今回は前回のオートチューニングで得られたパラメーターを当てはめてPIDしてみます。制御工学この記事中のプログラムコードには誤りがあります。詳細は下記記事を御覧ください。今回はJuliaでPIDパラメーターのオートチューニングを行ってみます。その12019-04-05JuliaでPID制御を実装する。前回、JuliaでPI制御を実装しましたので、今回はそれに微分ゲインを加えてPID制御にします。離散的であるプログラムでの実装では微分制御は以下の式で実装できます。\end...JuliaでPID制御を実装する。今回はプログラミング言語JuliaでPI制御を実装してみます。練習のため制御関係のパッケージは使用せず、一から実装してみました。例題は下記書籍のものになります。今回はJuliaの制御工学系パッケージ"ControlSystems.jl"で遊んでみたいと思います。下記Matlab用制御工学の参考書から例題を解いてみます。今回も前回に引き続き行列を使用して1階微分の微分方程式を解いてみます。ない場合について考えていきたいと思います。前回までの内容は以下になります。モデリング前回からの続きです。今回は計算にJulia言語を使用しています。この文章自体も一部Weave.jlを使って書いています。使用するパッケージを読み込んでおきます。notebookのように文章とプログラミングの演算結果を同時に掲載できるノートブックを作っていきます。"Weave.jl"と数値計算に強いJuliaを使えば簡単に技術資料を作ることができそうです。モデリング今回はJuliaのパッケージ"JuMP.jl"を使って最適化問題を解いてみます。Juliaの最適化計算パッケージです。特徴は下記のようですね。モデリング今回はJuliaのパッケージ"Modia.jl"を使用してModelicaのような非因果モデリングでシミュレーションしてみます。Modelicaについては以前の記事を参照してください。今回はプログラミング言語Modelicaを使用してプログラミングの練習をしていきたいと思います。Modelicaはオブジェクト指向のプログラミング言語ですが、その中でも物理モデリングに特化しています。他モデリング全般について。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。