ファイル:Inverted pendulum oscillatory base.svg

元のファイル ‎(SVG ファイル、1,070 × 431 ピクセル、ファイルサイズ: 60キロバイト)

このファイルは、ウィキメディア・コモンズから呼び出されています。
このファイルは、他のプロジェクトで使用されている可能性があります。
このファイルの解説やノートへの記入、履歴などの詳細の確認は、ウィキメディア・コモンズのファイルページ（ノート/履歴/ログ）を使用してください。

ウィキメディア・コモンズのファイルページにある説明を、以下に表示します。

概要

解説Inverted pendulum oscillatory base.svg	English: Plots illustrating the behaviour of an inverted pendulum mounted on an oscillatory base. The first plot shows the response of the pendulum on a slow oscillation ( $\omega =10$ ), the second the response on a fast oscillation ( $\omega =200$ ). The motion of the pendulum is given by the the following differential equation ${\ddot {\theta }}-{g \over \ell }\sin \theta =-{A \over \ell }\omega ^{2}\sin \omega t\sin \theta .$
日付	2016年7月4日
原典	投稿者自身による著作物
作者	Nicoguaro
SVG 開発 InfoField	このSVGのソースコードは正しい。このベクター画像はMatplotlibで作成されました。
ソースコード InfoField	Python code from __future__ import division import numpy as np from numpy import sin, cos, pi from scipy.integrate import odeint import matplotlib.pyplot as plt from matplotlib import rcParams rcParams['font.size'] = 14 rcParams['legend.handlelength'] = 0 def pend(x, t, grav=9.81, length=1, amplitude=0.05, omega=10): theta, dtheta = x dxdt = [dtheta, sin(theta)/length(grav - amplitudeomega*2 sin(omegat))] return dxdt plt.figure(figsize=(14, 5)) plt.subplot(121) time = np.linspace(0, 1.5, 101) sol = odeint(pend, [0.1, 0], time, args=(9.81, 1, 0.05, 10)) theta, _ = sol.T plt.plot(time, theta180/pi, color="#e41a1c", lw=2) plt.xlim(0, 1.5) plt.ylim(0, 100) plt.xlabel(r"$t$ (s)") plt.ylabel(r"$\theta$ (deg)") plt.legend([r"$\omega=10$"], numpoints=1, framealpha=0) plt.subplot(122) time = np.linspace(0, 3, 1001) sol = odeint(pend, [0.05, 0], time, args=(9.81, 1, 0.05, 200)) theta, omega = sol.T plt.plot(time, theta*180/pi, color="#e41a1c", lw=2) plt.xlim(0, 2.5) plt.xlabel(r"$t$ (s)") plt.ylabel(r"$\theta$ (deg)") plt.legend([r"$\omega=200$"], numpoints=1, framealpha=0) plt.savefig("inverted_pendulum_oscillatory_base.svg", bbox_inches="tight") plt.show()