File:Inverted pendulum oscillatory base.svg

此SVG文件的PNG预览的大小：800 × 322像素。其他分辨率：320 × 129像素 | 640 × 258像素 | 1,024 × 412像素 | 1,280 × 516像素 | 2,560 × 1,031像素 | 1,070 × 431像素。

原始文件 ‎（SVG文件，尺寸为1,070 × 431像素，文件大小：60 KB）

摘要

描述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 divisionimport numpy as npfrom numpy import sin, cos, pifrom scipy.integrate import odeintimport matplotlib.pyplot as pltfrom matplotlib import rcParamsrcParams['font.size'] = 14rcParams['legend.handlelength'] = 0def 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 dxdtplt.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.Tplt.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.Tplt.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()