基于随机减量法的结构自由衰减响应提取平台

资 源 简 介

本系统旨在实现随机减量法（Random Decrement Method, RDM）的核心算法，用于从结构在环境激励或未知平稳随机激励下的振动响应信号中提取出等效的自由衰减响应信号。该功能通过设定特定的触发条件（如位移触发、过零触发或斜率触发），在原始响应时间序列中搜索并截取满足条件的信号片段。通过对海量截取片段进行同步叠加和平滑平均处理，利用统计学原理使信号中的随机激励成分相互抵消为零，从而剥离出由系统初始状态决定的相干成分，即得到反映系统动力特性的随机减量签名信号。此过程无需获知具体的输入激励数据，能

详 情 说 明

随机减量法结构自由衰减响应提取平台

项目介绍

本平台是一个专注于从结构随机振动响应中提取等效自由衰减信号的算法实现系统。在实际工程中，大型桥梁、高层建筑等结构受环境激励（如风载、地波）的影响，其输入激励往往难以获得。本系统利用随机减量法（RDM），通过对响应信号进行统计平均处理，消除随机激励的成分，保留由系统初值决定的自由振动特征。提取出的“随机减量签名”可直接用于后续的模态频率识别和阻尼比评估，是结构健康监测领域的核心预处理工具。

功能特性

1. 动力响应数值建模：内置单自由度系统仿真功能，通过离散化系统的传递函数模拟结构在白噪声激励下的真实响应过程。
2. 智能化信号触发机制：支持位移阈值触发模式，能够自动捕捉信号中满足特定物理条件的起始时刻（如正向穿过阈值点）。
3. 统计平均去噪：通过对海量提取片段进行同步重叠平均，利用统计相关性原理抵消非相关的随机白噪声及测量噪声。
4. 阻尼性能自动评估：集成基于希尔伯特变换（Hilbert Transform）的响应包络提取算法，并通过对数衰减拟合估计结构的动力学阻尼比。
5. 多维度可视化展示：提供时域响应监测、RD签名信号对比以及频域功率谱密度分析的综合看板。

实现逻辑与算法流程

系统内核严格遵循随机减量法的数学定义，具体流程如下：

1. 结构模拟与信号合成

系统预设了动力学参数（频率、阻尼比），利用双线性变换将连续时间系统转换为离散时间数字滤波器。输入标准高斯白噪声，通过滤波器模拟生成带有物理特性的位移响应，并叠加高斯随机噪声以模拟传感器采集环境。

1. 触发准则判定

采用“位移阈值且斜率为正”的复合触发判定方式。系统计算全局信号的标准差，以标准差的1.2倍作为硬触发阈值。遍历时间序列，当检测到信号从低于阈值跳变至等于或高于阈值时（Up-crossing），判定为一个触发原点。

1. 同步累加平均处理

一旦触发点确立，系统将截取从该点开始的一个固定时间窗口（如4秒）的信号片段。所有被截取的片段将在时域进行同步对齐叠加。随着片段数量增加，服从均值为零分布的随机激励成分逐渐抵消，而具有确定相干性的系统自由衰减波形得以凸显。

1. 模态参数逆向推导

对生成的RD签名信号进行解析处理：

• 希尔伯特变换：计算信号的解析表示，提取其瞬时振幅包络线。
• 线性回归拟合：选取包络线的前端稳定部分，在对数域执行一阶多项式拟合。根据衰减常数与阻尼比的线性关系，计算出系统的估算阻尼比。

1. 频域验证

通过Welch法计算原始随机信号的功率谱，并与RD签名信号的幅频响应进行对比，验证系统在处理过程中是否保持了动力特性（频率一致性）的完整性。

关键技术分析

• c2d转换：在模拟阶段，使用了连续系统到离散系统的映射技术，确保了数值仿真中的频率响应与理论物理解完全匹配。
• 阈值触发逻辑：通过 (x(i) >= threshold) && (x(i-1) < threshold) 逻辑，精确捕捉穿过阈值的瞬时相位，这对于减小RD签名的初始相位偏差至关重要。
• Hilbert包络提取：相比于峰值搜索法，希尔伯特变换提供了更平滑的信号包络，增强了阻尼比计算在噪声环境下的鲁棒性。
• 统计量保证：系统设计了自动统计机制，若有效触发片段数量过少，系统将发出预警，确保由于样本量不足导致的计算误差在可控范围内。

系统要求

• 运行环境：MATLAB R2018b 或更高版本。
• 必备工具箱：Signal Processing Toolbox（用于滤波器设计、希尔伯特变换、功率谱计算等）。
• 硬件建议：4GB以上内存，以支持长序列信号的高速循环检索与海量矩阵平均运算。

使用方法

1. 启动环境：在MATLAB开发环境中配置好工作路径。
2. 配置参数：根据实际需要，在代码初始化区域修改结构频率、阻尼比、采样频率及触发阈值倍数。
3. 执行分析：运行主函数。系统将自动完成信号生成、片段提取、平均计算及后处理分析。
4. 查看报告：分析完成后，系统将自动弹出三栏联动图表，并在命令行窗口打印包含有效片段总数、理论值与估计值对比的结构动力特性报表。
5. 数据导出：处理后的RD签名信号及其对应的时间轴会自动保存至MATLAB基础工作区，变量名分别为rd_output_data和rd_time_axis，方便后续进行深度模态分析。