网球拍避震器到底“避”掉了什么——一个高中生研究案例

网球拍避震器对击球振动的抑制作用：基于声学模态分析的实验研究及安装位置效应

作者： 一名高中生（实验设计与数据采集；姓名不公开）；数据分析与写作在 AI（Claude）协助下完成
指导： Tian2
版本： 2026-09-24

摘要

目的 网球拍弦间避震器（vibration dampener）被广泛认为能"减震"，但其作用对象究竟是弦床振动还是拍框振动、效果多大、以及安装位置是否重要，在业余实验条件下缺少定量证据。本研究利用手机麦克风记录球拍受落球冲击后的声辐射，通过模态分辨的衰减分析来回答上述问题。

方法 球拍拍柄固定，网球自固定高度自由落至弦床中心，用 Audacity 记录冲击声。数据集 A（2026-04-19）：无避震器 3 次录音（33 次有效冲击）与有避震器 3 次录音（29 次有效冲击）；数据集 B1（2026-08-30）：避震器分别安装在自拍喉起第 1、4、7、10、13 根横弦处，每位置 5 次录音（共 264 次有效冲击）；数据集 B2（2026-09-13）：第 16、19 根横弦，每位置 5 次录音（共 102 次有效冲击）。自写解码器从 .aup3 工程还原 44.1 kHz 浮点音频；自动检测冲击并剔除二次回弹、低信噪比事件；在频谱上识别出四个物理频带（118 Hz 拍框一阶弯曲模态、~270 Hz、~435 Hz、549 Hz 弦床基频模态），对每次冲击计算噪声感知的窄带包络衰减、宽带 Schroeder 衰减、双指数包络模型参数等指标。统计分析同时报告轨级（录音为单位）置换检验、事件级线性混合模型（录音为随机效应）与分层 bootstrap 置信区间；用留一轨交叉验证检验结论的可推广性；位置效应用线性、二次与模态振型 sin² 模型拟合并以交叉验证比较；两个位置录音日之间先用"应与避震器位置无关"的装置指标检查可比性。

结果 (1) 避震器的作用高度模态选择性：549 Hz 弦床模态在无避震器时形成持续到 0.4 s 以上的慢衰减尾（尾段相对拍框模态电平 −14.9 ± 2.7 dB），有避震器时在 0.2 s 内降至峰值以下 60 dB（尾段相对电平 −23.9 ± 4.6 dB）；"弦床模态相对拍框模态的额外衰减"由 −2.0 ± 4.2 dB 升至 16.9 ± 3.3 dB（差值 18.9 dB，95% CI 16.3–21.4；Hedges' g = 4.9）。弦床模态的衰减率由 12.7 ± 6.2 s⁻¹ 升至 31.6 ± 4.2 s⁻¹（比值 2.5，95% CI 2.1–3.0；T60 由 0.62 s 缩短至 0.22 s，损耗因子 η 由 0.007 增至 0.019）。(2) 118 Hz 拍框模态几乎不受影响：峰后 100 ms 衰减量 23.2 vs 22.9 dB，衰减率 11.3 ± 4.6 vs 12.6 ± 3.4 s⁻¹（比值 1.11，CI 0.89–1.35）。(3) 宽带指标随之改变：50 ms 后残余能量占比由 0.231 降至 0.102（比值 0.44，CI 0.33–0.58），早期衰减时间 EDT 由 0.390 s 降至 0.292 s；而宽带双指数模型的慢分量（≈ 9 s⁻¹）在两条件下相同，对应的是拍框模态的长尾而非弦床模态。(4) 留一轨交叉验证下，单一特征"弦床模态额外衰减"即可以 98.4% 的准确率判别有无避震器，弦床模态衰减率单独为 93.5%，四特征逻辑回归为 93.5%。(5) 安装位置：在同一录音日（B1）内，随避震器由拍喉（第 1 根）移向第 13 根，残余能量占比总体下降（0.081 → 0.037，轨级 Spearman ρ = −0.91，p < 0.001），早期谱质心由 311 Hz 升至 478 Hz，冲击峰值升高；线性模型的交叉验证误差优于以第 10 根为中心的 sin² 振型模型。(6) 补充的 B2 数据（第 16、19 根）显示：549 Hz 弦床模态在全部七个位置都被抑制（各位置衰减率轨中位数 24–32 s⁻¹，远离对照组无避震器的 11 s⁻¹、接近有避震器的 32 s⁻¹；属跨录音日比较，仅作定性）；但两个录音日的装置状态明显不同（拍框模态衰减率约 11 vs 24–39 s⁻¹，~270 Hz 模态峰频 263 vs 291–296 Hz），跨日的宽带振铃指标不能直接比较。同一天内，第 19 根比第 16 根的弦床模态衰减率低（25.2 vs 31.8 s⁻¹，轨级精确置换 p = 0.008）、相对拍框的额外衰减少约 6 dB（p = 0.03）、EDT 更长（p = 0.02），即越过弦床中心后效果减弱，方向与 sin² 振型模型一致；允许场次偏移的探索性合并模型中，sin² 模型在残余能量占比、EDT、早期谱质心等部分终点上交叉验证误差最小，但它在本应与位置无关的拍框模态衰减率上同样"最优"，说明录制漂移本身就可能造成这种结果。由于 16、19 同样按顺序录制且其间拍框频率漂移约 3 Hz，这一证据只是提示性的。

结论 避震器通过对弦床基频模态施加强阻尼来消除击球后的"叮"声，而对拍框低频模态无可测影响，与 Brody (1989)、Stroede 等 (1999) 及 Li 等 (2004) 的力学测量一致。避震器装在第 1–19 根任一横弦处都能抑制弦床模态；在同一录音日内，装得越靠近弦床中心（也越靠近击球点），残余振铃越少、冲击声越"脆"，越过中心后（第 19 根 vs 第 16 根）效果有减弱迹象。由于位置实验缺少同日无避震器基线、位置与录制顺序混杂、且两个录音日装置状态不同，位置效应的形状与机制（振型阻尼 vs 改变冲击激励）尚不能完全确定，本文提出了能够桥接不同录音日的改进方案。

