4. 基于峰值拾取法的平台加速度频谱初步模态识别
任务对象与数据条件
本节使用某导管架式海洋平台在环境激励下的 6 个加速度测点同步记录（采样频率 50 Hz，单段时长 600 s，共 30000 点）。测点布置为：P1、P2 位于顶层甲板两对角，P3、P4 位于中层甲板，P5、P6 位于导管架腿柱水下段。数据已完成例 2 的去趋势与例 3 的 0.2–10 Hz 带通滤波，本节直接在其基础上做频域初步识别。
峰值拾取法的核心假设是：环境激励近似白噪声、结构阻尼较小、各阶模态在功率谱上表现为可分辨的窄峰。因此处理链为：分段平均周期图 → 功率谱密度（PSD）→ 峰值搜索 → 半功率带宽估阻尼 → 汇总模态清单。为保证频率分辨率，取每段 2048 点、50% 重叠（Hann 窗），则 Δf = 50/2048 ≈ 0.0244 Hz，单测点约 28 段平均，谱估计的随机误差可接受。
功率谱密度计算的关键设置
用 Welch 平均周期图法计算 PSD，参数选择直接影响峰值可分辨性与阻尼估计精度：
- 窗长 N=2048：Δf≈0.0244 Hz，能分辨间距约 0.05 Hz 以上的峰；若平台低阶模态密集（如 1.1 Hz 与 1.2 Hz），需增至 4096 点（Δf≈0.0122 Hz），代价是平均段数减半、谱更粗糙。
- 重叠 50%：在段数与统计独立性之间折中；重叠过高会使各段相关性增强，平均降噪效果被高估。
- Hann 窗：主瓣宽、旁瓣低，适合寻找弱峰；但主瓣展宽会使半功率带宽偏大，阻尼比估计偏高，需在报告中注明窗型影响。
- 单边谱与幅值标定：加速度单位 m/s²，PSD 单位为 (m/s²)²/Hz；若只关心峰值位置，标定不影响频率，但影响半功率点的绝对幅值判读。
峰值拾取与多测点一致性校核
单测点谱上的局部极大不一定是结构模态：可能是波浪低频能量泄漏、设备运转窄带干扰或谱估计随机起伏。因此对每个候选峰执行以下校核：
- 在 6 个测点的 PSD 上分别搜索 0.2–10 Hz 内的局部极大，记录频率与幅值；
- 以 ±1 个频率分辨率（±0.0244 Hz）为容差，统计同一频率附近出现峰的测点数；
- 保留出现在 ≥4 个测点上的候选频率，作为结构模态的初选；
- 对每个初选频率，检查其峰值是否显著高于左右各 5 个频点的均值（建议 ≥3 倍），剔除宽带隆起；
- 检查峰间距：若两峰间距小于 2Δf，视为同一峰的分裂或窗泄漏，合并处理。
按上述流程，本例 6 测点数据得到三个通过校核的低阶频率：f₁≈1.17 Hz、f₂≈2.44 Hz、f₃≈4.30 Hz。其中 f₁ 在 P1、P2、P3、P4 上均清晰，P5、P6 幅值较低但仍有峰，符合整体弯曲模态在甲板测点响应更大的预期；f₃ 仅在 4 个测点通过门限，标记为“待确认”，建议后续用频域分解法复核。
半功率带宽估计阻尼比
对每个确认的峰，在其 PSD 曲线上找到幅值等于峰值一半的两个频率点 f₁′、f₂′，则阻尼比
ζ = (f₂′ − f₁′) / (2 · f_pk)以 f₂ 为例：峰值 2.44 Hz，半功率点约 2.38 Hz 与 2.50 Hz，Δf=0.12 Hz，得 ζ≈2.46%。需要注意三点工程约束：
- 频率量化误差：Δf_res=0.0244 Hz，半功率点落在离散频点上，Δf 的相对误差可达 ±20%，因此 ζ 只宜作为量级估计，不宜作为精确设计值。
- 谱插值：对峰值附近频点做抛物线或对数线性插值，可把半功率点定位精度提高到约 0.005 Hz，显著改善 ζ 的稳定性。
- 窗与平均的影响：Hann 窗主瓣展宽会使 Δf 偏大、ζ 偏高；段数少时谱起伏大，半功率点难以稳定读取。建议同时用不同窗长（1024/2048/4096）计算，观察 ζ 的离散程度，作为不确定度参考。
初步模态参数清单
将上述结果整理为可直接交付后续识别的清单（示例数据，非实测）：
| 阶次 | 频率 f (Hz) | 阻尼比 ζ (%) | 通过测点 | 峰值幅值相对值 | 状态 |
|---|---|---|---|---|---|
| 1 | 1.17 | 1.8 | P1–P6（6/6） | 1.00 | 确认 |
| 2 | 2.44 | 2.5 | P1–P6（6/6） | 0.72 | 确认 |
| 3 | 4.30 | 3.1 | P1–P4（4/6） | 0.35 | 待确认 |
清单中“峰值幅值相对值”以 f₁ 在 P1 的峰值为 1.00 归一化，仅用于比较各阶在测点上的相对强弱，不代表真实振型幅值。阻尼比一栏应随清单附上计算条件（窗长、重叠、平均段数、是否插值），否则不同批次结果无法比较。