📚 课程目录第4章宽带噪声激励建模

4. 宽带噪声激励建模

在随机振动试验与仿真中，激励信号的质量直接决定了响应分析的可信度。宽带噪声是模拟随机振动环境最常用的激励形式，但“宽带”并非一个单一概念——白噪声、有色噪声和限带噪声在功率谱密度（PSD）形状、时域统计特性和工程适用场景上存在显著差异。理解它们的生成方法以及如何通过数字滤波实现目标谱形，是把随机振动环境“搬进”实验室或仿真模型的关键一步。

三类宽带噪声的本质区别

白噪声是最基础的宽带噪声模型，其特点是功率谱密度在整个频率轴上为常数，即各频率分量的功率均匀分布。理想白噪声总功率无穷大，工程中只能使用有限带宽的“带限白噪声”。白噪声的自相关函数为冲激函数，意味着任意两个不同时刻的采样值互不相关，这一特性使它成为许多随机过程建模的起点。

有色噪声则指功率谱密度随频率变化的噪声，其各频率分量功率不再均匀。实际随机振动环境几乎都是有色的：例如喷气噪声在低频段能量较高，路面激励在特定空间频率上有峰值。有色噪声的生成通常以白噪声为源，通过具有特定频率响应的滤波器整形得到。

限带噪声是频谱被限制在某一频率区间内的噪声，区间外功率近似为零。它可以是带限白噪声（通带内平坦），也可以是带限有色噪声（通带内按目标谱形变化）。在随机振动试验中，限带噪声用于模拟特定频段内的振动能量，例如 20 Hz～2000 Hz 的宽带随机振动。

_out(f) = |H(f)|² · S_in(f)，其中 S_in(f) = S₀（白噪声） • 均方根值：Grms = √(∫ S_out(f) df)，积分区间由限带范围决定 • 滤波器设计：由目标PSD反推 |H(f)| = √(S_target(f) / S₀)，再通过双线性变换得到 H(z)

白噪声的数字化生成

在数字域生成白噪声，最直接的方法是使用伪随机数发生器产生独立同分布的序列。若要求高斯白噪声，可对均匀分布随机数做 Box-Muller 变换或利用中心极限定理叠加多个均匀随机数。设采样率为 f_s，则离散白噪声序列 w[n] 的方差 σ² 与单边功率谱密度 S₀ 的关系为：

σ² = S₀ · f<sub>s</sub> / 2    （双边谱密度折算到单边）

例如：目标 S₀ = 0.01 g²/Hz，f<sub>s</sub> = 5120 Hz
则 σ² = 0.01 × 5120 / 2 = 25.6 g²
σ = 5.06 g

这意味着，若直接以该白噪声作为激励，其总均方根值约为 5.06 g。但实际随机振动试验谱通常不是平坦的，因此需要后续滤波整形。

数字滤波实现有色噪声

将白噪声变为有色噪声的核心是设计一个数字滤波器 H(z)，使其幅频特性满足 |H(f)| = √(S_target(f) / S₀)。工程中常用两类方法：

时域 IIR 滤波法：用双线性变换法设计巴特沃斯、切比雪夫等模拟原型滤波器，再离散化。例如要生成 20 Hz～500 Hz 带限白噪声，可设计一个带通滤波器，通带内增益为 1，阻带衰减不低于 40 dB。以二阶节级联形式实现时，每个二阶节的系数由中心频率和品质因数决定。

频域整形法：先产生白噪声序列，做 FFT 得到频谱，再乘以目标幅值谱 |H(f)|，最后 IFFT 回到时域。这种方法适合任意形状的目标 PSD，但需要注意频域相乘对应时域循环卷积，需通过补零避免混叠。对于 1024 点 FFT，频率分辨率为 f_s/1024，若 f_s=5120 Hz，则 Δf=5 Hz，足以刻画 20 Hz 以上的谱形变化。

限带噪声与随机振动环境的对应

随机振动试验谱通常以 PSD 包络形式给出，例如某机载设备振动谱：

频率范围	PSD 斜率	典型量值
20～100 Hz	+6 dB/oct	0.01→0.04 g²/Hz
100～500 Hz	平坦	0.04 g²/Hz
500～2000 Hz	-6 dB/oct	0.04→0.01 g²/Hz

要生成与该谱对应的时域激励，需构造一个分段整形的滤波器：低频段用一阶高通或二阶高通实现 +6 dB/oct 上升，中频段直通，高频段用一阶低通实现 -6 dB/oct 下降。滤波后信号的 PSD 应通过 Welch 法估计并与目标谱对比，容差通常控制在 ±3 dB 以内。同时计算 Grms 值，若目标谱的 Grms 为 6.5 g，则生成信号的 Grms 应在 6.5 g 附近，偏差过大说明滤波器增益或频带设置有问题。

实践提示：生成有色噪声时，先确认目标 PSD 是单边谱还是双边谱。数字仿真中 FFT 得到的是双边谱，而试验规范通常给单边谱。若混淆两者，幅值会差 √2 倍，导致 Grms 计算错误。

从生成到验证的闭环

完整的宽带噪声激励建模流程包括：确定目标 PSD 与频率范围 → 选择白噪声源并设定采样率 → 设计数字滤波器 → 生成时域序列 → 估计实际 PSD 与 Grms → 与目标对比并迭代修正滤波器参数。其中采样率至少为目标最高频率的 2.5 倍，通常取 2.56 倍以配合 FFT 整周期采样。例如目标最高 2000 Hz，则 f_s 取 5120 Hz。序列长度应保证在最低关心频率上有足够周期数，若最低频率为 20 Hz，则至少需要 10 个周期即 0.5 s，对应 2560 个采样点。

通过这一套方法，白噪声经过数字滤波后即可转化为与随机振动环境谱形一致的有色或限带噪声激励，既可用于振动台控制，也可用于有限元仿真的载荷输入。