关键词 网球拍；避震器；阻尼振动；模态分析；衰减时间；混合效应模型；交叉验证

1 引言

1.1 研究背景

网球击球时，球与弦床的接触约持续 4–5 ms，随后弦床与拍框以各自的固有模态自由振动，振动经拍喉、拍柄传至手臂。弦间避震器是一枚约 2 g 的橡胶纽扣，夹在两根中央主弦之间、通常紧贴最下方的横弦。厂商与球员普遍认为它"减震护臂"，但已有力学研究给出更精细的图景：Brody (1989) 用加速度计测得避震器显著缩短弦振动（约 500–600 Hz）的持续时间，却对拍框振动（100–200 Hz）没有影响；Stroede 等 (1999) 与 Li 等 (2004) 进一步证实避震器不能减少传至手柄和前臂的拍框振动。也就是说，避震器改变的主要是声音与手感，而非可能引发网球肘的低频拍框振动。

对中学生研究者而言，上述结论提示了一个可用低成本设备检验的清晰假设：如果避震器只作用于弦床模态，那么在击球声的频谱中，弦床模态峰应随避震器的安装而消失或加速衰减，而拍框模态峰应保持不变。手机麦克风加自由落球即可实现这一测量。

1.2 研究问题与假设

本研究围绕两个问题：

• Q1：避震器对击球后各振动模态的衰减有何影响？影响多大？
• Q2：避震器沿主弦安装的位置（自拍喉向弦床中心）是否改变其效果？

据此提出可检验的假设：

1.3 本研究的分析策略

与常见的"看波形衰减快慢"不同，本文强调三点：(a) 以物理模态为终点——分频带分析，使结论具有力学含义并自带阴性对照；(b) 增益无关指标优先——不同录音的绝对幅度不可比，事件内的相对量（如弦床模态电平相对拍框模态电平）对麦克风距离、录音增益等混杂免疫；(c) 诚实的统计单位——每个条件只有 3 次录音，轨内的 10 次冲击并非独立样本，因此同时给出轨级检验、混合模型与交叉验证。

2 理论基础

2.1 球拍—弦床系统的模态

冲击后的自由振动可写成若干阻尼简谐模态之和：

$$x(t)=\sum_k A_k\,e^{-\sigma_k t}\cos(2\pi f_k t+\varphi_k),$$

其中 $f_k$ 为第 $k$ 阶模态频率，$\sigma_k$ 为幅度衰减率（s⁻¹）。对手柄固定的球拍，主要模态包括：拍框一阶弯曲模态（约 100–150 Hz）、拍框高阶/环向模态（约 300–500 Hz）以及弦床模态（约 500–700 Hz；Brody, Cross & Lindsey, 2002）。弦床模态由整张弦床像膜一样整体振动形成，其基频近似为

$$f_1\approx\frac{1}{2L}\sqrt{\frac{T}{\mu_{\text{eff}}}},$$

$L$ 为主弦有效长度，$T$ 为张力，$\mu_{\text{eff}}$ 为考虑交织横弦后的等效线密度。

2.2 衰减率、T60 与损耗因子

单一模态的幅度按 $e^{-\sigma t}$ 衰减，能量按 $e^{-2\sigma t}$ 衰减。声学中常用"降 60 dB 所需时间"

$$T_{60}=\frac{\ln 10^{3}}{\sigma}=\frac{6.908}{\sigma},$$

模态损耗因子与品质因数为 $\eta=\sigma/(\pi f)$，$Q=1/\eta$。对含多个模态的宽带信号，Schroeder (1965) 的反向积分法给出能量衰减曲线（EDC）

$$\mathrm{EDC}(t)=\int_t^{\infty}x^2(\tau)\,d\tau,$$

由其 −5…−25 dB 段的斜率外推得 T60，由 0…−10 dB 段得早期衰减时间 EDT（ISO 3382-1）。

2.3 局部阻尼器：附加阻尼与振型平方成正比

在弦上 $x_d$ 处附加一个小的粘性阻尼器（阻尼系数 $c$），对第 $n$ 阶模态（振型 $\phi_n(x)=\sin(n\pi x/L)$，模态质量 $M_n=\mu L/2$）产生的附加衰减率为

$$\Delta\sigma_n=\frac{c\,\phi_n^2(x_d)}{2M_n}=\frac{c}{\mu L}\sin^2\!\left(\frac{n\pi x_d}{L}\right).$$

因此避震器安装在弦端（节点附近）时对基频几乎无效，安装在弦中点（波腹）时效果最大。若以横弦序号 $j$ 表示位置、横弦总数加一为 $N$，则 $x_d/L\approx j/N$。对 16×19 线床 $N\approx20$，预测效果在第 10 根横弦附近达到峰值并对称下降（图 1b）。

2.4 质量加载引起的频移

附加质量 $m$ 使模态频率按微扰理论降低

$$\frac{\Delta f_n}{f_n}\approx-\frac{m\,\phi_n^2(x_d)}{2M_n}.$$

取 $m=2$ g，弦床模态质量 4–8 g，避震器处 $\phi^2\approx0.05$–0.3，得 $\Delta f/f\approx-1\%$ 至 $-4\%$，即 549 Hz 附近下移 5–20 Hz。

2.5 可检验的预测

3 材料与方法

3.1 实验装置与材料

球拍拍柄以夹具固定于桌面（拍面水平、朝上），网球自固定高度经导管释放，落在弦床中心（图 1a；实物见图 2 照片）。避震器为圆形纽扣式橡胶避震器（直径约 25 mm，质量约 2 g），安装于两根中央主弦之间。声音由手机麦克风录制并存为 Audacity 工程（44.1 kHz，单声道，浮点）。数据集 B 中第 1 根线的首条音轨命名为"5cm Trial 1"，提示该批次落球高度为 5 cm；其余参数（落球高度、麦克风距离）在原始记录中未逐轨注明，这一点在第 5.4 节作为局限讨论。

