M063 · 振动 / 参数识别

阻尼振动与振铃测量

从拍下的振动峰值，什么时候能算出阻尼和劲度？ 运动欠阻尼、每个峰都看得清的时候。峰值取对数后随时间均匀下降，下降的快慢就是阻尼率；峰与峰的间隔给出阻尼振动频率。再知道质量，就能换算出阻尼系数和劲度。画面噪声大，峰挑得不准，两个结果都会偏。

释放带黏性阻尼器的弹簧振子，用相机拍下带噪声的画面，把能分辨出的正峰值取对数，再做线性回归。由观测估出衰减率和阻尼频率，再用已标定的质量算出阻尼系数和刚度。比较欠阻尼时的螺旋、临界阻尼时的回归和过阻尼时的松弛。

经过时间 $0\,\mathrm{s}$

相对平衡位置的位移 $0\,\mathrm{m}$

速度 $0\,\mathrm{m\,s^{-1}}$

无阻尼参考周期 $0\,\mathrm{s}$

相对平衡位置的机械能 $0\,\mathrm{J}$

积分耗散能量 $0\,\mathrm{J}$

积分驱动做功 $0\,\mathrm{J}$

能量收支误差 $0\,\mathrm{J}$

相对解析解的位移误差 $0\,\mathrm{m}$

步长减半后的终点变化 $0\,\mathrm{m}$

实际积分步长 $0\,\mathrm{ms}$

沿运动坐标的有向弹簧力 $0\,N$

沿运动坐标的加速度 $0\,\mathrm{m/s^2}$

执行器瞬时功率 $0\,W$

瞬时耗散功率 $0\,W$

阻尼区间 $\text{采样中}$

包络 / 慢极点时间尺度 $0\,\mathrm{s}$

标称阻尼频率 $0\,\mathrm{Hz}$

相机峰值拟合衰减率 $\text{采样中}$

相机峰值拟合阻尼系数 $\text{采样中}$

相机峰值拟合刚度 $\text{采样中}$

峰值对数残差均方根 $\text{采样中}$

拟合使用的峰值数 $0$

已采集相机帧数 $1$

拍下振铃，拟合阻尼和刚度

识别衰减系统

物理受力矢量弹簧力、驱动力、阻尼力、重力和轨道反力，每一时刻都按同一个比例画出。速度矢量蓝色箭头，比例与力分开。正方向是向右或向下。近期已采集轨迹只显示最近两个参考周期，不画未来的运动。解析轨迹参考隐藏解析轨迹，只看相机采到的数据。运动质量

m

1.00\,\mathrm{kg}

弹簧刚度

k

20.00\,\mathrm{N\,m^{-1}}

释放位移

x_0

0.40\,\mathrm{m}

初速度

v_0

0.00\,\mathrm{m\,s^{-1}}

阻尼比

\zeta

0.06\,

相机帧率

f_s

40.00\,\mathrm{Hz}

位置噪声标准差

\sigma_x

0.00\,\mathrm{m}

最大积分步长

\Delta t_{\max}

20.00\,\mathrm{ms}

检查实验时刻

t/t_{\mathrm{end}}

0.00\,\%

可以旋转三维视角查看模拟的弹簧和阻尼器。在坐标视图里拖动质量块，或用键盘调节滑块。切到能量视图不会推进记录。重新采集只换一组随机噪声。蓝点是已采到的相机帧，金线是峰值对数的拟合；拟合不会提前读取记录中后面的帧。

振子在水平轨道上运动，阻尼是线性黏性阻尼，弹簧是无质量的胡克弹簧。三维装置显示的是模拟运动，不是拟合结果，也不是带噪声的相机位置；长度按整段运动的范围缩放，硬件只是示意，不计算碰撞和弹簧丝形变。模拟相机按精确均匀的时间间隔拍摄，每个位置加上独立的高斯噪声。要拟合峰值，至少需要三个高于八倍噪声标准差的正峰值，且每个标称阻尼周期至少八帧。已知质量和标称刚度用来设定峰值间隔和采样门槛。挑选峰值、有限采样和噪声都会让拟合结果偏离；对数拟合的残差不是置信区间。运动不振荡，或峰值分辨不出来时，不给出峰值拟合。

