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 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()