CN113009828B - 一种复杂力学系统动力学参数的抗干扰标定方法 - Google Patents
一种复杂力学系统动力学参数的抗干扰标定方法 Download PDFInfo
- Publication number
- CN113009828B CN113009828B CN202110201532.1A CN202110201532A CN113009828B CN 113009828 B CN113009828 B CN 113009828B CN 202110201532 A CN202110201532 A CN 202110201532A CN 113009828 B CN113009828 B CN 113009828B
- Authority
- CN
- China
- Prior art keywords
- interference
- mechanical system
- complex mechanical
- model
- dynamic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000009795 derivation Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000012795 verification Methods 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 239000013598 vector Substances 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
Abstract
本发明公开一种复杂力学系统动力学参数的抗干扰标定方法,针对慢变干扰情况进行分析,在理想化的复杂力学系统动力学模型中加入集总力矩干扰项,从而使得所构建的模型能够更为准确地反映系统的实际状况;通过对预估干扰的复杂力学系统动力学模型中的已知的状态变量进行对应的时间域上的求导,将预估干扰的复杂力学系统动力学模型中的干扰化解,削弱干扰对系统参数标定的影响。本发明既考虑到外界的干扰,又保证了复杂力学系统模型的准确性,其能够在不增加外部设备和经济成本的条件下,有效且准确的标定复杂力学系统中的动力学参数。
Description
技术领域
本发明涉及自动控制技术领域,具体涉及一种复杂力学系统动力学参数的抗干扰标定方法。
背景技术
机械系统力学模型在自动控制领域中占据着关键的地位。在机械系统的描述和控制器的设计中,采用动力学模型对系统进行动态表征,设计相关控制器实现自动控制,这是一种常用的技术方法,其中力学系统模型中动力学参数是用于表征系统动态的重要组成部分。动力学参数的精确性对系统的描述及其自动控制具有不可或缺的影响。
力学系统的复杂程度决定着获取动力学参数的难易程度,而动力学参数的精确度决定着力学系统的准确性。一般情况下,力学系统越复杂其动力学参数的标定就越困难,常规标定系统动力学参数可采用实际测量物理参数、三维模型测量参数或结合离线分体测量参数的方式按一定的规则进行数值计算标定;也可以通过逆向参数辨识的方法标定。然而,上述提及的方式均在理想条件或不考虑干扰影响的情况下标定动力学参数,其标识结果与真实值误差较大,因此提高复杂力学系统动力学参数抗干扰标定的精确度,具有重要的理论意义和工程价值。
发明内容
本发明所要解决的是现有力学系统的动力学参数的标定方法的标识结果与真实值误差较大的问题,提供一种复杂力学系统动力学参数的抗干扰标定方法,能有效的计算复杂力学系统的动力学参数,同时也能削弱因外界客观产生的未知干扰对复杂力学系统动力学参数标定的影响。
为解决上述问题,本发明是通过以下技术方案实现的:
一种复杂力学系统动力学参数的抗干扰标定方法,包括步骤如下:
步骤1、针对复杂力学系统建立对应的复杂力学系统动力学模型:
步骤2、分析复杂力学系统在实际工况下可能受到的干扰情况,并以集总力矩干扰项的形式补偿于复杂力学系统动力学模型中,得到预估干扰的复杂力学系统动力学模型:
步骤3、先根据系统设计需求对预估干扰的复杂力学系统动力学模型进行变量分离,分离出其中的已知的状态变量和未知的动力学参数变量;再将已知的状态变量作为常数矩阵,将未知的动力学参数变量作为待求解矩阵,得到重组后的复杂力学系统动力学模型;
步骤4、对重组后的复杂力学系统动力学模型中的已知的状态变量进行对应的时间域上的求导,其中求导阶次与预估的干扰模型的集总力矩干扰项Δf(t)的阶次相同,以将重组后的复杂力学系统动力学模型转换为无干扰的复杂力学系统动力学模型;
步骤5、利用不同环境下的样本数据对无干扰的复杂力学系统动力学模型进行验证:如果满足验证要求,则无干扰的复杂力学系统动力学模型即为抗干扰标定的复杂力学系统动力学模型;否则,返回步骤2,对集总力矩干扰项进行调整;
上述步骤2中,需要确定预估干扰的复杂力学系统动力学模型中的集总力矩干扰项的阶次。
上述步骤5的验证过程如下:
步骤5.1、将样本数据中的已知的状态变量带入到无干扰的复杂力学系统动力学模型中,求解得到动力学参数变量;
步骤5.2、将求解得到的动力学参数变量与对应样本数据中的动力学参数变量进行比较:若求解得到的动力学参数变量与样本数据中的动力学参数变量相差在允许范围内时,认为该求解得到的动力学参数变量准确;否则,认为该求解得到的动力学参数变量不准确;
步骤5.3、统计所有求解得到的动力学参数变量的准确率,若准确率达到设定的阈值时,则认为无干扰的复杂力学系统动力学模型满足验证要求;否则,认为无干扰的复杂力学系统动力学模型不满足验证要求。
与现有技术相比,本发明具有如下特点:
1)考虑到复杂力学系统的动力学参数在标定时会受到外界干扰因素的影响,本发明针对慢变干扰情况进行分析,在理想化的复杂力学系统动力学模型中加入集总力矩干扰项,从而使得所构建的模型能够更为准确地反映系统的实际状况;
2)考虑到复杂力学系统的干扰无法直接获取,本发明通过对预估干扰的复杂力学系统动力学模型中的已知的状态变量进行对应的时间域上的求导,将预估干扰的复杂力学系统动力学模型中的干扰化解,削弱干扰对系统参数标定的影响;
3)本发明既考虑到外界的干扰,又保证了复杂力学系统模型的准确性,其能够在不增加外部设备和经济成本的条件下,有效且准确的标定复杂力学系统中的动力学参数。
附图说明
图1为一种复杂力学系统动力学参数的抗干扰标定方法的流程图。
图2为一阶导数下的车架横滚角状态图。
图3为一阶导数下的车把转角状态图。
图4为一阶导数下的车体线速度状态图。
图5为一阶导数下的力矩对比曲线图。
图6为一阶导数下的力矩误差曲线图。
图7为一阶导数下的不同工作环境下的车架横滚角状态图。
图8为一阶导数下的不同工作环境下的车把转角状态图。
图9为一阶导数下的不同工作环境下的车体线速度状态图。
图10为一阶导数下的力矩对比曲线图。
图11为一阶导数下的力矩误差曲线图。
图12为二阶导数下的车架横滚角状态图。
图13为二阶导数下的车把转角状态图。
图14为二阶导数下的车体线速度状态图。
图15为二阶导数下的力矩对比曲线图。
图16为二阶导数下的力矩误差曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实例,对本发明进一步详细说明。
本发明以复杂力学系统动力学参数的抗干扰标定技术为基础,其基本原理是预估干扰项阶次,结合系统动力学模型,通过微分求导,提高系统动态响应方程阶次,削弱干扰对系统参数标定的影响,实现对复杂力学系统动力学参数的标定。
针对标定某机械系统数学模型的动力学参数,本发明所设计的一种复杂力学系统动力学参数的抗干扰标定方法,如图1所示,其具体实施步骤如下:
1)建立复杂力学系统动力学模型:
针对复杂力学系统建立该系统对应的复杂力学系统动力学。复杂力学系统动力学模型的通式为:
2)预估计系统干扰项,得到预估干扰的复杂力学系统动力学模型:
复杂力学系统在不同的工作条件下,会受到不同程度的未知干扰,干扰项不易直接由传感器等装置直接测量获取,而且其组成成分复杂多变,可采用数学函数组合的形式进行表征。分析系统在实际工况下可能受到的干扰情况,用函数的形式对其做预估计,即可以通过简单的线性函数或非线性函数叠加对干扰项展开建模。例如采用不同幅值、不同相位的三角函数、衰减函数等进行叠加表征物理系统在实际运动过程中可能受到的干扰情况,其中记干扰项为Δf(t)。
由于干扰是不易测量的,但是在系统实际工作过程中,干扰真实存在,并且不可不免,将干扰项考虑进入模型中计算中是具有一定工程意义的,所以需要对系统的干扰项进行预估分析,将干扰考虑到动力学模型中,得到预估干扰的复杂力学系统动力学模型:
其中,Δf(t)是以集总力矩干扰项的形式补偿于动力学模型中,同时需要对预估的干扰模型进行分析,确定干扰阶次。
3)对预估干扰的复杂力学系统动力学模型进行变量分离与模型重组,得到重组后的复杂力学系统动力学模型:
先根据系统设计需求对预估干扰的复杂力学系统动力学模型进行变量分离,即将预估干扰的复杂力学系统动力学模型中的变量分为已知状态变量和未知的动力学参数变量,其中已知的状态变量可以通过复杂力学系统即物理实体运动状况获取。再将已知的状态变量作为常数矩阵,将未知的动力学参数变量作为待求解矩阵,对预估干扰的复杂力学系统动力学模型进行重组,得到重组后的复杂力学系统动力学模型。
4)对重组后的复杂力学系统动力学模型进行时间域上的求导,得到无干扰的复杂力学系统动力学模型:
对重组后的复杂力学系统动力学模型中的已知的状态变量进行对应的时间域上的一次或多次求导,通过模型对时间求导的方式,提升系统模型状态变量的阶次,消除预估的干扰项的影响,从而将重组后的复杂力学系统动力学模型转换为无干扰的复杂力学系统动力学模型。
求导阶次的依据是以模型中的干扰项Δf(t)的阶次所决定,即求导阶次与预估的干扰模型的集总力矩干扰项Δf(t)的阶次相同。根据干扰项的阶次,提高已知的状态变量的阶次,削弱干扰对系统参数标定的影响。
5)利用不同环境下的样本数据对无干扰的复杂力学系统动力学模型进行验证:
步骤5.1、将样本数据中的已知的状态变量带入到无干扰的复杂力学系统动力学模型中,求解得到动力学参数变量;
步骤5.2、将求解得到的动力学参数变量与对应样本数据中的动力学参数变量进行比较:若求解得到的动力学参数变量与样本数据中的动力学参数变量相差在允许范围内时,认为该求解得到的动力学参数变量准确;否则,认为该求解得到的动力学参数变量不准确;
步骤5.3、统计所有求解得到的动力学参数变量的准确率,若准确率达到设定的阈值时,则认为无干扰的复杂力学系统动力学模型满足验证要求,无干扰的复杂力学系统动力学模型即为抗干扰标定的复杂力学系统动力学模型;否则,认为无干扰的复杂力学系统动力学模型不满足验证要求,返回步骤2,对集总力矩干扰项进行调整,并执行后续步骤。
求解动力学参数变量需要获取大量的状态变量数据,一部分状态数据可通过采集复杂力学系统的关节或构件上的传感器数据获取,另一部分高阶状态数据不能通过传感器直接测量获得。高阶状态数据可采用常用的差分方,或者其他通用的数据的滤波的方法获取。
通过不同工作环境下获得的数据代入到确定动力学参数的动力学模型中,判断获取的动力学参数的准确性,从而实现对复杂力学系统动力学参数的抗干扰标定方法。
下面以复杂力学系统动力学模型中无人自行车的力学模型为例,对本发明进行进一步说明:
1)自行车模型的线性运动方程如下所示:
其中,S为与系统质量相关的对称矩阵,v(t)U1为系统阻尼系数相关的矩阵,gJ0+v(t)2J2为系统刚度和重力相关的矩阵,q(t)为车架横滚角和车把转角组成的向量,f(t)为力矩向量。
无人自行车由车架、车把、前轮、后轮、惯性测量单元、工控机、驱动器模块以及DSP模块等部分组成,其中前轮为驱动轮,IMU为系统惯性测量单元。
2)对预估的干扰模型进行分析,确定干扰阶次,在无人自行车的力学模型中加入预测干扰项。以无人自行车相应传感器测量的系统状态及通过滤波后得到的自行车车体速度为输入,以待识别的动力学参数为输出,利用变量分离技术提取模型中动力学参数变量的系数;动力学参数组成一个矢量,构造系统的识别模型。
3)采用变量分离技术提取模型中动力学参数变量的系数,将自平衡状态下的无人自行车通过传感器测量的系统状态,以及经过滤波处理后的自行车车体速度作为输入,将动力学模型中待识别的动力学参数为输出;构造系统参数的识别模型,其模型为:
AX=y
其中,A为一个2×15的矩阵,X为一个15×1的列向量,y为一个2×1的列向量。
采集传感器数据,得到动力学参数识别结果的初值。由于构造的系统参数的模型是一个由15个未知数、2个方程的欠定方程组。故采用按一定周期连续采集系统状态,扩展为一个新的矩阵,重新构造一个新的动力学参数识别模型。采取最小二乘法的方法,处理新的模型,得到最优解。
5)利用不同环境下的样本数据对无干扰的复杂力学系统动力学模型进行验证。
该模型在数学上一个含有15个未知数、2个方程的方程组,即方程数少于变量个数的欠定方程,可见通过上述步骤求得的X值是不唯一的,不能直接对变量进行求解。本发明采取一种有效的计算方法,即按一定周期连续采集系统状态,扩展为一个新的矩阵,重新构造一个新的动力学参数模型,然后通过最小二乘法求其最优解。
如图2和3所示,其车架横滚角不超过±0.06rad,而车把转角不超过±0.26rad,将实验采集的状态数据代入卡尔曼滤波算法中,得到如图4所示的实时速度值。将滤波后的车体速度及其他传感器采集到的状态变量代入到一阶导数后的模型中,得到系统动力学参数如表1所示:
表1动力学参数识别结果
符号 | 数值 | 符号 | 数值 | 符号 | 数值 |
M<sub>11</sub> | 0.0453 | M<sub>12</sub> | 0.2319 | M<sub>22</sub> | 1.5435 |
C<sub>11</sub> | 0.0093 | C<sub>12</sub> | 1.1814 | C<sub>21</sub> | -0.1458 |
C<sub>22</sub> | 11.7435 | K<sub>011</sub> | -0.3129 | K<sub>012</sub> | -0.1537 |
K<sub>021</sub> | 0.2115 | K<sub>022</sub> | 0.0176 | K<sub>211</sub> | 0 |
K<sub>212</sub> | 1.1457 | K<sub>221</sub> | 0 | K<sub>222</sub> | -0.0111 |
用传感器采集到相关数据,并将表1中的动力学参数代入到一阶导数后的模型中,计算得系统的驱动力矩和力矩误差,如图5和6所示。
6)对5)中识别动力学参数进行响应分析并计算误差。
为了验证识别到的动力学参数的可靠性,进行一组新的实验,即在不同工况下驱动无人自行车开展另外一组平衡行走实验,其车架横滚角、车把转角和车体线速度的状态图,如图7-9所示,对比该计算力矩与传感器的实测力矩的差异,结果如图10和11所示。实验结果表明,通过本发明提出的复杂力学系统动力学参数的抗干扰标定方法可以更好的反映自行车的动态,证明此方法是有效的。
对该模型进行二阶求导,其车架横滚角、车把转角和车体线速度的状态图,如图12-14所示,对比该计算力矩与传感器的实测力矩的差异,如图15和16所示。经实验结果分析,对于不同阶次干扰项的估计均可重构动力学模型,求解得到模型中的动力学参数,与此同时,预估不同阶次的干扰构建的模型求解出的动力学参数在精度上是有一定影响的,可根据实际情况对干扰阶次进行选取。
本发明以复杂力学系统动力学模型为基础,分析干扰的情况;预估干扰项的阶次,构造包含干扰项的力学模型;利用变量分离的方法使模型重组,将已知的状态变量与未知的动力学参数变量分离;根据干扰项的阶次,提高变量分离后已知的状态变量的阶次,即进行对应的时间域上的求导,削弱干扰对系统参数标定的影响,标定动力学参数。考虑到部分高阶状态数据不能通过传感器直接获得,使用差分方式会引入大量噪声,所以对数据进行高阶滤波处理,将处理后的数据代入到求导后的新模型中,求解得到动力学参数。上述方案的优点在于既将外界客观存在的干扰考虑在内又保证了复杂力学系统模型的准确性。本发明所述的抗干扰标定方法能有效的标定复杂力学系统中的动力学参数。
需要说明的是,尽管以上本发明所述的实施例是说明性的,但这并非是对本发明的限制,因此本发明并不局限于上述具体实施方式中。在不脱离本发明原理的情况下,凡是本领域技术人员在本发明的启示下获得的其它实施方式,均视为在本发明的保护之内。
Claims (3)
1.一种复杂力学系统动力学参数的抗干扰标定方法,该复杂力学系统为无人自行车,其特征是,包括步骤如下:
步骤1、针对复杂力学系统建立对应的复杂力学系统动力学模型:
步骤2、分析复杂力学系统在实际工况下可能受到的干扰情况,并以集总力矩干扰项的形式补偿于复杂力学系统动力学模型中,得到预估干扰的复杂力学系统动力学模型:
步骤3、先根据系统设计需求对预估干扰的复杂力学系统动力学模型进行变量分离,分离出其中的已知的状态变量和未知的动力学参数变量;再将已知的状态变量作为常数矩阵,将未知的动力学参数变量作为待求解矩阵,得到重组后的复杂力学系统动力学模型;
步骤4、对重组后的复杂力学系统动力学模型中的已知的状态变量进行对应的时间域上的求导,其中求导阶次与预估的干扰模型的集总力矩干扰项Δf(t)的阶次相同,以将重组后的复杂力学系统动力学模型转换为无干扰的复杂力学系统动力学模型;
步骤5、利用不同环境下的样本数据对无干扰的复杂力学系统动力学模型进行验证:如果满足验证要求,则无干扰的复杂力学系统动力学模型即为抗干扰标定的复杂力学系统动力学模型;否则,返回步骤2,对集总力矩干扰项进行调整;
2.根据权利要求1所述的一种复杂力学系统动力学参数的抗干扰标定方法,其特征是,步骤2中,需要确定预估干扰的复杂力学系统动力学模型中的集总力矩干扰项的阶次。
3.根据权利要求1所述的一种复杂力学系统动力学参数的抗干扰标定方法,其特征是,步骤5的验证过程如下:
步骤5.1、将样本数据中的已知的状态变量带入到无干扰的复杂力学系统动力学模型中,求解得到动力学参数变量;
步骤5.2、将求解得到的动力学参数变量与对应样本数据中的动力学参数变量进行比较:若求解得到的动力学参数变量与样本数据中的动力学参数变量相差在允许范围内时,认为该求解得到的动力学参数变量准确;否则,认为该求解得到的动力学参数变量不准确;
步骤5.3、统计所有求解得到的动力学参数变量的准确率,若准确率达到设定的阈值时,则认为无干扰的复杂力学系统动力学模型满足验证要求;否则,认为无干扰的复杂力学系统动力学模型不满足验证要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110201532.1A CN113009828B (zh) | 2021-02-23 | 2021-02-23 | 一种复杂力学系统动力学参数的抗干扰标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110201532.1A CN113009828B (zh) | 2021-02-23 | 2021-02-23 | 一种复杂力学系统动力学参数的抗干扰标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113009828A CN113009828A (zh) | 2021-06-22 |
CN113009828B true CN113009828B (zh) | 2022-07-22 |
Family
ID=76407633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110201532.1A Active CN113009828B (zh) | 2021-02-23 | 2021-02-23 | 一种复杂力学系统动力学参数的抗干扰标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113009828B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BG50580A1 (en) * | 1988-07-18 | 1992-09-15 | Univ Tekhnicheski | Method for technilogical objects control |
EP0856784A2 (de) * | 1997-01-31 | 1998-08-05 | Daimler-Benz Aktiengesellschaft | Verfahren und Vorrichtung zur bordautonomen Bestimmung der Position eines Satelliten |
CN108673502A (zh) * | 2018-05-18 | 2018-10-19 | 昆明理工大学 | 一种基于未知模型动态观测器的机械臂运动控制方法 |
CN111198570A (zh) * | 2020-02-05 | 2020-05-26 | 哈尔滨工业大学 | 一种基于固定时间微分器预测的抗时延高精度自抗扰姿态控制方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1249478A1 (ru) * | 1985-01-03 | 1986-08-07 | Ленинградский Ордена Ленина Электротехнический Институт Им.В.И.Ульянова/Ленина/ | Самонастраивающа с система управлени |
JP5011757B2 (ja) * | 2005-08-02 | 2012-08-29 | 日産自動車株式会社 | 車両用操舵装置 |
CN100468427C (zh) * | 2007-05-11 | 2009-03-11 | 山东大学 | 一种利用计算机模拟蛋白质相互作用的方法 |
US8380473B2 (en) * | 2009-06-13 | 2013-02-19 | Eric T. Falangas | Method of modeling dynamic characteristics of a flight vehicle |
CN102298390A (zh) * | 2011-06-24 | 2011-12-28 | 北京航空航天大学 | 一种抗干扰挠性航天器姿态和振动复合控制方法 |
JP6310722B2 (ja) * | 2013-11-08 | 2018-04-11 | 国立大学法人大阪大学 | システムのパラメータ同定方法 |
CN107357166B (zh) * | 2017-04-29 | 2020-11-10 | 天津大学 | 小型无人直升机的无模型自适应鲁棒控制方法 |
CN108180910B (zh) * | 2017-12-26 | 2019-01-08 | 北京航空航天大学 | 一种基于气动参数不确定的飞行器快速高精度制导方法 |
CN108710302B (zh) * | 2018-06-20 | 2021-07-06 | 天津大学 | 无源性全方位移动机器人轨迹跟踪自抗扰控制方法 |
CN111444623B (zh) * | 2020-03-31 | 2024-01-26 | 桂林电子科技大学 | 一种阻尼非线性商用车悬架动力学的协同优化方法及系统 |
CN111546336B (zh) * | 2020-04-29 | 2021-01-29 | 北京航空航天大学 | 一种用于机器人系统的灰箱模型参数辨识方法及系统 |
-
2021
- 2021-02-23 CN CN202110201532.1A patent/CN113009828B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BG50580A1 (en) * | 1988-07-18 | 1992-09-15 | Univ Tekhnicheski | Method for technilogical objects control |
EP0856784A2 (de) * | 1997-01-31 | 1998-08-05 | Daimler-Benz Aktiengesellschaft | Verfahren und Vorrichtung zur bordautonomen Bestimmung der Position eines Satelliten |
CN108673502A (zh) * | 2018-05-18 | 2018-10-19 | 昆明理工大学 | 一种基于未知模型动态观测器的机械臂运动控制方法 |
CN111198570A (zh) * | 2020-02-05 | 2020-05-26 | 哈尔滨工业大学 | 一种基于固定时间微分器预测的抗时延高精度自抗扰姿态控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113009828A (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109940622B (zh) | 一种基于电机电流的机器人机械臂无传感碰撞检测方法 | |
CN107000755B (zh) | 用于估计影响车辆动力学变量的方法和对应的虚拟传感器 | |
JP4724163B2 (ja) | 車体速度推定装置 | |
Venture et al. | Modeling and identification of passenger car dynamics using robotics formalism | |
CN106125548A (zh) | 工业机器人动力学模型参数辨识方法 | |
CN112577516B (zh) | 一种用于车辆轮速误差识别和补偿的方法和系统 | |
JP2009073466A (ja) | 車両姿勢角推定装置及びプログラム | |
CN111949036B (zh) | 一种轨迹跟踪控制方法、系统及两轮差速移动机器人 | |
CN103279675A (zh) | 轮胎-路面附着系数与轮胎侧偏角的估计方法 | |
CN108068882A (zh) | 控制转向系统的方法 | |
CN112270039A (zh) | 基于分布式异步融合的线控底盘车辆非线性状态估计方法 | |
Turnip et al. | Estimation of the wheel-ground contacttire forces using extended kalman filter | |
Carratù et al. | Velocity prediction from acceleration measurements in motorcycle suspensions | |
Kelly | Fast and easy systematic and stochastic odometry calibration | |
JP4568996B2 (ja) | 操舵装置と操舵装置の抗力算出装置 | |
JP2008094375A (ja) | 車両物理量推定装置及びプログラム | |
Tong | An approach for vehicle state estimation using extended Kalman filter | |
CN113009828B (zh) | 一种复杂力学系统动力学参数的抗干扰标定方法 | |
Thaller et al. | Fast determination of vehicle suspension parameters via continuous time system identification | |
Jiang et al. | Modeling of unbounded long-range drift in visual odometry | |
Fouka et al. | Motorcycle state estimation and tire cornering stiffness identification applied to road safety: Using observer-based identifiers | |
CN117091592A (zh) | 姿态解算方法、姿态解算装置以及计算机存储介质 | |
Khaknejad et al. | Identification of vehicle parameters using modified least square method in ADAMS/Car | |
CN114084199B (zh) | 一种基于递归图分析的列车稳定性评估方法及系统 | |
Samuel et al. | A reduced-order multisensor-based force observer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210622 Assignee: Guangxi Cuizhi Rongchuang Technology Co.,Ltd. Assignor: GUILIN University OF ELECTRONIC TECHNOLOGY Contract record no.: X2023980046599 Denomination of invention: An anti-interference calibration method for dynamic parameters of complex mechanical systems Granted publication date: 20220722 License type: Common License Record date: 20231108 |