Question

摆锤的质量为m，正电荷为q。电容器极板平行于地球表面。电容器内部的电场垂直向上指向，其大小通过摆运动进行调制，其中E（t）= E_0 * | sin（θ）（t），其中qE_0 / mg = <1，且θ是摆之间的角度和垂直线。初始条件为θ（0）= pi / 2 rad和dθ/ dt = 0 rad / s。

L = 1.0m，g = 9.8 m / s ^ 2

（a）首先，取= 0，估计数值获得的周期等于公式T = 2pi * sqrt（L / g）所预测的周期的初始角度θ（0），精度优于1 ％。

（b）查找并绘制该摆的振荡周期与参数的依赖关系。

（c）如果= 1，摆将发生什么？

所以对于a）这就是我所做的

#import
%pylab nbagg
import numpy as np
from scipy.integrate import odeint

#Solve for T first
#let the length L be equal to 1
L = 1
g=9.8
T = 2*np.pi*np.sqrt(g)
print(T)

#let delta = d
#if d = 0 then the ODE becomes

#Define the ODE

w_0 = np.sqrt(g/L)

def dy_dt(y,t):
  y1 = y
  y2 = -w_0**2*sin(y1)
  dydt = (y1,y2)
  return dydt

#Integration values and interval
t_0 = 0
t_f = T
nt = 10000
t = linspace(t_0, t_f, nt)

现在我不确定如何解决dθ/ dt时如何进行，但前提是dθ/ dt仅为0。

Answer 1

所以，我已经看了好几次了。看来这是一个非常有趣的问题，但我认为StackOverflow尚不是寻求帮助的正确位置。我们发现您是该网站的新手，因此，我鼓励您将其分解为多个问题，并将它们定向到适当的组。

由于此时您的代码实际上正在工作，因此StackOverflow并不是解决此问题的合适场所。我建议您在physics stack exchange上寻求有关此问题背后的物理问题的帮助（例如，如何建立所需的微分方程组），以及在computational science stack exchange网站上寻求帮助解决您通过数值获得的微分方程。

如果您在进一步解决问题时无法使代码本身正常工作，那么StackOverflow绝对是寻求帮助的最佳位置。

大平面平行平板电容器内部的简单摆

1 个答案: