2015IYPT及CUPT参考资料.docx

Cupt T13：磁力摆一.问题回顾：Magnetic PendulumMake a light pendulum with a small magnet at the free end. An adjacent electromagnet connected to an AC power source of a much higher frequency than the natural frequency of the pendulum can lead to undamped oscillations with various amplitudes. Study and explain the phenomenon.制作一个自由端有一个小磁体的轻摆。一个临近的通有频率远高于摆的固有频率交流电的电磁铁可以导致不同振幅的无阻尼振动。研究并解释这个现象。二.关键词解读：1.交流电频率、摆的固有频率；2.无阻尼振动；3.研究、解释。三：模型假设：1.小磁体看作质点；2.空气阻力不可忽略；3.绳的质量不计；4.摩擦不计。四：模型建立：受力分析：微分方程： ，， 令，写成标量式：，两式联立得： ，化简得的二阶非线性微分方程：.在磁场中，方程化为：令，方程化为一阶的非线性微分方程组：再令，得一阶的非线性微分方程组：令，，，方程化为：回到题目，题目要求，即，还要求不同振幅的无阻尼振动，无阻尼振动的定义即：振幅不变的振动。所以要求交流电场能提供的能量抵消阻力的损耗；现在的关键是，如何造成不同振幅的无阻尼运动。分三方面来考虑：1.力学条件：，应维持在与其他力同一量级上，故初步取，后面将看到，的最优值的确在1左右；2.频率要求：要远大于1，还是远小于1（题设有无问题，之前有同学提出的主要疑问）；3.实现的条件：不论是化为方程组还是看原方程，均存在非线性的敏感项，即，其敏感程度取决于参数。方程的求解借助于MATLAB，并不困难。但求解的可信度有多大？首先我们先按题设要求的远大于的情形来讨论。不妨取=10。若绳长为半米，则约为60，接近200。若达到此数值，则周期为0.033s左右，可以看到，若我们选定某一初始条件直接代入方程求解，其必对应于周期内某一点。但是如果偏差值达到0.01，则已错过约1/3周期，这意味着cos的值与代入值已相差甚远。若绳长进一步缩小，则对偏差值要求更小。如果按疑问中所提出的，远小于1，同样的，我们取=0.1。这时，周期3.3s。0.1s的偏差影响也不大。这样的分析虽然在理论上并无太大意义，因为周期再小，也有方法将结果逼近，但对于实验却非常重要。有同学提到，远大于1时，在实验上看不到所谓不同振幅的振动，其实这是非常正常的。T如此小的情况下，已经无法判断初始条件，甚至可以说，初始条件是随机的。而远小于1时，情况则大不相同。实验的初始条件是可测的，且误差值较小。所以说，对于远大于1的情形，从实验上判断时，仅研究一两组或者同一个位置并无太大的意义。这是因为随机性下研究逐个的点是无意义的。对于远小于1的情形，实验则简单得多，得到的现象也较明显。但这并不意味着说，远小于1时，摆一定能做不同振幅的无阻尼运动。对于初始条件仍有要求。但是条件较易满足。换言之，无论取何值，要达到所谓不同振幅的无阻尼运动，都是有办法的。只不过题设给出的情形较难满足。这也很容易理解，因为方程的敏感性并不取决于，而是取决于非线性性，这也是为何文献中改变除了的其他参量时，也得到了非常敏感的结果。解决了这个疑惑后，我们的着眼点在于，如何进行实验？由于理论和实验间的联系较少，故采用逼近的方法。由大数定律，设是独立同分布的随机变量序列，具有相同的有限数学期望，，则对于任意有：这意味着，只要样本值足够多，完全可以逼近理论值。对于一个点，这就是所谓的取平均值的方法，其效果是显然的。但前面也已提到过，研究一个点的意义并不大，那么我们反其道而行，取非常多的点，比较通常的取法是将圆周360度等分成50份或更多。由此得到随机变量的观察值，然后改变方程中的参量，得到另一组观察值，将所有这些观察值画在图上，会有什么效果？从理论上讲，就是所谓的解的歧化图，下面是改变一系列参量时的歧化图：可以看出，歧化图的歧化几乎是无规律的（虽然部分简单方程有周期性规律），那么我们得到了试验和理论的图象，剩下的只需要比较两图象的重合度即可。重合度虽然可以看重合的面积的大小，但这么做并不可取。主要应该看两方面：1.歧化点的位置；2.歧化点的数量。然后再比较面积，不过很可能偏差较大，因为样本数量即使有数百个，也不能保证说能较均匀地分布在可能的点上，所以其实比较面积意义不大。从以上分析可以看出，理论和实验联系的关键是大样本实验。换言之，实验可以独立于理论之外，为了提高理论的精确性，再考虑考虑理论中唯一不确定的阻