CN114397816B - 基于状态反馈x-LMS算法的发动机主动悬置控制方法 - Google Patents
基于状态反馈x-LMS算法的发动机主动悬置控制方法 Download PDFInfo
- Publication number
- CN114397816B CN114397816B CN202111537727.XA CN202111537727A CN114397816B CN 114397816 B CN114397816 B CN 114397816B CN 202111537727 A CN202111537727 A CN 202111537727A CN 114397816 B CN114397816 B CN 114397816B
- Authority
- CN
- China
- Prior art keywords
- cycle
- matrix
- kth
- kth cycle
- formula
- 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
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明公开了一种基于状态反馈x‑LMS算法的发动机主动悬置控制方法,其步骤包括:1基于Sage‑Husa卡尔曼滤波器的系统状态估计;2基于状态反馈的主动控制力的计算;3基于残余振动信号均方差的抽头权向量计算。本发明能够在不依赖激励信号的条件下,以系统状态为参考信号,实现一种具备全频带控制能力的主动悬置x‑LMS状态反馈控制方法。
Description
技术领域
本发明涉及机械系统动力学建模及振动控制领域,更具体地说是一种基于状态反馈x-LMS算法的发动机主动悬置控制方法。
背景技术
发动机作为汽车振动的主要振源之一,直接影响汽车的NVH性能。发动机悬置是减小发动机振动对汽车NVH性能的影响的最有效的手段。在不同频率条件下,发动机悬置性能的要求是不同的:在低频区大振幅振动时,为了限制动力总成的振动幅度并快速衰减振动,需要悬置元件具有较大的刚度和阻尼;在高频小振幅振动时,为了获得较低的振动传递率,需要悬置元件具有较小的刚度和阻尼。因此,理想的悬置需要能够根据发动机的运行工况,自适应调节自身参数以满足发动机不同的工作需求。发动机主动悬置在相应的控制方法作用下,能够在低频时利用被动悬置的隔振性能获得大刚度、大阻尼,来抑制发动机的位移幅值,而在高频时利用作动器的作动力,降低动刚度来有效隔离发动机振动,降低车内噪音。因此,在发动机主动振动控制技术中,主动振动控制方法直接决定了主动振动控制的效果,是主动振动控制的核心技术之一。
在发动机主动悬置控制技术方面,已经有大量的控制方法得到了应用PID方法、LQR方法、模型参考控制方法、H2和H∞两种鲁棒控制方法、滑膜控制方法、线性矩阵不等式控制方法等。但是这些主动振动控制方法所擅长的频率范围是不同,比如基于随机干扰假设的最优控制,在原理上存在不同程度的低频相对位移控制和高频隔振性能的折中处理,不能适应宽频带振动控制。
发明内容
本发明为解决上述现有技术所存在的不足,提出一种具有计算量小、易于实现且自适应能力强的基于状态反馈x-LMS算法的发动机主动悬置控制方法,以期能够实现整个频率范围内的最优控制,从而提高发动机主动悬置的控制性能,并能克服传统x-LMS算法以激励信号作为参考信号的缺点。
本发明为解决技术问题采用如下技术方案:
本发明一种基于状态反馈x-LMS算法的发动机主动悬置控制方法的特点在于,按如下步骤进行:
步骤1:采集振动数据:
步骤1.1:设定循环次数为k,并且初始化k=1;
步骤1.2:定义第k次循环的状态向量为Xk=[x1.k x2.k x3.k x4.k]T,其中,x1.k为第k次循环的发动机振动位移,x2.k为第k次循环的发动机振动速度,x3.k为第k次循环的簧载质量振动位移,x4.k为第k次循环的簧载质量振动速度;
步骤1.3:由位移传感器测量第k次循环的发动机与簧载质量之间的相对位移xr.k,由力传感器测量第k次循环的发动机振动时通过主动悬置传递给簧载质量的动态力Fd.k;
步骤2:根据发动机、主动悬置、簧载质量、悬架组成的发动机主动悬置系统的离散状态方程,建立如式(1)-式(9)所示的Sage-Husa卡尔曼滤波器,用于估计状态变量值:
步骤2.1、利用式(1)计算第k次循环的一步预测均方误差Pk,k-1:
步骤2.3、利用式(3)计算第k次循环的残差vk:
式(3)中,e(k)表示第k次循环的残余振动,且e(k)=[xr.k,Fd.k]T,C为第一输出系数矩阵,D为第二输出系数矩阵,rk为第k次循环的测量噪声ηk的均值;
步骤2.4、利用式(4)计算第k次循环的增益矩阵Kk:
Kk=Pk,k-1CT[CPk,k-1CT+Rk]-1 (4)
式(4)中,CT为第一输出系数矩阵C的转置矩阵,Rk为第k次循环的测量噪声ηk的协方差;
步骤2.5、利用式(5)计算第k次循环的均方误差Pk:
Pk=(I-KkC)Pk(I-KkC)T+KkRkKk T (5)
式(5)中,I为单位矩阵,Kk T为增益矩阵Kk的转置矩阵;
步骤2.7、利用式(7)计算第k次循环的加权系数dk:
dk=(1-b)/(1-b(k+1)) (7)
式(7)中,b为遗忘因子;
步骤2.8、利用式(8)计算第k次循环的第二输入变量Fin.k的均值:
步骤2.9、利用式(9)计算第k次循环的第二输入变量Fin.k的协方差矩阵:
步骤3:利用式(10)计算第k次循环的第一输入变量Fc.k:
式(10)中,为第k次循环的第j个状态变量xj.k的估计值对应的N×1阶矢量,且其中,/>为第k-(N-1)次循环第j个状态变量xj.k-(N-1)的估计值,N为LMS滤波器的阶数,/>为/>的转置向量,wj.k为第k次循环的N×1阶抽头权向量,且wj.k=[wj0.k wj1.k wj2.k…wjN-1.k]T,wjN-1.k为第k次循环第N-1个抽头权值;
步骤5:利用式(12)得到第k+1次循环的抽头向量wj.k+1:
式(12)中,L为目标权值矩阵;λ固定收敛因子,γ为控制因子;
步骤6:所述主动悬置根据所述第k次循环的第一输入变量Fc.k输出相应主动控制力;
步骤7:将k+1赋值给k后,返回步骤1.2执行,直至由位移传感器测量得到发动机与簧载质量之间的相对位移xr.k以及力传感器测量得到发动机振动通过主动悬置传递给簧载质量的动态力Fd.k的均方差值最小为止。
与已有技术相比,本发明有益效果体现在:
1、本发明基于状态反馈x-LMS算法,以发动机主动悬置系统的状态作为参考信号,克服了传统x-LMS算法依赖激励信号作为参考信号的缺点,并具有传统x-LMS算法的优点,适应于发动机主动悬置的宽频带控制。
2、本发明基于状态反馈x-LMS算法针对发动机主动悬置系统的状态变量并不是全部可测的,采用步骤2所述的Sage-Husa卡尔曼滤波器来实现发动机主动悬置系统状态的估计,控制方法适应能力强。
附图说明
图1为本发明所涉的发动机主动悬置系统的示意图;
图2为本发明所涉的基于状态反馈控制x-LMS算法的发动机主动悬置系统的控制原理。
具体实施方式
本实施例中,一种基于状态反馈x-LMS算法的发动机主动悬置控制方法,是先基于Sage-Husa卡尔曼滤波器的系统状态估计;再基于状态反馈的主动控制力的计算;最后基于残余振动信号均方差的抽头权向量计算。具体的说,是按如下步骤进行:
步骤1:如图1所示,利用式(1)建立发动机主动悬置系统动力学模型:
式(1)中,X为状态向量且X=[x1 x2 x3 x4]T,为X的导数,x1为发动机振动位移,x2为发动机振动速度且/> 为x1的导数,x3为簧载质量振动位移,x4为簧载质量振动速度且/> 为x3的导数,A为状态矩阵,B1为主动控制力输入矩阵,B2为激振力输入矩阵,e为残余振动,即发动机与簧载质量之间的相对位移和发动机振动时通过主动悬置传递给簧载质量的动态力,C为观测矩阵,D为传递矩阵,Fc为主动控制力,Fin为发动机垂向激振力,η为测量噪声,me和mb分别为四分之一发动机质量及簧载质量,ks和bs分别为悬架刚度和阻尼,km和bm分别为主动悬置刚度和阻尼。
步骤2:对式(1)表示的发动机主动悬置系统动力学模型离散化:
式(2)中,G=expm(AT),H1=A-1·(G-I)·B1,H2=A-1·(G-I)·B2,Xk=[x1.k x2.kx3.k x4.k]T为第k次循环的状态向量,其中,x1.k为第k次发动机振动位移,x2.k为第k次循环的发动机振动速度,x3.k为第k次循环的簧载质量振动位移,x4.k为第k次循环的簧载质量振动速度。Xk-1为第k-1次循环的状态向量,Fc.k-1为第k-1次循环的主动控制力,Fin.k-1为第k-1次循环的发动机垂向激振力,Fc.k为第k次循环的主动控制力,ηk为第k次循环的测量噪声,ek为第k次循环的残余振动。
步骤3:采集振动数据:
步骤3.1:设定循环次数为k,并初始化k=1;
步骤3.3:由位移传感器测量第k次循环的发动机与簧载质量之间的相对位移xr.k,由力传感器测量第k次循环的发动机振动时通过主动悬置传递给簧载质量的动态力Fd.k;
步骤4:根据式(2)所示发动机主动悬置系统离散状态方程,建立如式(3)-式(11)所示的Sage-Husa卡尔曼滤波器,用于估计状态变量值。
步骤4.1、利用式(3)计算第k次循环的一步预测均方误差Pk,k-1:
步骤4.3、利用式(5)计算第k次循环的残差vk:
式(5)中,e(k)=[xr.k,Fd.k]T,rk为第k次循环的测量噪声ηk的均值。
步骤4.4、利用式(6)计算第k次循环的增益矩阵Kk:
Kk=Pk,k-1CT[CPk,k-1CT+Rk]-1 (6)
式(6)中,CT为观测矩阵C的转置矩阵,Rk为第k次循环的测量噪声ηk的协方差。
步骤4.5、利用式(7)计算第k次循环的均方误差Pk:
Pk=(I-KkC)Pk(I-KkC)T+KkRkKk T (7)
式(7)中,I为单位矩阵,Kk T为增益矩阵Kk的转置矩阵。
步骤4.7、利用式(9)计算第k次循环的加权系数dk:
dk=(1-b)/(1-b(k+1)) (9)
式(9)中,b为遗忘因子。
步骤4.8、利用式(10)计算第k次循环的发动机垂向激振力Fin.k的均值:
步骤4.9、利用式(11)计算第k次循环的发动机垂向激振力Fin.k协方差矩阵:
步骤5:利用式(12)计算第k次循环的主动控制力Fc.k:
式(12)中,为第k次循环的第j个状态变量xj.k的估计值对应的N×1阶矢量,其中,/>为第k-(N-1)次循环第j个状态变量xj.k-(N-1)的估计值,N为LMS滤波器阶数,/>为/>的转置向量,wj.k为第k次循环的N×1阶抽头权向量,wj.k=[wj0.k wj1.k wj2.k…wjN-1.k]T,wjN-1.k为第k次循环第N-1个抽头权值。
步骤7:利用式(14)得到第k+1次循环的抽头向量wj.k+1:
式(14)中,L为目标权值矩阵;λ固定收敛因子,γ为避免fj(k)Tfj(k)过小的控制因子。
步骤8:主动悬置输出第k次循环的主动控制力Fc.k;
步骤9:k+1赋值给k后,返回步骤3.2执行,直至由位移传感器测量得到发动机与簧载质量之间的相对位移xr.k以及力传感器测量得到发动机振动通过主动悬置传递给簧载质量的动态力Fd.k的均方差值最小为止。
Claims (1)
1.一种基于状态反馈x-LMS算法的发动机主动悬置控制方法,其特征在于,按如下步骤进行:
步骤1:采集振动数据:
步骤1.1:设定循环次数为k,并且初始化k=1;
步骤1.2:定义第k次循环的状态向量为Xk=[x1.k x2.k x3.k x4.k]T,其中,x1.k为第k次循环的发动机振动位移,x2.k为第k次循环的发动机振动速度,x3.k为第k次循环的簧载质量振动位移,x4.k为第k次循环的簧载质量振动速度;
步骤1.3:由位移传感器测量第k次循环的发动机与簧载质量之间的相对位移xr.k,由力传感器测量第k次循环的发动机振动时通过主动悬置传递给簧载质量的动态力Fd.k;
步骤2:根据发动机、主动悬置、簧载质量、悬架组成的发动机主动悬置系统的离散状态方程,建立如式(1)-式(9)所示的Sage-Husa卡尔曼滤波器,用于估计状态变量值:
步骤2.1、利用式(1)计算第k次循环的一步预测均方误差Pk,k-1:
步骤2.3、利用式(3)计算第k次循环的残差vk:
式(3)中,e(k)表示第k次循环的残余振动,且e(k)=[xr.k,Fd.k]T,C为第一输出系数矩阵,D为第二输出系数矩阵,rk为第k次循环的测量噪声ηk的均值;
步骤2.4、利用式(4)计算第k次循环的增益矩阵Kk:
Kk=Pk,k-1CT[CPk,k-1CT+Rk]-1 (4)
式(4)中,CT为第一输出系数矩阵C的转置矩阵,Rk为第k次循环的测量噪声ηk的协方差;
步骤2.5、利用式(5)计算第k次循环的均方误差Pk:
Pk=(I-KkC)Pk(I-KkC)T+KkRkKk T (5)
式(5)中,I为单位矩阵,Kk T为增益矩阵Kk的转置矩阵;
步骤2.7、利用式(7)计算第k次循环的加权系数dk:
dk=(1-b)/(1-b(k+1)) (7)
式(7)中,b为遗忘因子;
步骤2.8、利用式(8)计算第k次循环的第二输入变量Fin.k的均值:
步骤2.9、利用式(9)计算第k次循环的第二输入变量Fin.k的协方差矩阵:
步骤3:利用式(10)计算第k次循环的第一输入变量Fc.k:
式(10)中,为第k次循环的第j个状态变量xj.k的估计值对应的N×1阶矢量,且其中,/>为第k-(N-1)次循环第j个状态变量xj.k-(N-1)的估计值,N为LMS滤波器的阶数,/>为/>的转置向量,wj.k为第k次循环的N×1阶抽头权向量,且wj.k=[wj0.k wj1.k wj2.k…wjN-1.k]T,wjN-1.k为第k次循环第N-1个抽头权值;
步骤5:利用式(12)得到第k+1次循环的抽头向量wj.k+1:
式(12)中,L为目标权值矩阵;λ固定收敛因子,γ为控制因子;
步骤6:所述主动悬置根据所述第k次循环的第一输入变量Fc.k输出相应主动控制力;
步骤7:将k+1赋值给k后,返回步骤1.2执行,直至由位移传感器测量得到发动机与簧载质量之间的相对位移xr.k以及力传感器测量得到发动机振动通过主动悬置传递给簧载质量的动态力Fd.k的均方差值最小为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111537727.XA CN114397816B (zh) | 2021-12-15 | 2021-12-15 | 基于状态反馈x-LMS算法的发动机主动悬置控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111537727.XA CN114397816B (zh) | 2021-12-15 | 2021-12-15 | 基于状态反馈x-LMS算法的发动机主动悬置控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114397816A CN114397816A (zh) | 2022-04-26 |
CN114397816B true CN114397816B (zh) | 2023-06-27 |
Family
ID=81226304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111537727.XA Active CN114397816B (zh) | 2021-12-15 | 2021-12-15 | 基于状态反馈x-LMS算法的发动机主动悬置控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114397816B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416845A (en) * | 1993-04-27 | 1995-05-16 | Noise Cancellation Technologies, Inc. | Single and multiple channel block adaptive methods and apparatus for active sound and vibration control |
CN105738817A (zh) * | 2016-01-29 | 2016-07-06 | 深圳市沃特玛电池有限公司 | 一种基于aekf的电池荷电状态估计方法及估计系统 |
WO2020245902A1 (ja) * | 2019-06-04 | 2020-12-10 | 三菱電機株式会社 | 方位角推定装置、ジャイロシステム、方位角推定方法及びプログラム |
-
2021
- 2021-12-15 CN CN202111537727.XA patent/CN114397816B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416845A (en) * | 1993-04-27 | 1995-05-16 | Noise Cancellation Technologies, Inc. | Single and multiple channel block adaptive methods and apparatus for active sound and vibration control |
CN105738817A (zh) * | 2016-01-29 | 2016-07-06 | 深圳市沃特玛电池有限公司 | 一种基于aekf的电池荷电状态估计方法及估计系统 |
WO2020245902A1 (ja) * | 2019-06-04 | 2020-12-10 | 三菱電機株式会社 | 方位角推定装置、ジャイロシステム、方位角推定方法及びプログラム |
Non-Patent Citations (1)
Title |
---|
基于卡尔曼滤波算法的车辆振动状态估计与最优控制研究;钟孝伟;陈双;张不扬;;汽车技术(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114397816A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pinte et al. | A piezo-based bearing for the active structural acoustic control of rotating machinery | |
JP3451891B2 (ja) | 能動型振動制御装置 | |
US7974769B2 (en) | Integrated open and closed-loop control method for active engine mounts | |
JP4857897B2 (ja) | 騒音制御方法および騒音制御装置 | |
CN104993766B (zh) | 一种二质量系统谐振抑制方法 | |
CN109426150B (zh) | 基于扩张状态观测器的负载模拟器反步控制方法 | |
Svaricek et al. | Automotive applications of active vibration control | |
CN110829885A (zh) | 一种磁致伸缩精密驱动装置的机械阻抗匹配控制方法 | |
CN114397816B (zh) | 基于状态反馈x-LMS算法的发动机主动悬置控制方法 | |
CN107102658B (zh) | 一种离线次级通道的辨识方法及装置 | |
CN110045603B (zh) | 压电陶瓷驱动部件变载荷环境下的鲁棒自适应控制方法 | |
CN113467236B (zh) | 一种对误差信号进行时滞补偿的方法 | |
CN112731814B (zh) | 基于自适应谐波识别频响修正的直升机振动主动控制方法 | |
JP4857907B2 (ja) | 騒音制御装置および騒音制御方法 | |
CN110928184B (zh) | 面向军用计算机应用的主动减振控制方法及装置 | |
Zhang et al. | Simulation and experiment on active vibration isolation with an adaptive method | |
CN107807532B (zh) | 一种用于超磁致伸缩隔振平台的自适应逆隔振控制方法 | |
CN114035627B (zh) | 基于单自由度隔振平台的主动复合控制系统及方法 | |
Qiu et al. | Adaptive resonant vibration control of a piezoelectric flexible plate implementing Filtered-X LMS algorithm | |
JP4857928B2 (ja) | 騒音制御装置および騒音制御方法 | |
Ján et al. | Dynamic Properties Modeling Analysis of the Rubber-Metal Elements for Electric Drive | |
Riley et al. | An adaptive strategy for vehicle vibration and noise cancellation | |
Li et al. | DOB-based piezoelectric vibration control for stiffened plate considering accelerometer measurement noise | |
Hausberg et al. | Incorporation of adaptive grid-based look-up tables in adaptive feedforward algorithms for active engine mounts | |
Beijen et al. | Filtered-error recursive least squares optimization for disturbance feedforward control in active vibration isolation |
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 |