图 1 (a) 装置与避震器安装位置示意（横弦自拍喉向拍头编号）；(b) 局部阻尼器理论：附加阻尼正比于该处模态振型的平方，N≈20 时基频模态在第 10 根横弦处最大。

图 2（球拍结构示意与避震器实物照片）引用了第三方图片，出于版权考虑不在网页中展示。

3.2 数据集

数据集 B1 按 1→4→7→10→13 的顺序连续录制（16:36–17:11），且没有同日的无避震器基线。数据集 B2 于两周后录制，文件保存时间为 17:17（第 16 根）与 17:21（第 19 根），同样按顺序录制、没有无避震器基线；16×19 线床共 19 根横弦，第 19 根为最靠拍头的一根。B2 随 2026-09-22 的压缩包提供，包内另含 B1 的五个文件（SHA-256 与原文件完全一致）以及每个位置文件最后一条音轨的 16 位 WAV 导出；WAV 与对应 .aup3 音轨逐样本相关系数为 1.000，但在 |x| > 1 处被削波，因此不使用。

3.3 原始数据解码与完整性

.aup3 为 SQLite 数据库，样本以小端 float32 分块存储，工程结构以二进制 tokenised XML 记录。自写解码器（s01_extract_aup3.py）按 waveblock.start 拼接样本，校验块内极值与数据库摘要一致，并记录源文件 SHA-256。全部 41 条音轨无削波（无平台段，最大值处样本数 ≤ 2）；样本绝对值超过 1 说明学生在 Audacity 中做过放大，浮点格式下不损失信息。

3.4 冲击事件检测与质控

冲击起点定义为：|x| 峰值 ≥ 噪声底（50 ms 窗 RMS 的 10% 分位）的 100 倍，事件间隔 ≥ 1.5 s，自峰值向前回溯到包络连续 2 ms 低于峰值 5% 处。每个事件取起点前 250–150 ms 作为噪声参考（起点前 150 ms 内会受零相位窄带滤波的预振铃污染，见 3.5）、起点后至多 0.6 s 作为分析窗，并在下列情况截断或剔除（图 3）：

• 二次回弹：20 ms 平滑包络相对前 60 ms 的最小值上升 ≥ 10 dB 且高于峰值 −40 dB（最早可检测到 0.10 s）。回弹在 0.25 s 内的事件剔除，否则在回弹前 10 ms 截断；截断后分析窗短于 0.3 s 的事件亦剔除；
• 峰值出现在起点 30 ms 之后（起点误判）；峰值低于该轨中位峰值的 50%；信噪比 < 35 dB。

共检测到 449 个冲击；其中 1 个（第 19 根第 5 轨的最后一次冲击）距录音结束仅 0.30 s，分析窗不完整，未计算指标；其余 448 个中 428 个通过质控进入分析（B2 计算了指标的 102 个全部通过）。

3.5 指标定义

宽带指标在 60–8000 Hz 带通后的信号上计算，模态指标在 40 Hz 高通后再分频带计算；除峰值幅度外均与录音增益无关。

宽带指标 - $t_{-20}$、$t_{-30}$：5 ms RMS 包络自峰值降至 −20/−30 dB 的时间； - Schroeder T60（−5…−25 dB 外推）与 EDT（0…−10 dB）； - 残余能量占比：冲击后 50 ms 至窗末的能量 / 总能量； - 能量/峰值²：等效持续时间； - 双指数包络模型：$\hat A(t)=A_f e^{-\sigma_f t}+A_s e^{-\sigma_s t}+n$（$n$ 为噪声底），在对数域最小二乘拟合 10 ms–窗末；与单指数模型比较 AICc，$\Delta\mathrm{AICc}>0$ 表示需要慢分量。

模态指标（四个固定频带，4 阶零相位 Butterworth） - 频带：拍框模态 100–135 Hz、~270 Hz（240–300）、~435 Hz（410–460）、弦床模态 530–570 Hz； - 2.5 ms RMS 包络；频带噪声底取起点前 250–150 ms 窗的 RMS——零相位（filtfilt）窄带滤波会把冲击能量向前"拖"约 100 ms（预振铃），若把起点前 150 ms 内也计入噪声窗，噪声底会被高估 40 dB 以上、拟合被过早截断（本研究早期版本即因此几乎无法拟合弦床模态）；峰后 10 ms 起、仅使用高于该频带噪声底 10 dB 的样本做 $\ln A=\ln A_0-\sigma t$ 拟合（≥ 8 个点），得 $\sigma$、T60、$\eta$； - 峰后 100/250 ms 衰减量（dB）；尾段（150–300 ms）电平； - 相对拍框模态电平：弦床频带尾段电平 − 拍框频带尾段电平； - 额外衰减：(早期 10–80 ms 相对电平) − (尾段相对电平)，即弦床模态比拍框模态多衰减了多少 dB，是完全事件内的增益无关量； - 模态频率：20–150 ms 窗零填充 FFT 的抛物线插值峰位。

3.6 模态频带的识别依据

对数据集 A 六条录音分别计算振铃窗（30–300 ms）与尾段窗（100–400 ms）的逐事件中位 PSD（Welch，4096 点，s05b_track_psd.py）。两组在 60–2000 Hz 的逐频点中位电平差仅 0.3–0.6 dB，说明录音增益一致；差异集中在少数峰：549 Hz（尾段无避震器高约 16 dB）、250–290 Hz（尾段中位高约 7 dB）与 430–450 Hz（高 2–4 dB），而 118 Hz 峰在两组几乎重合（图 5）。据此固定上述四个频带。118 Hz 与 Brody 等报告的固定球拍拍框一阶弯曲模态一致，549 Hz 与弦床基频范围一致。

3.7 统计分析

分析单位 每条件 3 次录音、每次约 10 次冲击。轨内冲击共享装置状态，不能视为独立。因此对每个终点同时报告： 1. 轨级：每轨中位数（n = 3 vs 3）的 Welch t 检验与精确置换检验（20 种分配，最小可达 p = 0.10）； 2. 事件级线性混合模型：$y\sim\text{condition}+(1\,|\,\text{track})$（REML；正值指标取对数，估计以比值报告），给出估计值、95% CI 与 p； 3. 分层 bootstrap：先重抽录音、再重抽录音内事件（4000 次），给出均值差与均值比的 95% 百分位区间； 4. 事件级 Hedges' g 作为效应量。

