同济大学《现代计算物理方法》讲义设计chap.pptVIP

* 第七章 常微分方程：牛顿运动方程求解 采用Matlab中的ode函数求解常微分初值问题 m文件：f.m ode函数 微分方程 求解区间 初始值 例一 例二：求单摆的运动方程 首先做单位变换 L 例二：求单摆的运动方程 摆幅较小时解析近似较好 在符合牛顿定律的运动系统中，存在一些不随时间改变的 运动常数： （1）系统的能量： （2）系统的动量： （3）系统的角动量： 通过运动常数来检验数值解的准确性 通过运动常数来检验数值解的准确性 单摆: 对于单摆系统，由于与外界支点的相互作用，只有系统能量是运动常数 L 通过运动常数来检验数值解的准确性 单摆: L 通过运动常数来检验数值解的准确性 能量随时间而减小 通过运动常数来检验数值解的准确性 摆动角度和角速度的幅度随时间而减小 通过运动常数来检验数值解的准确性 摆动角度和角速度的幅度随时间而减小 通过运动常数来检验数值解的准确性 在相空间中，运动轨迹不断向中心移动， 同样说明摆动角度和角速度的人为衰减。 采用Verlet算法求解牛顿方程 采用Verlet算法求解牛顿方程 Verlet蛙跳法 Euler算法 采用Verlet算法求解单摆运动方程 蛙跳法 采用Verlet算法求解单摆运动方程 摆动角度和角速度的幅度随时间保持恒定 采用Verlet算法求解单摆运动方程 摆动角度和角速度的幅度随时间保持恒定 采用Verlet算法求