CN105404609A - 一种新型多目标动力系统参数辨识方法 - Google Patents
一种新型多目标动力系统参数辨识方法 Download PDFInfo
- Publication number
- CN105404609A CN105404609A CN201510687739.9A CN201510687739A CN105404609A CN 105404609 A CN105404609 A CN 105404609A CN 201510687739 A CN201510687739 A CN 201510687739A CN 105404609 A CN105404609 A CN 105404609A
- Authority
- CN
- China
- Prior art keywords
- identified
- power system
- value
- evaluation index
- power systems
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 238000011156 evaluation Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 230000033001 locomotion Effects 0.000 claims abstract description 13
- 230000010354 integration Effects 0.000 abstract description 8
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000000638 solvent extraction Methods 0.000 abstract description 2
- 238000009826 distribution Methods 0.000 description 14
- 238000013016 damping Methods 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000009467 reduction Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 239000013598 vector Substances 0.000 description 5
- 239000006096 absorbing agent Substances 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 230000035939 shock Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种新型多目标动力系统参数辨识方法,包括以下步骤:(1)针对待辨识动力系统建立运动微分方程;(2)选取多个待辨识动力系统输出响应的特征参数作为评价指标,并通过试验获得评价指标的值及待辨识动力系统的输入数据;(3)应用改进的分块增维精细积分法和待辨识动力系统的输入数据,通过待辨识动力系统的运动微分方程计算评价指标的值,将计算获得的评价指标的值与试验获得的评价指标的值相减获得差值,以差值的平方和作为待辨识动力系统的目标函数;(4)应用NSGA-II算法对目标函数进行多次寻优计算,直至输出结果稳定,获得最终的非劣解集;(5)选择结果。该方法效率高,适用面广,稳定性高。
Description
技术领域
本发明涉及系统参数辨识领域,更具体的说,是涉及一种新型多目标动力系统参数辨识方法。
背景技术
参数辨识是根据实验数据和建立的模型来确定一组参数值,使得由模型计算得到的数值结果能最好的拟合测试数据,从而可以为生产过程进行预测,提供一定的理论指导,在工程实际中拥有非常广泛的应用。传统参数辨识一般采用单目标辨识,或者将多目标辨识进行线性加权转化成单目标辨识。单目标辨识方法简便,但准确性较差,由于目标函数少,辨识结果很难全面反映被辨识模型的全部特性,而线性加权的多目标参数辨识方法在使用过程中必须首先确定每个目标函数的权重值,而权重值的确定十分困难,因此在实际使用过程中通常采用经验值,而这又会极大地降低辨识结果的准确性。所以,开发一种不需要线性加权的新型多目标参数辨识方法是十分必要的。
发明内容
本发明的目的是为了克服现有技术中的不足,提供一种新型多目标动力系统参数辨识方法,该方法从功能上可以分为两部分,第一部分是计算目标函数的值,第二部分是对目标函数进行寻优计算以获得最优解。为了实现上述功能,本发明方法主要由两种算法构成,分别为改进的分块增维精细积分法和NSGA-II算法;其中,改进的分块增维精细积分法主要用来计算目标函数的值,NSGA-II算法主要用来对目标函数进行寻优计算。
本发明的目的是通过以下技术方案实现的:
一种新型多目标动力系统参数辨识方法,包括以下步骤:
(1)选定要进行参数识别的待辨识动力系统,根据待辨识动力系统的结构建立运动微分方程,并给出运动微分方程中已知参数的值和待辨识参数的范围及约束条件;
(2)选取多个待辨识动力系统输出响应的特征参数作为评价指标,并通过试验获得评价指标的值及待辨识动力系统的输入数据;
(3)应用改进的分块增维精细积分法和待辨识动力系统的输入数据,通过待辨识动力系统的运动微分方程计算评价指标的值,将计算获得的评价指标的值与试验获得的评价指标的值相减获得差值,以差值的平方和作为待辨识动力系统的目标函数;
(4)寻优计算,对NSGA-II算法进行相应的参数设定,应用NSGA-II算法对目标函数进行多次寻优计算,直至输出结果稳定,获得最终的非劣解集;
(5)选择结果,结合工程实际,在非劣解集中选取合适的点为最优点,以最优点对应的参数为最终辨识结果。
步骤(5)中所述最优点的选取原则是最优点使不同的目标函数值都达到最小。
与现有技术相比,本发明的技术方案所带来的有益效果是:
1.本发明方法中同时采用了改进的分块增维精细积分法和NSGA-II算法,其中改进的分块增维精细积分法是在精细积分法的基础之上发展而来的,它不但具有精细积分法求解精度高,稳定性和收敛性好的特点,而且克服了精细积分法不能求解非线性系统的缺点,提高了求解效率;NSGA-II算法具有计算量小、搜索效率高的特点,在保持种群多样性和防止优秀个体流失方面表现突出,而且不需要设定初值。
2.本发明方法的优点在于目标函数的个数不受限制,求解时不需要对多个目标函数进行线性加权,并且由于改进的分块增维精细积分法可以方便地求解线性/非线性偏微分方程任意精度的数值解,使得本发明的多目标动力系统参数辨识方法适用范围非常广泛,只要系统的模型可以转化为偏微分方程,均可以使用本方法进行系统辨识,尤其适合动力系统的参数辨识,无论是线性还是非线性动力系统,无论有无外部激励,无论激励是谐波激励还是随机激励,本方法都能很方便的应用。
附图说明
图1是履带急救车非线性减振系统结构示意图。
图2是履带急救车非线性减振系统四自由度振动模型。
图3是本发明方法的流程图。
图4是本发明方法的计算机操作流程图。
图5是应用本发明方法对履带急救车非线性减振系统进行参数辨识所获得的第一前端个体分布示意图。
图6是由辨识结果计算获得的担架台垂向振动加速度功率谱密度和由试验数据计算获得的担架台垂向振动加速度功率谱密度的对比示意图。
图7是由辨识结果计算获得的担架台垂向振动加速度概率分布和由试验数据计算获得的担架台垂向振动加速度概率分布的对比示意图。
附图标记:1-大板车厢2-担架台3-零刚度减振器4-橡胶阻尼减振器5-履带底盘M1-担架抬质量M2-大板车厢质量Ks-零刚度减振器的线性刚度Kz-零刚度减振器的三次非线性刚度C1-零刚度减振器的阻尼K2-橡胶阻尼减振器的刚度C2-橡胶阻尼减振器的阻尼J1、J2-担架台和大板车厢的质量在几何中心处绕z轴的转动惯量l1、l2-担架台和大板车厢的几何中心距离各自减振器沿y轴的距离q1、q2-振动输入fun1-目标函数1fun2-目标函数2
具体实施方式
为了更加清楚的叙述本发明的实施过程,下面结合工程实例及说明书附图对本发明作进一步的详细描述:
图3为本发明方法的流程图,应用本发明所提出的新型多目标动力系统参数辨识方法对某履带急救车非线性减振系统进行参数辨识,履带急救车非线性减振系统结构示意图如图1所示。忽略其他两个方向,仅保留对成员舒适性影响最大的垂直方向,建立减振系统四自由度振动模型如图2所示,其中Ks、Kz、C1、K2和C2作为待辨识参数,M1、M2、J1、J2、l1和l2为已知参数,取值如下:M1=180kg、M2=2000kg、J1=58kg·m2、J2=3665kg·m2、l1=0.8m、l2=1.539m。
本发明方法中所提出的新型多目标动力系统参数辨识方法由改进的分块增维精细积分法和NSGA-II算法构成,具体使用时需要通过计算机编程来实现,本实施例采用MATLAB编程,NSGA-II算法采用MATLAB内置的gamultiobj函数来实现,gamultiobj函数所采用的是NSGA-II算法的一种变形。而改进的分块增维精细积分法需要自行编程实现,下面结合履带急救车非线性减振系统的振动模型对改进的分块增维精细积分法的计算过程进行介绍。
对图2所示振动模型建立运动微分方程如下:
其中,
对运动微分方程中的Ks、Kz、C1、K2和C2进行参数辨识,首先需要将式(1)改写如下形式
式中 In为n阶单位矩阵。
引入新的变量X2n+1=1,则式(2)增加一维后变成齐次方程形式,即
记时间步长为τ。令Δt=τ/m,其中m=2L,L=20。经过推导可得如下计算过程
当循环(5)结束后,可得
Y(kτ)=Tk-1Y[(k-1)τ](k=1,2,…)(7)
对矩阵B和矩阵F进行如下分块:
由于矩阵F2n×1的前n行为0,分块以后得到所以矩阵G的计算公式变为
考虑增维后H的特点,Y的最后一行实际上不会对所求结果产生影响,则式(7)可变成如下形式
即
公式(2)-(11)即为改进的分块增维精细积分法的计算过程,具体应用时需要将公式(2)-(11)转换为计算机程序。
以上即为构成新型多目标动力系统参数辨识方法的两种基本算法的实现方式。
在参数辨识之前需要对履带急救车进行野外行驶振动试验以获得履带底盘5、大板车厢1、担架台3等部位的振动加速度数据信号,为本方法提供数据基础。
如图4所示,是利用本发明的实施方法编制出的计算机流程图,具体实施步骤包括如下:
(1)应用MATLAB的load函数导入由履带急救车野外行驶振动试验获得的非线性减振系统的振动数据,并存储在相应的向量内;
(2)对履带急救车非线性减振系统的已知参数、改进的分块增维精细积分法的积分步长、计算过程中所需的向量空间等进行设定;
(3)应用MATLAB的pwelch函数和ksdensity函数计算由试验获得的振动数据中担架台垂向振动加速度的功率谱密度pxx和概率分布hx,并存储在相应的向量之中;
(4)根据式(2)的形式建立履带急救车非线性减振系统振动模型的运动微分方程,在计算程序中只需要给出矩阵A和矩阵F即可;
(5)将公式(2)-(11)编写成MATLAB程序段,并应用该程序段对履带急救车非线性减振系统振动模型的运动微分方程进行求解,获得履带急救车非线性减振系统的振动数据,并存储在相应的向量内;
(6)应用MATLAB的pwelch函数和ksdensity函数计算由仿真获得的振动数据中担架台垂向振动加速度的功率谱密度pxx1和概率分布hx1,并存储在相应的向量之中;
(7)以pxx与pxx1、hx与hx1差值的平方和建立目标函数1和目标函数2,分别表示为fun1和fun2;
(8)应用gaoptimset函数对NSGA-II算法的最优前端个体系数、种群大小、最大进化代数、停止代数、适应度函数值偏差等进行设置;
(9)给出待辨识参数的范围以及约束条件;
(10)应用gamultiobj函数对目标函数进行优化计算,获得第一前端个体分布以及第一前端个体分布所对应的非劣解集;
(11)从非劣解集中选取最优点,并以最优点对应数据为最终辨识结果。
以上即为本发明方法的计算机实现过程,具体编程时,第(1)-(7)步编写在一个M文件里,第(8)-(10)步编写在另一个M文件里,第(11)步为选择参数辨识结果,选择的原则是尽可能使两个目标函数值都达到最小。
应用本发明所提出的新型多目标动力系统参数辨识方法对工程实例中的履带急救车非线性减振系统进行参数辨识,第一前端个体分布(ParetoFront)如图5所示,第一前端个体分布对应的非劣解集如表1所示,其中加粗栏为参数辨识结果。
表1
将参数辨识结果重新带入图2所示模型的运动微分方程中计算担架台垂向振动加速度的功率谱密度和概率分布(以“仿真结果”来表示),并与由试验数据计算获得的功率谱密度和概率分布(以“试验结果”来表示)进行对比,如图6和图7所示;观察图6和图7可以发现,试验结果与仿真结果基本吻合,尤其是在概率分布方面,吻合程度较高,这充分证明了本发明所提出的新型多目标动力系统参数辨识方法的有效性。
Claims (2)
1.一种新型多目标动力系统参数辨识方法,其特征在于,包括以下步骤:
(1)选定要进行参数识别的待辨识动力系统,根据待辨识动力系统的结构建立运动微分方程,并给出运动微分方程中已知参数的值和待辨识参数的范围及约束条件;
(2)选取多个待辨识动力系统输出响应的特征参数作为评价指标,并通过试验获得评价指标的值及待辨识动力系统的输入数据;
(3)应用改进的分块增维精细积分法和待辨识动力系统的输入数据,通过待辨识动力系统的运动微分方程计算评价指标的值,将计算获得的评价指标的值与试验获得的评价指标的值相减获得差值,以差值的平方和作为待辨识动力系统的目标函数;
(4)寻优计算,对NSGA-II算法进行相应的参数设定,应用NSGA-II算法对目标函数进行多次寻优计算,直至输出结果稳定,获得最终的非劣解集;
(5)选择结果,结合工程实际,在非劣解集中选取合适的点为最优点,以最优点对应的参数为最终辨识结果。
2.根据权利要求1所述的一种新型多目标动力系统参数辨识方法,其特征在于,步骤(5)中所述最优点的选取原则是最优点使不同的目标函数值都达到最小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510687739.9A CN105404609A (zh) | 2015-10-21 | 2015-10-21 | 一种新型多目标动力系统参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510687739.9A CN105404609A (zh) | 2015-10-21 | 2015-10-21 | 一种新型多目标动力系统参数辨识方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105404609A true CN105404609A (zh) | 2016-03-16 |
Family
ID=55470104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510687739.9A Pending CN105404609A (zh) | 2015-10-21 | 2015-10-21 | 一种新型多目标动力系统参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105404609A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908595A (zh) * | 2017-10-30 | 2018-04-13 | 燕山大学 | 基于卷积神经网络和数据驱动的非线性动态系统辨识方法 |
CN109271652A (zh) * | 2018-07-12 | 2019-01-25 | 哈尔滨工程大学 | 基于多目标进化算法的船舶横向三自由度运动参数辨识方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425811A (zh) * | 2012-11-05 | 2013-12-04 | 上海理工大学 | 数控机床进给系统的等效惯量和等效阻尼的辨识方法 |
-
2015
- 2015-10-21 CN CN201510687739.9A patent/CN105404609A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425811A (zh) * | 2012-11-05 | 2013-12-04 | 上海理工大学 | 数控机床进给系统的等效惯量和等效阻尼的辨识方法 |
Non-Patent Citations (2)
Title |
---|
PACHECO R P EL AL.: "On the identification of non-linear mechanical systems using orthogonal functions", 《INTERNATIONAL JOURNAL NON-LINEAR MECHANICS》 * |
杨猛 等: "履带急救车非线性减振系统动态性能优化", 《振动与冲击》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908595A (zh) * | 2017-10-30 | 2018-04-13 | 燕山大学 | 基于卷积神经网络和数据驱动的非线性动态系统辨识方法 |
CN107908595B (zh) * | 2017-10-30 | 2021-02-23 | 燕山大学 | 基于卷积神经网络和数据驱动的非线性动态系统辨识方法 |
CN109271652A (zh) * | 2018-07-12 | 2019-01-25 | 哈尔滨工程大学 | 基于多目标进化算法的船舶横向三自由度运动参数辨识方法 |
CN109271652B (zh) * | 2018-07-12 | 2022-05-31 | 哈尔滨工程大学 | 基于多目标进化算法的船舶横向三自由度运动参数辨识方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Du et al. | Fuzzy control for nonlinear uncertain electrohydraulic active suspensions with input constraint | |
Shook et al. | GA-optimized fuzzy logic control of a large-scale building for seismic loads | |
CN107357951B (zh) | 基于参数识别技术的桥梁冲击系数优化方法 | |
Liu et al. | Adaptive finite-time control for half-vehicle active suspension systems with uncertain dynamics | |
Jin et al. | Stochastic optimal active control of a half-car nonlinear suspension under random road excitation | |
Sibielak et al. | Optimal control of slow-active vehicle suspension—results of experimental data | |
Rong et al. | Transfer matrix method for dynamics modeling and independent modal space vibration control design of linear hybrid multibody system | |
CN114676598B (zh) | 汽车车身系统的整车道路耐久试验的加速方法和装置 | |
CN112949131B (zh) | 连续桥梁集群损伤诊断的概率损伤定位向量法 | |
CN102968665A (zh) | 并联机器人的正运动学求解方法 | |
Pang et al. | Stability analysis and fuzzy smith compensation control for semi-active suspension systems with time delay | |
Kaldas et al. | Rule optimized fuzzy logic controller for full vehicle semi-active suspension | |
CN105404609A (zh) | 一种新型多目标动力系统参数辨识方法 | |
Zhang et al. | Novel model‐free optimal active vibration control strategy based on deep reinforcement learning | |
Dessort et al. | Explicit model predictive control of semi-active suspension systems using Artificial Neural Networks (ANN) | |
Kaleemullah et al. | Comparative analysis of LQR and robust controller for active suspension | |
Wu et al. | HIS‐Based Semiactive Suspension Dual‐Frequency‐Range Switching Control to Improve Ride Comfort and Antiroll Performance | |
CN106125549A (zh) | 一种顶置武器站行进间连发射击炮口振动分步抑制方法 | |
Félix-Herrán et al. | Disturbance Rejection in a One‐Half Semiactive Vehicle Suspension by means of a Fuzzy-H∞ Controller | |
Wang | Active fault tolerance control for active suspension with control input time-delay and actuator gain variation fault | |
CN107655463A (zh) | 基于模拟退火算法的电子罗盘校准方法 | |
Likaj et al. | Optimisation and control of vehicle suspension using linear quadratic gaussian control | |
Vemireddy et al. | Development of a driving dynamics-oriented suspension design during the early concept phase | |
Qi et al. | Symmetrical valve controlled asymmetrical cylinder based on wavelet neural network | |
CN114722676B (zh) | 一种乘用车悬置系统垂向安装位置功率流优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160316 |
|
WD01 | Invention patent application deemed withdrawn after publication |