位置效应：趋势统计只在同一录音日（B1）内进行——轨级中位数（每位置 5 个）的 Kruskal–Wallis 检验与 Spearman 趋势；事件级 Spearman 及聚类 bootstrap CI；混合模型 $y\sim\text{position}+(1\,|\,\text{track})$。B2 只有两个位置，用 5 vs 5 轨中位数的精确置换检验（全部 252 种分组，最小可达 p = 0.008）、Hedges' g 与混合模型比较第 16 与第 19 根。合并两个录音日之前，先比较理论上不受避震器位置影响的装置指标（拍框模态衰减率与频率、~270 Hz 模态峰频、噪声底）。显著性水平 0.05，未做多重比较校正，但按预先假设区分主终点（H1–H3）与探索性终点。

3.8 建模与交叉验证

• 判别交叉验证（数据集 A）：留一轨（LOTO）逻辑回归，特征为四个增益无关量（弦床模态 100 ms 衰减量、额外衰减、残余能量占比、EDT，标准化）；同时评估单特征阈值规则的 LOTO 准确率。
• 位置模型（数据集 B1）：在轨级中位数上拟合常数、线性、二次、$a+b\sin^2(\pi j/20)$、$a+b\sin^2(\pi j/N)$（$N$ 自由，14–30）五个模型，用留一轨与留一位置交叉验证的 RMSE 及 AICc 比较。
• 七位置合并模型（B1 + B2，探索性）：每个模型都加一个 B2 场次偏移项（常数 + 场次、线性 + 场次、$\sin^2(\pi j/20)$ + 场次），同样用 AICc、留一轨与留一位置 CV 比较。由于 B2 只有第 16、19 两个位置，场次偏移只能借助假定的曲线形状来估计，结果只作参考。

3.9 软件与可复现性

Python 3.12，NumPy、SciPy、pandas、statsmodels、scikit-learn、Matplotlib。全部分析可以用一条命令从原始录音文件完整重建，两次重建的结果逐字节一致。

4 结果

4.1 数据质量

41 条音轨共检测到 449 个冲击，448 个计算了指标（1 个位于录音末尾），428 个通过质控（图 3）。有效冲击信噪比 36–70 dB（中位 55 dB；B2 为 53–65 dB）。数据集 B1 中检测到 8 次早于 0.25 s 的二次回弹并剔除，另有 7 个事件因截断后分析窗短于 0.3 s 被剔除；A 与 B1 共 5 个事件（A 组 3 个、B1 组 2 个）因峰值过小被剔除；B2 中检测到 4 次二次回弹（0.33–1.37 s，均晚于 0.25 s，仅截断），计算了指标的事件无一被剔除。对照组两条件的噪声底相近，但有避震器录音的冲击峰值约为无避震器的 1.9 倍（6.29 ± 1.84 vs 3.31 ± 0.58 a.u.）；由于拍框模态电平在两组间仅相差 3–5 dB（远小于峰值的 5.6 dB 差异），此差异更可能反映冲击瞬态本身的变化或麦克风相对位置的差异，而非增益，绝对幅度因此只作为次要指标。

图 3 数据质控：(a) 每轨检测到与通过质控的冲击数；(b) 有效冲击的信噪比；(c) 检测到的球二次回弹时间分布。

图 4 (a,b) 对照组原始录音（竖线为检测到的冲击起点）；(c,d) 单次冲击归一化波形；(e,f) 频谱图。无避震器时 118、435、549 Hz 处可见持续的水平亮线，有避震器时 549 Hz 亮线消失。

4.2 避震器敏感模态的识别

图 5 对照组每条录音的中位功率谱：(a) 振铃窗 30–300 ms；(b) 尾段窗 100–400 ms。灰色带为四个分析频带。

图 5 显示，六条录音在 118 Hz 拍框模态处电平一致（尾段 −44 至 −51 dB），而 549 Hz 处形成鲜明对比：在 549 Hz 频点，无避震器三条录音的尾段谱为 −60 至 −63 dB 的尖峰，有避震器三条录音仅为 −75 至 −78 dB 的背景（三轨中位数相差约 16 dB）。250–290 Hz 区域无避震器尾段中位高约 7 dB（有避震器的一条录音例外），430–450 Hz 高 2–4 dB；振铃窗中 330 Hz 附近则有避震器高约 7 dB，可能对应避震器与弦的局部模态。这一模式在三条重复录音间高度一致。

4.3 弦床模态的衰减（H1）与拍框模态的不变性（H2）

图 6 549 Hz 弦床频带（左）与 118 Hz 拍框频带（右）的逐事件包络（dB re 频带峰值）。上排无避震器，下排有避震器。

图 6 是本研究的核心证据。无避震器时（a），549 Hz 频带在冲击后先快速下降约 25 dB，随后以约 −6 dB / 0.1 s 的斜率缓慢衰减，0.45 s 时仍高于 −50 dB；对峰后 10 ms 至噪声底 +10 dB（平均 0.33 s）的整段作指数拟合，衰减率 σ = 12.7 ± 6.2 s⁻¹（各轨中位数 9.2、12.5、11.2 s⁻¹），T60 = 0.62 ± 0.18 s，损耗因子 η = 0.007。有避震器时（c），同一频带在 0.2 s 内即降至 −60 dB 以下，拟合段（平均 0.19 s，跨度 54 dB）给出 σ = 31.6 ± 4.2 s⁻¹（各轨中位数 32.2、33.0、32.5 s⁻¹），T60 = 0.22 ± 0.03 s，η = 0.019——衰减率提高 2.5 倍（bootstrap 95% CI 2.1–3.0）。相反，118 Hz 拍框频带（b, d）的包络在两条件下几乎重合：峰后 100/250 ms 的衰减量为 23.2/35.2 dB（无）与 22.9/35.7 dB（有），整段拟合的 σ 为 11.3 ± 4.6 与 12.6 ± 3.4 s⁻¹。

