C++向けの行列計算ライブラリEigenを使って姿勢運動を扱う話． 初回はウォーミングアップということで，Quaternionクラスの使い方の練習をする． はじめに 僕はC++をメインで使い、行列計算ライブラリとしてEigenを使っている．C++の高速性とMatlab的な行列…

octaveでisstable( )する話． 追記：octave-forgeのcontrolコンポーネントの中にisstable()関数がありました．そちらを使うことを推奨します． 線形システムの内部安定性 フライトシミュレーションの記事でも扱ったけれど，線形なシステムは一般にという形で…

6DOF flight simulation編の最終回． 飛行機の絵を書こう 飛行運動の可視化には様々なスタイルがあるけれど，やはり3次元的な飛行軌道の可視化が最も良い。octaveにもそういう機能はある(comet3関数)．ただこれは質点運動をアニメーションしてくれるのみで，…

ここからはポスト処理について解説する。 軌道データを生成する これまでに、定常飛行状態からの微小擾乱の時間発展を計算してきた。 そこから飛行軌道を生成する方法は5章にサラッと書かれている。1章に基づいて軌道計算すればよい。 1章では，慣性座標系と…

渦格子法で3次元空力計算を行う． Vortex Lattice Method + High AoA Delta Wing 揚力面理論の一種である渦格子法（VLM）を用いて高AoAデルタ翼の空力計算をした． 使用PC MacBookAir2011 使用ツール Xcode C++ Eigen gnuplot LEVも含めて解いている．数値計…

本編．前回の記事で示したスクリプトの中身の解説． 概要 Matlab互換のフリーソフトoctaveを使って、簡単な6自由度フライトシミュレーションを行う。6自由度というのは，xyzの3軸並進運動とピッチ/ロール/ヨーの3軸姿勢運動の自由度を合わせたもの。要するに…

はじめに 日本の飛行力学の教科書と言えば「航空機力学入門」。 一般向けの航空工学本と比べると内容が難しいので、学生の間では鬼門扱いされるとか。 反面、勘所(第3章)さえ押さえれば、一気に取っつき易くなる本でもある。 1~2章は前置き，3章が力学系の構…