倒立振子のソースを表示

{{翻訳直後|[[:en:Special:Redirect/revision/749658225|en:Inverted pendulum]]|date=2016年11月}}
[[ファイル:Balancer_with_wine_3.JPG|サムネイル|単純な[[ロボット工学|ロボティクス]]系、台車駆動型倒立振子。1976年。]]
'''倒立振子'''（とうりつしんし {{lang-en-short|inverted pendulum}}）とは、[[支点]]よりも[[重心]]が高い位置にある[[振り子]]をいう。写真に示すように、支点を台車に載せて実装する、台車駆動型倒立振子がよく知られている<ref>C.A. Hamilton Union College Senior Project 1966</ref>。ほとんどの応用例において振り子はある[[自転|回転軸]]まわりにのみ運動するよう固定されており、[[自由度]]は1に制限されている。振り子は吊り下げられた状態が安定であり、したがって倒立振子は本質的に不安定であるため倒立状態を保つためには常に能動的に制御する必要がある。このためには支軸に直接トルクを加えたり、支点を水平方向もしくは鉛直方向に移動させることにより振り子の錘の回転速度を変化させることにより総トルクを変化させたりして[[フィードバック]]系を構築する。支点を動かす型の倒立振子の一番簡単な例は、箒を手の上に立たせる系である。倒立振子は[[動力学]]と[[制御理論|制御論]]における古典的な問題であり、制御戦略の試験用ベンチマークとして用いられる。また、[[パーソナルモビリティ]]分野において二輪倒立振子（例: [[セグウェイ]]）や[[自立安定一輪車]]（例: [[ホンダ・U3-X]]）などの形式で応用されつつある。

