一种基于细化谱分析的加速度计动态模型结构参数辨识的方法
【专利摘要】本发明公开了一种基于细化谱分析的加速度计动态模型结构参数辨识的方法，属于信息技术信号处理及测控领域。通过加速度计绝对法冲击校准实验获取加速度计输入输出数据及数据截取；对所获得的输入输出数据进行离散傅里叶变换，确定加速度计频率响应估计；对频率响应估计进行分析，得到谐振频率一步估计值；利用频谱细化方法对谐振频率及零频附近区段细化，得到谐振幅值谐振频率及零频幅值二步估计值；利用模型参数与特征点坐标二步估计值关系，确定加速度计动态模型结构参数估计，实现基于细化谱分析的加速度计动态模型结构参数辨识。本方法简化了加速度计二阶动态模型参数辨识的过程，提高了频域内对加速度计进行建模和辨识的精度及计算效率。
【专利说明】
一种基于细化谱分析的加速度计动态模型结构参数辨识的 方法
技术领域
[0001] 本发明涉及一种传感器动态模型参数辨识方法，属于信息技术信号处理及测控领 域，尤其涉及一种基于细化谱分析的加速度计动态模型结构参数辨识的方法。
【背景技术】
[0002] 加速度计是一种被广泛应用于测量载体(如机械设备或建筑结构等)所受冲击或 振动的传感器，为振动分析及故障检测的数据来源。
[0003]在线性动态范围内，加速度计通常等效为一个单自由度质量弹簧阻尼系统，质量 块内置于加速度计壳体内被弹簧及阻尼器支撑。质量块作用在加速度计压电元件上，压电 元件产生正比于质量块相对于壳体位移的电荷量输出。对加速度计输出加速度的测定，是 将该加速度计与电荷放大器相连，加速度计的输出作为电荷放大器的输入，测量信号实际 为被放大的加速度计的输出。
[0004] 利用加速度计的输入输出信号及离散傅里叶变换(DFT)方法所获得的加速度计频 率响应非参数模型，具有一定精度，但谐振频率及对应的谐振幅值受频率分辨率限制，误差 较大。当采样频率和采样数据长度一定时，采用上述方法无法通过提高频率分辨率减小谐 振频率及谐振幅值的误差。
【发明内容】
[0005] 本发明以提高加速度计模型参数辨识精度及计算效率为目的。在直接傅里叶变换 求取加速度计频率响应估计的基础上，通过细化频响曲线谐振频率及频率为零的频率点 (以下简称为零频)附近区域，仅通过两个特征点坐标值，经简单计算，便可得到加速度计动 态模型结构参数，实现基于单自由度质量弹簧阻尼系统的加速度计动态模型参数辨识，具 有较强抗噪性能，且计算过程快速、稳定。
[0006] -种基于细化谱分析的加速度计动态模型结构参数辨识的方法，包括以下步骤：
[0007] 步骤一:通过加速度计绝对法冲击校准实验获取加速度计输入输出数据及数据截 取；
[0008] 步骤二:对步骤一所获得的输入输出数据进行离散傅里叶变换，确定加速度计频 率响应估计。对频率响应估计进行分析，得到谐振频率一步估计值；
[0009] 步骤三:利用频谱细化方法对谐振频率及零频附近区段细化，得到谐振幅值谐振 频率及零频幅值二步估计值；
[0010] 步骤四：利用模型参数与特征点坐标二步估计值关系，确定加速度计动态模型结 构参数估计，实现基于细化谱分析的加速度计动态模型结构参数辨识。
[0011] 所述步骤一，具体包括：
[0012] 基于加速度计绝对法冲击校准实验，利用数据采集软硬件采集冲击激励加速度及 加速度计响应的时间序列。加速度计绝对法冲击校准实验的装置是依据IS016063-13绝对 法冲击校准标准建立的冲击激励系统，要求冲击激励加速度信号频率响应范围能够覆盖加 速度计谐振频率。数据采集系统的采样频率为fs。
[0013] 将得到的时间序列进行截取，去掉数据中激励作用之前的部分与响应结束之后的 部分，以保留信噪比较大的数据部分，即保留冲击激励有效作用的数据部分，用于后续处理 计算。
[0014] 所述步骤二，具体包括：
[0015] 利用步骤一所获得的输入加速度时间序列和加速度计输出时间序列，分别记为u (k)和7(1〇沽=0，1，2，...少-14为时间序列长度汰为离散数据序号。利用公式（1-1)确定 输入加速度离散频谱U(1 )和加速度计输出离散频谱Y( 1 )，1为离散谱线序号，
。利用公式(1-2)确定加速度计频率响应估计(^/q)。
[0019] 其中，反为U的共辄。
[0020] 利用公式（1-3)确定A f，确定横坐标序列
；确定幅频特性纵坐标序列 10 x )))，确定相频特性纵坐标序列§/#〇,)) x 180/;r，其中，))为 频率响应估计的模，为频率响应估计的相角，所述相角单位为弧度。绘出加速 度计频率响应估计曲线，确定幅值峰值点横坐标值，记为氣，则4即为加速度计谐振频率一 步估计值。
(1-3)
[0022]所述步骤三，具体包括：
[0023] 以4为中心频率，确定谐振频率附近频率响应特征区域+/L4/1，A为调 整因子，取值为2~5;利用公式(1-3)确定A f。
[0024] 确定细化段数为L，确定细化离散频率序列
[0026]利用公式（1-4)分别确定细化后加速度计输入输出序列的离散傅里叶变换Xx(f)、 Yx(f)。利用公式(1-5)确定谐振频率附近区段细化后的加速度计频率响应估计)。
[0029] 其中，R为Ux的共辄。
[0030] 利用l〇x丨确定谐振频率附近区段频率响应幅值序列，搜索其中最大 值以及该最大值所对应的频率值，即为谐振幅值及谐振频率二步估计值，分别记为im与成
[0031] 然后，采样相同方法对零频附近区段进行细化，确定零频附近区段频率响应幅值 序列，搜索零频所对应的频率响应幅值，即为零频幅值二步估计值，记为歲0_)、
[0032] 所述步骤四，具体包括：
[0033] 根据步骤三获得的谐振幅值、谐振频率和零频幅值的二步估计值，利用模型参数 与特征点坐标值二步估计值关系公式(1-6)确定对应二阶传递函数模型参数。
(1-6)
[0035]将参数K、con、S代入公式（1-7)，得到等效为单自由度质量弹簧阻尼系统的加速度 计二阶传递函数模型。
(1-7)
[0037] 其中，K、《"和5分别表示加速度计二阶模型静态传递系数、固有频率和阻尼比。
[0038] 本发明的优点：利用细化谱分析细化加速度计频响曲线谐振频率及零频附近区 域，得到了精度较高的两个特征点坐标值，提高了谐振频率、谐振幅值及零频幅值的估计精 度，降低了加速度计动态模型参数的估计误差，提高了参数辨识精度。同时，仅利用两个特 征点坐标值参与计算，提高了辨识过程抗噪性能，计算过程简单、稳定。
【附图说明】
[0039] 图1是本发明所述的一种基于细化谱分析的加速度计动态模型参数辨识的方法的 流程图；
[0040] 图2是本发明所述加速度计绝对法冲击校准实验中典型冲击激励加速度波形图； [0041 ]图3是【具体实施方式】所述加速度计频率响应；
[0042]图4是相同激励下本发明所述方法得到的模型输出与实际输出的波形图；
[0043]图5是图4波形的局部放大。
【具体实施方式】
[0044] 下面结合实例及附图对本发明作进一步的描述，需要说明的是，实施例并不限定 本发明要求保护的范围。
【具体实施方式】 [0045] ：利用依据ISO 16063-13绝对法冲击校准标准建立的冲击激励装 置，该装置典型输入加速度信号波形如图2所示。本实施例以10MHz采样频率对冲击激励加 速度信号和加速度计响应信号进行同步采集，每个通道共采集数据样本1.2 X 105个。冲击 激励加速度峰值约为3.83X 102g，其中，g = 9.81m/s2。为提高信噪比，去掉数据中激励作用 之前的部分与响应结束之后的部分，冲击激励有效作用的数据共51000个点。
[0046] 按照公式（1-1)分别对冲击激励加速度时间序列和加速度计响应时间序列进行 51000个点的快速傅里叶变换，将零频分量移至频谱中心，按照公式（1-2)确定加速度计频 率响应估计(》(_/以.）。利用公式（1-3)确定A f = 1.9608 X 102Hz，分别确定幅频响应及相频响 应的横纵坐标序列，绘出加速度计频率响应估计曲线，如图3所示。
[0047]相频曲线随着频率的增加，相角由0°快速降低至-180°，在这一变化过程对应的频 率范围内，确定幅频特性曲线的幅值峰值点即为谐振频率点，该特征点对应谐振频率的一 步估计值为4.4118X104Hz。
[0048] 选取谐振频率附近区段[4-I4f，4+Mf]，其中彳取为4,即将[4. 3333 X 104， 4.4902 X 104]作为细化区域，细分段数为100，细化后确定谐振频率的二步估计值为4.4180 X 104Hz，谐振幅值的二步估计值为3.6496 X 101;
[0049] 选取零频附近区段作为细化区域，细分段数为100,细化后确定零频幅值二步估计 值为 4 ? 1687X10-、
[0050] 利用公式(1-6)确定对应二阶传递函数模型参数，计算结果见表1。
[0051] 表1二阶传递函数模型参数计算结果
[0053]由对应二阶传递函数模型参数依据公式(1-7)构建加速度计传递函数模型，为：
[0055]以实际加速度作为输入，由模型（1-8)得到的输出与实际测得输出进行比较，如图 4、图5所示，图5为图4的局部放大。易知，所得传递函数模型能够很好的描述加速度计输入 输出特性，由模型得到的输出信号与实际测得输出信号差别不大，部分拟合效果不好的原 因可能是由量测噪声及环境噪声扰动造成的。
【主权项】
1. 一种基于细化谱分析的加速度计动态模型结构参数辨识的方法，其特征在于:该方 法包括以下步骤， 步骤一:通过加速度计绝对法冲击校准实验获取加速度计输入输出数据及数据截取； 步骤二:对步骤一所获得的输入输出数据进行离散傅里叶变换，确定加速度计频率响 应估计;对频率响应估计进行分析，得到谐振频率一步估计值； 步骤三:利用频谱细化方法对谐振频率及零频附近区段细化，得到谐振幅值谐振频率 及零频幅值二步估计值； 步骤四：利用模型参数与特征点坐标二步估计值关系，确定加速度计动态模型结构参 数估计，实现基于细化谱分析的加速度计动态模型结构参数辨识。2. 根据权利要求1所述的一种基于细化谱分析的加速度计动态模型结构参数辨识的方 法，其特征在于:所述步骤一，具体包括， 基于加速度计绝对法冲击校准实验，利用数据采集软硬件采集冲击激励加速度及加速 度计响应的时间序列；加速度计绝对法冲击校准实验的装置是依据IS016063-13绝对法冲 击校准标准建立的冲击激励系统，要求冲击激励加速度信号频率响应范围能够覆盖加速度 计谐振频率;数据采集系统的采样频率为f s; 将得到的时间序列进行截取，去掉数据中激励作用之前的部分与响应结束之后的部 分，以保留信噪比较大的数据部分，即保留冲击激励有效作用的数据部分，用于后续处理计 算。3. 根据权利要求1所述的一种基于细化谱分析的加速度计动态模型结构参数辨识的方 法，其特征在于:所述步骤二，具体包括： 利用步骤一所获得的输入加速度时间序列和加速度计输出时间序列，分别记为u(k)和 y (k)，k = 0，1，2，. . .，N-1，N为时间序列长度，k为离散数据序号;利用公式（1-1)确定输入加 速度离散频谱U( 1)和加速度计输出离散频谱Y( 1 )，1为离散谱线序号，；利 用公式(1-2)确定加速度计频率响应估计其中，疗为U的共辄； 利用公式（1 - 3 )确定△ f，确定横坐标序列；确定幅频特性纵坐标序列 10 X丨加;)))，确定相频特性纵坐标序列o?g/e( ?ο;)) X180/;r,其中，j·%))为 频率响应估计的模，/q))为频率响应估计的相角，所述相角单位为弧度;绘出加速 度计频率响应估计曲线，确定幅值峰值点横坐标值，记为成.，则即为加速度计谐振频率一 步估计值；4. 根据权利要求1所述的一种基于细化谱分析的加速度计动态模型结构参数辨识的方 法，其特征在于:所述步骤三，具体包括： 以4为中心频率，确定谐振频率附近频率响应特征区域|；氣-乂^,4+从/；^为调整因 子，取值为2~5;利用公式(1-3)确定Af; 确定细化段数为L，确定细化离散频率序列利用公式（1-4)分别确定细化后加速度计输入输出序列的离散傅里叶变换Xx(f)、Yx (f);利用公式(1-5)确定谐振频率附近区段细化后的加速度计频率响应估计0乂/%);其中，t为Ux的共辄； 利用l〇x )))确定谐振频率附近区段频率响应幅值序列，搜索其中最大值以 及该最大值所对应的频率值，即为谐振幅值及谐振频率二步估计值，分别记为与成； 然后，采样相同方法对零频附近区段进行细化，确定零频附近区段频率响应幅值序列， 搜索零频所对应的频率响应幅值，即为零频幅值二步估计值，记为5. 根据权利要求1所述的一种基于细化谱分析的加速度计动态模型结构参数辨识的方 法，其特征在于:所述步骤四，具体包括： 根据步骤三获得的谐振幅值、谐振频率和零频幅值的二步估计值，利用模型参数与特 征点坐标值二步估计值关系公式(1-6)确定对应二阶传递函数模型参数；将参数Κ、ωη、δ代入公式（1-7)，得到等效为单自由度质量弹簧阻尼系统的加速度计二 阶传递函数模型；其中，Κ、ω4Ρδ分别表示加速度计二阶模型静态传递系数、固有频率和阻尼比。
【文档编号】G01P15/09GK105929201SQ201610232361
【公开日】2016年9月7日
【申请日】2016年4月14日
【发明人】王建林, 魏青轩, 赵利强, 于涛, 王颖
【申请人】北京化工大学