表 1 汇总主要终点（完整表见附录 A）。

表 1 对照组主要终点（均值 ± SD；混合模型估计；分层 bootstrap 95% CI）

注：轨级置换检验在 3 vs 3 设计下最小可达 p = 0.10，达到该值表示三条有避震器录音的中位数全部落在三条无避震器录音之外。

图 7 主要终点指标的事件级分布与各轨中位数。(a) 弦床模态额外衰减；(b) 尾段弦床模态相对拍框电平；(c) 残余能量占比（对数轴）；(d) EDT；(e) 双指数慢分量衰减率（两条件相同）；(f) 拍框模态衰减率（阴性对照）。

H1 得到强支持：三种独立定义的弦床模态指标一致，且三条有避震器录音的中位数全部落在三条无避震器录音的范围之外（图 7a–b），轨内离散远小于组间差。H2 亦得到支持：拍框模态的固定窗衰减量相同，整段拟合衰减率的比值 CI（0.89–1.35）虽略偏向有避震器更快（混合模型 p = 0.04），但与弦床模态 2.5 倍的变化相比可以忽略。

4.4 宽带衰减与双指数模型（H3）

图 8 宽带双指数包络模型。(a) 无避震器示例：快分量 σ_f ≈ 33 s⁻¹ 与慢分量 σ_s ≈ 10 s⁻¹ 叠加；(b) 有避震器示例：同样存在 σ_s ≈ 7–10 s⁻¹ 的慢分量；(c) 全部有效冲击的 (σ_s, 慢分量占比) 分布，两条件重叠。

宽带包络在两条件下都需要双指数描述（ΔAICc > 0 的事件占 94% 与 86%），且慢分量的衰减率相同（8.6 vs 10.1 s⁻¹，比值 1.17，CI 0.83–1.60）。对照图 6b/d 可知，这一慢分量是 118 Hz 拍框模态在 −35 dB 附近的长尾，避震器不改变它；549 Hz 弦床模态虽然在无避震器时也持续到 0.4 s，但其尾段电平比拍框模态低约 15 dB，宽带模型分不出这两个衰减率相近（≈ 9 与 ≈ 13 s⁻¹）的成分。换言之，避震器的效应在宽带包络的形状上不明显，必须分频带才能看到。尽管如此，弦床模态在 50–300 ms 内的能量足以改变宽带能量指标：残余能量占比减半、EDT 缩短 25%、t₋₂₀ 与 t₋₃₀ 缩短约 20%。值得注意的是，Schroeder T60（−5…−25 dB 段）对避震器不敏感（比值 0.95，CI 含 1），因为该段被快速衰减的冲击瞬态主导——这说明在冲击类信号中，EDT 与残余能量比 T60 更适合作为"振铃"指标。

4.5 频率偏移（H4）

弦床模态频率有避震器时低 7.3 Hz（−1.3%），方向与质量加载预测一致，但 bootstrap CI（−18.8 至 +4.1 Hz）包含零；且有避震器时该频带内已无清晰模态峰，频率估计主要反映宽带背景，故 H4 只能视为"不矛盾"。另一方面，拍框模态频率在两组间相差 1.6 Hz（1.4%，混合模型 p ≈ 1×10⁻⁸）。2 g 的避震器不可能使拍框频率升高，这一差异更可能来自两组录音之间的夹持状态差异，是一个提示"两组录音装置不完全相同"的内部指标（见 5.4）。

4.6 交叉验证

图 9 (a) 留一轨交叉验证：以其余 5 条录音训练的逻辑回归在留出录音上的逐冲击判别准确率；(b) 单特征阈值规则的留一轨准确率。

四特征逻辑回归的留一轨合并准确率为 93.5%（58/62），六条留出录音的准确率分别为 9/9、10/10、8/10、11/11、10/11、10/11。单一特征"弦床模态额外衰减"即达 98.4%（61/62），弦床模态衰减率 σ 为 93.5%，均远高于宽带指标（EDT 80.6%、残余能量占比 79.0%、T60 75.8%）。这表明：(i) 避震器效应可推广到未见过的录音而非某一录音的特例；(ii) 模态分辨指标的判别力显著优于宽带指标。

4.7 安装位置效应（B1：2026-08-30，第 1–13 根）

图 10 安装位置（自拍喉起的横弦序号）与各终点的关系。点=单次冲击，黑线=各轨中位数（每位置 5 轨），实线=线性拟合，虚线=以第 10 根为中心的 sin² 振型模型，两者都只用 B1（第 1–13 根）拟合，浅灰线为其向第 13 根以外的外推；标题给出 B1 的轨级 Spearman ρ、Kruskal–Wallis p 与混合模型斜率 p。紫色底纹为另一录音日（B2）的第 16、19 根，右上角为两者之间的轨级精确置换 p。B2 与 B1 的绝对水平差异主要反映装置状态不同（见 4.9），不应按横轴连成一条曲线来读。

表 2 位置组 B1 各终点的位置均值（每位置 5 个轨中位数的均值）与趋势统计（完整表见附录 C）

位置效应呈现清晰而有选择性的模式：

1. 振铃减少：残余能量占比从第 1 根到第 13 根下降 54%（0.081 → 0.037），第 7、10、13 根之间五条轨中位数完全不重叠，第 1–7 根的相邻位置之间有重叠（图 10a）；EDT 下降 19%。
2. 冲击变"脆"：早期（10–80 ms）谱质心升高 54%，700–2000 Hz 能量占比升高 2.5 倍，60–250 Hz 占比减半，冲击峰值升高 53%（图 10c、d 及图 11）。这些是激励谱的变化，而非衰减率的变化。
3. 宽带衰减率不变，但拍框模态漂移：宽带包络衰减率与 t₋₃₀ 无趋势；拍框模态的衰减率却由 22 升至 39 s⁻¹（图 10f），与其频率的 +2.3 Hz 漂移同步——夹持得更紧会同时提高拍框频率和阻尼，这更像录制过程中夹具状态的变化，而非避震器位置的作用（2 g 的避震器无法改变拍框模态）。~435 Hz 模态的衰减率随位置由 29 降至 16 s⁻¹，与其尾段相对电平的升高一致。
4. 435 Hz 模态尾段相对电平随位置升高并趋于饱和（−10.4 → −1.5 dB，图 10e）。

