《计算机模拟物理第二次作业报告》.pdfVIP

计算机模拟物理作业2 向左鲜（2008011833 ） 姚浩田（2008011823 ） 赵鲁宝（2008011812 ） 目录 一 问题描述2 二 计算机模拟结果2 2.1  3, 0.5, A 1, 2 的情形2 0 0 2.2  1和 4 时的情形7 2.3  4 时的情形10 0 2.4 A 0, 3 时的临界阻尼系数11 0 0 2.5 （补充内容）其他驱动力的条件下的位移曲线12 附 录16 一 振动问题的解析解16 1.1 无强迫力即f (t ) 0 的情形16 1.2 有强迫力的情形17 二 几种方法的收敛性分析20 2.1 Euler 法的收敛性分析20 2.2 Euler_Cromer 法的收敛性分析22 2.3 Euler_richardson 法的收敛性分析23 2.4 几种方法收敛性总结25 n i 三 n 阶常微分方程a (t) d y a (t) y f (t) 数值解法分析.25 i i 0 i 1 dt 参考文献27 一 问题描述 描述阻尼线性振子的方程为 d 2 x 2 dx  x  A cost dt2 0 dt 0 选取您认为合适的模拟方法，进行下面模拟。 1、 3, 0.5, A 1, 2 ，这是一个小阻尼振子，绘出初始条件为x 1,v 0 0 0 0 0 时的t~x(t) 图。 a、x(t) 的行为与无扰动时有何不同？ b 、从初始开始数次振动后，x(t) 的周期和角频率。 c、x(t) 的极限性质是否与初始条件相关。 d、x(t) 的短期行为是否取决于初始条件。划分x(t) 的与初始条件相关部分和与初始条件无 关的稳态部分。同时在本题和下面的模拟中考察能量守恒情况。 2 、分别计算 1和 4 时的x (t ) 。每种情况下的周期和角频率是多少？ 3、 4 计算时的x (t ) 。确定稳态运动时的角频率。根据1~3的计算结果，解释哪些参 0 数决定稳态频率。 4 A 0, 3 确定阻尼系数的临界值（临界阻尼系数） 0 0 在模拟时，由于方程为二阶方程，因此将方程化为两个一阶方程计算：  dx v(t)  dt  dv 2  x v A cost f (x,v,t ) dt 0 0 在计算时，先计算