一种基于增强响应灵敏度的多体系统参数识别方法技术方案

本发明专利技术公开了一种基于增强响应灵敏度的多体系统参数识别方法，该方法包括：构建多体系统的动力学方程和待识别参数最小二乘优化方程；载入测量数据求解多体系统的动力学方程；求取多体系统的响应灵敏度；建立迭代参数的最小二乘方程；引入Tikhonov正则化计算迭代参数；对待识别参数进行更新；再次求解多体系统的动力学方程；引入置信域；判断到更新后的待识别参数符合预设条件，输出更新后的待识别参数。通过使用本发明专利技术，可以准确识别多体系统的参数，并且具有良好的抗噪性。本发明专利技术可广泛应用于参数识别领域。应用于参数识别领域。应用于参数识别领域。

【技术实现步骤摘要】
一种基于增强响应灵敏度的多体系统参数识别方法


[0001]本专利技术涉及参数识别领域，尤其涉及一种基于增强响应灵敏度的多体系统参数识别方法。

技术介绍

[0002]典型的多体系统包括运载器、航天器、机器人、各种机械装置等。多体系统广泛存在于工业生产与日常生活中，对于多体系统的模型建立、优化设计、测试诊断过程等，识别其系统参数都显得尤为重要。多体系统的参数识别主要是根据系统的响应识别出系统的几何参数和惯性参数。传统的参数识别方法主要缺陷在于形成的是有偏估计，识别准确度不高。

技术实现思路