图 11 (a,b) 各位置的中位归一化 PSD（以 0–50 ms 冲击能量为参考）；(c,d) 逐频点电平与位置的 Spearman 秩相关。80–300 Hz 成分随位置向中心移动而降低，550–600 Hz 及 700 Hz 以上成分升高。

4.8 位置模型比较（H5，B1 场次内）

表 3 轨级中位数上的模型比较（LOTO = 留一轨 CV-RMSE，LOPO = 留一位置 CV-RMSE；粗体为该列最小值）

对振铃类指标（残余能量占比、EDT），以第 10 根横弦为中心的 sin² 模型交叉验证误差明显大于线性模型；令 N 自由时拟合收敛到搜索边界 N = 30，即数据在 1–13 范围内单调、无饱和，与 H5 预测的"中心处饱和、对称下降"不符。仅 435 Hz 模态相对电平与 60–250 Hz 占比表现出饱和特征，sin² 模型与二次模型相当。需要注意，1–13 只覆盖了 sin² 曲线的上升段和顶点附近，线性与 sin² 两种形状在这一段的区别主要落在第 10、13 根；真正能区分二者的是第 10 根之后的位置——这正是 B2 补充第 16、19 根的意义（4.9）。

4.9 第二个位置场次（B2：2026-09-13，第 16、19 根）

(1) 两个录音日不能直接合并。 表 4 列出了理论上与避震器安装位置无关的装置指标。B2 的拍框模态衰减率约 11 s⁻¹，只有 B1 的一半以下（24–39 s⁻¹），而与对照组（2026-04-19）的 10–12 s⁻¹ 相同；~270 Hz 模态峰频在 B2 的中位数为 262–264 Hz（第 19 根有一条录音为 295 Hz），比 B1（291–296 Hz）和对照组（296–303 Hz）低约 30 Hz；B2 的噪声底也低于 B1。2 g 的避震器既不能让拍框阻尼减半，也不能让一个模态下移 30 Hz，这些差异只能来自装置状态（夹持方式与力度、弦的张力或状态、麦克风位置等）的改变。因此，B2 与 B1 在残余能量占比、EDT、60–250 Hz 占比等宽带指标上的巨大差异（图 10 紫色区域）主要是拍框模态在 B2 振得更久的结果，而不是第 16、19 根位置本身的效应。

表 4 各录音组的装置一致性指标（各轨中位数的中位数；方括号为轨间范围）

(2) 在所有七个位置，弦床模态都被抑制。 弦床模态自身的衰减率（只看 530–570 Hz 频带，不以拍框为参照）在第 1–19 根各位置的轨中位数为 24–32 s⁻¹，全部远离对照组无避震器的 11 s⁻¹、接近有避震器的 32 s⁻¹（图 12c；这是跨录音日比较，只能定性地看——各位置之间仍有差异，例如同日第 19 根显著低于第 16 根，见下文）；峰后 250 ms 衰减量 54–62 dB，也与有避震器（67 dB）接近、远大于无避震器（35 dB）。也就是说，无论装在哪根横弦处，避震器都能消除 549 Hz 的"叮"声；位置只改变次级效应。

(3) 同一天内，第 19 根的效果弱于第 16 根。 表 5 比较 B2 内两个位置（各 5 条录音）。第 19 根的弦床模态衰减率更低（25.2 vs 31.8 s⁻¹，g = −5.0，精确置换 p = 0.008，已是 5 vs 5 设计能达到的最小值）、峰后 100 ms 衰减量少 4 dB、相对拍框的额外衰减少约 6 dB（17.8 → 11.6 dB），宽带上 EDT 更长（0.33 → 0.42 s）、早期谱质心更低（226 → 170 Hz），尾段 ~435 Hz 相对电平更低。按局部阻尼器理论，第 16 与第 19 根处的振型平方 sin²(πj/20) 分别为 0.35 与 0.02，第 19 根几乎在弦的端点，预测附加阻尼明显减小——观测方向与此一致，而与"越往拍头越好"的线性外推相反。

表 5 B2 场次内第 19 根 vs 第 16 根（完整表见附录 E）

但这一对比同样有顺序混杂：两个文件相隔 4 分钟按 16 → 19 的顺序录制，其间拍框模态频率下降了约 3 Hz（轨中位数 115.4 → 112.4 Hz，p = 0.008），说明装置在同一天内也在漂移。冲击峰值在 B1（1 → 13）和 B2（16 → 19）中都随录制先后上升（B2：1.40 → 1.67，p = 0.016），很可能主要反映录制过程中的漂移，而不是避震器位置。

(4) 含场次偏移的七位置模型（探索性）。 表 6 把 B1 与 B2 合在一起，给每个模型加一个"B2 场次偏移"。在残余能量占比、EDT、早期谱质心、~435 Hz 尾段相对电平与 60–250 Hz 占比上，sin²(πj/20) + 场次模型的 AICc、留一轨与留一位置 CV 误差都是三个模型中最小的——与只用 B1 时线性模型占优的结论相反；而冲击峰值、700–2000 Hz 占比和尾段弦床模态相对拍框电平上线性 + 场次模型较好，弦床模态额外衰减则没有一致的赢家（AICc 与留一轨偏向线性 + 场次，留一位置偏向常数 + 场次）。一个重要的反证：理论上与避震器位置无关的拍框模态衰减率，在合并模型中同样是 sin² + 场次"最优"（AICc 100.6 vs 线性 + 场次 106.5）——B1 中拍框阻尼随录制先升后降的漂移恰好呈"中间高"的形状。也就是说，录制漂移本身就能让 sin² 模型胜出。再加上场次偏移只能靠假定的形状来估计（B2 只有两个位置），这一结果只能说明"数据与振型模型不矛盾，且不支持效果随位置单调增强"，不能作为 H5 成立的证据。

