受激结构振动阻尼系数测定：V(t)曲线蓝色异常区域排查

一、数值计算逻辑排查 #

• 验证Adams-Bashford初始步生成：Adams-Bashford是多步法，必须用欧拉法/龙格-库塔法生成前1~2步的初始位移、速度值，确认初始值符合实验设定，初始步计算无错误。
• 检查微分方程离散化：受激振动的牛顿方程为 m*x'' + c*x' + k*x = F(t)，需转化为一阶方程组：

  x' = v
  v' = (F(t) - c*v - k*x)/m
  

  核对Adams-Bashford迭代公式是否正确应用到该方程组，重点检查阻尼项、刚度项的符号是否与物理模型一致。
• 校验步长合理性：过大的步长会引发数值不稳定，尝试缩小计算步长，观察异常区域是否消失。
• 检查激励函数实现：确认F(t)（如正弦激励）的频率、幅值、起始时间是否匹配实验参数，排查是否存在函数突变导致数值跳变。
• 输出中间结果：打印前20步的x、v计算值，手动代入离散公式验证，定位数值首次出现异常的步骤。

二、绘图代码排查 #

• 匹配数据维度：确认时间数组t与速度数组v的长度完全一致，维度不匹配极易导致绘图填充异常。
• 检查绘图函数调用：若使用matplotlib，排查是否误调用fill_between等填充函数，或plot函数的参数（如linewidth、color）设置错误，导致折线被渲染为填充区域。
• 排查异常数值：检查v数组是否存在NaN、无穷大或突增的极值，可通过print(np.min(v), np.max(v))快速查看数值范围，这类异常值会直接导致绘图显示异常。
• 调整坐标轴范围：手动设置纵轴范围（如plt.ylim([-10, 10])，根据实际数据调整），避免异常值挤压正常曲线显示空间。

内容的提问来源于stack exchange，提问作者Sonpipe