先抓重点

核心问题: 靠相机拍到的峰值，能求出两个物理系数吗？什么时候这个方法应该不给答案？
一句话直觉: 峰值取对数后下降多快，告诉你阻尼；峰与峰隔多远，告诉你阻尼频率。有了已知质量，就能换算成系数，但从噪声里挑峰值仍可能带来偏差。

核心数学模型

三种阻尼情形

m\ddot x+c\dot x+kx=0,\quad\gamma=\frac c{2m},\quad\zeta=\frac c{2\sqrt{mk}},\quad\omega_d=\sqrt{\omega_0^2-\gamma^2}

只有欠阻尼时，阻尼频率才是实数。临界阻尼时两个衰减极点重合；过阻尼时没有振荡，只剩两个实数衰减率。

由可分辨的峰值求衰减率

\log x_{\mathrm{peak},i}=a-\gamma t_i+r_i,\qquad\widehat\omega_d=\frac{2\pi(N-1)}{t_N-t_1}

对相机峰值的对数拟合直线，再量相邻正峰值的间隔，就能估出两个参数。每个采样峰值的位置用二次插值确定，离噪声底太近的峰值会被探测器舍弃。

求出刚度和阻尼

\widehat c=2m\widehat\gamma,\qquad\widehat k=m(\widehat\omega_d^2+\widehat\gamma^2)

质量需要单独标定。标称质量和刚度只用来设定峰值间隔和采样门槛，不会代替拟合出的衰减率和频率。从噪声里挑峰值，仍可能让结果偏离。

耗散能量单独积分

Q(t)=\int_0^t c\dot x^2\,\mathrm dt,\qquad E(t)-E(0)+Q(t)=0

耗散能量由阻尼功率积分得到，而不是拿少掉的机械能倒推。剩下的不平衡量用来检验数值解和步间插值。

常见难点

临界阻尼不是振得快

容易误解临界阻尼的振荡周期最短。

正确理解临界阻尼根本不振荡。量峰值间隔的方法用不上，但模拟本身仍然有效。

带噪声的峰值越多，偏差可能越大

容易误解每一个看起来像峰值的起伏都会改善拟合。

正确理解后期的峰值可能只是噪声。探测器要求峰值高过噪声底、彼此间隔不小于最短距离；这样挑选本身，也可能让对数振幅偏离。

采样与积分是两个时钟

容易误解求解器步长小，相机就一定能分辨振荡。

正确理解刚度高、相机帧率低时，采样门槛可能通不过，而解析参考和能量检验照样通过。

动手实验

01
先在无噪声时测量
把位置噪声标准差设为零，点“采集新的振铃”，把拟合出的刚度和设定值比一比。
你会看到： 由可分辨峰值求出的阻尼和刚度，与设定值只差有限采样带来的误差。
02
换一种阻尼情形
比较“可分辨振铃”“临界阻尼”“过阻尼回归”三个预设。绕着阻尼器旋转视角，暂停后用“精细运动单步”对照阻尼力和速度的方向。
你会看到： 螺旋消失了；运动不再振荡，也就不给出峰值拟合。
03
让相机更难测
选“不理想相机”，多采集几次，再调高刚度、调低相机帧率。
你会看到： 噪声会改变估计值。采样太粗时，实验干脆不给拟合，不会报出一个虚假的系数。
04
追踪损失的能量
切到“能量与功率”，拖到振铃后期，比较储存的能量、积分得到的耗散和剩余误差。三维运动显示的是模拟状态，拟合系数仍然只用相机采到的峰值。
你会看到： 耗散的能量正好对上损失的机械能，数值误差一直很小。