表 6 七个位置、含场次偏移的模型比较（轨级中位数；粗体为该列最小值；最后一列为只用 B1 时的对照；完整表见附录 F）

图 12 三个录音日的装置一致性检查与 B2 场次内的位置对比。(a) 拍框模态衰减率与 (b) ~270 Hz 模态峰频理论上与避震器位置无关，却在录音日之间明显不同；(c) 弦床模态衰减率在所有位置都与对照组"有避震器"相当；(d–f) 第 16 与第 19 根的对比，右上角为轨级精确置换 p。点=单次冲击，黑线=轨中位数。

5 讨论

5.1 主要发现与理论、文献的对照

本研究以声学方法重现了力学测量的经典结论：避震器强烈抑制弦床模态、不改变拍框模态（Brody, 1989；Stroede 等, 1999；Li 等, 2004）。三条相互独立的证据链支持这一点——(i) 频谱上 549 Hz 峰的消失（图 5）；(ii) 窄带包络显示弦床模态衰减率由 12.7 升至 31.6 s⁻¹、T60 由 0.62 s 缩短至 0.22 s（图 6）；(iii) 事件内相对量"弦床模态相对拍框模态的额外衰减"增加 18.9 dB（图 7）。留一轨交叉验证表明该效应在每一条未见过的录音上都能被识别，不依赖任何单条录音。

与实验方案中"T60 由约 1.1 s 缩短至 0.15 s（7 倍）"的预期相比，本研究测得的弦床模态 T60 由 0.62 s 缩短至 0.22 s（2.8 倍），损耗因子由 0.0074 增至 0.0187（2.5 倍）；方向一致，倍数小于预期。差异来源包括：本实验为拍柄固定的落球激励而非拨弦，且以整体弦床模态而非单根弦为对象。

拍框模态频率 113–117 Hz、衰减率在对照组两条件下基本一致（11.3 vs 12.6 s⁻¹，T60 ≈ 0.6 s，η ≈ 0.03；固定窗衰减量相同），说明夹持球拍的一阶弯曲模态由夹具与拍框材料共同决定，避震器（≈ 拍框质量的 0.7%）对其没有实质影响（衰减率比值 1.11，CI 0.89–1.35；混合模型 p = 0.04，效应远小于弦床模态）。位置组中拍框模态的频率与衰减率随录制顺序同步漂移（5.4 节），进一步说明拍框模态受夹持状态而非避震器支配。这与"避震器不能减少传至手臂的低频振动"的力学结论一致，因此避震器的价值主要在声音与主观手感，而不应期望它预防网球肘。

5.2 为什么宽带 T60 不敏感而 EDT 敏感

Schroeder T60 由 EDC 的 −5…−25 dB 段外推，而冲击类信号的前 25 dB 能量衰减几乎完全由 0–50 ms 的冲击瞬态决定，两条件下的快分量衰减率相近（≈ 30 s⁻¹）。弦床频带（450–700 Hz）的贡献在无避震器时仅占总能量约 7%，只影响 EDC 的后段。EDT 与残余能量占比恰好捕捉这一后段，因此更适用于评价"振铃"。这一方法学结论对同类学生实验有直接指导意义：不要只报告一个 T60。

5.3 安装位置效应的两种机制

B1 数据支持"避震器越靠近弦床中心，击球后残余振铃越少、冲击声越明亮"，B2 数据进一步提示越过中心后效果减弱；但对形状与机制的判断都需要谨慎：

• 振型阻尼机制（H5） 预测效果正比于 sin²(πj/N)，在中心处饱和。这与 435 Hz 模态相对电平和低频占比的饱和趋势相符，却与残余能量占比、EDT 的单调下降不符（表 3）。
• 激励改变机制 落球点在弦床中心，避震器移向中心即移向冲击点。橡胶块处于接触区附近会改变局部接触刚度与有效质量，使冲击更短、更尖锐——这正是早期谱质心、高频占比与冲击峰值随位置升高的含义。振铃能量占比的下降可能主要是分母（冲击能量）增大而非分子减小。宽带包络衰减率与拍框衰减率均不随位置变化，也支持"变化在激励而不在阻尼"。

两种机制并不互斥；在 549 Hz 弦床模态已被彻底抑制（七个位置下弦床模态衰减率都与对照组有避震器相当，4.9）的前提下，位置实验实际观察到的是次级效应。

B2 的第 16、19 根带来一条新信息，但它区分的不是这两种机制：落球点在弦床中心附近，振型阻尼（峰在第 10 根）和激励改变（峰在击球点）都预测"越过中心后效果减弱"，第 19 根弱于第 16 根与二者都相容。它真正排除的是第三种解释——B1 中的单调趋势只是录制顺序上的装置漂移或"越往拍头越好"的单调关系：若果真如此，B2 中按顺序录制的第 19 根应当延续 B1 的方向（振铃更少），实际却相反。要区分振型阻尼与激励改变，仍需要在同一场次中改变落球点，并加入无避震器基线。

5.4 局限性与混杂因素

