弹簧振子运动方程推导

    昨日准备一份物理家教,关于弹簧振子的。有种推导其运动方程的冲动。下面的推导也不难,但大学物理教材中貌似没有,还是写下来吧,免得忘了。

推导过程

考虑一个简谐振子,描述一个质量在弹簧上做简谐运动的情况。该系统可以用以下常微分 方程组表示:    dx dt = v dv dt = −ω 2x 其中: • x 是质量的位置。 • v 是质量的速度。 • ω 是角频率,与弹簧的劲度系数 k 和质量 m 有关,即 ω = √ k m 。 给定初始条件 x(0) = x0 和 v(0) = v0,该方程组的解析解为:    x(t) = x0 cos(ωt) + v0 ω sin(ωt) v(t) = −x0ω sin(ωt) + v0 cos(ωt) 3ω = √ k m 假设我们有一个弹簧振子,其劲度系数 k = 4 N/m,质量 m = 1 kg,初始位置 x0 = 1 m,初始速度 v0 = 0 m/s。我们可以计算角频率 ω 和位置 x(t)、速度 v(t) 的表达式。 1. 写出这个问题的解析解表达式, 并画出 t ∈ [0, 20] 时间段的解的函数图像. 2. 分别利用向前、向后 Euler 格式、Euler 预估-校正格式、四级四阶 Runge-Kutta 求解这 个问题,分别取同样的时间步长为 ∆t = 0.1, 0.05, 0.01, 0.001,分析它们的数值误差。将 这四种方法的数值解与解析解画在同一张图中. 再把它们每一个时间步的误差也画在同 一张图中. 3. 你会发现随着时间的增长,数值解的振幅会逐渐增大或衰减. 从 n 到 n + 1 时间层的数 值解的增长矩阵 定义如下: ( x n+1 v n+1 ) = G ( x n v n ) 分析每个方法的行列式,是大于 1 还是小于 1,亦或等于 1?分别对应于数值解的振幅增 长还是衰减的? 代码均为matlab代码
最新发布
06-26
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值