[0003]为了解决上述技术问题，本专利技术的目的是提供一种基于增强响应灵敏度的多体系统参数识别方法，准确识别多体系统的参数，并且具有良好的抗噪性。
[0004]本专利技术所采用的技术方案是：一种基于增强响应灵敏度的多体系统参数识别方法，包括以下步骤：
[0005]S1、基于绝对坐标方法构建多体系统的动力学方程，并建立待识别参数最小二乘优化方程；
[0006]S2、载入测量数据并基于增广法求解多体系统的动力学方程，计算得到系统响应和响应残差；
[0007]S3、基于中心差分求取多体系统的响应灵敏度；
[0008]S4、基于待识别参数最小二乘优化方程建立迭代参数的最小二乘方程；
[0009]S5、引入Tikhonov正则化对迭代参数的最小二乘方程进行处理，结合响应残差和相应灵敏度计算得到迭代参数；
[0010]S6、根据迭代参数对待识别参数进行更新，得到更新后的待识别参数；
[0011]S7、基于更新后的待识别参数重新求解多体系统的动力学方程，计算得到新系统响应和新响应残差。
[0012]S8、引入置信域计算一致性指标并与预设置信因子比较，调整迭代参数；
[0013]S9、判断到迭代参数与更新后的待识别参数的比值超出误差容限内，输出更新后的待识别参数，否则返回步骤S2。
[0014]优选地，所述多体系统的动力学方程公式如下：
[0015][0016]上式中，b∈R
k
表示一个系统参数组成的列向量，分别是n
×
1广义坐标系下的位移、速度和加速度，Φ∈R
m
表示完整约束方程，是Jacobian矩阵，与Lagrange乘子λ∈R
m
一同表示由约束方程引起的系统的反作用力，M∈R
n
×
n
表示系统的广义
质量矩阵，而Q∈R
n
表示系统受到的广义外力矢量。
[0017]对于位置φ∈R
n
‑
m
和速度φ
′
∈R
n
‑
m
的系统的初始条件可以表示为如下形式：
[0018][0019]利用传统增广法求解方程。由于式是指标3的微分代数方程，先对中的约束方程关于时间t求二阶导，将其化为指标1的微分代数方程，并写成如下的等价形式：
[0020][0021]其中，
[0022]由式中后两式直接增广计算，得到
[0023][0024][0025]由的第一式和即构成常规的一阶常微分方程，
[0026][0027]上式可以利用Matlab的ode进行求解。
[0028]然而，由于上述方法只考虑了加速度形式的约束方程，所得到的响应未能严格符合位移或者速度的约束方程，会引起约束违约的问题。为了进行违约修正，通过引入人工阻尼将不稳定的约束方程转化为
[0029][0030]从而，的第三式变为
[0031][0032]再利用Matlab求解即可。
[0033]对于阻尼参数的选取，提出了参数α和β的一种自动选取方法，能得到很好的数值稳定性，即，
[0034][0035]其中，h为求解微分方程的离散差分格式的步长。
[0036]优选地，在典型参数识别流程中，对于一组实际测量得到的数据待识别的系统参数为b∈R
k
应该是如下非线性最小二乘优化问题的解，所述待识别参数最小二乘优化方程的公式如下：
[0037][0038]上式中，表示待识别参数b的可行域，表示测量数据，R
i
(b)表示通过求解动力
学方程得到的理论解，为长度是pl的列向量的形式，符号表示关于权重矩阵W的加权范数。
[0039]方程中的R(b)通常是关于b的非线性隐函数，于是需要利用迭代法求解如上的非线性最小二乘问题。接下来的关键就在于如何从已知的参数确定一个合适的迭代更新Δb。通常将在点处线性化，保留一阶项，即
[0040][0041]优选地，所述响应灵敏度的形式如下：
[0042][0043]上式中，表示响应灵敏度矩阵，表示对应的响应灵敏度。
[0044]具体地，基于梯度法解决反问题，灵敏度分析是不可或缺的。由于直接利用传统增广法求解多体系统的灵敏度复杂耗时，故采用有限差分法求取。本文中心差分求取响应灵敏度，即，
[0045][0046]其中，ε是一个充分小的正数；而q(t,b
i
+ε)和q(t,b
i
‑
ε)可以通过求解正问题得到。本文ε取0.001。
[0047]优选地，至此，即可得到关于Δb的一个线性最小二乘问题，所述迭代参数的最小二乘方程公式如下：
[0048][0049]上式中，Δb表示迭代参数，表示响应残差。
[0050]然而，问题通常是不适定的，因此需要引入Tikhonov正则化进行处理，所述引入Tikhonov正则化对迭代参数的最小二乘方程进行处理，公式表示如下：
[0051][0052]上式中，μ≥0为正则化参数，||
·
||2表示L2范数，从而可得，
[0053][0054]上式中，I为单位矩阵。正则化参数μ的值通常通过L曲线法选取，将其记为
[0055]优选地，通过线性化得到近似的线性最小二乘问题。要使线性化后得到的能较好地近似原非线性函数，需要在参数的邻域内，亦即置信域[]内。那么，对于式加入的置信域约束可以表示为，存在足够小的δ＞0，使得||Δb||＜δ。
由此可见，解得的Δb应该是δ的一个函数。为了选取合适的δ值，需要引入如下的一致性指标，所述一致性指标的公式如下：
[0056][0057]良好的一致性指标应该使得比某一特定的正数大，亦即，
[0058][0059]从而需要选取合适的δ使得一致性条件满足。
[0060]满足置信域的条件下，响应R(b)通常是关于b的光滑函数，且至少满足因此在满足的条件下，应有
[0061][0062]以及
[0063][0064]从上可知，通过适当的增加正则化参数μ，可以使迭代更新Δb
μ
变小，并且满足置信域条件。
[0065]最后，更新待识别的参数b
k+1
＝b
k
+Δ本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于增强响应灵敏度的多体系统参数识别方法，其特征在于，包括：S1、基于绝对坐标方法构建多体系统的动力学方程，并建立待识别参数最小二乘优化方程；S2、载入测量数据并基于增广法求解多体系统的动力学方程，计算得到系统响应和响应残差；S3、基于中心差分求取多体系统的响应灵敏度；S4、基于待识别参数最小二乘优化方程建立迭代参数的最小二乘方程；S5、引入Tikhonov正则化对迭代参数的最小二乘方程进行处理，结合响应残差和相应灵敏度计算得到迭代参数；S6、根据迭代参数对待识别参数进行更新，得到更新后的待识别参数；S7、基于更新后的待识别参数重新求解多体系统的动力学方程，计算得到新系统响应和新响应残差；S8、引入置信域计算一致性指标并与预设置信因子比较，调整迭代参数；S9、判断到迭代参数与更新后的待识别参数的比值超出误差容限内，输出更新后的待识别参数，否则返回步骤S2。2.根据权利要求1所述一种基于增强响应灵敏度的多体系统参数识别方法，其特征在于，所述多体系统的动力学方程公式如下：上式中，b∈R
k
表示一个系统参数组成的列向量，分别是n
×
1广义坐标系下的位移、速度和加速度，Φ∈R
m
表示完整约束方程，是Jacobian矩阵，与Lagrange乘子λ∈R
m
一同表示由约束方程引起的系统的反作用力，M∈R
n
×
n
表示系统的广义质...

【专利技术属性】
技术研发人员：吕中荣，陈敏，刘广，汪利，
申请(专利权)人：中山大学，
类型：发明
国别省市：