1. 非交错的组间比较。数据集 A 的两条件为分别连续录制的 3 条录音，任何在两组之间发生的装置变化（夹持、麦克风位置、落球高度）都与避震器混杂。拍框模态频率的 1.4% 差异与冲击峰值的 1.9 倍差异提示确有此类差异。所幸主终点是事件内的相对量（弦床模态相对拍框模态），对增益与距离免疫；且拍框模态衰减率的不变性表明夹持阻尼没有系统变化。
2. 位置与录制顺序完全混杂。位置组按 1→13 的顺序录制，期间噪声底自第 7 根第 4 轨起下降、至第 13 根前三轨又回到原水平，拍框模态频率净漂移 +2.3 Hz（非严格单调）。虽然主要位置指标为增益无关量，仍不能排除落球高度或麦克风位置的缓慢漂移贡献了部分趋势。
3. 无同日基线、无重复安装。位置组缺少无避震器基线，且每个位置只安装一次，避震器按压深度等安装变异无法评估。
4. 样本层级。每条件 3 条录音使轨级检验的最小 p 为 0.10；事件级混合模型的极小 p 值以录音内冲击可交换为前提，应作为效应方向与大小的度量而非严格的假设检验证据。
5. 声学测量的间接性。麦克风记录的是声辐射，不同模态的辐射效率不同；本文不对模态间的绝对能量做比较，只比较同一模态在不同条件下的衰减。
6. 实验参数记录不完整。落球高度、麦克风距离与球拍型号/张力未逐轨记录，限制了与文献的定量对照。
7. 不同录音日的装置状态不同。三个录音日的拍框模态衰减率（约 11、24–39、11 s⁻¹）与 ~270 Hz 模态峰频（约 300、293、263 Hz）明显不同，而这些量与避震器位置无关；不同录音日之间也没有共同的参考条件（同一位置或无避震器基线）可以"桥接"。因此 B1 与 B2 只能分别在场次内比较，跨日合并的模型只有探索意义。B2 仅两个位置，同样按顺序录制，场次内也有约 3 Hz 的拍框频率漂移。

5.5 对后续实验的建议

• 采用 ABAB 交错设计：无→有→无→有各 ≥ 3 条录音，位置实验在每个位置前后各录一条无避震器基线，位置顺序随机化。
• 每个录音日都录同一组参考条件（例如无避震器 + 第 10 根），用来把不同录音日的数据桥接到同一标尺上；本研究 B1 与 B2 之间正是缺少这一环，才无法合并。
• 位置覆盖整根主弦（例如第 1、4、7、10、13、16、19 根）并在同一天内随机顺序完成，直接检验 sin² 曲线越过中心后的下降段。
• 固定并记录落球高度（导管）、麦克风距离（用尺）、拍框夹持力矩；每条录音前录 5 s 静音用于噪声底校准。
• 增加拍柄加速度计（手机 phyphox 或 MPU-6050）与声学同步，直接验证"拍框振动不变"。
• 位置实验同时改变落球点（中心 vs 偏喉 vs 偏头），以区分振型阻尼与激励改变两种机制。
• 分析时沿用本文的模态分辨指标（弦床模态额外衰减、尾段相对电平）作为主终点。

6 结论

1. 弦间避震器使弦床基频模态（549 Hz）的衰减率由 12.7 升至 31.6 s⁻¹（比值 2.5，95% CI 2.1–3.0；T60 由 0.62 s 缩短至 0.22 s），其相对拍框模态的额外衰减增加 18.9 dB（95% CI 16.3–21.4），效应量极大（g ≈ 3.5–4.9），并在留一轨交叉验证中以 93–98% 的准确率区分有无避震器。
2. 拍框一阶弯曲模态（118 Hz）的衰减基本不受避震器影响（固定窗衰减量相同；衰减率比值 1.11，CI 0.89–1.35），证实避震器不改变传向手臂的低频拍框振动；两组间 1.6 Hz 的拍框频率差异方向与质量加载相反，归因于两组录音夹持状态不同（5.4 节）。
3. 宽带上，避震器使 50 ms 后残余能量占比减半、EDT 缩短 25%；而 −5…−25 dB 段的 Schroeder T60 与宽带双指数模型的慢分量（由拍框长尾主导）几乎不变，说明评价击球振铃应使用 EDT、残余能量占比，最好是模态分辨指标。
4. 避震器装在第 1 到第 19 根横弦的任一处，弦床模态的衰减率（24–32 s⁻¹）都远离对照组无避震器的 11 s⁻¹、接近有避震器的 32 s⁻¹（跨录音日比较，定性），即安装位置不决定"有没有效果"，只影响次级效应的大小。
5. 同一录音日内，避震器由拍喉移向第 13 根时残余振铃总体减少（第 13 根较第 1 根低 54%，轨级 ρ = −0.91），冲击声更明亮，在 1–13 范围内线性模型优于振型 sin² 模型；另一录音日的第 19 根则比第 16 根弦床模态衰减慢、额外衰减少约 6 dB、EDT 更长，即越过弦床中心后效果减弱；探索性的合并模型在部分终点上偏向 sin²，但同样的模型也会被录制漂移"骗"出 sin² 形状。由于两个录音日的装置状态明显不同、每个场次内位置都与录制顺序混杂，位置效应的准确形状以及它来自振型阻尼还是冲击激励改变，仍需在同一天、随机顺序、带参考条件的实验中确认。
6. 方法学上，本研究建立了一条从 Audacity 工程解码、事件检测与质控、噪声感知的模态衰减估计、到分层统计与交叉验证的完整可复现流程，适合作为中学生声学实验的分析范式。

参考文献

1. Brody H. (1989). Vibration damping of tennis rackets. International Journal of Sport Biomechanics, 5(4), 451–456.
2. Stroede C. L., Noble L., Walker H. S. (1999). The effect of tennis racket string vibration dampers on racket handle vibrations and discomfort following impacts. Journal of Sports Sciences, 17(5), 379–385.
3. Li F.-X., Fewtrell D., Jenkins M. (2004). String vibration dampers do not reduce racket frame vibration transfer to the forearm. Journal of Sports Sciences, 22(11–12), 1041–1052.
4. Brody H., Cross R., Lindsey C. (2002). The Physics and Technology of Tennis. Racquet Tech Publishing.
5. Cross R. (1998). The sweet spots of a tennis racket. Sports Engineering, 1(2), 63–78.
6. Schroeder M. R. (1965). New method of measuring reverberation time. Journal of the Acoustical Society of America, 37(3), 409–412.
7. ISO 3382-1:2009. Acoustics — Measurement of room acoustic parameters — Part 1: Performance spaces.
8. Ewins D. J. (2000). Modal Testing: Theory, Practice and Application (2nd ed.). Research Studies Press.
9. Hedges L. V. (1981). Distribution theory for Glass's estimator of effect size and related estimators. Journal of Educational Statistics, 6(2), 107–128.
10. Seabold S., Perktold J. (2010). statsmodels: Econometric and statistical modeling with Python. Proceedings of the 9th Python in Science Conference, 92–96.