擅长:python、mysql、java
<p>这里的问题是术语<code>np.cos(Wd*i)</code>。应该是<code>np.cos(Wd*i*dT)</code>,也就是说,自从<code>t = i*dT</code>以来,<code>dT</code>已经被添加到正确的方程中。</p>
<p>如果进行了此校正,则模拟看起来是合理的。这是一个带有<code>F0=0.001</code>的版本。注意,在阻尼条件下的持续振荡中,驱动力是明显的。</p>
<p><img src="https://i.stack.imgur.com/DF2JF.png" alt="enter image description here"/></p>
<p>原始方程的问题是<code>np.cos(Wd*i)</code>只是在圆周围随机跳跃,而不是在圆周围平滑移动,最终不会产生净效应。这可以通过直接绘制来最好地看到,但是最简单的事情是使用非常大的<code>F0</code>运行原始表单。下面是<code>F0 = 10</code>(即,1000x是正确方程式中使用的值),但使用了不正确的方程式形式,很明显,这里的驱动力只是在它随机绕圆移动时添加噪声。</p>
<p><img src="https://i.stack.imgur.com/n5LSh.png" alt="enter image description here"/></p>