別の種類の倒立振子として、高層建築に用いられる傾斜計が挙げられる。建築物の基部にワイヤの一端を固定し、他端につけた浮きを建築物の最上部に油で浮かべ、浮きの平衡位置の変化を測定することで傾斜を検出する<ref>{{cite web|url=http://soil.co.uk/products/inclination/inverted-hanging-pendulum-systems/|title=Inverted & Hanging Pendulum Systems|publisher=Soil Instruments|accessdate=2016-11-19}}</ref><ref>{{Cite conference|和書|url=http://library.jsce.or.jp/jsce/open/00035/2003/58-6/58-6-0178.pdf|title=浮き式傾斜計の鋼管傾斜計測への適用について|conference=土木学会第58回年次学術講演会|date=September 2003|format=PDF}}</ref>。

== 概要 ==
倒立振子は[[動力学]]および[[制御理論|制御論]]における古典的な問題であり、制御アルゴリズム（[[PID制御]]、[[状態空間 (制御理論)|状態空間表示]]、[[ニューラルネットワーク]]、[[ファジィ制御]]、[[遺伝的アルゴリズム]]など）の試験用ベンチマークとして広く用いられている。この問題は[[リンク機構|リンク]]の追加や台車への移動指令、台車をシーソーの上に載せるなどの様々な拡張が可能である。また、ロケットやミサイルは空気抵抗の力の中心が重心よりも前側にあるため[[空気力学]]的に不安定であり、その[[ミサイルの誘導方式|誘導制御問題]]と倒立振子問題は関連している<ref>http://exploration.grc.nasa.gov/education/rocket/rktstab.html</ref>。このような問題を理解するため、単純なロボット台車の上に載せた倒立振子が用いられることが多い。手の上に箒を立てる遊びや、自律平衡型の輸送機械である[[セグウェイ]]などもこの問題の解である。

フィードバックに頼らず倒立振子を安定化する別の方法として、支点を上下に高速に振動させる方法がある。このような装置を{{仮リンク|カピッツァの振り子|en|Kapitza's pendulum}}と呼ぶ。支点を（加速度と振幅の意味で）十分に強く振動させることにより、直感には強く反するが倒立振子に加えられた擾乱から回復することが可能となる。支点が[[自由振動|単振動]]する場合は、振り子の動きは[[マシュー函数|マシュー方程式]]により記述することができる。

== 運動方程式 ==
倒立振子の[[運動方程式]]は、倒立振子に課せられる拘束によって変化する。倒立振子の構成には様々なものがありうるため、それらを記述する運動方程式も数多く存在する。

=== 支点固定型 ===
振り子の支点が空間的に固定されている場合、運動方程式は{{仮リンク|label=非倒立振子|振り子 (数学)|en|Pendulum (mathematics)}}のものと同一となる。次に示す方程式は運動に抗する摩擦力がなく、棒は質量のない剛体と見做すことができ、[[2次元|平面内]]に運動が制限されているという仮定をおいた場合の運動方程式である。
: <math>\ddot \theta - {g \over \ell} \sin \theta = 0</math>
ここで、{{Mvar|{{Ddot|&theta;}}}} は[[角加速度]]を、{{mvar|g}} は[[重力加速度]]を、{{mvar|ℓ}} は振り子の長さを、{{mvar|&theta;}} は平衡位置からの変位角を表わす。

片々足し合わせることにより角加速度項と重力項が同符号となるように変形することもできる。
: <math>\ddot \theta = {g \over \ell} \sin \theta</math>
したがって、倒立振子は不安定平衡点からみた初期位置の方向に向かって、長さに反比例した加速度で遠ざかることがわかる。したがって、振り子は長い方が短い振り子よりも遅く落ちる。

'''トルクと慣性モーメントを用いた導出'''

質量 {{mvar|m}} の質点を、長さ {{mvar|ℓ}} の質量のない剛体棒の支点と逆側に固定した構成の振り子を考える。

系の総[[トルク]]は[[慣性モーメント]]と角加速度をかけた値に等しくなる。
: <math>\boldsymbol{\tau}_{\mathrm{net}}=I \ddot \theta</math>
重力による総トルクは以下のように書ける。
: <math>\boldsymbol{\tau}_{\mathrm{net}}= m g \ell \sin \theta\,\!</math>
ここで、{{Mvar|&theta;}} は倒立平衡点から測った変位角である。

したがって、次の運動方程式を得る。
: <math> I \ddot \theta= m g \ell \sin \theta\,\!</math>
質点の慣性モーメントは次のように得られる。
: <math>I = m R^2</math>
倒立振子の場合、動径は棒の長さ {{mvar|ℓ}} に等しい。

{{Math|1=''I'' = ''mℓ''<sup>2</sup>}} を代入すると以下が得られる。
: <math> m \ell ^2 \ddot \theta= m g \ell \sin \theta\,\!</math>
{{Math|''mℓ''<sup>2</sup>}} で辺々割って次を得る。
: <math>\ddot \theta = {g \over \ell} \sin \theta</math>

=== 台車駆動型倒立振子 ===
[[ファイル:Cart-pendulum.svg|サムネイル|330x330ピクセル| 台車上の倒立振子を表わす模式図。棒は質量を持たないものとする。台車の質量と棒の端点の質量はそれぞれ {{mvar|M}} および {{mvar|m}} とする。棒の長さは {{mvar|ℓ}} とする。]]
台車駆動型倒立振子は図のように、質量 {{mvar|m}} の質点を長さ {{mvar|ℓ}} の棒の先に付け、逆側を水平に動ける台につけた構成である。台車は{{仮リンク|線形運動|en|Linear motion}}しかしないものとし、運動を起こさせたり妨げたりするような力がはたらくものとする。

=== 安定化の要点 ===
倒立振子を安定化するための要点は定性的に三つのステップにまとめることができる。
# もし傾き角 {{mvar|&theta;}} が右向きならば台車は右に加速する必要があり、逆もまたなりたつ。
# 軌条中心からの相対台車位置 {{mvar|x}} を安定させるには、ヌル角度（制御系がゼロにしようとする角度誤差）を位置で変調すればよい。具体的にはヌル角度を {{Math|''&theta;'' + ''kx''}} とする。ただし、{{mvar|k}} は小さくとる。これにより、台車が軌条中心から外れているなら棒は中心に向けてわずかに傾こうとするので、台車は中心向きに加速する。軌条中心では棒が鉛直に立つことで安定となる。傾きセンサのオフセットや軌条の傾きなどの不安定要素は安定位置のオフセットとして表われる。さらにオフセットを追加すると位置制御ができる。
# クレーンで荷を運ぶときのように、運動する支点に吊り下げられた非倒立振子の応答には振り子の角振動数 {{Math|1=''&omega;''<sub>p</sub> = {{sqrt|''g''/''ℓ''}}}} にピークがある。制御不能な振動を防ぐため、{{Math|''&omega;''<sub>p</sub>}} に近い支点運動の周波数スペクトルを抑制する必要がある。倒立振子でも、安定化のために同じようなフィルタが必要となる。
ヌル角度変調戦略の結果として、位置フィードバックは正のフィードバックとなり、突然右に動くよう指令すると台車はまず左へ動いてから右へ傾いた振り子を再び平衡にするため右へと移動する。倒立振子の不安定性と、正の位置フィードバックによる不安定性との相互作用により系は安定となっており、これを数学的に解析することは興味深く奥の深い問題となっている。
[[ファイル:Inverted_Pendulum_control_essentials.JPG|サムネイル|center|600x600ピクセル| ワイングラスを載せた台車の単純な安定化制御系]]{{clear}}
=== ラグランジュ方程式 ===
[[ラグランジュ力学|ラグランジュ方程式]]から運動方程式を導出することもできる。前掲の図のように、{{Math|''&theta;''(''t'')}} を長さ {{mvar|ℓ}} の振り子の直立位置からの変位角とし、作用する力は重力および {{mvar|x}} 方向への外力 {{mvar|F}} とする。 {{Math|''x''(''t'')}} を台車の位置と定義すると、系の[[ラグランジュ力学|ラグランジアン]] {{Math|1=''L'' = ''T'' &minus; ''V''}} は以下のように書ける。
: <math>
L = \frac{1}{2} M v_1^2  + \frac{1}{2} m v_2^2 - m g \ell\cos\theta
</math>
ここで {{math|''v''<sub>1</sub>}} は台車の速度、 {{math|''v''<sub>2</sub>}} は質量 {{mvar|m}} の質点の速度とする。{{math|''v''<sub>1</sub>}} および {{math|''v''<sub>2</sub>}} は {{mvar|x}} と {{mvar|&theta;}} の導関数を用いて次のように書ける。
: <math>
v_1^2=\dot x^2
</math>
: <math>
v_2^2=\left({\frac{\mathrm d}{\mathrm dt}}{\left(x- \ell\sin\theta\right)}\right)^2 + \left({\frac{\mathrm d}{\mathrm dt}}{\left( \ell\cos\theta \right)}\right)^2
</math>
{{math|''v''<sub>2</sub>}} の表式を展開すると以下のようになる。
: <math>
v_2^2= \dot x^2 -2 \ell \dot x \dot \theta\cos \theta + \ell^2\dot \theta^2

</math>
すると、ラグランジアンは次のように書ける。
: <math>
L = \frac{1}{2} \left(M+m \right ) \dot x^2 -m \ell \dot x \dot\theta\cos\theta + \frac{1}{2} m \ell^2 \dot \theta^2-m g \ell\cos \theta
</math>
ここで、[[オイラー＝ラグランジュ方程式|オイラー・ラグランジュの運動方程式]]は次の表式である。
: <math>
\frac{\mathrm{d}}{\mathrm{d}t}{\partial{L}\over \partial{\dot x}} - {\partial{L}\over \partial x} = F
</math>

: <math>
\frac{\mathrm{d}}{\mathrm{d}t}{\partial{L}\over \partial{\dot \theta}} - {\partial{L}\over \partial \theta} = 0
</math>
この方程式系における {{mvar|L}} に上の表式を代入すると、倒立振子の運動を記述する方程式系が次のように得られる。
: <math>
\left ( M + m \right ) \ddot x - m \ell \ddot \theta \cos \theta + m \ell \dot \theta^2 \sin \theta = F
</math>
: <math>
\ell \ddot \theta - g \sin \theta = \ddot x \cos \theta
</math>
この方程式系は非線形であるが、制御系の目標は振り子を直立に保つことなので、{{Math|''&theta;'' &asymp; 0}} 近傍で線形化することが多い。

=== ニュートンの第二法則 ===
[[ラグランジュ力学|ラグランジュ方程式]]の代わりに、[[運動の第2法則|ニュートンの第二法則]]を用いる方が利点が多いことが多い。ニュートン方程式を用いれば、振り子と台車の間の結節点における反力が得られる。すると、それぞれの物体について {{mvar|x}} 方向と {{mvar|y}} 方向のあわせて二つの方程式が得られる。台車の運動方程式は次のようになる。ただし、左辺は合力、右辺は加速度を表わす。
: <math>
F-R_x = M \ddot x
</math>
: <math>
F_N - R_y - M g = 0
</math>
ここで、{{Math|''R''<sub>''x''</sub>}} および {{Math|''R''<sub>''y''</sub>}} は結節点における反力を表わす。{{Math|''F''<sub>''N''</sub>}} は台車にかかる垂直抗力である。二つ目の方程式は鉛直方向の反力にしか依存しないので、垂直抗力を解くのに用いることができる。一つ目の方程式は水平方向の反力を解くのに用いることができる。この運動方程式系を完成させるには、振り子に固定された質点の加速度を計算する必要がある。慣性座標系における質点の位置は次のように表わすことができる。
: <math>
\vec r_P = (x-\ell \sin \theta) \hat x_I + \ell \cos \theta \hat y_I
</math>
二階微分をとることにより慣性系における加速度ベクトルを得る。
: <math>
\vec a_{P/I} = (\ddot x + \ell \dot \theta^2 \sin \theta - \ell \ddot \theta \cos \theta ) \hat x_I + (-\ell \dot \theta^2 \cos \theta - \ell \ddot \theta \sin \theta) \hat y_I
</math>
第二運動法則を用いて、 {{mvar|x}} 方向と {{mvar|y}} 方向の二つの方程式を得る。質点にかかる反力は台車にかかる反力と逆符号になることに注意されたい。これは[[運動の第3法則|ニュートンの第三法則]]から帰結する。
: <math>
R_x = m(\ddot x + \ell \dot \theta^2 \sin \theta - \ell \ddot \theta \cos \theta )
</math>
: <math>
R_y - m g = m (-\ell \dot \theta^2 \cos \theta - \ell \ddot \theta \sin \theta)
</math>
一つ目の方程式から、外力 {{mvar|F}} によらず水平方向の反力を計算する方法が得られる。 二つ目の方程式は鉛直方向の反力の計算に用いることができる。一つ目の運動方程式は <math>F-R_x = M \ddot x</math> を <math>R_x = m(\ddot x + \ell \dot \theta^2 \sin \theta - \ell \ddot \theta \cos \theta )</math> に代入して次のように得られる。
: <math>
\left (M+m \right) \ddot x - m \ell \ddot \theta \cos \theta  + m \ell \dot \theta^2 \sin \theta  = F
</math>
この方程式をよく見れば、ラグランジュ法による方程式と同じであることがわかる。二つ目の方程式を得るには、振り子の運動方程式と、常に振り子に直交する単位ベクトルとのドット積をとる必要がある。このベクトルは物体系の {{mvar|x}} 座標と呼ばれることが多い。慣性座標系ではこのベクトルは単純な二次元座標変換を用いて以下のように表わされる。
: <math>
\hat x_B = \cos \theta \hat x_I + \sin \theta \hat y_I
</math>
ベクトル表記で書いた振り子の運動方程式は <math>\sum \vec F = m \vec a_{P/I}</math> のようになる。この両辺と {{Math|{{Hat|''x''<sub>B</sub>}}}} とのドット積を取ると、求める運動方程式の左辺は次のようになる（転置をとってかけるとドット積になることに注意）。
: <math>
(\hat x_B)^T\sum \vec F = (\hat x_B)^T (R_x \hat x_I + R_y \hat y_I - m g \hat y_I) = (\hat x_B)^T(R_p \hat y_B - m g \hat y_I) = -m g \sin \theta
</math>
上の式展開においては反力の物体系成分と慣性系成分との関係式が用いられている。質点と台車を繋ぐ棒が質量を持たないという仮定から、この棒は自身に垂直な負荷をまったく伝えないと考えられる。したがって、反力の物体系成分は単に {{Math|''R''<sub>p</sub> {{Hat|''y''<sub>B</sub>}}}} のように書け、棒に沿った成分しか存在しないことがいえる。このことから、棒の張力に関して次の方程式が得られる。
: <math>
R_p = \sqrt{ R_x^2 + R_y^2}
</math>
同様に、先に示した加速度 {{Math|{{Vec|''a''}}<sub>P/I</sub>}} の式と {{Math|{{hat|''x''<sub>B</sub>}}}} とのドット積を取ることで、求める運動方程式の右辺も計算することができる。その結果は（少々整理すると）以下のようになる。
: <math>
m(\hat x_B)^T(\vec a_{P/I}) = m(\ddot x \cos \theta - \ell \ddot \theta)
</math>
左辺と右辺を繋いで {{mvar|m}} で割ると以下を得る。
: <math>
\ell \ddot \theta - g \sin \theta = \ddot x \cos \theta
</math>
これもまた、ラグランジュ法による方程式と同一である。ニュートン法を用いる利点は、全ての反力が明らかであり、損壊がないことを保証できる点にある。

=== 支点の振動する振り子 ===
質量の無い、振動する台に繋がれた振り子の運動方程式は台車に繋がれた振り子と同じように導出することができる。質点の位置はこの場合は次のように表わされる。
: <math>\left( -\ell \sin \theta , y + \ell \cos \theta    \right)</math>
そして、位置の一階微分をとることにより速度が得られる。
: <math>v^2=\dot y^2-2 \ell \dot y \dot \theta \sin \theta  + \ell^2\dot \theta ^2</math>
この系の[[ラグランジュ力学|ラグランジアン]]は次のように書ける。
: <math>
 L = \frac{1}{2} m \left ( \dot y^2-2 \ell \dot y \dot \theta \sin \theta  + \ell^2\dot \theta ^2   \right) - m g \left( y + \ell \cos \theta  \right )
</math>
そして、オイラー・ラグランジュ方程式
: <math>
{\mathrm{d} \over \mathrm{d}t}{\partial{L}\over \partial{\dot \theta}} - {\partial{L}\over \partial \theta} = 0
</math>
から次の運動方程式が得られる。
: <math>
\ell \ddot \theta - \ddot y \sin \theta = g \sin \theta
</math>
ここで、{{mvar|y}} が[[自由振動|単振動]] {{Math|1=''y'' = ''A'' sin ''&omega;t''}} で表わされる場合、次の[[微分方程式]]を得る。
: <math>
\ddot \theta - {g \over \ell} \sin \theta = -{A \over \ell} \omega^2 \sin \omega t \sin \theta
</math>
この方程式は閉形式の基本解を持たないが、様々な方法で調べることができる。たとえば、振幅が小さい場合は[[マシュー函数|マシュー方程式]]により精度よく近似することができる。解析により、振動が速ければ振り子が直立を保つことがわかる。下図のプロットの一つ目は {{mvar|y}} の振動が遅い場合には振り子を直立状態からずらすとすぐに倒れてしまうことを示している。短時間で {{mvar|&theta;}} は 90°を超えており、つまり振り子は倒れてしまっている。{{mvar|y}} の振動が速ければ振り子は直立位置の周りで安定に保たれる。二つ目のプロットは振り子が直立位置 ({{Math|1=''&theta;'' = 0}}) からずれてもその近傍で振動することを示している。直立位置からのずれは小さくたもたれ、倒れてしまうことはない。
[[ファイル:Inverted_pendulum_oscillatory_base.svg|center|サムネイル|500x500ピクセル|振動する台に載せられた倒立振子のプロット。左のプロットは遅い振動に対する振り子の応答を、二つ目のプロットは速い振動に対する応答を示す。]]

== 倒立振子の種類 ==
倒立振子を安定させることは広く取り組まれている工学的課題である<ref>http://robotics.ee.uwa.edu.au/theses/2003-Balance-Ooi.pdf</ref>。台車駆動型倒立振子にも、棒を載せただけのものから振り子が複数の節に分かれているものまで様々な変種がある。その他にも、倒立振子の棒もしくは節分けされた棒を回転部品の端にとりつける種類のものがある。台車型と回転型のいずれの場合でも、倒立振子は平面内でしか運動しない。平衡位置にある倒立振子を維持するのが課題である場合もあれば、自律的に平衡位置を達成することが求められる場合もある。また、二輪倒立振子という別のプラットフォームもある。二輪倒立振子では一点での回転が可能であり、格段に高い操作性を実現できる<ref name="#1">http://csuchico-dspace.calstate.edu/bitstream/handle/10211.4/145/4%2022%2009%20Jose%20Miranda.pdf?sequence=1 Resource not found. {{404|date=2024-05}}</ref>。その他にも、一点上での平衡をとるものもある。[[独楽|コマ]]や [[一輪車]]、ボール上の倒立振子はすべて一点上で平衡をとる。上で説明したように、鉛直に振動する台によっても倒立振子を安定化することができる。

== 倒立振子の例 ==
倒立振子には人造のものも自然のものも含めて様々な例がある。

最も身近な倒立振子は[[ヒト]]であるといわれる。直立した体を持つヒトは、常に平衡を保つための調整を行なっていなければ立つことも歩くことも走ることもできない。

単純な例として、箒や定規を手の上で直立させる遊びが挙げられる。

倒立振子は様々な装置に組込まれており、また特異的な工学的問題として研究されている<ref name="#1"/>。倒立振子は本質的に不安定なので擾乱に対して計測可能な応答を示すため、初期の[[地震計]]の設計における中心的要素として採用されていた<ref>http://earthquake.usgs.gov/learn/topics/seismology/history/part12.php Not Found0. {{404|date=2024-05}}</ref>。

倒立振子模型はいくつかの[[パーソナルモビリティ]]にも採用されている。二輪車椅子やその他の二輪モーター車両により高い移動性を確保することができる。

== 関連項目 ==
* [[振り子]]
* [[二重振り子]]
* [[振動]]
* [[ゼロモーメントポイント]]
* {{仮リンク|二重倒立振子|en|Double inverted pendulum}}
* {{仮リンク|慣性車輪振子|en|Inertia wheel pendulum}}
* {{仮リンク|古田の振子|en|Furuta pendulum}}
* [[一輪車]]
* [[自立安定一輪車]]
* [[セグウェイ]]
* [[ホンダ・U3-X]]
* [[ハロー! ミップ]]
* [[球体ロボット]]
* {{仮リンク|iBOT|en|iBOT}}
* [[二足歩行ロボット]]
* {{仮リンク|人型ロボット|en|Humanoid robot}}

== 出典 ==
{{Reflist}}

== 関連文献 ==
*{{Cite book|last=Franklin|author=|first=Gene|title=Feedback control of dynamic systems|edition=5|date=|year=2005|publisher=Prentice Hall|ISBN=0-13-149930-0|last2=Powell|last3=Emami-Naeini|first2=J. D.|first3=Abbas}}
* {{Cite journal|和書|author=|year=1993|title=現代制御理論を使った倒立振り子の実験|journal=[[トランジスタ技術]]|volume=|issue=5月号|publisher=[[CQ出版]]|pages=315-322}}
* {{Cite journal|和書|author=|year=1993|title=現代制御理論を使った倒立振り子の実験|journal=[[トランジスタ技術]]|volume=|issue=6月号|publisher=[[CQ出版]]|pages=367-373}}
* {{Cite journal|和書|author=|year=1993|title=現代制御理論を使った倒立振り子の実験|journal=[[トランジスタ技術]]|volume=|issue=7月号|publisher=[[CQ出版]]|pages=363-370}}

== 外部リンク ==
* [https://www.youtube.com/watch?v=-vpc-mt_wPg YouTube - Inverted Pendulum - Demo #3]
* [https://www.youtube.com/watch?v=MWJHcI7UcuE YouTube - inverted pendulum]
* [https://www.youtube.com/watch?v=B6vr1x6KDaY&feature=relmfu YouTube - Double Pendulum on a Cart]
* [https://www.youtube.com/watch?v=cyN-CRNrb3E&feature=related YouTube - Triple Pendulum on a Cart]
* [http://mw.concord.org/modeler1.3/mirror/mechanics/inversependulum.html A dynamical simulation of an inverse pendulum on an oscillatory base]
* [http://www.library.cmu.edu/ctms/ctms/examples/pend/invpen.htm Inverted pendulum modeling with several control systems in Matlab]
* [http://www.engr.usask.ca/classes/EE/480/Inverted%20Pendulum.pdf Inverted Pendulum: Analysis, Design, and Implementation]
* [http://lars.mec.ua.pt/public/LAR%20Projects/Humanoid/2006_MiltonRuas-Projecto/Simuladores/Pendulo%20Invertido/bugeja.pdf Non-Linear Swing-Up and Stabilizing Control of an Inverted Pendulum System]
* [http://neuron.tuke.sk/~vascak/predmety/FSR/Eseje/Seliga%20-%20Stabilization%20fuzzy%20control%20of%20inverted%20pendulum%20systems.pdf Stabilization fuzzy control of inverted pendulum systems]
* [http://mech.u-fukui.ac.jp/~Kawa-Lab/kenkyu/pend/pend.htm 倒立振子系に対する安定化制御]
* [http://kinno-homepage.sakura.ne.jp/pendulum/report.pdf ワンチップマイコンを応用した自立型倒立振り子]
* {{YouTube|J4CP9I88y_g|Hello! MiP}}
* {{YouTube|AfFBa4JG1n4|倒立2輪 新タイプ}}

{{DEFAULTSORT:とうりつしんし}}
[[Category:動力学]]
[[Category:制御工学]]
[[Category:振り子]]
[[Category:機械要素]]
[[Category:ロボット工学]]